搜索

文章查询

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

2 μm人眼安全波段太阳光抽运激光器的理论研究

林学彤 杨苏辉 王欣 李卓 张金英

2 μm人眼安全波段太阳光抽运激光器的理论研究

林学彤, 杨苏辉, 王欣, 李卓, 张金英
PDF
HTML
导出引用
导出核心图
  • 太阳光直接抽运激光器在空间光通信、遥感等领域有着重要的潜在应用, 但是一直以来人们对太阳光抽运激光器的研究局限于以掺Nd3+粒子为增益介质的1 μm波段. 通过对现有固体激光工作物质的吸收谱进行分析, 发现掺Tm3+离子在太阳辐射较强的可见光波段具有强的吸收峰, 使2 μm人眼安全波段实现太阳光直接抽运激光输出成为可能. 本文对Tm:YAG和Tm:YAP两种常见晶体的吸收谱与太阳光谱匹配度进行了分析计算, 得出两种材料用于太阳光抽运激光器的阈值抽运功率密度分别为1.14和1.434 kW/cm3. 选择与抽运阈值功率密度低的Tm:YAG晶体作为增益介质, 使用TracePro软件建立太阳光抽运激光器的二级抽运模型, 并对模型进行优化, 得到了锥形腔窗口与菲涅耳透镜的最佳距离、晶体棒的最佳长度以及锥形腔最佳锥度. 本文的工作为实现太阳光直接抽运2 μm激光输出做了理论上的准备.
      通信作者: 杨苏辉, suhuiyang@bit.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 61835001, 61875011)资助的课题
    [1]

    Shi J, Wang H, Qian J, He X 2016 Opt. Commun. 363 21

    [2]

    Zhu X, Lu Z, Wang Y 2015 Laser Part. Beams 33 11

    [3]

    Kavaya M J, Beyon J Y, Koch G J, Petros M, Petzar P, Singh U N, Trieu B C, Yu J 2014 J. Atmos. Ocean. Tech. 31 826

    [4]

    Wagener T J, Demma N, Kemetec J D, Kubo T S 1995 IEEE Aero. El. Sys. Mag. 10 23

    [5]

    Marano M, Galzerano G, Svelto C, Laporta P 2004 IEEE T. Instrum. Meas. 53 571

    [6]

    Koch G J, Beyon J Y, Petzar P, Petros M, Yu J, Trieu B C, Kavaya M J, Singh U N, Modlin E A, Bames B W, Demoz B B 2010 J. Appl. Remote Sens. 4 043512

    [7]

    Koch G J, Beyon J Y, Bames B W, Petros M, Yu J, Amzajerdian F, Kavaya M J, Singh U N 2007 Opt. Eng. 46 16201

    [8]

    Geisthoff U W, Zenk J, Steinhart H, Iro H 2001 HNO 49 194

    [9]

    Ma Q L, Bo Y, Zong N, Pan Y B, Peng Q J, Cui D F, Xu Z Y 2011 Opt. Commun. 284 1645

    [10]

    Young C G 1966 Appl. Opt. 5 993

    [11]

    Thompson G A, Krupkin V, Yogev A 1992 Opt. Eng. 31 2644

    [12]

    Yabe T, Ohkubo T, Uchida S, Yoshida K, Nakatsuka M, Funatsu T, Mabuti A, Oyama A, Nakagawa K, Oishi T, Daito K, Behgol B, Naayama Y, Yoshida M, Motokoshi S, Sato Y, Baasandash C 2007 Appl. Phys. Lett. 90 261120

    [13]

    Saiki T, Funahashi K, Motokoshi S, Imasaki K, Fujioka K, Fujita H, Nakatsuka M, Yamanaka C 2009 Opt. Commun. 282 614

    [14]

    Saiki T, Motokoshi S, Imasaki K, Nakatsuka M, Yamanaka C, Fujioka K, Fujita H 2009 Opt. Commun. 282 936

    [15]

    杨扬 2007 博士毕业论文 (上海: 上海交通大学)

    Yang Y 2007 Ph. D. Dissertation (Shanghai: Shanghai Jiao Tong University) (in Chinese)

    [16]

    O'Hare J M, Donlan V L 1976 Pyhs. Rev. B 14 3732

    [17]

    Beyatli E, Sumpf B, Demirbas U 2019 Appl. Opt. 58 2973

    [18]

    赵彬, 赵长明, 何建伟, 杨苏辉 2007 光学学报 27 1797

    Zhao B, Zhao C M, He Z W, Yang S H 2007 Acta Opt. Sin. 27 1797

    [19]

    方容川 2001 固体光谱学 (合肥: 中国科学技术大学出版社) 第4页

    Fang R C 2001 Solid State Spectroscopy (Vol. 1) (Hefei: Press of University of Science and Technology of China) p4 (in Chinese)

    [20]

    赵立伟 2010 硕士毕业论文 (北京: 北京理工大学)

    Zhao L W 2010 M. S. Dissertation (Beijing: Beijing Institute of Technology) (in Chinese)

    [21]

    徐鹏 2019 博士毕业论文 (北京: 北京理工大学)

    Xu P 2019 Ph. D. Dissertation (Beijing: Beijing Institute of Technology) (in Chinese)

  • 图 1  太阳光谱曲线

    Fig. 1.  Curve of solar spectrum.

    图 2  太阳光谱与Tm:YAP, Tm:YAG吸收谱 (a) Tm:YAG; (b) Tm:YAP

    Fig. 2.  Matching curve of crystals and solar spectrum: (a) Tm:YAG; (b) Tm:YAP.

    图 3  Tm3+离子能级跃迁示意图

    Fig. 3.  Schematic diagram of Tm3+ ion level transition.

    图 4  (a) TracePro软件建立的太阳光抽运激光器二级抽运模型; (b)锥形腔结构图

    Fig. 4.  (a) Two-stage pumping model; (b) structure diagram of conical cavity.

    图 5  锥形腔窗口位置与接受光功率关系图

    Fig. 5.  Curve of the relationship between the position of conical cavity and the received solar power.

    图 6  100 mm晶体棒侧面抽运光分布图

    Fig. 6.  Distribution map of side pump power on 100 mm-length crystal rod.

    图 7  不同晶体长度下晶体棒轴向光辐照度分布

    Fig. 7.  Axial irradiance distribution of different-length crystals.

    图 8  不同锥度下晶体棒轴向光辐照度分布图

    Fig. 8.  Axial irradiance distribution of different-taper crystals.

    表 1  晶体光谱匹配分析结果

    Table 1.  Spectral matching analysis results of crystals.

    Active
    medium
    Doping density/cm3Absorption
    band/nm
    Irradiance in absorption
    band/W·m–2
    Percentage of solar
    radiance/%
    Tm:YAG1.261×1020 (1 at.%)360—41021.9104122.9
    456—48038.78578
    656—72084.80076
    747—81269.35269
    Tm:YAP1.965×1020 (1 at.%)360—39419.5293829.8
    450—49468.23487
    643—726108.88686
    744—83695.81626
    下载: 导出CSV

    表 2  晶体参数

    Table 2.  Crystal parameters.

    Tm:YAGTm:YAP
    Doping density/cm31.26×1020
    (1 at.%)
    1.965×1020
    (1 at.%)
    Upper level lifetime/ms10.54.4
    Boltzmann factor
    in upper level
    0.460.29
    Boltzmann factor
    in lower level
    0.0170.015
    Emission cross section/cm22.5×10–213.81×10–21
    Quantum efficiency1.81.9
    Calculation results of
    absorption coefficient
    curve: $\sum\nolimits_i { {\eta _i}\overline { {a_i} } {\lambda _i} } $/cm
    46.0069×10–751.5739×10–7
    Refractive index1.821.91
    下载: 导出CSV
  • [1]

    Shi J, Wang H, Qian J, He X 2016 Opt. Commun. 363 21

    [2]

    Zhu X, Lu Z, Wang Y 2015 Laser Part. Beams 33 11

    [3]

    Kavaya M J, Beyon J Y, Koch G J, Petros M, Petzar P, Singh U N, Trieu B C, Yu J 2014 J. Atmos. Ocean. Tech. 31 826

    [4]

    Wagener T J, Demma N, Kemetec J D, Kubo T S 1995 IEEE Aero. El. Sys. Mag. 10 23

    [5]

    Marano M, Galzerano G, Svelto C, Laporta P 2004 IEEE T. Instrum. Meas. 53 571

    [6]

    Koch G J, Beyon J Y, Petzar P, Petros M, Yu J, Trieu B C, Kavaya M J, Singh U N, Modlin E A, Bames B W, Demoz B B 2010 J. Appl. Remote Sens. 4 043512

    [7]

    Koch G J, Beyon J Y, Bames B W, Petros M, Yu J, Amzajerdian F, Kavaya M J, Singh U N 2007 Opt. Eng. 46 16201

    [8]

    Geisthoff U W, Zenk J, Steinhart H, Iro H 2001 HNO 49 194

    [9]

    Ma Q L, Bo Y, Zong N, Pan Y B, Peng Q J, Cui D F, Xu Z Y 2011 Opt. Commun. 284 1645

    [10]

    Young C G 1966 Appl. Opt. 5 993

    [11]

    Thompson G A, Krupkin V, Yogev A 1992 Opt. Eng. 31 2644

    [12]

    Yabe T, Ohkubo T, Uchida S, Yoshida K, Nakatsuka M, Funatsu T, Mabuti A, Oyama A, Nakagawa K, Oishi T, Daito K, Behgol B, Naayama Y, Yoshida M, Motokoshi S, Sato Y, Baasandash C 2007 Appl. Phys. Lett. 90 261120

    [13]

    Saiki T, Funahashi K, Motokoshi S, Imasaki K, Fujioka K, Fujita H, Nakatsuka M, Yamanaka C 2009 Opt. Commun. 282 614

    [14]

    Saiki T, Motokoshi S, Imasaki K, Nakatsuka M, Yamanaka C, Fujioka K, Fujita H 2009 Opt. Commun. 282 936

    [15]

    杨扬 2007 博士毕业论文 (上海: 上海交通大学)

    Yang Y 2007 Ph. D. Dissertation (Shanghai: Shanghai Jiao Tong University) (in Chinese)

    [16]

    O'Hare J M, Donlan V L 1976 Pyhs. Rev. B 14 3732

    [17]

    Beyatli E, Sumpf B, Demirbas U 2019 Appl. Opt. 58 2973

    [18]

    赵彬, 赵长明, 何建伟, 杨苏辉 2007 光学学报 27 1797

    Zhao B, Zhao C M, He Z W, Yang S H 2007 Acta Opt. Sin. 27 1797

    [19]

    方容川 2001 固体光谱学 (合肥: 中国科学技术大学出版社) 第4页

    Fang R C 2001 Solid State Spectroscopy (Vol. 1) (Hefei: Press of University of Science and Technology of China) p4 (in Chinese)

    [20]

    赵立伟 2010 硕士毕业论文 (北京: 北京理工大学)

    Zhao L W 2010 M. S. Dissertation (Beijing: Beijing Institute of Technology) (in Chinese)

    [21]

    徐鹏 2019 博士毕业论文 (北京: 北京理工大学)

    Xu P 2019 Ph. D. Dissertation (Beijing: Beijing Institute of Technology) (in Chinese)

  • [1] 吕志伟, 丁迎春, 何伟明. 抽运功率密度对布里渊放大的影响. 物理学报, 2002, 51(6): 1286-1290. doi: 10.7498/aps.51.1286
    [2] 林志锋, 张云山, 高春清, 高明伟. LD抽运Cr,Tm,Ho∶YAG微片激光器单纵模运转特性的研究. 物理学报, 2009, 58(3): 1689-1693. doi: 10.7498/aps.58.1689
    [3] 詹敏杰, 邹育婉, 林清峰, 王兆华, 韩海年, 吕亮, 魏志义, 章建, 唐定远. 钛宝石激光抽运的被动锁模Tm:YAG陶瓷激光实验研究. 物理学报, 2014, 63(1): 014205. doi: 10.7498/aps.63.014205
    [4] 柳 强, 巩马理, 李 晨, 宫武鹏, 陆富源, 陈 刚. 角抽运Yb:YAG激光器. 物理学报, 2005, 54(2): 721-725. doi: 10.7498/aps.54.721
    [5] 张新陆, 李 立, 崔金辉, 王月珠, 鞠有伦. 端面抽运Tm,Ho∶YLF连续激光器的参数优化与实验研究. 物理学报, 2008, 57(6): 3519-3524. doi: 10.7498/aps.57.3519
    [6] 张新陆, 李立, 姜波, 王月珠, 鞠有伦. 端面抽运Tm,Ho:YLF激光器双稳特性的理论分析与实验研究. 物理学报, 2009, 58(2): 964-969. doi: 10.7498/aps.58.964
    [7] 延凤平, 魏淮, 傅永军, 王琳, 郑凯, 毛向桥, 刘鹏, 彭健, 刘利松, 简水生. 石英基掺Tm3+包层抽运光纤激光器. 物理学报, 2009, 58(9): 6300-6303. doi: 10.7498/aps.58.6300
    [8] 王月珠, 史洪峰, 张新陆. 激光二极管端面抽运室温Tm,Ho:YLF连续固体激光器. 物理学报, 2006, 55(4): 1787-1792. doi: 10.7498/aps.55.1787
    [9] 张新陆, 李 立, 王月珠, 鞠有伦. 激光二极管端面抽运Tm,Ho:YLF激光器双稳特性研究. 物理学报, 2008, 57(3): 1699-1703. doi: 10.7498/aps.57.1699
    [10] 柳 强, 巩马理, 潘圆圆, 李 晨. 边缘抽运复合Yb:YAG/YAG薄片激光器设计与功率扩展. 物理学报, 2004, 53(7): 2159-2164. doi: 10.7498/aps.53.2159
    [11] 李 立, 张新陆, 王月珠, 鞠有伦. 端面抽运Tm, Ho:YLF激光器热转换系数及热透镜效应的研究. 物理学报, 2007, 56(4): 2196-2201. doi: 10.7498/aps.56.2196
    [12] 刘宏, 王正平, 张怀金, 王继扬, 于浩海, 吕耀辉, 桑元华, 许心光, 兰瑞君, 任诠. LD抽运的Nd:YAG陶瓷/KTP绿光激光器. 物理学报, 2009, 58(11): 7686-7689. doi: 10.7498/aps.58.7686
    [13] 汪超, 韦辉, 王江峰, 姜有恩, 范薇, 李学春. 激光二极管抽运的高重频高平均功率Nd:YAG激光器. 物理学报, 2014, 63(22): 224204. doi: 10.7498/aps.63.224204
    [14] 王月珠, 鞠有伦, 张新陆. 能量传递上转换对Tm,Ho:YLF激光器阈值的影响. 物理学报, 2005, 54(1): 117-122. doi: 10.7498/aps.54.117
    [15] 李哲, 江海河, 王礼, 杨经纬, 吴先友. 2 m Cr,Tm,Ho:YAG激光热退偏效应的数值模拟及实验研究. 物理学报, 2012, 61(4): 044205. doi: 10.7498/aps.61.044205
    [16] 朱洪涛, 楼祺洪, 漆云凤, 马海霞, 董景星, 魏运荣. 钛宝石激光器端面抽运Nd:YAG陶瓷激光器热沉积理论和实验研究. 物理学报, 2005, 54(12): 5648-5653. doi: 10.7498/aps.54.5648
    [17] 刘 欢, 姚建铨, 郑芳华, 路 洋, 王 鹏. LD端面抽运Nd:YAG 1319/1338nm双波长激光器研究. 物理学报, 2008, 57(1): 230-237. doi: 10.7498/aps.57.230
    [18] 张玉萍, 张会云, 王鹏, 李喜福, 姚建铨, 何志红. 36 W侧面抽运腔内倍频Nd:YAG/KTP连续绿光激光器. 物理学报, 2009, 58(7): 4647-4651. doi: 10.7498/aps.58.4647
    [19] 刘欢, 王巍, 巩马理. 角抽运Nd:YAG复合板条946 nm连续运转激光器 . 物理学报, 2013, 62(14): 144205. doi: 10.7498/aps.62.144205
    [20] 张新陆, 陈历学, 李 立. 648nm激光雪崩抽运掺Tm晶体的本征光学双稳特性研究. 物理学报, 2008, 57(1): 278-284. doi: 10.7498/aps.57.278
  • 引用本文:
    Citation:
计量
  • 文章访问数:  150
  • PDF下载量:  3
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-12-27
  • 修回日期:  2020-01-27
  • 刊出日期:  2020-05-01

2 μm人眼安全波段太阳光抽运激光器的理论研究

  • 1. 北京理工大学光电学院, 北京 100081
  • 2. 精密光电测试仪器及技术北京市重点实验室, 北京 100081
  • 3. 信息光子技术工业和信息化部重点实验室, 北京 100081
  • 通信作者: 杨苏辉, suhuiyang@bit.edu.cn
    基金项目: 国家自然科学基金(批准号: 61835001, 61875011)资助的课题

摘要: 太阳光直接抽运激光器在空间光通信、遥感等领域有着重要的潜在应用, 但是一直以来人们对太阳光抽运激光器的研究局限于以掺Nd3+粒子为增益介质的1 μm波段. 通过对现有固体激光工作物质的吸收谱进行分析, 发现掺Tm3+离子在太阳辐射较强的可见光波段具有强的吸收峰, 使2 μm人眼安全波段实现太阳光直接抽运激光输出成为可能. 本文对Tm:YAG和Tm:YAP两种常见晶体的吸收谱与太阳光谱匹配度进行了分析计算, 得出两种材料用于太阳光抽运激光器的阈值抽运功率密度分别为1.14和1.434 kW/cm3. 选择与抽运阈值功率密度低的Tm:YAG晶体作为增益介质, 使用TracePro软件建立太阳光抽运激光器的二级抽运模型, 并对模型进行优化, 得到了锥形腔窗口与菲涅耳透镜的最佳距离、晶体棒的最佳长度以及锥形腔最佳锥度. 本文的工作为实现太阳光直接抽运2 μm激光输出做了理论上的准备.

English Abstract

    • 2 μm人眼安全波段激光在遥感[1,2]、激光雷达[3-7]、激光医疗[8]、空间光通信等领域有着巨大的应用前景, 目前国内外已有多家研究单位在2 μm波段固体激光器的研究上取得重要进展[9]. 太阳光抽运激光器作为一种新型能量利用装置, 能够将太阳光辐射作为抽运源, 将太阳光直接转化为激光, 具有能量转换环节少、转换效率高、使用寿命长等多项优点, 应用前景广泛, 特别是在星载遥感及空间光通信等领域有重要的应用潜力[10].

      目前, 国内外关于太阳光抽运激光器的研究局限于以掺Nd3+离子为增益介质的1 μm波段, 常用的晶体材料为单掺杂的Nd:YAG或双掺杂的Cr:Nd:YAG, Cr:Nd:GSGG晶体和陶瓷材料等[11-14]. 那么太阳光直接抽运是否可以在其他波段实现输出呢? 通过对现有固体激光工作物质的吸收谱进行分析, 发现掺铥离子在太阳辐射较强的可见光波段具有强吸收峰[15], 这使得2 μm人眼安全波段实现太阳光直接抽运激光输出成为可能. 2 μm太阳光直接抽运激光器在星载测风雷达、星载测量CO2温室气体分布等遥感领域有着重要的应用前景.

      掺Tm3+离子的常用抽运波长为793 nm, 由于交叉弛豫效应, 其量子效率可以大于1[16]. 而掺Tm3+离子在可见光波段的吸收是否会为2 μm波段的辐射做出贡献这个问题, 最近也有了答案. Beyatli等[17]的研究结果表明, 对于Tm:YAG晶体, 其吸收谱中位于可见光吸收波段的681 nm吸收峰是793 nm吸收系数的3倍. 通过实验, 采用波长681 nm的激光二极管抽运, 最终实现了波长为2019 nm, 功率为250 mW的连续波输出与超过50%的斜率效率. 由此证明, 掺Tm3+离子晶体在太阳辐射强的可见光波段内能够利用的波带较宽, 在太阳光抽运激光器领域具有巨大发展前景.

      本文基于材料的吸收谱与太阳光谱匹配度计算得到太阳光抽运Tm:YAG/Tm:YAP激光器输出的阈值抽运功率密度. 采用阈值抽运密度较低的晶体作为激光增益介质, 以菲涅耳透镜作为第一级汇聚系统, 锥形聚光腔作为第二级汇聚系统. 使用TracePro软件建立太阳光抽运激光器抽运系统模型, 对锥形腔结构进行优化, 得到了锥形腔最佳位置、晶体最佳长度与锥形腔的最佳锥度等优化参数, 为实验研究做准备.

    • 对于太阳光谱而言, “太阳常数”是指在日地平均距离处垂直于太阳光线的平面上, 单位时间内单位面积能够接收到的辐照度[18]. 通常, 人们所使用的“太阳常数”数值为世界气象组织所确立的1367 W/m2. 但是, 实际在地面上进行测量时, 往往测得的辐照度小于1000 W/m2. 图1为ASTM G173-03 (2012)文件标定在纬度37°下测量得到的太阳光谱曲线. 在对图1所示的光谱曲线原始数据建立插值模型后转化为波长间隔1 nm的数据点, 计算得到该情形下地表接收到的“太阳常量”大小为981.7808 W/(m2·nm).

      图  1  太阳光谱曲线

      Figure 1.  Curve of solar spectrum.

      对于太阳光抽运激光器而言, 确定晶体的吸收光谱是确定激光器抽运波长和偏振方向的依据. 本文结合所确立的太阳光谱模型对Tm:YAG和Tm:YAP晶体做匹配分析发现: Tm:YAG晶体为各向同性介质, 无偏振特性; Tm:YAP为双轴晶体, 取吸收系数最高的b轴偏振方向作为测量方向. 在已知太阳在地面的辐照光谱与晶体吸收光谱后, 假设一定波长下太阳光能够全部照射到材料上, 其入射辐照度为${W_0}(\lambda )$, 则可以推断出晶体对该波长的光辐照度吸收值为[19]

      $W(\lambda ) = {W_0}(\lambda ) \cdot \left( {1 - {{\rm{e}}^{ - l \cdot a}}} \right),$

      式中, l为晶体长度, a为晶体吸收系数. Tm: YAG晶体共有四条主要吸收带: 360—410 nm, 456—480 nm, 656—720 nm, 747—812 nm. Tm: YAP晶体同样有四条主要吸收带: 360—394 nm, 450—494 nm, 643—726 nm, 744—836 nm. 两种晶体的光谱曲线如图2所示.

      图  2  太阳光谱与Tm:YAP, Tm:YAG吸收谱 (a) Tm:YAG; (b) Tm:YAP

      Figure 2.  Matching curve of crystals and solar spectrum: (a) Tm:YAG; (b) Tm:YAP.

      将两种晶体吸收光谱应用于太阳光谱模型后, 对图2中Tm:YAG和Tm:YAP晶体的吸收谱进行插值计算, 分别得到两种晶体各吸收带能够吸收的辐照度能量之和:

      $I = \sum\nolimits_i {M \cdot \Delta \lambda }, $

      式中, M代表晶体吸收对应波长间隔内的平均辐照度密度, $\Delta \lambda $代表插值后的对应波长间隔. 将各能量带吸收能量和与太阳光发射总的辐照度能量相除即得到各能量带的光谱匹配度, 具体分析结果如表1所列.

      Active
      medium
      Doping density/cm3Absorption
      band/nm
      Irradiance in absorption
      band/W·m–2
      Percentage of solar
      radiance/%
      Tm:YAG1.261×1020 (1 at.%)360—41021.9104122.9
      456—48038.78578
      656—72084.80076
      747—81269.35269
      Tm:YAP1.965×1020 (1 at.%)360—39419.5293829.8
      450—49468.23487
      643—726108.88686
      744—83695.81626

      表 1  晶体光谱匹配分析结果

      Table 1.  Spectral matching analysis results of crystals.

      上述计算表明, 两种晶体中Tm:YAP吸收谱与太阳的辐射谱有更好的匹配度, 但是光谱匹配度并不是选择晶体唯一需要考虑的因素, 还要考虑晶体导热、对振荡光的再吸收等其他因素, 下面分别计算以两种晶体为增益介质的太阳光直接抽运的阈值抽运密度.

    • 图3显示了掺Tm3+晶体的能级结构图, 分析能级跃迁过程后, 建立上能级粒子数密度的速率方程及光子数密度的速率方程:

      图  3  Tm3+离子能级跃迁示意图

      Figure 3.  Schematic diagram of Tm3+ ion level transition.

      $\begin{split} & {{\partial {N_2}(x,y,z)}}/{{\partial t}}\\ =\;& {R_4}{r_0}(x,y,z) - {{{N_2}(x,y,z)}}/{{{\tau _2}}} - {K_{\rm{m}}}N_2^2(x,y,z) \\ & - \sigma \frac{c}{n}({f_{\rm{u}}}{N_2}(x,y,z) - {f_1}{N_1})\varPhi {s_0}(x,y,z),\\[-12pt] \end{split} $

      $\begin{split}\frac{{\partial \varPhi }}{{\partial t}} =\; & \frac{{\sigma c}}{n}\iiint\limits_V \left[({f_u} + {f_1}){N_2}(x,y,z) - {f_1}{N_{{\rm{Tm}}}}\right]\\ & \times\varPhi {s_0}(x,y,z){\rm{d}}V - {\varPhi }/{{{\tau _{\rm{c}}}}},\end{split}$

      式中, ${R_4}$为(下同)抽运速率, ${N_2}(x, y, z)$代表上能级的粒子数密度, $\sigma $代表晶体的受激发射面积, ${K_{\rm{m}}}$代表上转换损耗系数, ${\tau _2}$代表上能级寿命, $n$代表晶体折射率, $\varPhi $代表腔内的光子数, ${f_{\rm{u}}}$代表处于上能级的粒子玻尔兹曼分布因子, ${f_1}$代表处于下能级的粒子玻尔兹曼分布因子, ${N_{{\rm{Tm}}}}$代表总的粒子数浓度, ${\tau _{\rm{c}}}$代表光子数平均寿命, ${r_0}(x, y, z)$${s_0}(x, y, z)$分别代表归一化的抽运速率分布及腔内光子数分布.

      $ \begin{split} & \iiint\limits_V {r_0}(x,y,z){\rm{d}}V = \iiint\limits_V {{s_0}(x,y,z)}{\rm{d}}V = 1. \end{split}$

      当激光器在稳态下工作时, 上能级粒子数及腔内光子数均不随时间变化,

      $\frac{{\partial {N_2}}}{{\partial t}} = \frac{{\partial \varPhi }}{{\partial t}} = 0.$

      且当激光器处于阈值状态时, 腔内光子数较少, 可以认为$\varPhi \approx 0$. 将各项条件代入速率方程后上能级粒子数密度${N_2}(x, y, z)$与抽运速率${R_4}$可以表示为

      ${N_2}(x,y,z) = \frac{{ - \dfrac{1}{{{\tau _2}}} + \sqrt {\dfrac{1}{{\tau _2^2}} + 4{R_4}{r_0}(x,y,z){K_{\rm{m}}}} }}{{2{K_{\rm{m}}}}},$

      ${R_4} = \dfrac{{\dfrac{n}{{{\tau _{\rm{c}}}\sigma c}} + {f_1}{N_{{\rm{Tm}}}}}}{{({f_{\rm{u}}} + {f_1}){\tau _2}}},$

      式中, ${\tau _{\rm{c}}} ={{2 L}}/{{c\delta }}$, $\delta $代表谐振腔往返损耗, $c$代表光速. 最终得到阈值条件下抽运速率的表达式为

      ${R_4} = \dfrac{{\dfrac{{n\delta }}{{2L\sigma }} + {f_1}{N_{{\rm{Tm}}}}}}{{({f_{\rm{u}}} + {f_1}){\tau _2}}}.$

    • 与传统激光器抽运阈值不同, 对于太阳光抽运激光器, 可以将晶体材料各吸收带吸收的太阳光辐射功率用来表示激光器的抽运速率[20], 即

      ${R_4} = \sum\nolimits_i {{N_i}} .$

      ${N_i}$(i = 1, 2, 3, ···)代表各吸收带吸收的太阳辐射功率对抽运速率的贡献. 对于激光介质的各个吸收带, 均有${P_i} {{\bar a_i}} {\eta _{\rm{u}}} = {N_i}h{\upsilon _{\rm{L}}}$. 式中, ${\eta _{\rm{u}}}$代表吸收的太阳光转化为激光发射功率的效率, ${\eta _{\rm{u}}} = {\eta _{\rm{Q}}}{\eta _{\rm{S}}}$, ${\eta _{\rm{Q}}}$代表量子效率, ${\eta _{\rm{S}}}$为斯托克斯因子, 表示激光跃迁时发射的光子能量与抽运光子能量之比; $ {{\bar a_i}} $代表整个吸收带内的平均吸收系数, ${P_i}$代表第i个吸收带内的太阳辐射功率, ${P_i} = {P_{{\rm{solar}}}} \cdot {\eta _i}$(${\eta _i}$表示各吸收带内的太阳辐射功率占整个太阳光辐射功率的比例), ${\nu _{\rm{L}}}$代表激光跃迁频率.

      最终, 将(9)式代入(10)式后得到太阳光抽运激光器的阈值抽运功率密度为

      ${P_{{\rm{solar}}{\rm{.th}}}} = \dfrac{{hc}}{{{\eta _Q}{\tau _2}}} \cdot \dfrac{{\dfrac{{n\delta }}{{2L\sigma }} + {f_1}{N_{{\rm{Tm}}}}}}{{({f_1} + {f_{\rm{u}}})}} \cdot \dfrac{1}{{\sum\nolimits_i {{\eta _i}\overline {{a_i}} {\lambda _i}} }}.$

      表2列出了Tm:YAG与Tm:YAP两种晶体材料的晶体参数:

      Tm:YAGTm:YAP
      Doping density/cm31.26×1020
      (1 at.%)
      1.965×1020
      (1 at.%)
      Upper level lifetime/ms10.54.4
      Boltzmann factor
      in upper level
      0.460.29
      Boltzmann factor
      in lower level
      0.0170.015
      Emission cross section/cm22.5×10–213.81×10–21
      Quantum efficiency1.81.9
      Calculation results of
      absorption coefficient
      curve: $\sum\nolimits_i { {\eta _i}\overline { {a_i} } {\lambda _i} } $/cm
      46.0069×10–751.5739×10–7
      Refractive index1.821.91

      表 2  晶体参数

      Table 2.  Crystal parameters.

      此外, 在计算阈值的过程中用到的其他参数有: 晶体棒长度$l$取70 mm, 谐振腔光学长度$L = l \cdot n + x$, $x$为腔内除晶体棒外其余部分长度, $n$为晶体棒折射率, 最终取谐振腔光学长度为150 mm, 谐振腔损耗$\delta $ = 0.01. 最终结算结果为: Tm:YAG晶体的抽运阈值密度为1.14 kW/cm3, Tm:YAP晶体的抽运阈值密度为1.434 kW/cm3.

    • 基于上文中计算的阈值参数以及菲涅耳透镜、锥形抽运腔的二级抽运系统, 在TracePro软件中建立太阳光抽运Tm:YAG激光器的抽运模型. 图4为TracePro软件中建立的太阳光抽运Tm:YAG激光器系统二级抽运系统模型与模型中使用的锥形腔结构图. 为了完整构建太阳光抽运Tm:YAG激光器的系统模型, 需要在模型中分别对太阳光源、菲涅耳透镜、锥形抽运腔、Tm:YAG晶体棒与冷却液体忽略冷却液流动对抽运光分布带来的影响以及锥形腔入射窗口对入射光的吸收.

      图  4  (a) TracePro软件建立的太阳光抽运激光器二级抽运模型; (b)锥形腔结构图

      Figure 4.  (a) Two-stage pumping model; (b) structure diagram of conical cavity.

      1)软件中设置太阳光源功率密度与上文中计算参数一致, 为981.7808 W/m2, 太阳辐射的发散角设为0.27°, 由于距离较远, 将光源设置为表面光源, 定义其一面为发光面.

      2)菲涅耳透镜使用PMMA (聚甲基丙烯酸甲酯)材料, 焦距为1500 mm, 环距为0.33 mm, 厚度为5 mm, 通光面面型为直径2 m的圆.

      3)锥形腔材料设为紫铜结构, 内表面使用陶瓷反射面, 腔内冷却液设置为水.

      4)晶体选择Tm:YAG材料, 其侧面设置为磨砂面.

    • PMMA材质的菲涅耳透镜作为一种投射式聚光器, 其材料的色散特性会导致不同波长的光经过菲涅耳透镜后汇聚焦点不同[21]. 因此, 对于太阳光抽运激光器而言, 不能简单地将锥形腔窗口设置在菲涅耳透镜焦点处, 需要寻找锥形腔窗口的最佳位置. 在软件中, 设置模拟追迹光线条数为50000条, 锥形腔窗口尺寸设置为Ф40 mm, 将锥形腔在菲涅耳透镜的焦点附近来回移动, 改变腔与菲涅耳透镜之间的距离, 对锥形腔窗口位置进行优化. 图5为模拟的锥形腔在不同位置下接收入射光的能量分布.

      图  5  锥形腔窗口位置与接受光功率关系图

      Figure 5.  Curve of the relationship between the position of conical cavity and the received solar power.

      图5的横坐标代表的是锥形腔窗口距菲涅耳透镜的距离, 当锥形腔窗口位置设置在距离菲涅耳透镜1506 mm处时, 接收到的光强最大. 因此, 锥形腔的最佳位置应设置在距离菲涅耳透镜1506 mm处.

    • 由于锥形腔的存在, 由菲涅耳透镜汇聚的太阳光能够以端面抽运与侧面抽运结合的方式混合抽运激光晶体, 由于晶体长度较长, 晶体棒接收到的侧面抽运光通量远远大于端面抽运. 因此计算晶体的侧面抽运阈值辐照度对于晶体长度的选择具有重要意义. 设晶体棒半径为$r$, 由(11)式得晶体抽运阈值密度为${P_{{\rm{solar}}{\rm{.th}}}}$, 晶体长度为$l$, 晶体侧面抽运阈值辐照度为${E_0}$. 仅当满足(12)式时, 晶体能够出光.

      ${E_0}{\text{π}2}r{\rm{d}}l \geqslant {P_{{\rm{soalr}}{\rm{.th}}}}{\text{π}}{r^2}{\rm{d}}l.$

      最终计算得到半径为r的晶体棒侧面抽运阈值辐照度为

      ${E_0} = {{{P_{{\rm{solar}}{\rm{.th}}}} \cdot r}}/{2}.$

      在TracePro软件中将Tm:YAG晶体长度设置为100 mm, 锥形腔入射窗口尺寸为Ф40 mm, 出射窗口尺寸为Ф6 mm. 如图6所示, 将锥形腔位置放置于距离菲涅耳透镜1506 mm处, 进行光线追迹, 得到晶体棒侧面入射光通量分布.

      图  6  100 mm晶体棒侧面抽运光分布图

      Figure 6.  Distribution map of side pump power on 100 mm-length crystal rod.

      根据(13)式, 对于尺寸为Ф6 mm的Tm:YAG晶体, 仅考虑侧面抽运情况下, 计算得到侧面抽运的阈值辐照度为1.71 × 106 W/m2. 从图6可以地明显看出, 尽管晶体棒侧面接收到的入射总光通量够大, 但是位于锥形腔后端的晶体侧面接收到的光很少, 无法达到阈值, 晶体未能起振, 造成了抽运光的浪费. 因此, 需要对锥形腔与晶体的长度进行优化. 将晶体长度分别减少为90, 80, 70, 60和50 mm进行光线追迹, 观测晶体侧面入射光分布并进行对比, 得到了图7所示的晶体棒轴向抽运光辐照度分布图.

      图  7  不同晶体长度下晶体棒轴向光辐照度分布

      Figure 7.  Axial irradiance distribution of different-length crystals.

      图7可知, 当逐渐减小晶体棒长时, 晶体棒前端接收到的抽运光辐照度逐步增大, 证明原本晶体后端并未吸收的抽运光重新分布在晶体前端, 增大了前端晶体的抽运功率密度, 使得腔内抽运光的利用率得到了有效提升. 当长度减小到70 mm时, 晶体棒侧面接收到的光辐照度值均大于Ф6 mm晶体的抽运阈值辐照度1.71 × 106 W/m2, 证明该长度下整个晶体均能起振, 输出激光. 继续减小晶体棒长度, 可以发现晶体棒接收到的抽运光平均辐照度继续增大, 但抽运总功率变小. 为确定晶体棒的最佳长度, 对介质增益系数进行计算. 对于激光介质, 小信号增益系数可表示为

      $\begin{split}{g_0} \; &= \Delta N \cdot \sigma = \left( {{f_{\rm{u}}}{N_2} - {f_1}{N_1}} \right) \cdot \sigma \\ & = \left[ {({f_{\rm{u}}} + {f_1}){N_2} - {f_1}{N_{{\rm{Tm}}}}} \right] \cdot \sigma ,\end{split}$

      式中, $\Delta N$表示反转粒子数密度, $\sigma $表示发射截面面积, 将(8)式与(13)式代入(14)式后得到

      $\begin{split}{g_0} =\; & \Bigg\{\frac1{{2{K_{\rm{m}}}}}\bigg[{{ - \dfrac{1}{{{\tau _2}}} + \sqrt {\dfrac{1}{{{\tau _2}}} + 4{r_0}(x,y,z){K_{\rm{m}}} \cdot \dfrac{{{P_{{\rm{solar}}}}{\eta _{\rm{u}}}}}{{h{\upsilon _L}}} \cdot \sum\limits_i {{\eta _i} {{\bar a_i}} } } }} \bigg] ({f_{\rm{u}}} + {f_1}) - {f_1}{N_{{\rm{Tm}}}}\Bigg\}\sigma . \end{split}$

      从(15)式可得小信号增益系数随抽运太阳光功率增大而增大, 因此在辐照度达到阈值、晶体利用率最大的情况下, 70 mm长晶体受到的抽运光功率最大, 其增益系数最大, 因此最终选择晶体长度为70 mm.

    • 使用软件对锥形腔锥度进行优化设计, 晶体长度设置为70 mm, 锥形腔有效聚光长度同样设置为70 mm, 将锥形腔入射窗口尺寸设置为Φ40 mm, 改变锥形腔出射窗口尺寸, 即改变锥形腔的锥度, 让锥形腔出射窗口尺寸分别为Φ6 mm, Φ12 mm, Φ18 mm, Φ24 mm, Φ30 mm和Φ36 mm时, 晶体棒轴向光辐照度分布如图8所示.

      图  8  不同锥度下晶体棒轴向光辐照度分布图

      Figure 8.  Axial irradiance distribution of different-taper crystals.

      图8可以发现, 当锥形腔锥度较大时, 其出射窗口尺寸小, 晶体棒接收到的侧面抽运总光通量大, 但是在晶体后端接收到的光通量却较少, 无法达到阈值, 晶体无法起振, 造成了晶体材料的浪费. 当逐渐减少锥形腔锥度, 即增大锥形腔出射窗口尺寸时, 晶体接收到的总光通量变少, 但晶体棒各部分接收到的光通量更加均匀, 晶体利用率得到提高. 但当锥形腔锥度减小到一定度数后, 锥形腔的二次汇聚作用减小, 入射抽运光在腔内反射次数过少, 溢出锥形腔的光过多, 入射光利用率过低, 同样也会造成晶体后端无法起振现象, 使晶体利用率降低. 综合考虑光通量与抽运光均匀度后, 确定锥形腔最佳锥度, 出射窗口尺寸选择为Φ12 mm, 此时入射总光通量最大, 光照均匀度好, 且晶体棒利用率高.

    • 本文提出了一种采用掺Tm3+离子晶体作为增益介质的新型太阳光抽运激光器的设计方法, 能够发出2 μm的人眼安全波段激光. 对两种常见掺Tm3+离子晶体Tm:YAG和Tm:YAP的吸收光谱与光谱匹配度进行了计算, 并计算出以两种晶体用于太阳光抽运激光器的阈值抽运功率密度分别为1.14 和1.434 kW/cm3. 采用阈值较小的Tm:YAG晶体, 运用TracePro软件建立了太阳光抽运激光器的两级抽运模型, 并确定了其锥形腔窗口与菲涅耳透镜的最佳距离为1506 mm, 锥形腔的最佳长度为70 mm, 锥形腔最佳锥度, 即出射窗口尺寸选择为Φ12 mm.

      在未来工作方面有一点特别需要注意的地方: 以Tm:YAG或Tm:YAP晶体作为增益介质的太阳光抽运激光器的抽运阈值功率密度比传统Nd:YAG激光器大, 为达到工作阈值, 需要汇聚更多的入射光, 从而导致锥形腔内热量增加, 且掺Tm:YAG与Tm:YAP晶体为准三能级结构, 其工作时热效应更加严重. 对于本太阳光抽运激光器而言, 腔内晶体的冷却工作是未来工作重点, 需要寻找新型漫反射冷却液或采用热键合材质晶体以缓解热效应.

参考文献 (21)

目录

    /

    返回文章
    返回