Er³⁺掺杂锗酸盐激光玻璃猝灭浓度的计算与预测 增强出版

Er 3+ 掺杂锗酸盐激光玻璃猝灭浓度的计算与预测 Calculation and prediction of quenching concentration of Er 3+ -doped germanate glass first-author 万杰（1996－），女，重庆开州人，博士研究生，2019年于重庆交通大学获得学士学位，主要从事稀土掺杂激光玻璃的研究。E-mail： wanjie_12345@163.com http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40845456&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40845458&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40845457&type= 21.06447792 31.51574707 万杰 （1996－），女，重庆开州人，博士研究生，2019年于重庆交通大学获得学士学位，主要从事稀土掺杂激光玻璃的研究。E-mail： wanjie_12345@163.com wanjie_12345@163.com corresp E-mail： wangweichao@scut.edu.cn E-mail： wangweichao@scut.edu.cn 王伟超（1988－），男，陕西咸阳人，博士研究生，副教授，2017年于华南理工大学获得博士学位，主要从事特种激光玻璃及光纤激光器的研究。E-mail： wangweichao@scut.edu.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40845459&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40845461&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40845460&type= 21.42521667 31.21496582 王伟超 （1988－），男，陕西咸阳人，博士研究生，副教授，2017年于华南理工大学获得博士学位，主要从事特种激光玻璃及光纤激光器的研究。E-mail： wangweichao@scut.edu.cn wangweichao@scut.edu.cn 华南理工大学，物理与光电学院，发光材料与器件国家重点实验室， 广东省光纤激光材料与应用技术重点实验室，　广东 广州510640 华南理工大学，物理与光电学院，发光材料与器件国家重点实验室， 广东省光纤激光材料与应用技术重点实验室，　广东 广州510640 State Key Laboratory of Luminescent Materials and Devices， Guangdong Provincial Key Laboratory ofFiber Laser Materials and Applied Techniques， School of Physics and Optoelectronics， South China University of Technology， Guangzhou 510640， China State Key Laboratory of Luminescent Materials and Devices， Guangdong Provincial Key Laboratory ofFiber Laser Materials and Applied Techniques， School of Physics and Optoelectronics， South China University of Technology， Guangzhou 510640， China 稀土掺杂激光玻璃光纤是光纤激光器的核心增益介质，其中稀土离子掺杂浓度是决定激光玻璃增益特性的重要参数，如何快速有效地确定最佳掺杂浓度是稀土掺杂激光玻璃光纤的关键科学问题之一。本文以Er 3+ 掺杂锗酸盐激光玻璃为例，利用Er 3+ ： 4 I 13/2 → 4 I 15/2 跃迁的自发辐射跃迁寿命和实测寿命，预测了锗酸盐激光玻璃的猝灭浓度。研究表明，在该玻璃体系中，理论预测与实际值的最大绝对误差小于0.4 wt.%。对比唯象模型和有限扩散模型分别拟合多个样品发光强度和实测寿命确定猝灭浓度的方法，本文提出的方法仅需通过少量样品的测试参数即可确定激光玻璃猝灭浓度，简单快捷且计算误差小，对高增益激光玻璃与光纤研究具有指导意义。 Rare-earth doped laser glass fiber is the core gain medium of fiber laser.The concentration of the rare-earth ions is one of the important parameters to determine the gain characteristics of laser glass. How to determine the quenching concentration（QC） of laser glass quickly and effectively is a key scientific problem. In this paper， the QC of Er 3+ -doped germanateglass is predicted by using the spontaneous emission lifetime and the measured lifetime of the Er 3+ ： 4 I 13/2 → 4 I 15/2 transition. The results show that the maximum absolute error between the theoretical prediction and the actual value is less than 0.4 wt.%.Compared with the phenomenological model and finite diffusion model to determine the QC by fitting the luminescent intensity and measured lifetime of multiple experimental samples， our proposed method can determine the QC of laser glass by the test parameters of one or two samples， which is simple， fast， and has little calculation error. It has a guiding significance for the research of high-gain laser glass and optical fiber. 锗酸盐玻璃 稀土离子 Er 3+ 猝灭浓度 理论计算 germanate glass rare-earth ions Er 3+ quenching concentration theoretical calculation 1　引　　言 1 光纤激光器是大功率激光、卫星激光通信、引力波探测等国家安全与科学前沿领域发展的迫切与重大需求 [ 1 1 - 3 3 1-3 ] 。稀土掺杂激光玻璃光纤是光纤激光器的核心增益介质，稀土离子掺杂浓度是决定激光玻璃增益特性的重要参数。提高激光玻璃与光纤的稀土离子掺杂浓度，有利于增加激发态上能级粒子数，从而提高激光输出功率和斜率效率，同时可以缩短激光腔长而有利于器件的小型化和集约化。然而，如何快速有效地确定稀土最佳掺杂浓度是稀土掺杂激光玻璃光纤的关键科学问题之一 [ 4 4 - 5 5 4-5 ] 。 掺杂离子发光猝灭主要由多声子弛豫和浓度猝灭所致。其中，多声子弛豫发生几率与掺杂浓度相关性较小，更多地取决于玻璃基质的最大声子能量。浓度猝灭作为限制发光效率的重要因素，主要由掺杂浓度引起。当掺杂离子浓度增加到一定值时，离子间过强的相互作用一方面会增加交叉弛豫等能量传递过程. 另一方面，也会增加能量传递到杂质（包括羟基、过渡金属离子、其他稀土离子）的几率，从而导致非辐射跃迁作用增强，使得激发态粒子布居数降低，最终导致上能级寿命和发光强度等降低 [ 6 6 ] 。当前，研究者们一般通过实验，基于不同稀土离子浓度掺杂，判断荧光强度最大值对应的浓度为猝灭浓度。该方法依赖一系列实验判定，准确性较高。然而，由于缺乏理论指导，实际工作中实验周期长、效率低下、成本高。探寻简洁快速的理论预测研究方法，对研发新材料具有重要的理论和实际意义。 美国贝尔实验室Van Uitert [ 7 7 ] 建立了一种唯象模型描述了发光强度与掺杂离子浓度之间的关系，但并未涉及猝灭浓度计算；Auzel等 [ 8 8 - 9 9 8-9 ] 在稀土掺杂单晶中利用有限扩散模型提出了荧光寿命与掺杂浓度之间的定量关系公式，一些学者也将其应用在激光玻璃中计算理论猝灭浓度 [ 10 10 - 12 12 10-12 ] 。黄莉蕾等 [ 13 13 - 14 14 13-14 ] 基于Burshtein模型描述的荧光衰减规律推导了激发态粒子布居数与掺杂浓度的关系，利用少量实验数据计算了不同稀土离子单掺YAG激光晶体的最佳掺杂浓度，为理论预测猝灭浓度提供了一种新的方法 [ 15 15 ] 。徐淑君等 [ 16 16 ] 利用效应面优化模型，优化了红色荧光粉中两种离子的最佳掺杂浓度，提供了一种直接寻找荧光粉最佳掺杂浓度的方法。上述关于猝灭浓度理论预测的研究主要集中于晶体材料，鲜有关于激光玻璃的报道。不同于晶体，玻璃中稀土离子处于复杂的网络结构中，难以建立玻璃结构与发光性质之间的定量关系，这也是本工作的难点之一。目前激光玻璃猝灭浓度的理论预测研究尚处于初步阶段，建立简单快捷准确的猝灭浓度理论预测方法仍需要大量的研究探索。 锗酸盐玻璃具有声子能量低、红外透过率高、玻璃形成能力好、热稳定性和化学耐久性好等特点。 [ 17 17 ] Er 3+ ： 4 I 13/2 → 4 I 15/2 跃迁产生的1.5 μm发光波长位于第三通信窗口，在通信和传感等领域具有重要应用 [ 18 18 ] 。因此，本文以Er 3+ 掺杂锗酸盐激光玻璃为例，利用Er 3+ ： 4 I 13/2 → 4 I 15/2 跃迁的自发辐射跃迁寿命和实测寿命计算了猝灭浓度，通过实验值进行了验证与分析对比，并讨论了该方法与其他不同方法计算猝灭浓度的简易性与准确性。 2　计算猝灭浓度理论基础 1 2.1　荧光强度衰减规律 2 考虑稀土离子的无辐射跃迁作用 W ，荧光强度衰减规律为 [ 19 19 - 22 22 19-22 ] ： <math display="block" id="M1"><mi>I</mi><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo><mi>I</mi><mo stretchy="false">(</mo><mn mathvariant="normal">0</mn><mo stretchy="false">)</mo><msup><mrow><mi>e</mi></mrow><mrow><mo>-</mo><mfrac><mrow><mi>t</mi></mrow><mrow><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow></mfrac><mo>-</mo><mi>W</mi><mi>t</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40845464&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40845462&type= 26.07733345 6.01133299 其中， τ 0 是Judd-Ofelt(简记为J-O)理论计算的自发辐射跃迁寿命， I （0）和 I （ t ）分别表示时间为0和 t 时的荧光强度。在不考虑其他稀土离子与过渡金属离子引起的非辐射作用时， W 主要由稀土与基质晶格振动之间相互作用的多声子弛豫几率 W MP 、稀土离子间能量传递几率 W ET 、稀土离子和OH - 基团间的能量传递几率 W OH - 构成，即 W=W MP +W ET +W OH - 。若只考虑无辐射跃迁作用中稀土离子间的能量传递，Inokuti-Hirayama（简记为I-H）和Burshtein（简记为BN）分别假设稀土离子间主要通过供体-受体、供体-供体进行能量传递和能量迁移，提出了荧光强度衰减规律的I-H模型（ 公式（2） 公式（2） ）和BN模型（ 公式（3） 公式（3） ） [ 23 23 - 24 24 23-24 ] ： <math display="block" id="M2"><mi>I</mi><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo><mi>I</mi><mo stretchy="false">(</mo><mn mathvariant="normal">0</mn><mo stretchy="false">)</mo><msup><mrow><mi>e</mi></mrow><mrow><mo>-</mo><mfrac><mrow><mi>t</mi></mrow><mrow><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow></mfrac><mo>-</mo><mi>γ</mi><mroot><mrow><mi>t</mi></mrow><mrow/></mroot></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40845466&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40845465&type= 27.68600082 6.01133299 <math display="block" id="M3"><mi>I</mi><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo><mi>I</mi><mo stretchy="false">(</mo><mn mathvariant="normal">0</mn><mo stretchy="false">)</mo><msup><mrow><mi>e</mi></mrow><mrow><mo>-</mo><mfrac><mrow><mi>t</mi></mrow><mrow><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow></mfrac><mo>-</mo><mi>γ</mi><mroot><mrow><mi>t</mi></mrow><mrow/></mroot><mo>-</mo><mi>ω</mi><mi>t</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40845468&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40845467&type= 31.91933441 6.01133299 其中， γ 为稀土离子间的偶极-偶极相互作用， γ=γ 0 x A ， x A 是受体离子的相对浓度； ω 为稀土离子间的迁移作用， ω=ω 0 x A x D ， x D 是施主离子的相对浓度， γ 0 和 ω 0 均为不随掺杂离子浓度变化的常数。在单掺情况下 ω=ω 0 x 2 ， γ=γ 0 x ， x 为稀土离子掺杂浓度。 由于激发态布居数 N （ t ）与荧光强度 I （ t ）成正比， N （ t ）可写成 公式（3） 公式（3） 的形式。在连续抽运情况下， t 时刻激发态布居数可看作0- t 时刻无数个 δ 脉冲在 t 时刻抽运的累积结果，即： <math display="block" id="M4"><mi>N</mi><mo>∝</mo><mstyle displaystyle="true"><munderover><mo>∫</mo><mrow><mn mathvariant="normal">0</mn></mrow><mrow><mi>∞</mi></mrow></munderover></mstyle><mrow><mi>N</mi><mo stretchy="false">(</mo></mrow><mi>t</mi><mo stretchy="false">)</mo><mi mathvariant="normal">d</mi><mi mathvariant="normal">t</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40845471&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40845470&type= 19.98133469 9.90600014 故得到在单掺时激发态粒子布居数 N 与掺杂浓度 x 之间的关系 [ 13 13 ] : <math display="block" id="M5"><mi>N</mi><mo>∝</mo><mfrac><mrow><mi>x</mi><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mfenced open="[" close="]" separators="|"><mrow><mn mathvariant="normal">2</mn><mo>+</mo><mn mathvariant="normal">2</mn><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msub><mrow><mi>ω</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msup><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>-</mo><mroot><mrow><mi>π</mi></mrow><mrow/></mroot><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msub><mrow><mi>γ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mi>x</mi><mroot><mrow><msub><mrow><mi>ω</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msup><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow></mfrac></mrow><mrow/></mroot><mi>e</mi><mi>r</mi><mi>f</mi><mi>c</mi><mo stretchy="false">(</mo><mfrac><mrow><msub><mrow><mi>γ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mi>x</mi></mrow><mrow><mn mathvariant="normal">2</mn><mroot><mrow><msub><mrow><mi>ω</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msup><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow></mfrac></mrow><mrow/></mroot></mrow></mfrac><mo stretchy="false">)</mo><msup><mrow><mi>e</mi></mrow><mrow><mfrac><mrow><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msubsup><mrow><mi>γ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><msup><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow><mn mathvariant="normal">4</mn><mo>+</mo><mn mathvariant="normal">4</mn><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msub><mrow><mi>ω</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msup><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac></mrow></msup></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo>+</mo><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msub><mrow><mi>ω</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msup><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><msup><mrow><mo stretchy="false">)</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40845473&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40845472&type= 116.07800293 26.58533096 其中erfc( x )为互补误差函数。根据荧光衰减规律模型拟合荧光寿命曲线得到参数（ ω 、 γ、ω 0 和 γ 0 ），再建立 N 与 x 之间的关系曲线，可以在曲线极值处得到理论猝灭浓度 X q 。需要指出的是： 公式（5） 公式（5） 只考虑了激发态粒子布居数，未考虑激光下能级，因此该方法主要适用于激光下能级为基态或不考虑激光下能级的情况。 2.2　理论猝灭浓度近似解 2 进一步简化 公式（5） 公式（5） 得到三种特殊情况：（1）当 ω =0， γ =0时， N ∝ x ，即没有离子间的相互作用。激发态粒子布居数与掺杂浓度成正比，这与发光效率和掺杂浓度成正比的规律一致。（2）当 ω =0， γ ≠0时，离子间相互作用以偶极-偶极相互作用为主，黄莉蕾等得到猝灭浓度近似解 [ 14 14 ] 。（3）当 γ =0， ω ≠0时，离子间相互作用以迁移作用为主，本文作者推导猝灭浓度 X q 近似解： <math display="block" id="M6"><msub><mrow><mi>X</mi></mrow><mrow><mi>q</mi></mrow></msub><mo>=</mo><mroot><mrow><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><msub><mrow><mi>ω</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow></mfrac></mrow><mrow/></mroot><mo>=</mo><mroot><mrow><mfrac><mrow><msubsup><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></mrow><mrow><mi>ω</mi><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow></mfrac></mrow><mrow/></mroot></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40845475&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40845474&type= 32.76599884 9.90600014 当荧光寿命符合单指数衰减时（ γ =0），结合 公式（1） 公式（1） 和（3）得到 W = ω ，进一步推导 X q ： <math display="block" id="M7"><msub><mrow><mi>X</mi></mrow><mrow><mi>q</mi></mrow></msub><mo>=</mo><mroot><mrow><mfrac><mrow><msubsup><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></mrow><mrow><mfrac><mrow><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mrow><msub><mrow><mi>τ</mi></mrow><mrow><mi>m</mi></mrow></msub></mrow></mfrac><mo>-</mo><mn mathvariant="normal">1</mn></mrow></mfrac></mrow><mrow/></mroot></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40845477&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=40845476&type= 22.52133369 14.56266594 相比于 公式（6） 公式（6） 和 公式（7） 公式（7） 的 X q 可根据掺杂浓度 x 0 以及相应的自发辐射跃迁寿命 τ 0 和实测寿命 τ m 计算。值得注意的是，BN模型忽略了多声子弛豫以及OH - 引起的非辐射跃迁，只考虑了稀土离子间的迁移作用，所以基于BN模型推导的 公式（6） 公式（6） 和（7）存在一定的理论近似和误差。为了直观描述计算结果的准确性，本文采用绝对误差而非相对误差。由于稀土离子浓度变化间隔一般取0.5 wt.%或1 wt.%，因此以绝对误差在±0.5 wt.%范围内判定为精确计算，±1 wt.%范围内为准确计算，大于±1 wt.%为不可以计算。 3　实　　验 1 3.1　样品制备 2 采用熔融-冷却法制备Er 3+ 掺杂锗酸盐玻璃，其玻璃摩尔组成为70GeO 2 -10Ga 2 O 3 -20BaO（简记为GGB玻璃），外掺1、2、3、4、5wt.%Er 2 O 3 。原料包括GeO 2 （99.99%，Aladdin），Ga 2 O 3 （99.99%，Aladdin），BaCO 3 （99.99%，Aladdin）和Er 2 O 3 （99.99%，Aladdin）。所有玻璃样品按如下步骤制备：根据玻璃组分准确称取原料15 g，将称好的原料转移至玛瑙研钵中均匀混合，置于刚玉坩埚，放入1400℃的玻璃熔炉中保温30 min后得到锗酸盐熔体，浇注于预热石墨模具上成型；随后于600℃的马弗炉中退火2 h。所有样品切割研磨至约1 mm厚并抛光以备测试表征。 3.2　样品表征 2 样品的密度根据阿基米德原理确定，浸泡液为蒸馏水，测量精度为±0.01g/cm 3 。折射率采用Metricon Model 2010棱镜耦合仪测量，准确度为±0.0005。吸收光谱由Perkin-Elmer Lambda 900紫外/可见/近红外双光束分光光度计（Waltham，MA）测定，分辨率为1 nm。荧光光谱通过iHR320光谱仪（JobinYvon Corp.，Horiba Scientific）和液氮冷却的PbSe探测器在980nm LD激发下测得。荧光衰减曲线用数字示波器（Tektronix TDS3012C）记录。红外透过光谱用傅里叶红外光谱仪（Vector-33, Bruker, 瑞士）测量。所有测试在室温下进行。 4　结果与讨论 1 4.1　实验结果 2 图1 图1 是不同浓度Er 2 O 3 掺杂GGB玻璃的光谱性质。根据 图1 图1 (a)吸收光谱的366、378、407、489、521、546、652、801、979 nm吸收峰，结合J-O理论采用最小二乘法拟合得到Er 3+ 离子的J-O强度参数 Ω t （t=2、4、6），拟合得到的均方差 δ 仅为1.22×10 -6 ，表明拟合结果可以接受，采用的约化矩阵元见参考文献 [ 25 25 ] 。 Ω t 反映了稀土离子周围的局域结构和键合特性， Ω 2 越大，稀土离子周围局域环境的不对称性越小， Ω 6 越大，稀土离子与配体阴离子之间的共价性越弱 [ 26 26 , 27 27 ] 。 Ω t 还可以进一步计算用于拟合荧光寿命曲线及预测猝灭浓度的 τ 0 ，计算结果见 表1 表1 。从表中可以看出，J-O强度参数和自发辐射跃迁寿命随掺杂浓度增加变化很小，表明玻璃中稀土离子周围配位场受稀土掺杂浓度的影响较小，自发辐射跃迁寿命可认为不随掺杂浓度变化。 x 0 /（wt.%） Ω 2 /（10 –20 cm 2 ） Ω 4 /（10 –20 cm 2 ） Ω 6 /（10 –20 cm 2 ） τ 0 /（ms） τ m /（ms） R 2 - 1 7.75 1.69 0.74 6.70 5.68 0.999 2 7.20 1.86 0.79 6.63 4.20 0.999 3 8.18 2.09 0.96 6.30 3.26 0.999 4 7.01 1.99 0.76 6.50 2.58 0.999 5 6.98 2.01 0.89 6.43 2.11 0.999 Er 3+ 的1.5 μm发光属于三能级系统，其发光强度由上能级（ 4 I 13/2 ）粒子布居数决定，这满足 公式（5） 公式（5） 的前提条件。 图1 图1 (b)是980 nmLD泵浦Er 3+ 的荧光发射光谱，随Er 2 O 3 浓度增加，源自 4 I 13/2 → 4 I 15/2 跃迁的1540 nm发光强度先上升后下降。发光材料的猝灭浓度可根据发光强度或荧光寿命这两个指标的变化确定，在稀土掺杂激光玻璃体系中，因为随着稀土掺杂浓度的增加，无辐射跃迁几率增加，荧光寿命大多单调降低，因此通常以发光强度开始下降时的浓度作为猝灭浓度。据此可知，GGB玻璃发生明显浓度猝灭现象，实验猝灭浓度为3 wt.%。 不同掺杂浓度Er 3+ ： 4 I 13/2 能级的实测寿命（ τ m ）均呈单指数衰减，如 图1 图1 (c)所示，结果列于 表1 表1 。 τ m 随Er 2 O 3 浓度增加而逐渐降低，归因于无辐射跃迁几率不断增加。 图1 图1 (d)是Er 2 O 3 掺杂3 wt.%的红外透过光谱，根据红外透过光谱2600 nm和3000 nm处透过率 [ 24 24 ] 计算该玻璃体系OH吸收系数约为4-5 cm -1 。以上实验结果表明，Er 3+ 掺杂GGB玻璃发生明显浓度猝灭效应，基于Er 3+ 三能级跃迁特点及荧光寿命呈单指数衰减特性可利用公式（ 4 4 ）-（ 6 6 ）计算猝灭浓度。 4.2　基于I-H模型拟合荧光寿命曲线计算猝灭浓度 2 图2 图2 (a)-(e)是基于I-H模型拟合不同浓度Er 2 O 3 掺杂GGB玻璃荧光寿命曲线及拟合参数 γ 。拟合公式为 y =exp(- t / τ 0 - γt 1/2 )，固定 τ 0 为1 wt.%时的计算值（6.7 ms）， 忽略了不同浓度 τ 0 的细微差异。由图可知随着Er 2 O 3 浓度增加，拟合参数 γ 和 γ 0 增大，拟合度减小。这表明稀土离子浓度增加，Er 3+ 离子间的偶极-偶极相互作用增强，然而在高浓度时偶极相互作用不是导致荧光强度衰减的主要原因。Er 3+ ： 4 I 13/2 激发态粒子布居数与掺杂浓度之间的关系，如 图2 图2 (f)所示。所有曲线都呈现出随Er 2 O 3 浓度增加，激发态布居数先增加后降低的趋势。4 wt.%与5 wt.%两条曲线由于参数 γ 0 相近，两条曲线重合。曲线极大值处对应猝灭浓度 X q 随Er 2 O 3 浓度（1-5 wt.%）增大而减小，所有 X q 在6.4-12.6 wt.%范围内与实验猝灭浓度3 wt.%相差较大。尽管泵浦功率增大，锗酸盐玻璃中Er 3+ 有较强的上转换发光，会影响 4 I 13/2 能级上粒子数布居，但总体上1.5 μm发光强度会随泵浦功率增大而变强， 4 I 13/2 能级粒子数布居增多。然而，本文的激发态布居数与掺杂浓度的关系基于荧光寿命曲线得到，而实际上泵浦功率对荧光寿命的测量以及猝灭浓度值影响不明显。 4.3　基于BN模型拟合荧光寿命曲线计算猝灭浓度 2 图3 图3 (a)-(e)是BN模型拟合不同Er 2 O 3 浓度掺杂GGB玻璃荧光寿命曲线及拟合参数 γ 和 ω ，拟合公式为 y =exp(- t /0.0067- γt 1/2 - ωt )。对比 图2 图2 和 图3 图3 可知，相比于I-H模型，BN模型拟合度更好。不同Er 2 O 3 浓度的拟合参数 γ 都为0，而拟合参数 ω 随浓度增大而增大，这表明离子间的迁移作用大于离子间偶极-偶极相互作用，且迁移作用随Er 2 O 3 浓度增加而增强。Er 3+ ： 4 I 13/2 激发态粒子布居数与掺杂浓度之间的关系见 图3 图3 (f)，所有曲线仍呈现出随Er 2 O 3 浓度增加，激发态布居数先增加后降低的趋势。随Er 2 O 3 浓度增加， ω 0 减小， X q 增大。这与 X q 与 ω 0 成反比结论一致，见 公式（6） 公式（6） 。 X q 在2.6-3.6 wt.%范围内与实验猝灭浓度3 wt.%相差较小。由于不同Er 2 O 3 浓度计算自发辐射跃迁寿命 τ 0 存在些许差异，探究了自发辐射跃迁寿命 τ 0 及Er 2 O 3 掺杂浓度 x 0 对BN模型拟合参数 γ ， ω 及计算猝灭浓度 X q 的影响。计算了 τ 0 为6.3 ms、6.5 ms和6.7 ms时GGB玻璃的猝灭浓度，结果见 表2 表2 。由表可知，固定 τ 0 时， γ 都等于0， ω 随 x 0 增大而增大， X q 在2.6-3.7 wt.%范围内。当 τ 0 =6.3与6.5 ms时，随 x 0 增大, X q 先减小后增大，这是因为当掺杂浓度在1-2wt.%时， ω 与 x 0 2 的 x 0 /（wt.%） τ 0 /（ms） γ /（s -1 ） ω /（s -1 ） γ 0 /（s -1 ） ω 0 /（s -1 ） X q /（wt.%） 1 5 6.3 0 12.27 0 1.23×10 5 3.6 2 0 66.39 0 1.66×10 5 3.1 3 0 126.20 0 1.40×10 5 3.4 4 0 212.42 0 1.33×10 5 3.5 5 0 289.41 0 1.16×10 5 3.7 1 5 6.5 0 17.15 0 1.72×10 5 3 2 0 71.27 0 1.78×10 5 2.9 3 0 131.08 0 1.46×10 5 3.2 4 0 217.30 0 1.36×10 5 3.4 5 0 294.30 0 1.18×10 5 3.6 1 5 6.7 0 21.85 0 2.19×10 5 2.6 2 0 77.59 0 1.94×10 5 2.8 3 0 136.10 0 1.51×10 5 3.1 4 0 221.89 0 1.39×10 5 3.3 5 0 298.97 0 1.20×10 5 3.5 斜率分别大于 k 1 和 k 2 ，而在2-5wt.%时斜率小于 k 1 和 k 2 ，见 图4 图4 (a)和 图4 图4 (b)。当 τ 0 =6.7ms时，随 x 0 增大， X q 增大， ω 与 x 0 2 的斜率一直小于 k 3 ，见 图4 图4 (c)。 k 1 , k 2 , k 3 分别是不同 τ 0 时，掺杂浓度为1 wt.%时 ω 与 x 0 2 的斜率。 X q 随 x 0 增大而增大的主要原因在于 x 0 越大， ω 0 越小，而在理论 公式(5) 公式(5) 中假设参数 ω 0 为常数。此外，在 图4 图4 (c)中，随掺杂浓度增大， ω 与 x 0 2 的斜率与 k 3 的偏离度越大。偏离度增大的原因可能是Er 2 O 3 浓度增加，Er 3+ 离子之间的能量迁移到OH基团的几率变大。由于拟合参数 ω 近似等于无辐射跃迁几率 W 时包含了 W OH ，而 W OH 主要与 x 0 成正比，不与 x 0 2 成正比，则当 W OH 增大，易使 ω 与 x 0 2 的关系偏离线性关系。以上结果表明，当 τ 0 一定时， X q 的变化取决于 ω 0 ，即 ω 与 x 0 2 的斜率。而当 x 0 一定时， X q 随 τ 0 增加而减小，这从 表2 表2 和 公式(6) 公式(6) 都可看出。 4.4　利用自发辐射跃迁寿命和实测寿命计算猝灭浓度 2 本文方法（ 公式（7） 公式（7） ）计算的不同 τ 0 ， x 0 时非辐射跃迁几率 W 以及猝灭浓度 X q ，见 表3 表3 。对比 表2 表2 和 表3 表3 可以看出， ω 与 W 两者相近， X q 亦相近。 图5 图5 (a)和 图5 图5 (b)分别对比了基于BN模型和利用自发辐射跃迁寿命与实测寿命计算的 X q 与实验猝灭浓度。由图可知，这两种方法计算 X q 均随掺杂浓度增加而增大，在掺杂浓度2-3wt.%时计算结果最准确，预测误差均小于0.4 wt.%。以上两种理论计算猝灭浓度的方法都做了必要的近似，存在一定误差。如 γ =0即忽略了离子间的偶极-偶极相互作用以及参数 W、ω 实际上包含了多声子弛豫和OH引起的非辐射跃迁，理论上忽略了 W OH 和 W MP 的影响。根据文献 [ 1 1 ] 可知锗酸盐玻璃的最大声子能量约为900 cm -1 ，计算OH吸收系数约为4-5 cm -1 ，这两部分将引起一定误差。此外，利用BN模型计算猝灭浓度的其他误差来源还包括：J-O理论计算 τ 0 的计算误差，如不同浓度下的计算值存在一定差异；名义掺杂浓度和实际浓度之间的差异；由于荧光寿命曲线精度有限，拟合参数过程中的拟合误差. 相较而言，直接利用自发辐射跃迁寿命和实测寿命计算猝灭浓度的方法不需要拟合荧光寿命曲线。 x 0 /（wt.%） τ 0 /（ms） τ m /（ms） W /（s -1 ） X q /（wt.%） 1 5 6.3 5.68 17.33 3.0 2 4.2 79.37 2.8 3 3.26 148.02 3.1 4 2.58 228.87 3.3 5 2.11 315.20 3.5 1 5 6.5 5.68 22.21 2.6 2 4.2 84.25 2.7 3 3.26 152.90 3.0 4 2.58 233.75 3.2 5 2.11 320.09 3.5 1 5 6.7 5.68 26.80 2.4 2 4.2 88.84 2.6 3 3.26 157.49 2.9 4 2.58 238.34 3.2 5 2.11 324.68 3.4 4.5　不同方法计算理论猝灭浓度简易性与准确性 2 除本文的方法外，根据唯象模型 [ 7 7 ] 和有限扩散模型 [ 8 8 ] 分别可以计算Er 3+ 掺杂GGB玻璃的理论猝灭浓度。本文通过对比不同方法得到的Er 2 O 3 掺杂GGB玻璃的猝灭浓度理论值与实验值，进一步对不同方法的简易性与准确性进行了探讨： （1）根据唯象模型计算的猝灭浓度理论值为2.67 wt.%，绝对误差较小，仅为0.33 wt.%。表明该方法可以精确计算猝灭浓度，但该方法基于多个实验值拟合发光强度，原理同实验，没有减少实验量，且影响发光强度的因素较多，易使不同浓度样品的发光强度值偏离唯象模型，难以拟合参数。 （2）根据有限扩散模型计算的猝灭浓度理论值为3.41 wt.%，绝对误差较小，为0.41 wt.%。表明该方法也可以精确计算猝灭浓度，但拟合所需参数也基于多个实验值，且存在寿命变化趋势拟合度差等问题。 （3）基于I-H和BN模型拟合荧光寿命曲线得到拟合参数并建立激发态布居数与浓度的关系，可利用一个或两个浓度样品的光谱数据（吸收光谱、荧光寿命）计算猝灭浓度，减少实验量。其中I-H模型预测猝灭浓度的绝对误差范围为3.4-9.6 wt.%，误差过大，不可以计算；BN模型绝对误差范围＜0.7 wt.%，可以准确计算。对比两种模型可知，I-H模型没有参数 ω ，BN模型没有参数 γ （ γ ＝0时），这表明相比于 γ ， ω 对计算猝灭浓度准确性具有重要影响，BN模型也因此计算准确高。利用不同浓度下的光谱数据根据BN模型可以计算猝灭浓度，为获得较小的预测误差，需选择相对较低浓度下（2-3wt.%）的光谱数据。 （4）利用自发辐射跃迁寿命和实测寿命计算预测猝灭浓度的方法，也只需要一个或两个浓度样品的光谱数据，且采用掺杂浓度2-3wt.%的样品参数，误差小于0.4 wt.%，可精确预测。相比实验试错法，该方法需要样品数少，预测误差小；相比基于BN模型拟合荧光寿命曲线的方法，该方法不需要拟合荧光寿命曲线，更简单快捷。 对于该方法的适用性和准确性，仍需进行更多探究。我们利用同样的方法对Yb 3+ 、Er 3+ 和Nd 3+ 掺杂磷酸盐玻璃也做了初步研究。结果发现，磷酸盐玻璃体系由于OH - 含量较高，计算得到猝灭浓度误差较大，如果采用低浓度下的实测寿命代替自发辐射跃迁寿命，计算较准确。因此，该方法对不同的玻璃体系和稀土离子具有一定的适用性，通过适当优化后其适用范围更广。 5　结　　论 1 本文提出了直接利用Er 3+ ： 4 I 13/2 → 4 I 15/2 跃迁的自发辐射跃迁寿命和实测寿命预测Er 3+ 掺杂锗酸盐激光玻璃猝灭浓度的方法。研究表明，利用低浓度下Er 2 O 3 掺杂样品参数，理论预测与实际值的最大绝对误差小于0.4 wt.%。与以往利用唯象模型和有限扩散模型拟合不同浓度下的荧光寿命和BN模型拟合荧光寿命曲线的方法相比较，本方法仅需利用一个或两个浓度样品的参数，极大地减少了实验工作量，同时保持了较高的预测准确性。本工作对新型高增益激光玻璃与光纤研究具有指导意义。 参考文献 WANG W C ， ZHOU B ， XU S H ， et al . Recent advances in soft optical glass fiber and fiber lasers ［J］. Prog. Mater. Sci. ， 2019 ， 101 ： 90 - 171 . doi: 10.1016/j.pmatsci.2018.11.003 http://dx.doi.org/10.1016/j.pmatsci.2018.11.003 Recent advances in soft optical glass fiber and fiber lasers ZERVAS M N ， CODEMARD C A . High power fiber lasers： areview ［J］. IEEE J. Sel. Top. Quant. ， 2014 ， 20 ： 219 - 241 . doi: 10.1109/jstqe.2014.2321279 http://dx.doi.org/10.1109/jstqe.2014.2321279 High power fiber lasers： areview 张炳涛 ， 陈月娥 ， 赵兹罡 ， 等 . 有源光纤的进展与应用 ［J］. 应用物理 ， 2018 ， 8 ： 256 - 268 . doi: 10.12677/app.2018.85032 http://dx.doi.org/10.12677/app.2018.85032 有源光纤的进展与应用 ZHANG B T ， CHEN Y E ， ZHAO Z G ， et al . Progress and applications of active optical Fibers ［J］. Appl. Phys. ， 2018 ， 8 ： 256 - 268 . （in Chinese） . doi: 10.12677/app.2018.85032 http://dx.doi.org/10.12677/app.2018.85032 Progress and applications of active optical Fibers SNOEKS E ， KIK P G ， POLMAN A . Concentration quenching in erbium implanted alkali silicate glasses ［J］. Opt. Mater. ， 1996 ， 5 （ 3 ）： 159 - 167 . doi: 10.1016/0925-3467(95)00063-1 http://dx.doi.org/10.1016/0925-3467(95)00063-1 Concentration quenching in erbium implanted alkali silicate glasses MYSLINSKI P ， NGUYEN D ， CHROSTOWSKI J . Effects of concentration on the performance of erbium-doped fiber amplifiers ［J］. J.Lightwave Technol. ， 1997 ， 15 ： 112 - 120 . doi: 10.1109/50.552118 http://dx.doi.org/10.1109/50.552118 Effects of concentration on the performance of erbium-doped fiber amplifiers HOUDE-WALTER S N ， PETERS P M ， STEBBINS J F ， et al . Hydroxyl-contents and hydroxyl-related concentration quenching in erbium-doped aluminophosphate， aluminosilicate and fluorosilicate glasses ［J］. J. Non.Cryst.Solids ， 2001 ， 286 ： 118 - 131 . doi: 10.1016/s0022-3093(00)00445-2 http://dx.doi.org/10.1016/s0022-3093(00)00445-2 Hydroxyl-contents and hydroxyl-related concentration quenching in erbium-doped aluminophosphate， aluminosilicate and fluorosilicate glasses UITERT V L G . Characterization of energy transfer interactions between rare earth ions ［J］. J. Electrochem. Soc. ， 1967 ， 114 ： 1048 . doi: 10.1149/1.2424184 http://dx.doi.org/10.1149/1.2424184 Characterization of energy transfer interactions between rare earth ions AUZEL F ， BONFIGLI F ， GAGLIARI S ， et al . The interplay of self-trapping and self-quenching for resonant transitions in solids： role of a cavity ［J］. J. Lumin. ， 2001 ， 94 ： 293 - 297 . doi: 10.1016/s0022-2313(01)00308-8 http://dx.doi.org/10.1016/s0022-2313(01)00308-8 The interplay of self-trapping and self-quenching for resonant transitions in solids： role of a cavity AUZEL F ， BALDACCHINI G ， LAVERSENNE L ， et al . Radiation trapping and self-quenching analysis in Yb 3+ ， Er 3+ ， and Ho 3+ doped Y 2 O 3 ［J］. Opt. Mater. ， 2003 ， 24 ： 103 - 109 . doi: 10.1016/s0925-3467(03)00112-5 http://dx.doi.org/10.1016/s0925-3467(03)00112-5 Radiation trapping and self-quenching analysis in Yb3+， Er3+， and Ho3+ doped Y2O3 PUGLIESE D ， BOETTI N G ， LOUSTEAU J ， et al . Concentration quenching in an Er-doped phosphate glass for compact optical lasers and amplifiers ［J］. J. Alloys Compd. ， 2016 ， 657 ： 678 - 683 . doi: 10.1016/j.jallcom.2015.10.126 http://dx.doi.org/10.1016/j.jallcom.2015.10.126 Concentration quenching in an Er-doped phosphate glass for compact optical lasers and amplifiers GEBAVI H ， MILANESE D ， LIAO G H ， et al . Spectroscopic investigation and optical characterization of novel highly thulium doped tellurite glasses ［J］. J. Non.Cryst.Solids ， 2009 ， 355 ： 548 - 555 . doi: 10.1016/j.jnoncrysol.2009.01.016 http://dx.doi.org/10.1016/j.jnoncrysol.2009.01.016 Spectroscopic investigation and optical characterization of novel highly thulium doped tellurite glasses RAVAGL A ， BOETTI N G ， CRUZ F A G ， et al . Structural and spectral characterisation of Er 3+ and Nd 3+ doped Ga-La-S-Se glasses ［J］. Rsc. Adv. ， 2018 ， 8 ： 27556 - 27564 . Structural and spectral characterisation of Er3+ and Nd3+ doped Ga-La-S-Se glasses 黄莉蕾 ， 陈继勤 ， 赵渭忠 ， 等 . RE：YAG晶体中掺杂离子（Ho，Nd，Tm）的最佳浓度 ［J］. 量子电子学 ， 1995 ， 12 ： 227 - 231 . RE：YAG晶体中掺杂离子（Ho，Nd，Tm）的最佳浓度 HUANG L L ， CHEN J Q ， ZHAO W Z ， et al . The optimal concentration of doped rare-earth ions Nd， Ho， Tm in RE：YAG crystal ［J］. China J. Quantum. Electron. ， 1995 ， 12 ： 227 - 231 . （in Chinese） The optimal concentration of doped rare-earth ions Nd， Ho， Tm in RE：YAG crystal 黄莉蕾 ， 陈晓竹 . 激光晶体YAG中Er 3+ 离子最佳掺杂浓度 ［J］. 中国计量学院学报 ， 1996 ， 2 ： 15 - 20 . 激光晶体YAG中Er3+离子最佳掺杂浓度 HUANG L L ， CHEN X Z . Optimal concentrations doped Er 3+ ions in YAG laser crystal ［J］. Chin. Jiliang Uni. ， 1996 ， 2 ： 15 - 20 . （in Chinese） Optimal concentrations doped Er3+ions in YAG laser crystal 胡晓 ， 洪方煜 . 四能级和准四能级激活离子的最佳掺杂浓度 ［J］. 物理学报 ， 2002 ， 51 ： 2002 - 2010 . doi: 10.7498/aps.51.2002 http://dx.doi.org/10.7498/aps.51.2002 四能级和准四能级激活离子的最佳掺杂浓度 HU X ， HONG F Y . Optimal concentration of active particles in four-level and quasi-four-level laser systems ［J］. Acta. Phys. Sin. ， 2002 ， 51 ： 2002 - 2010 . （in Chinese） . doi: 10.7498/aps.51.2002 http://dx.doi.org/10.7498/aps.51.2002 Optimal concentration of active particles in four-level and quasi-four-level laser systems XU S J ， CHEN J ， CHEN L Y ， et al . Optimal doping content of red emitting Y 4 GeO 8 ： Bi 3+ ， Eu 3+ phosphor designed by response surface methodology ［J］. Chin. J. Lumin. ， 2022 ， 43 （ 05 ）： 633 - 641 . doi: 10.37188/cjl.20220021 http://dx.doi.org/10.37188/cjl.20220021 Optimal doping content of red emitting Y4GeO8： Bi3+， Eu3+phosphor designed by response surface methodology WEN X ， TANG G W ， YANG Q ， et al . Highly Tm 3+ doped germanate glass and its single mode fiber for 2.0μm laser ［J］. Sci. Rep. ， 2016 ， 6 ： 20344 . doi: 10.1038/srep20344 http://dx.doi.org/10.1038/srep20344 Highly Tm3+ doped germanate glass and its single mode fiber for 2.0μm laser KRUSTINA L ， FRANCESCA C ， JEAN-LUC AA. ， et al . Quench-free enhanced emission in cluster-free Er-doped heavy metal oxide glasses ［J］. Opt. Mater. Express ， 2019 ， 9 （ 3 ）： 1072 - 1084 . doi: 10.1364/ome.9.001072 http://dx.doi.org/10.1364/ome.9.001072 Quench-free enhanced emission in cluster-free Er-doped heavy metal oxide glasses 房喻 ， 王辉 . 荧光寿命测定的现代方法与应用 ［J］. 化学通报 ， 2001 ， 10 ： 631 - 636 . doi: 10.3969/j.issn.0441-3776.2001.10.006 http://dx.doi.org/10.3969/j.issn.0441-3776.2001.10.006 荧光寿命测定的现代方法与应用 FANG Y ， WANG H . Modern methods and applications for fluorescence lifetime determination ［J］. Chem. Bull. ， 2001 ， 10 ： 631 - 636 . （in Chinese） . doi: 10.3969/j.issn.0441-3776.2001.10.006 http://dx.doi.org/10.3969/j.issn.0441-3776.2001.10.006 Modern methods and applications for fluorescence lifetime determination CAMPBELL J H ， SURATWALA T I . Nd-doped phosphate glasses for high-energy/high-peak-power lasers ［J］. J. Non.Cryst.Solids ， 2000 ， 263 ： 318 - 341 . doi: 10.1016/s0022-3093(99)00645-6 http://dx.doi.org/10.1016/s0022-3093(99)00645-6 Nd-doped phosphate glasses for high-energy/high-peak-power lasers 刘永皓 ， 徐权 . 不同Er 3+ 离子掺杂浓度下硼铋酸盐玻璃荧光浓度猝灭研究 ［J］. 齐齐哈尔大学学报（自然科学版） ， 2010 ， 26 ： 41 - 45 . 不同Er3+离子掺杂浓度下硼铋酸盐玻璃荧光浓度猝灭研究 LIU Y H ， XU Q . The effect of Er 3+ concentrations on the fluorescence concentration quenching in bismuth borate glasses ［J］. Chin. J. Qiqihar Uni. ， 2010 ， 26 ： 41 - 45 . （in Chinese） The effect of Er3+ concentrations on the fluorescence concentration quenching in bismuth borate glasses RISEBERG L A ， MOOS H W . Multiphonon orbit-lattice relaxation of excited states of rare-earth ions in crystals ［J］. Phys. Rev. ， 1968 ， 174 ： 429 - 438 . doi: 10.1103/physrev.174.429 http://dx.doi.org/10.1103/physrev.174.429 Multiphonon orbit-lattice relaxation of excited states of rare-earth ions in crystals INOKUTI M ， HIRAYAMA F . Influence of energy transfer by exchange mechanism on donor luminescence ［J］. J. Chem. Phys. ， 1965 ， 43 ： 1978 - 1989 . doi: 10.1063/1.1697063 http://dx.doi.org/10.1063/1.1697063 Influence of energy transfer by exchange mechanism on donor luminescence BURSHTEIN A I . Hopping mechanism of energy transfer ［J］. J. Exp. Theor. Phys. ， 1972 ， 35 （ 5 ）： 882 - 885 . Hopping mechanism of energy transfer CARNALL W T ， FIELDS P R ， WYBOURNE B G . Spectral Intensities of the trivalent lanthanides and actinides in solution. I. Pr 3+ ， Nd 3+ ， Er 3+ ， Tm 3+ ， and Yb 3+ ［J］. J. Chem. Phys. ， 1965 ， 42 （ 11 ）： 3797 - 3806 . doi: 10.1063/1.1695840 http://dx.doi.org/10.1063/1.1695840 Spectral Intensities of the trivalent lanthanides and actinides in solution. I. Pr3+， Nd3+， Er3+， Tm3+， and Yb3+ KESAVULU C R ， SREEDHAR V B ， JAYASANKAR C K ， et al . Structural， thermal and spectroscopic properties of highly Er 3+ -dopednovel oxyfluoride glasses for photonic application ［J］. Mater. Res. Bull. ， 2014 ， 51 ： 336 - 344 . doi: 10.1016/j.materresbull.2013.12.023 http://dx.doi.org/10.1016/j.materresbull.2013.12.023 Structural， thermal and spectroscopic properties of highly Er3+-dopednovel oxyfluoride glasses for photonic application ZHAO Y ， SHI D . Effect of alkali metal oxides R 2 O （R=Na， K） on 1.53 μm luminescence of Er 3+ -doped Ga 2 O 3 -GeO 2 glasses for optical amplification ［J］. J. Rare Earths ， 2013 ， 31 ： 857 – 863 . doi: 10.1016/s1002-0721(12)60370-6 http://dx.doi.org/10.1016/s1002-0721(12)60370-6 Effect of alkali metal oxides R2O （R=Na， K） on 1.53 μm luminescence of Er3+-doped Ga2O3-GeO2glasses for optical amplification 10.37188/CJL.20230014 O482.31 A 1000-7032（XXXX）XX-0001-10 万杰,欧阳莎,姬瑶等.Er 3+ 掺杂锗酸盐激光玻璃猝灭浓度的计算与预测[J].发光学报,DOI:10.37188/CJL.20230014 万杰,欧阳莎,姬瑶等.Er3+掺杂锗酸盐激光玻璃猝灭浓度的计算与预测[J].发光学报,DOI:10.37188/CJL.20230014 WAN Jie,OUYANG Sha,JI Yao,et al.Calculation and prediction of quenching concentration of Er 3+ -doped germanate glass[J].Chinese Journal of Luminescence,DOI:10.37188/CJL.20230014 WAN Jie,OUYANG Sha,JI Yao,et al.Calculation and prediction of quenching concentration of Er3+-doped germanate glass[J].Chinese Journal of Luminescence,DOI:10.37188/CJL.20230014 国家自然科学基金（52172003）资助项目 Supported by National Natural Science Foundation ofChina （52172003） 2023-03-09 2023-03-09 发光学报