-
高功率掺镱光纤激光器凭借其结构紧凑、时序逻辑简单、热管理方便、功率拓展性强等优势, 在工业加工、材料处理、生物医疗、航空航天等领域得到了广泛应用[1–3]. 然而随着输出功率的不断提高, 光纤中的非线性效应, 如受激拉曼散射(stimulated Raman scattering, SRS), 受激布里渊散射(stimulated Brillouin scattering, SBS), 以及横向模式不稳定效应(transverse mode instability, TMI)等成为近衍射极限光纤激光器功率提升的主要限制因素[4–6]. 为获得更高功率输出, 通过光谱扩展实现光谱合成是提升光纤激光器输出功率的重要方法之一[7,8]. 对于掺镱光纤来说, 工作在1060—1080 nm的常规波长光纤激光器已经实现了数十千瓦近单模激光输出[9,10], 而波长小于1060 nm的短波长光纤激光器由于重吸收作用受到严重的放大自发辐射效应(amplified spontaneous emission, ASE)限制, 功率提升相对缓慢.
近年来, 国内外多家研究机构对短波长光纤激光器开展了一定的研究. 2011年, Schmidt等[11]在激光波长为1030 nm处实现了最高输出功率为697 W的激光输出, 有效地抑制了放大自发辐射. 2016年, Roman等[12]报道了一种基于相位调制的单频激光器, 在1030—1070 nm的工作波长内实现平均输出功率均超过1.5 kW, 在1032 nm波段输出功率达1.8 kW, 仅受限于泵浦功率, 而在波长1044 nm和1061 nm处输出功率仅为1.6 kW和1.48 kW. 同年, 孙殷宏等[13]基于主振荡功率放大(master oscillator power amplifier, MOPA)结构搭建了中心波长为1030 nm的光纤激光放大器, 通过对结构的优化实现了最高功率为1.01 kW的激光输出. 2020年, Chu等[14]基于光谱操控、放大级系统参数优化和模式控制等关键技术, 搭建了1030 nm窄线宽光纤放大器, 最高输出功率达到3 kW, 信噪比达到了37 dB. 2021年, Xu等[15]基于MOPA结构实现2.4 kW的1045 nm的近单模窄线宽光纤激光功率输出, 拉曼信噪比约33 dB, 结构通过优化增益光纤长度来抑制SRS和TMI. 2022年, Zheng等[16]提出了基于泵浦共享结构的1050 nm窄线宽光纤激光器, 即种子和主放大级之间没有包层光滤除装置, 振荡器和放大器之间可以共用泵浦光, 利用此结构可以有效抑制ASE, 在最大输出功率3.1 kW处, 信噪比达到45.5 dB, 光束质量因子
$M^2 \sim 1.33 $ .上述报道的短波长光纤激光器都是采用单端输出形式, 由于短波长光纤激光器具有较宽的吸收截面, 更容易受到ASE和SRS效应影响[17,18], 一定程度上限制了短波长光纤激光器的发展. 而双端输出光纤激光器功率分布更均衡, 减小了腔内的功率密度, 有利于提升非线性阈值. 同时, 双端输出的光纤激光器仅需一个谐振腔就能实现两路激光输出, 与单端输出的光纤激光器相比, 减少了器件数量, 简化了冷却装置和控制系统, 在降低系统成本的同时压缩了系统体积和重量. 根据其输出特性可以广泛应用于工业加工与制造, 材料处理和3D打印等领域, 可实现高功率并行加工, 提高效率. 在对光束质量要求不高时, 采用双端输出光纤激光器通过功率合束等技术, 能够实现更高的功率输出[9], 具有更加广阔的工业应用前景. 2004年, Liu等[19]基于纤芯为20 μm的双包层掺镱光纤(YDF), 实现了最高输出功率为810 W的空间结构光纤激光振荡器, 掺镱光纤的两端切平角, 分别提供3.5%的菲涅尔反射, 形成振荡腔, 由于没有高反和低反之分, 实现了掺镱光纤两端激光输出. 2018年, 王小林等[20]提出了一种双端输出的线性腔全光纤激光振荡器结构, 该结构利用低反射率光栅代替高反射率光栅, 从而实现全光纤激光器的双端功率输出. 2022年, Zhong等[21]基于该线性腔双端输出光纤激光振荡器结构, 首次实现中心波长为1080 nm的2×2 kW双端激光输出. 2023年, Liu等[22]基于纤芯/包层直径为30/400 μm的双包层掺镱光纤实现了中心波长为1070 nm的2×4 kW双端激光输出, 最高输出功率时两端输出激光的光束质量
$M^2 $ 分别约为 2.13和2.36, 没有观察到SRS和TMI现象. 然而当前已有研究中, 双端输出光纤激光器都是集中在常规波段, 短波长双端输出相关研究尚未开展, 双端输出体制下的短波长光纤激光器输出特性和非线性效应演化仍然值得探索.本文对中心波长为1050 nm的双端输出光纤激光振荡器输出特性进行研究, 首先基于光纤速率方程构建了1050 nm双端输出光纤激光器理论模型, 理论分析了1050 nm双端输出光纤振荡器中增益光纤长度与输出特性和非线性效应之间的影响规律. 然后, 基于纤芯/内包层直径为20/400 μm 的双包层掺镱光纤搭建了1050 nm双端输出光纤激光振荡器实验装置, 研究了不同泵浦方式下1050 nm双端输出光纤激光器的输出特性. 实验上首次实现了总输出功率约为4470 W的双端1050 nm近单模激光输出, 两端光束质量因子
$M^2$ 约为1.27和1.31. 最后优化谐振腔增益光纤长度, 进一步抑制了激光器的ASE和SRS效应, 在最大输出功率4.5 kW时, A端和B端的SRS抑制比分别提升约6.6 dB和8.1 dB. -
为研究1050 nm双端输出光纤激光振荡器的输出特性, 基于稳态速率方程构建双端输出光纤激光振荡器的仿真模型, 连续输出的光纤激光器在稳态条件时, 考虑放大自发辐射和受激拉曼散射的影响, 在增益光纤中, 上能级粒子数、泵浦光和信号光功率的速率方程可以表示如下, 方程中主要参数的物理含义在表1列出.
其中
$ {N_0} = {N_1} + {N_2} $ ,图1为双端输出光纤振荡器的简化结构图, 该结构的谐振腔取一段长度为L 的增益光纤, 定义最左端 z = 0, 最右端 z = L, 从左向右传输的泵浦和信号光为正向传输, 从右向左传输的泵浦和信号光为反向传输. 在谐振腔内传输正向信号光
$ P_n^{{\text{s}} + }\left( {\lambda _n^{\text{s}}, 0} \right) $ 和反向信号光$ P_n^{{\text{s}} - }\left( {\lambda _n^{\text{s}}, L} \right) $ 满足速率方程(3). 其中光纤光栅(FBG)反射率$ {R_{\text{a}}}\left( {\lambda _n^{\text{s}}} \right) $ 和$ {R_{\text{b}}}\left( {\lambda _n^{\text{s}}} \right) $ 是需要设置的参数, 如将两光栅的反射率均设为小于1的值, 则谐振腔两端都有激光透射, 此时可实现双端激光输出, 结合边界条件(4)式—(7)式选择合适的数值仿真算法可进行求解. -
本节主要研究了1050 nm双端输出光纤激光振荡器中增益光纤长度对输出功率、效率和SRS强度的影响. 表2列出了数值仿真中使用的各器件参数.
仿真结果如图2所示, 在两端泵浦功率均为3000 W, 最高总泵浦功率为6000 W时, 不同增益光纤长度对输出功率及效率的影响如图2(a)所示, 由于增益光纤长度越短, 泵浦吸收越少, 两端输 出功率、总输出功率以及光光转换效率均随光纤 长度减小而降低. 整体来说, 两端输出功率差异较小, 光光转换效率均在85%以上. 由图2(a)可知, 增益光纤长度为14.6 m和16.6 m时激光器输出功率和效率相近, 考虑非线性效应等影响, 光纤长度选定为14.6 m进行研究. 在增益光纤长度为14.6 m时采用双向泵浦方式下A, B两端输出光谱如图2(b)所示, 两光谱几乎重合. 不同增益光纤长度对受激拉曼散射强度的影响如图2(c)所示, 受激拉曼散射强度随增益光纤长度减小而降低, 在截短增益光纤2 m, 4 m和6 m的情况下拉曼抑制比分别提升了约6.3 dB, 12.5 dB和18.9 dB. 在相同条件下, 通过截短增益光纤长度对SRS有一定抑制效果.
-
图3为双向泵浦的1050 nm双端输出光纤激光振荡器实验结构示意图, 该结构由一对光纤布拉格光栅和一段增益光纤组成谐振腔. 实验中采用的增益光纤为20/400 μm的双包层掺镱光纤, 光纤长度约为14.6 m, 该光纤的泵浦吸收系数为0.44 dB/m@915 nm, 数值孔径约为0.061, 为了实现有效的模式控制, 将增益光纤弯曲盘绕在跑道型水冷板上, 光纤盘绕的最小弯曲直径约为10.8 cm, 最大弯曲直径约为13 cm. 泵浦源采用中心波长为976 nm的稳波长半导体激光器(LD), 每组LD可提供最大泵浦功率约为600 W. 泵浦光通过前向和后向(6+1)×1泵浦信号合束器(PSC)耦合进谐振腔, 合束器的泵浦注入纤和信号输出纤的纤芯/包层直径分别为220/242 μm和20/250 μm, 信号注入纤的纤芯/包层直径为20/400 μm. 前/后向合束器的5根泵浦尾纤分别与LD连接, 剩余1根泵浦纤切平角由600 W功率计监测回光功率. 谐振腔腔体由一对低反光纤光栅构成, 低反光纤光栅分别与增益光纤的两端连接, 光纤光栅的纤芯/包层直径为20/400 μm, 其中心波长约为1050 nm, 反射率约为10%, 3 dB带宽约为2 nm. 谐振腔内产生的激光经过两端低反射率光栅的透射后输出, 经过一段自制的长度约为15 cm的包层光滤除器(CLS)滤除残余泵浦光后, 由纤芯/包层直径为25/400 μm的光纤端帽(QBH)进行激光输出. 经QBH输出的激光扩束后进入测试系统, 测量其输出功率, 输出光谱, 光束质量和时频特性等. 实验测试系统由功率计、光谱仪、示波器和光束质量因子测试仪组成, 实验过程中水冷温度均设置为20 ℃.
-
在1050 nm双端输出光纤振荡器实验中, 首先采用单向泵浦方式对激光器的输出特性进行研究, 在A端泵浦方式下, 当注入的泵浦功率为2396 W时, 两端总输出功率达到1943 W, 其中A端输出功率为515 W, B端输出功率为1428 W, 此时光光转换效率为81.1%, 如图4(a)所示. 实验中单向泵浦注入方式下, 两端输出功率差异主要与器件差异和熔接损耗相关, 其中两个低反光栅的反射率虽然相近, 但B端低反光栅测得的反射谱较窄, 使谐振腔内更多的信号激光从B端透射出去, 从而导致B端的输出功率更高, 两端的熔接损耗也是造成功率差异的原因. 当输出功率接近1.94 kW时, 随着进一步增大A端的泵浦功率, A, B两端输出功率出现明显的滞涨, 对应的光光转换效率出现下降. 在此功率水平下, 测得A端的时域信号和对应的傅里叶变化频谱如图4(b)所示, 可以看出, 此时时域信号出现轻微波动, 同时傅里叶频谱图上在0.785 kHz出现了频率尖峰. 结合输出功率的滞涨和傅里叶频域信号的特征峰, 可以判断此时该激光器出现了TMI效应, 且TMI阈值约为1.94 kW.
考虑在双端输出光纤振荡器中, 相同泵浦波长下双向泵浦方式与单向泵浦相比可以显著提升TMI阈值[23], 因此采用双向泵浦方式来进一步提升输出功率, 双向泵浦实验结果如图5所示. 此时两端输出功率和总输出功率随着泵浦功率的增加呈线性增长, 且B端输出功率仍然高于A端, 这与单向泵浦时候的情况一致. 随着泵浦功率增大, 激光器光光效率整体保持在80%—85%之间. 在双端最高总泵浦功率为5262 W, A端的输出功率为1419 W, B端的输出功率为3051 W, 此时双端输出总功率为4470 W, 如图5(a)所示, 光光转换效率为84.9%. 在最高输出功率测得的输出激光时域信号和对应的傅里叶频域如图5(b)所示, 可以看到A, B两端的时域信号保持稳定, 对应的傅里叶频谱没有观察到特征尖峰, 表明此时该激光器没有发生TMI效应. 在最高输出功率下, 采集A, B两端输出光谱图如图5(c)所示, 可以看到信号激光中心波长在1050 nm, 由于镱离子对1050 nm短波长激光有一定的重吸收作用, 在光谱上出现了较强的ASE效应, 且光谱存在严重展宽. 另外从输出光谱图上可以看到在1100 nm附近出现明显的SRS效应, 在最大输出功率时, A端输出光谱的拉曼抑制比约为29.1 dB, B端拉曼抑制比约为24.6 dB. 在输出功率为1458 W和2733 W时测得的光束质量和光斑轮廓如图5(d)所示, 光束质量M2因子
$M_{\text{A}}^{2}$ 约为1.27,$M_{\text{B}}^{2}$ 约为1.31, 束腰处的光斑轮廓近似高斯形态, 表明该激光器具有良好的光束质量, 光束质量差异主要由两端熔接差异引起.为了进一步抑制非线性效应, 根据仿真结果综合考虑增益光纤长度对输出功率、效率和SRS效应的影响, 对实验结构进行优化. 由于输出功率及效率随光纤长度减小而降低, 为实现更高功率输出并达到良好的非线性效应抑制效果, 选择将YDF截短2 m, 光纤长度由14.6 m变为12.6 m, 其余实验条件保持不变. 在同样采用双向泵浦方式, 总泵浦功率为5262 W时, 双端总输出功率为4468 W, 与结构优化前相比两端输出功率及整体光光转换效率差异较小, 如图6(a)所示, 两端光束质量M2因子分别为1.34和1.32, 远场光斑与光束质量良好. 结构优化前后两端光谱对比如图6(b)所示, 由两端光谱可以看出, 与优化增益光纤长度前相比, 激光器SRS得到显著抑制, 由ASE引起的光谱展宽现象也有一定的改善, 在最高输出功率处A端拉曼抑制比提升约6.6 dB, B端拉曼抑制比提升约8.1 dB. 由此可得结论, 缩短增益光纤长度可以有效抑制双端输出短波长光纤振荡器中的拉曼效应. 虽然SRS得到一定抑制, 但由于选取的20/400 μm掺镱光纤有效模场面积较小, 在双端输出最高功率时光谱中仍存在明显的SRS现象, 限制了功率的进一步提升.
-
基于1050 nm双端输出光纤激光振荡器结构, 在总泵浦功率为5262 W时, 实现了双端总功率为4470 W的近单模输出, 两端功率分别为1419 W和3051 W, 光光转换效率为84.9%. 在输出功率为1458 W和2733 W时两端光束质量因子M2分别为1.27和1.31. 在A端泵浦功率为2386 W时出现TMI现象, 在相同泵浦功率下采用双向泵浦方式时域频域稳定, 未产生TMI效应. 优化增益光纤长度后实现了ASE和SRS的有效抑制, 在总泵浦功率为5262 W时, A端和B端拉曼抑制比分别提升约6.6 dB和8.1 dB. 双端输出功率主要受限泵浦功率和SRS效应, 后续通过增加泵浦和优化激光器结构有望实现更高的功率输出.
4.5 kW, 1050 nm双端输出近单模全光纤激光振荡器
4.5 kW, 1050 nm bidirectional output near-single-mode all-fiber laser oscillator
-
摘要: 随着工业应用对激光器性能要求的不断提高, 基于单一谐振腔结构实现两路激光同步输出的双端输出光纤激光器具有广阔的应用前景. 本文首先基于光纤稳态速率方程建立了1050 nm双端输出振荡器理论模型, 仿真分析增益光纤长度与输出功率、效率和受激拉曼散射效应强度间的关系. 实验上搭建中心波长为1050 nm的高功率双端输出全光纤激光振荡器, 详细研究了不同泵浦方式下(单向泵浦、双向泵浦) 1050 nm双端输出光纤激光器的输出特性. 在总泵浦功率为5262 W时, 首次实现了A端输出功率1419 W, B端输出功率3051 W, 总输出功率为4470 W的1050 nm近单模双端激光输出, 激光器光光转换效率达到84.9%, A端和B端测得的光束质量因子$M^2$分别为1.27和1.31. 进一步优化增益光纤长度, 有效抑制了放大自发辐射和受激拉曼散射效应, 最大输出功率下A端和B端的拉曼抑制比分别提升约6.6 dB和8.1 dB. 实验结果为设计和实现高功率高光束质量短波长双端输出光纤激光器提供参考.Abstract: High-power fiber laser oscillators have been widely used in industrial processing, material processing, biomedical and other fields due to their compact structure, simple logic and strong power scalability. With the increasing demand for laser performance in industrial applications, bidirectional output fiber laser based on a single resonator structure has a broad application prospect. In this work, we first establish a theoretical model for a 1050-nm bidirectional output fiber laser oscillator based on the steady-state rate equation, and simulate the relationship between the length of the gain fiber and output power, efficiency, and the intensity of stimulated Raman scattering (SRS). A high-power bidirectional output fiber laser with a central wavelength of 1050 nm is built using an ytterbium-doped fiber with a core/cladding diameter of 20/400 μm. The output characteristics of the 1050-nm bidirectional output fiber laser oscillator under different pump methods (unidirectional pump, bidirectional pump) are experimentally studied in detail. With a total pump power of 5262 W, A-end output power reaches 1419 W and B-end output power reaches 3051 W. Therefore, a total output power of 4470 W with an optical-to-optical conversion efficiency of 84.9% is achieved. The corresponding beam qualities (M2 factor) of both ends are 1.27 and 1.31 when the output powers reach 1458 W and 2733 W, respectively. By further optimizing the length of the gain fiber, the amplified spontaneous emission (ASE) and SRS are effectively suppressed. With a total pump power of 5262 W, the Raman suppression ratios at A-end and B-end are enhanced by ~6.6 dB and ~8.1 dB, respectively. It is expected that higher output power can be achieved by increasing the pump power and optimizing the laser structure in the future.
-
Key words:
- bidirectional output /
- fiber laser oscillator /
- near-single-mode /
- stimulated Raman scattering .
-
-
表 1 速率方程主要参数
Table 1. Main parameters of the rate equation.
物理量 物理意义 物理量 物理意义 R 反射率 ${N_2}$ 激发态粒子数 Z 增益光纤纵向坐标 ${N_1}$ 基态粒子数 $m$ 泵浦光波长序数 ${N_0}$ 掺杂离子浓度 $n$ 信号光波长序数 ${g_{\text{R}}}$ 拉曼增益系数 ${A_{{\text{eff}}}}$ 纤芯有效面积 $\sigma _n^{{\text{es}}}$ 第n个信号光吸收截面 ${\varGamma _{\text{s}}}$ 信号光填充因子 $\sigma _n^{{\text{as}}}$ 第n个信号光发射截面 ${\varGamma _{\text{p}}}$ 泵浦光填充因子 $\alpha _n^{\text{p}}(\lambda _n^{\text{p}})$ 信号光损耗系数 表 2 仿真主要参数
Table 2. Simulation parameter.
主要参数 值 信号光中心波长/nm 1050 泵浦光中心波长/nm 976 总泵浦功率/W 6000 纤芯/内包层/μm 20/400 泵浦吸收系数 0.44 dB/m@915 nm 增益光纤长度/m 16.6/14.6/12.6/10.6/8.6 FBG-A的反射率/% 10 FBG-B的反射率/% 10 FBG-A的半高全宽/nm 2 FBG-B的半高全宽/nm 2 -
[1] Richardson D J, Nilsson J, Clarkson W A 2010 J. Opt. Soc. Am. B 27 B63 doi: 10.1364/JOSAB.27.000B63 [2] Zervas, Michalis N 2014 Int. J. Mod. Phys. B 28 1442009 doi: 10.1142/S0217979214420090 [3] 王小林, 张汉伟, 杨保来, 奚小明, 王鹏, 史尘, 王泽锋, 周朴, 许晓军, 陈金宝 2021 中国激光 48 0401004 doi: 10.3788/CJL202148.0401004 Wang X L, Zhang H W, Yang B L, Xi X M, Wang P, Shi C, Wang Z F, Zhou P, Xu X J, Chen J B 2021 Chin. J. Lasers 48 0401004 doi: 10.3788/CJL202148.0401004 [4] Zhu J J, Zhou P, Ma Y X, Xu X J, Liu Z J 2011 Opt. Express 19 18645 doi: 10.1364/OE.19.018645 [5] Jauregui C, Limpert J, Tünnermann A 2013 Nat. Photonics 7 861 doi: 10.1038/nphoton.2013.273 [6] Zervas M N 2019 Opt. Express 27 19019 doi: 10.1364/OE.27.019019 [7] Augst S J, Ranka J K, Fan T Y, Sanchez A 2007 J. Opt. Soc. Am. B 24 1707 doi: 10.1364/JOSAB.24.001707 [8] 辛国锋, 皮浩洋, 沈力, 瞿荣辉, 蔡海文, 方祖捷, 陈高庭 2010 激光与光电子学进展 47 17 doi: 10.3788/LOP47.101404 Xin G F, Pi H Y, Shen L, Ju R H, Cai H W, Fang Z J, Chen G T 2010 Laser Optoelectron. Prog. 47 17 doi: 10.3788/LOP47.101404 [9] 王小林, 曾令筏, 叶云, 刘佳琪, 吴函烁, 王鹏, 杨保来, 奚小明, 张汉伟, 史尘, 习锋杰, 王泽锋, 周朴, 许晓军, 陈金宝 2024 中国激光 51 223 doi: 10.3788/CJL240948 Wang X L, Zeng L F, Ye Y, Liu J Q, Wu H S, Wang P, Yang B L, Xi X M, Zhang H W, Shi C, Xi F J, Wang Z F, Zhou P, Xu X J, Chen J B 2024 Chin. J. Laser 51 223 doi: 10.3788/CJL240948 [10] Zeng L F, Ding X Y, Liu J Q, Wang X L, Ye Y, Wu H S, Wang P, Xi X M, Zhang H W, Shi C, Xi F J, Xu X J 2024 Micromachines-Basel 15 153 doi: 10.3390/mi15010153 [11] Schmidt O, Wirth C, Rhein S, Rekas M, Kliner A, Schreiber T, Tünnermann R E, Andreas 2011 The European Conference on Lasers and Electro-Optics Munich, Germany, May 22–26, 2011 p1 [12] Roman Y, Nikolai P, Alexander Y, Valentin P G 2016 Proc. SPIE San Francisco, March 9, 2016 p972807 [13] 孙殷宏, 柯伟伟, 冯昱骏, 王岩山, 彭万敬, 马毅, 李腾龙, 王小军, 唐淳, 张凯 2016 中国激光 43 0601003 doi: 10.3788/CJL201643.0601003 Sun Y H, Ke W W, Feng Y J, Wang Y S, Peng W J, Ma Y, Li T L, Wang X J, Tang C, Zhang K 2016 Chin. J. Laser 43 0601003 doi: 10.3788/CJL201643.0601003 [14] Chu Q H, Shu Q, Liu Y 2020 Opt. Lett. 45 6502 doi: 10.1364/OL.405386 [15] Xu Y, Sheng Q, Wang P 2021 Appl. Opt. 60 3740 doi: 10.1364/AO.420708 [16] Zheng Y H, Han Z G, Li Y L 2022 Opt. Express 30 12670 doi: 10.1364/OE.456856 [17] Liu Z J, Ma P F, Tao R M, Wang X L, Zhou P 2015 Ieee J. Quantum Elect. 51 1 doi: 10.1109/JQE.2015.2442760 [18] Silva A, Boller K, Lindsay I D 2011 Opt. Express 19 10511 doi: 10.1364/OE.19.010511 [19] Liu C H, Galvanauskas A, Ehlers B, Doerfel F, Heinemann S, Carter A, Tankala K, Farroni J 2004 Advanced Solid-State Photonics Santa Fe, New Mexico, February 1–4, 2004 p17 [20] 王小林, 叶云, 奚小明, 史尘, 张汉伟, 韩凯, 王泽锋, 许晓军, 周朴, 司磊, 陈金宝 2018 中国专利201821644646.3 Wang X L, Ye Y, Xi X M, Shi C, Zhang H W, Han K, Wang Z F, Xu X J, Zhou P, Si L, Chen J B 2018 CN Patent 201821644646.3 [21] Zhong P L, Wang L, Yang B L 2022 Opt. Lett. 47 2806 doi: 10.1364/OL.458581 [22] Liu J Q, Zeng L F, Wang X L, Shi C, Wu H S, Wang P, Xi X M, Zhang H W, Ning Y, Xi F J 2024 Opt. Laser Technol. 169 110031 doi: 10.1016/j.optlastec.2023.110031 [23] Li F C, Ding X Y, Wang P, Yang B L, Xi X M, Zhang H W, Wang X L, Chen J B 2023 Photonics-Basel 10 912 doi: 10.3390/photonics10080912 -