搜索

x

留言板

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

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

基于抽运-探测法的皮秒反斯托克斯拉曼频移器的理论研究

王聪 吕冬翔

引用本文:
Citation:

基于抽运-探测法的皮秒反斯托克斯拉曼频移器的理论研究

王聪, 吕冬翔

Theoretical study of picosecond anti-Stokes Raman frequency converter based on pump-probe method

Wang Cong, Lü Dong-Xiang
PDF
HTML
导出引用
  • 采用拉曼频移器在晶体介质中利用相干反斯托克斯散射效应可以获得超短脉冲(皮秒)反斯托克斯激光. 基于抽运-探测法的晶体拉曼频移器可以实现相干反斯托克斯散射的共线相互作用, 从而可以有效提高反斯托克斯光的转化效率. 本文在平面波近似下建立了基于抽运-探测法的皮秒反斯托克斯拉曼频移器的耦合波方程组, 引入归一化参量对方程组进行了归一化处理. 通过数值计算, 得到了描述皮秒反斯托克斯拉曼频移器运行的一组普适理论曲线, 分析了归一化拉曼增益系数G、归一化相位失配参数ΔK以及探测光脉冲能量占基频光总能量的比值rprobe三个变量对反斯托克斯拉曼频移器性能的影响, 确定了实现高效反斯托克斯转化时各归一化变量的合理取值. 采用实验数据对该理论模型的正确性进行了验证, 反斯托克斯转化效率的理论值与文献数据基本一致.
    Ultra-short pulse (picosecond) anti-Stokes laser can be obtained by using Raman frequency converter in a crystal medium by the coherent anti-Stokes Raman scattering effect. The crystalline Raman frequency converter based on the pump-probe method can realize the collinear interaction of coherent anti-Stokes Raman scattering, thus effectively improving the conversion efficiency of the anti-Stokes light. Theoretical simulation is an important means to study laser operation. Coupled wave equation is widely used to study the characteristics of Raman laser and anti-Stokes laser. Although the coupling wave theory of anti-Stokes Raman frequency shifter reported previously can reflect the operation law of the frequency shifter, the optimization of the frequency shifter and the influence of the frequency shifter parameters on the output characteristics of anti-Stokes laser have not been reported so far. In this paper, the picosecond anti-Stokes Raman frequency converter based on the pump-probe method is studied theoretically. Considering the generation of the first Stokes light in the probe channel and the second Stokes light in the pump channel, the coupled wave equation of the collinear picosecond anti-Stokes Raman frequency converter is established under the plane wave approximation. Without loss of generality, four dimensionless comprehensive parameters are introduced to normalize the equations. A set of universal theoretical curves describing the operation of the Raman frequency converter is obtained. The numerical solutions of the equations show that the performance of the Raman frequency converter mainly depends on three parameters: the normalized phase mismatch parameter ΔK, the normalized Raman gain coefficient G, and the energy ratio of the probe light to the fundamental light rprobe. The reasonable values of normalized variables are determined when the high efficiency anti-Stokes conversion is realized. Experimental data are used to verify the correctness of the theoretical model. The theoretical value of the anti-Stokes conversion efficiency is basically consistent with the literature data. The normalized coupled wave theory proposed in this paper is helpful in understanding the operation law of the picosecond anti-Stokes Raman frequency shifter, and has guiding significance for the design of the frequency converter.
      通信作者: 王聪, wangc.sd@163.com
    • 基金项目: 国家自然科学基金(批准号: 11504268)和天津市教委科研计划(批准号: 2014KJ240)资助的课题
      Corresponding author: Wang Cong, wangc.sd@163.com
    • Funds: Project supported by the National Natural Science Foundation of China (Grant No. 11504268) and the Scientific Research Program of Tianjin Municipal Education Commission, China (Grant No. 2017KJ240)
    [1]

    孙瑛璐, 段延敏, 程梦瑶, 袁先漳, 张立, 张栋, 朱海永 2020 物理学报 69 124201Google Scholar

    Sun Y L, Duan Y M, Cheng M Y, Yuan X Z, Zhang L, Zhang D, Zhu H Y 2020 Acta Phys. Sin. 69 124201Google Scholar

    [2]

    张蕴川, 樊莉, 魏晨飞, 顾晓敏, 任思贤 2018 物理学报 67 024206Google Scholar

    Zhang Y C, Fan L, Wei C F, Gu X M, Ren S X 2018 Acta Phys. Sin. 67 024206Google Scholar

    [3]

    Zhou Q Q, Shi S C, Chen S M, Duan Y M, Zhang X M, Guo J, Zhao B, Zhu H Y 2019 Chin. Phys. Lett. 36 014205Google Scholar

    [4]

    Sheng Q, Li R, Lee A J, Spence D J, Pask H M 2019 Opt. Express 27 8540Google Scholar

    [5]

    Wei L J, Chen M T, Zhu S Q, Dai S B, Yin H, Chen Z Q 2018 Laser Phys. Lett. 15 125001Google Scholar

    [6]

    Duan Y M, Sun Y L, Zhu H Y, Mao T W, Zhang L, Chen X 2020 Opt. Lett. 45 2564Google Scholar

    [7]

    Liu J, Ding X, Jiang P B, Sheng Q, Yu X Y, Sun B, Wang J B, Shi R, Zhao L, Bai Y T, Zhang G Z, Wu L, Yao J Q 2018 Appl. Opt. 57 3154Google Scholar

    [8]

    Ding S H, Zhang X Y, Wang Q P, Jia P, Zhang C, Liu B 2006 Opt. Commun. 267 480Google Scholar

    [9]

    Smetanin S N, Doroshenko M E, Ivleva L I, Jelínek M, V. Kubeček, Jelínková H 2014 Appl. Phys. B 117 225Google Scholar

    [10]

    Ding S H, Zhang X Y, Wang Q P, Zhang J, Wang S M, Liu Y R, Zhang X H 2007 J. Phys. D 40 2736

    [11]

    Smetanin S N, Jelínek M, Tereshchenko D P, Kubeček V 2018 Opt. Express 26 22637Google Scholar

    [12]

    Smetanin S N, Jelínek M, Kubeček V 2017 Appl. Phys. B 123 203

    [13]

    Vermeulen N, Debaes C, Fotiadi A A, Panajotov K, Thienpont H 2006 IEEE J. Quantum Electron. 42 1144Google Scholar

    [14]

    Grasiuk A Z, Kurbasov S V, Losev L L 2004 Opt. Commun. 240 239Google Scholar

    [15]

    Mildren R P, Coutts D W, Spence D J 2009 Opt. Express 17 810Google Scholar

    [16]

    Wang C, Zhang X Y, Wang Q P, Cong Z H, Liu Z J, Wei W, Wang W T, Wu Z G, Zhang Y G, Li L, Chen X H, Li P, Zhang H J, Ding S H 2013 Opt. Express 21 26014Google Scholar

    [17]

    Wang C, Cong Z H, Qin Z G, Zhang X Y, Wei W, Wang W T, Zhang Y G, Zhang H J, Yu H H 2014 Opt. Commun. 322 44Google Scholar

    [18]

    Wei W, Zhang X Y, Wang Q P, Wang C, Cong Z H, Chen X H, Liu Z J, Wang W T, Wu Z G, Ding S H, Tu C Y, Li Y F, Cheng W Y 2014 Appl. Phys. B 116 561

    [19]

    Smetanin S N, Jelínek M, Tereshchenko D P, Shukshin V E, Konyukhov M V, Papashvili A G, Voronina I S, Ivleva L I, Kubeček V 2020 Opt. Express 28 22919Google Scholar

    [20]

    Shen Y R, Bloembergen N 1965 Phys. Rev. 137 1787Google Scholar

  • 图 1  抽运通道和探测通道中的拉曼散射能级图 (a) 抽运通道的SSRS能级图; (b) 探测通道的CARS能级图; (c) 探测通道的SSRS能级图

    Fig. 1.  Raman scattering energy levels in pump and probe channels: (a) SSRS in pump channel; (b) CARS in probe channel; (c) SSRS in probe channel.

    图 2  G = 90, Wp = 0.2, rprobe取不同值时(a) ηa, (b) η1s和(c) η2s随ΔK的变化

    Fig. 2.  (a) ηa, (b) η1s and (c) η2s versus ΔK for different rprobe with G = 90 and Wp = 0.2.

    图 3  rprobe = 0.3, G = 90, Wp = 0.2时, 抽运光、探测光、一阶斯托克斯光和反斯托克斯光归一化光强随ζ的空间演化 (a) |ΔK| = 0; (b) |ΔK| = 4; (c) |ΔK| = 8

    Fig. 3.  Plots of the spatial evolution of pump, probe, first Stokes, and anti-Stokes normalized intensities with rprobe = 0.3, G = 90 and Wp = 0.2: (a) |ΔK| = 0; (b) |ΔK| = 4; and (c) |ΔK| = 8.

    图 4  rprobe = 0.39, G = 90, Wp = 0.2时, 抽运光、探测光、一阶斯托克斯光和反斯托克斯光归一化光强随ζ的空间演化 (a) |ΔK| = 0; (b) |ΔK| = 4; (c) |ΔK| = 8

    Fig. 4.  Plots of the spatial evolution of pump, probe, first Stokes, and anti-Stokes normalized intensities with rprobe = 0.39, G = 90 and Wp = 0.2: (a) |ΔK| = 0; (b) |ΔK| = 4; and (c) |ΔK| = 8.

    图 5  G = 90, |ΔK| = 0, Wp = 0.2时抽运光、探测光、一阶斯托克斯光、二阶斯托克斯光和反斯托克斯光的脉冲形状 (a) rprobe = 0.3; (b) rprobe = 0.39

    Fig. 5.  Temporal profiles of the pump, probe, first Stokes, second Stokes and anti-Stokes pulses with G = 90, |ΔK| = 0 and Wp = 0.2: (a) rprobe = 0.3; (b) rprobe = 0.39.

    图 6  K| = 0, Wp = 0.2, rprobe取不同值时(a) ηa, (b) η1 s和(c) η2 sG的变化

    Fig. 6.  (a) ηa, (b) η1s and (c) η2s versus G for different rprobe with |ΔK| = 0 and Wp = 0.2.

    图 7  K| = 0, Wp = 0.2时, roptηamaxG的变化

    Fig. 7.  ropt and ηamax versus G with |ΔK| = 0 and Wp = 0.2

    图 8  K| = 0, Wp = 0.2时抽运光、探测光、一阶斯托克斯光、二阶斯托克斯光和反斯托克斯光的脉冲形状 (a) ropt = 0.270, G = 60; (b) ropt = 0.373, G = 110; (c) ropt = 0.414, G = 160

    Fig. 8.  Temporal profiles of the pump, probe, first Stokes, second Stokes, and anti-Stokes pulses with |ΔK |= 0 and Wp = 0.2: (a) ropt = 0.270, G = 60; (b) ropt = 0.373, G = 110; (c) ropt = 0.414, G = 160.

    图 9  K| = 0, Wp = 0.2时, Goptηamaxrprobe的变化

    Fig. 9.  Gopt and ηamax versus rprobe with |ΔK| = 0 and Wp = 0.2

    表 1  参考文献[12]中的参数

    Table 1.  Parameters in Ref. [12].

    参数参数
    νa/νp1.06Ap/cm22.83 × 10–4
    ν1s/νp0.94lR/cm3.2
    ν2s/νp0.88g/(cm·GW–1)13
    wp/ps20Δk0
    下载: 导出CSV

    表 2  不同情况下反斯托克斯转化效率的理论值与实验数据的对比结果

    Table 2.  Comparisons of theoretical and experimental results of anti-Stokes conversion efficiency under different conditions.

    抽运脉冲
    能量/μJ
    探测脉冲
    能量/μJ
    Grprobe实验
    值/%
    理论
    值/%
    30252060.452.802.54
    3021180.0631.881.86
    26121400.323.504.13
    下载: 导出CSV

    表 3  高斯近似时反斯托克斯转化效率的理论值

    Table 3.  Theoretical values of anti-Stokes conversion efficiency for Gaussian approximation.

    抽运脉冲能量/μJ探测脉冲能量/μJ$ \eta _{\rm{a}}^{\rm{g}}$/%
    30252.63
    3021.91
    26124.09
    下载: 导出CSV
  • [1]

    孙瑛璐, 段延敏, 程梦瑶, 袁先漳, 张立, 张栋, 朱海永 2020 物理学报 69 124201Google Scholar

    Sun Y L, Duan Y M, Cheng M Y, Yuan X Z, Zhang L, Zhang D, Zhu H Y 2020 Acta Phys. Sin. 69 124201Google Scholar

    [2]

    张蕴川, 樊莉, 魏晨飞, 顾晓敏, 任思贤 2018 物理学报 67 024206Google Scholar

    Zhang Y C, Fan L, Wei C F, Gu X M, Ren S X 2018 Acta Phys. Sin. 67 024206Google Scholar

    [3]

    Zhou Q Q, Shi S C, Chen S M, Duan Y M, Zhang X M, Guo J, Zhao B, Zhu H Y 2019 Chin. Phys. Lett. 36 014205Google Scholar

    [4]

    Sheng Q, Li R, Lee A J, Spence D J, Pask H M 2019 Opt. Express 27 8540Google Scholar

    [5]

    Wei L J, Chen M T, Zhu S Q, Dai S B, Yin H, Chen Z Q 2018 Laser Phys. Lett. 15 125001Google Scholar

    [6]

    Duan Y M, Sun Y L, Zhu H Y, Mao T W, Zhang L, Chen X 2020 Opt. Lett. 45 2564Google Scholar

    [7]

    Liu J, Ding X, Jiang P B, Sheng Q, Yu X Y, Sun B, Wang J B, Shi R, Zhao L, Bai Y T, Zhang G Z, Wu L, Yao J Q 2018 Appl. Opt. 57 3154Google Scholar

    [8]

    Ding S H, Zhang X Y, Wang Q P, Jia P, Zhang C, Liu B 2006 Opt. Commun. 267 480Google Scholar

    [9]

    Smetanin S N, Doroshenko M E, Ivleva L I, Jelínek M, V. Kubeček, Jelínková H 2014 Appl. Phys. B 117 225Google Scholar

    [10]

    Ding S H, Zhang X Y, Wang Q P, Zhang J, Wang S M, Liu Y R, Zhang X H 2007 J. Phys. D 40 2736

    [11]

    Smetanin S N, Jelínek M, Tereshchenko D P, Kubeček V 2018 Opt. Express 26 22637Google Scholar

    [12]

    Smetanin S N, Jelínek M, Kubeček V 2017 Appl. Phys. B 123 203

    [13]

    Vermeulen N, Debaes C, Fotiadi A A, Panajotov K, Thienpont H 2006 IEEE J. Quantum Electron. 42 1144Google Scholar

    [14]

    Grasiuk A Z, Kurbasov S V, Losev L L 2004 Opt. Commun. 240 239Google Scholar

    [15]

    Mildren R P, Coutts D W, Spence D J 2009 Opt. Express 17 810Google Scholar

    [16]

    Wang C, Zhang X Y, Wang Q P, Cong Z H, Liu Z J, Wei W, Wang W T, Wu Z G, Zhang Y G, Li L, Chen X H, Li P, Zhang H J, Ding S H 2013 Opt. Express 21 26014Google Scholar

    [17]

    Wang C, Cong Z H, Qin Z G, Zhang X Y, Wei W, Wang W T, Zhang Y G, Zhang H J, Yu H H 2014 Opt. Commun. 322 44Google Scholar

    [18]

    Wei W, Zhang X Y, Wang Q P, Wang C, Cong Z H, Chen X H, Liu Z J, Wang W T, Wu Z G, Ding S H, Tu C Y, Li Y F, Cheng W Y 2014 Appl. Phys. B 116 561

    [19]

    Smetanin S N, Jelínek M, Tereshchenko D P, Shukshin V E, Konyukhov M V, Papashvili A G, Voronina I S, Ivleva L I, Kubeček V 2020 Opt. Express 28 22919Google Scholar

    [20]

    Shen Y R, Bloembergen N 1965 Phys. Rev. 137 1787Google Scholar

  • [1] 田子阳, 赵会杰, 尉昊赟, 李岩. 基于混合飞秒/皮秒相干反斯托克斯拉曼散射的动态高温燃烧场温度测量. 物理学报, 2021, 70(21): 214203. doi: 10.7498/aps.70.20211144
    [2] 李健康, 李睿. 利用数值模拟研究表面增强相干反斯托克斯拉曼散射增强基底. 物理学报, 2021, 70(10): 104207. doi: 10.7498/aps.70.20201773
    [3] 彭亚晶, 孙爽, 宋云飞, 杨延强. 液相硝基甲烷分子振动特性的相干反斯托克斯拉曼散射光谱. 物理学报, 2018, 67(2): 024208. doi: 10.7498/aps.67.20171828
    [4] 侯国辉, 罗腾, 陈秉灵, 刘杰, 林子扬, 陈丹妮, 屈军乐. 双光子荧光与相干反斯托克斯拉曼散射显微成像技术的实验研究. 物理学报, 2017, 66(10): 104204. doi: 10.7498/aps.66.104204
    [5] 郑娟娟, 姚保利, 邵晓鹏. 基于光强传输方程相位成像的宽场相干反斯托克斯拉曼散射显微背景抑制. 物理学报, 2017, 66(11): 114206. doi: 10.7498/aps.66.114206
    [6] 刘双龙, 刘伟, 陈丹妮, 屈军乐, 牛憨笨. 相干反斯托克斯拉曼散射显微成像技术研究. 物理学报, 2016, 65(6): 064204. doi: 10.7498/aps.65.064204
    [7] 张赛文, 陈丹妮, 刘双龙, 刘伟, 牛憨笨. 纳米分辨相干反斯托克斯拉曼散射显微成像. 物理学报, 2015, 64(22): 223301. doi: 10.7498/aps.64.223301
    [8] 李亚晖, 梁闰富, 邱俊鹏, 林子扬, 屈军乐, 刘立新, 尹君, 牛憨笨. 紧聚焦条件下相干反斯托克斯拉曼散射信号场的矢量分析. 物理学报, 2014, 63(23): 233301. doi: 10.7498/aps.63.233301
    [9] 刘双龙, 刘伟, 陈丹妮, 牛憨笨. 超衍射极限相干反斯托克斯拉曼散射显微成像技术中空心光束的形成. 物理学报, 2014, 63(21): 214601. doi: 10.7498/aps.63.214601
    [10] 尹君, 余锋, 侯国辉, 梁闰富, 田宇亮, 林子扬, 牛憨笨. 多色宽带相干反斯托克斯拉曼散射过程的理论与实验研究. 物理学报, 2014, 63(7): 073301. doi: 10.7498/aps.63.073301
    [11] 刘伟, 陈丹妮, 刘双龙, 牛憨笨. 超衍射极限相干反斯托克斯拉曼散射显微成像技术及其探测极限分析. 物理学报, 2013, 62(16): 164202. doi: 10.7498/aps.62.164202
    [12] 于凌尧, 尹君, 万辉, 刘星, 屈军乐, 牛憨笨, 林子扬. 基于超连续光谱激发的时间分辨相干反斯托克斯拉曼散射方法与实验研究. 物理学报, 2010, 59(8): 5406-5411. doi: 10.7498/aps.59.5406
    [13] 汪 莎, 陈 军, 童立新, 高清松, 刘 崇, 唐 淳. 熔石英棒-光纤构成的新型复合相位共轭镜的实验和理论研究. 物理学报, 2008, 57(3): 1719-1724. doi: 10.7498/aps.57.1719
    [14] 王峨锋, 李宏福, 李 浩, 喻 胜, 牛新建, 刘迎辉. 螺旋波纹波导研究. 物理学报, 2005, 54(11): 5339-5343. doi: 10.7498/aps.54.5339
    [15] 李文博, 李克轩. 开普勒径向方程的赝角动量解法及其归一化本征态和相干态. 物理学报, 2004, 53(9): 2964-2969. doi: 10.7498/aps.53.2964
    [16] 牛新建, 李宏福, 喻胜, 谢仲怜, 杨仕文. 8mm高功率过模弯曲圆波导TE01—TM11模式变换. 物理学报, 2002, 51(10): 2291-2295. doi: 10.7498/aps.51.2291
    [17] 武鹏飞, 巩雄, 张桂兰, 汤国庆, 陈文驹. 饱和与反饱和吸收介质的简并四波混频和全息随时间变化规律. 物理学报, 1996, 45(3): 380-388. doi: 10.7498/aps.45.380
    [18] 张瑞华, 米辛, 周海天, 姜谦, 叶佩弦. 非相干光时延四波混频的多能级理论——非均匀加宽情形. 物理学报, 1991, 40(3): 414-423. doi: 10.7498/aps.40.414
    [19] 丁宏玉, 吴百诗, 殷大钧, 李锦泉. 简并四波混频的速率方程理论. 物理学报, 1988, 37(3): 408-415. doi: 10.7498/aps.37.408
    [20] 范俊颖, 吴存恺, 王志英. 共振吸收增强的四波混频理论. 物理学报, 1980, 29(7): 897-904. doi: 10.7498/aps.29.897
计量
  • 文章访问数:  1469
  • PDF下载量:  21
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-08-17
  • 修回日期:  2021-02-17
  • 上网日期:  2021-04-19
  • 刊出日期:  2021-05-05

基于抽运-探测法的皮秒反斯托克斯拉曼频移器的理论研究

  • 1. 天津理工大学理学院, 天津 300384
  • 2. 中国电子科技集团公司第十八研究所, 天津 300384
  • 通信作者: 王聪, wangc.sd@163.com
    基金项目: 国家自然科学基金(批准号: 11504268)和天津市教委科研计划(批准号: 2014KJ240)资助的课题

摘要: 采用拉曼频移器在晶体介质中利用相干反斯托克斯散射效应可以获得超短脉冲(皮秒)反斯托克斯激光. 基于抽运-探测法的晶体拉曼频移器可以实现相干反斯托克斯散射的共线相互作用, 从而可以有效提高反斯托克斯光的转化效率. 本文在平面波近似下建立了基于抽运-探测法的皮秒反斯托克斯拉曼频移器的耦合波方程组, 引入归一化参量对方程组进行了归一化处理. 通过数值计算, 得到了描述皮秒反斯托克斯拉曼频移器运行的一组普适理论曲线, 分析了归一化拉曼增益系数G、归一化相位失配参数ΔK以及探测光脉冲能量占基频光总能量的比值rprobe三个变量对反斯托克斯拉曼频移器性能的影响, 确定了实现高效反斯托克斯转化时各归一化变量的合理取值. 采用实验数据对该理论模型的正确性进行了验证, 反斯托克斯转化效率的理论值与文献数据基本一致.

English Abstract

    • 随着晶体材料研究的不断深入, 基于受激斯托克斯拉曼散射(stimulated Stokes Raman scattering, SSRS)效应的晶体拉曼激光器作为扩展激光波长范围的重要手段之一, 已成为固体激光器领域研究的热点[1-10]. 利用相干反斯托克斯拉曼散射(coherent anti-Stokes Raman scattering, CARS)效应, 可以在拉曼晶体中实现频率上转换, 从而获得反斯托克斯光[11-19]. 它可以进一步扩大相干光谱范围, 产生具有重要应用价值的相干光. 例如, 晶体中的CARS效应能够将532 nm激光转化为在彩色全息、光学对抗、激光显示等方面有重要应用的蓝光输出. 另外, 皮秒激光以其超高的峰值功率在许多领域得到了广泛的应用. 因此, CARS效应与皮秒激光的结合为反斯托克斯激光的应用开辟了更多的可能性[12,14].

      拉曼频移器是产生超短脉冲反斯托克斯激光的有效方法. 2004年, Grasiuk等[14]采用两级拉曼频移器在KGd(WO4)2晶体中实现了532 nm激光的频率上转换, 当第一级晶体产生的一阶斯托克斯种子光与抽运光在第二级晶体中的传播方向满足相位匹配条件时, 获得了转化效率为4%的511 nm反斯托克斯光输出. 然而, 第二级晶体中的非共线相位匹配导致抽运光、一阶斯托克斯光和反斯托克斯光的光束不能完全重合, 从而限制了反斯托克斯光的转化效率. 2017年, Smetanin等[12]设计了一种新的实验方案, 采用抽运-探测法实现了CaCO3晶体中共线相位匹配的皮秒反斯托克斯拉曼频移器, 由探测光向反斯托克斯光的转化效率高达30%. 基频光光源发出的基频光由分光装置分为不同偏振态的抽运光和探测光, 经光学延时系统进行时间和空间同步后入射到拉曼晶体中. 当抽运光和探测光束以一定的角度入射到拉曼晶体中时, 可以实现抽运光、探测光、一阶斯托克斯光和反斯托克斯光的共线相位匹配. 抽运光通过SSRS效应产生一阶斯托克斯光, 反斯托克斯光则由探测光经CARS效应产生.

      理论模拟是研究激光运转的重要手段. Shen和Bloembergen[20]采用耦合波方程解释了受激拉曼散射中高阶斯托克斯和反斯托克斯光的产生. 此后, 耦合波方程被广泛用来研究拉曼激光器[8,9]和反斯托克斯激光器特性[11-13]. Smetanin等[12]采用物质方程和耦合波方程分析了探测光和抽运光光强三种占比情况下CaCO3皮秒反斯托克斯拉曼频移器的转化效率. 然而, 以往报道的反斯托克斯拉曼频移器理论虽能反映频移器的运转规律, 但未有报道研究频移器的最优化问题, 也未有报道给出频移器参量对反斯托克斯激光输出特性的影响. 本文采用耦合波理论对基于抽运-探测法的皮秒反斯托克斯拉曼频移器进行了理论研究. 考虑探测通道中一阶斯托克斯光和抽运通道中二阶斯托克斯光的产生, 在平面波近似下, 建立了皮秒反斯托克斯拉曼频移器的耦合波方程. 对耦合波方程进行了归一化处理和数值求解, 得到了一组反映归一化参数对反斯托克斯拉曼频移器性能影响的曲线, 分析了归一化参量对反斯托克斯光转化效率的影响. 本文提出的归一化耦合波理论有助于了解皮秒反斯托克斯拉曼频移器的运转规律, 而且对频移器的设计具有指导意义.

    • 考虑探测通道中一阶斯托克斯分量的产生和抽运通道中二阶斯托克斯分量的产生, 忽略高阶斯托克斯光和反斯托克斯光的产生, 各分量之间的相互作用原理如图1所示. 图1(a)为抽运通道SSRS过程的能级图. 频率为ωpump的抽运光入射到拉曼晶体中, 与物质分子相互作用, 一个抽运光光子转化成一个频率为ω1s的一阶斯托克斯光光子, 当一阶斯托克斯光强度大于二阶斯托克斯光阈值时, 作为抽运光产生频率为ω2s的二阶斯托克斯光. 探测通道的CARS过程是抽运光、探测光、一阶斯托克斯光和反斯托克斯光的四波混频过程, 如图1(b)所示. 在这个过程中, 产生一个频率为ωa的反斯托克斯光子的同时, 消耗一个一阶斯托克斯光子和一个频率为ωprobe的探测光光子, 产生一个抽运光光子. 当反斯托克斯光强度足够大时, 还可以作为抽运光通过SSRS向探测光转化, 如图1(c)所示.

      图  1  抽运通道和探测通道中的拉曼散射能级图 (a) 抽运通道的SSRS能级图; (b) 探测通道的CARS能级图; (c) 探测通道的SSRS能级图

      Figure 1.  Raman scattering energy levels in pump and probe channels: (a) SSRS in pump channel; (b) CARS in probe channel; (c) SSRS in probe channel.

    • 在平面波近似下, 基于抽运-探测法的反斯托克斯拉曼频移器的耦合波方程为

      $\begin{split} &\left( {\frac{\partial }{{\partial z}} + \frac{{{n_{\rm{a}}}}}{c}\frac{\partial }{{\partial t}}} \right){E_{\rm{a}}} \\ =\;& - \frac{g}{2}\frac{{{\nu _{\rm{a}}}}}{{{\nu _{\rm{p}}}}}\left( {{{\left| {{E_{{\rm{probe}}}}} \right|}^2}{E_{\rm{a}}} + {E_{{\rm{probe}}}}{E_{{\rm{pump}}}}E_{1{\rm{s}}}^*{{\rm{e}}^{{\rm{i}}\Delta kz}}} \right), \end{split} \tag{1a}$

      $\begin{split} &\left( {\frac{\partial }{{\partial z}} + \frac{{{n_{{\rm{pump}}}}}}{c}\frac{\partial }{{\partial t}}} \right){E_{{\rm{pump}}}} \\ =\;& - \frac{g}{2}\left( {{{\left| {{E_{1{\rm{s}}}}} \right|}^2}{E_{{\rm{pump}}}} + {E_{1{\rm{s}}}}{E_{\rm{a}}}E_{{\rm{probe}}}^*{{\rm{e}}^{ - {\rm{i}}\Delta kz}}} \right), \end{split}\tag{1b}$

      $\begin{split} &\left( {\frac{\partial }{{\partial z}} + \frac{{{n_{{\rm{probe}}}}}}{c}\frac{\partial }{{\partial t}}} \right){E_{{\rm{probe}}}}\\ =\;& \frac{g}{2}\left( {{{\left| {{E_{\rm{a}}}} \right|}^2}{E_{{\rm{probe}}}} + {E_{1{\rm{s}}}}{E_{\rm{a}}}E_{{\rm{pump}}}^*{{\rm{e}}^{ - {\rm{i}}\Delta kz}}} \right), \end{split}\tag{1c}$

      $\begin{split} \left( {\frac{\partial }{{\partial z}} + \frac{{{n_{{\rm{1s}}}}}}{c}\frac{\partial }{{\partial t}}} \right){E_{{\rm{1s}}}}=\;& \frac{g}{2}\frac{{{\nu _{{\rm{1s}}}}}}{{{\nu _{\rm{p}}}}}\Big[ \left( {{{\left| {{E_{{\rm{pump}}}}} \right|}^2} - {{\left| {{E_{{\rm{2s}}}}} \right|}^2}} \right){E_{{\rm{1s}}}} \\ &+ {E_{{\rm{probe}}}}{E_{{\rm{pump}}}}E_{\rm{a}}^*{{\rm{e}}^{{\rm{i}}\Delta kz}} \Big],\\[-12pt] \end{split}\tag{1d}$

      $\left( {\frac{\partial }{{\partial z}} + \frac{{{n_{2{\rm{s}}}}}}{c}\frac{\partial }{{\partial t}}} \right){E_{2{\rm{s}}}} = \frac{g}{2}\frac{{{\nu _{{\rm{2s}}}}}}{{{\nu _{\rm{p}}}}}{\left| {{E_{1{\rm{s}}}}} \right|^2}{E_{{\rm{2s}}}}, \tag{1e}$

      式中, Ej (j = a, pump, probe, 1s, 2s)分别为沿z轴传播的反斯托克斯光、抽运光、探测光、一阶斯托克斯光和二阶斯托克斯光的缓变振幅, nj为拉曼晶体中各辐射分量的折射率, c为真空中的光速, νm (m = p, a, 1s, 2s)分别为基频光、反斯托克斯光、一阶斯托克斯光和二阶斯托克斯光的频率, g为拉曼晶体对基频光的拉曼增益系数, Δk = k1skpumpkprobe + ka为四波混频的相位失配参量. Ejtz的函数, 即Ej = Ej (t, z), 为方便起见, 这里采用简化形式.

      抽运光和探测光单程通过拉曼晶体, (1)式的初始条件为

      $\begin{split} &{E_{{\rm{probe}}}}\left( {t,0} \right) = \sqrt {{r_{{\rm{probe}}}}} {E_{\rm{p}}}\left( t \right),\;\\ &{E_{{\rm{pump}}}}\left( {t,0} \right) = \sqrt {1 - {r_{{\rm{probe}}}}} {E_{\rm{p}}}\left( t \right),\end{split} $

      式中, Ep(t)为从光源发出的基频光的缓变振幅, rprobe为探测光能量占基频光总能量的比例.

      为使耦合波方程具有一般性, 引入归一化空间坐标ζ、归一化时间τ、归一化缓变振幅Φj、归一化拉曼增益系数G和归一化相位失配参量$ \Delta K $:

      $\begin{split} &\zeta = \frac{z}{{{l_{\rm{R}}}}}, \;\tau = \frac{t}{{{t_{\rm{R}}}}},\; {\varPhi _j} = \frac{{{E_j}}}{{{E_{{\rm{pmax}}}}}},\\ &G = \frac{g}{2}{\left| {{E_{{\rm{p}}\max }}} \right|^2}{l_{\rm{R}}},\; \Delta K = \Delta k{l_{\rm{R}}},\end{split} $

      式中, lR为拉曼晶体的长度, tR为光在拉曼晶体中单程通过所需要的时间, Epmax为光源产生的基频光的最大振幅.

      假设基频光脉冲的强度在时间上为高斯分布, 脉冲宽度为wp, 则基频光的归一化振幅Φp与归一化时间τ之间的关系可表示为

      ${\varPhi _{\rm{p}}}\left( \tau \right) = \exp \left( { - \frac{{\tau - {\tau _{{\rm{pm}}}}}}{{{W_{\rm{p}}}/2\sqrt {\ln (2)} }}} \right)\exp \left( {{\rm{i}}{\varphi _{\rm{p}}}} \right), $

      式中, φp为随机相位, τpm为脉冲峰值对应的归一化时间, Wp = wp/tR为基频光的归一化脉冲宽度. 若基频光脉冲宽度wp = 20 ps, 拉曼晶体长度lR = 1.5 cm, 拉曼晶体折射率n = 2, 则估算出的Wp = 0.2.

      将(3)式代入(1)式和(2)式中, 可以得到归一化耦合波方程组为

      $\begin{split} &\left( {\frac{\partial }{{\partial \zeta }} + \frac{\partial }{{\partial \tau }}} \right){\varPhi _{\rm{a}}} \\ =\;& - \frac{{{\nu _{\rm{a}}}}}{{{\nu _{\rm{p}}}}}G\left( {{{\left| {{\varPhi _{{\rm{probe}}}}} \right|}^2}{\varPhi _{\rm{a}}} + {\varPhi _{{\rm{probe}}}}{\varPhi _{{\rm{pump}}}}\varPhi _{1{\rm{s}}}^*{{\rm{e}}^{{\rm{i}}\Delta K\varsigma }}} \right), \end{split}\tag{5a}$

      $\begin{split} &\left( {\frac{\partial }{{\partial \zeta }} + \frac{\partial }{{\partial \tau }}} \right){\varPhi _{{\rm{pump}}}} \\ =\;& - G\left( {{{\left| {{\varPhi _{1{\rm{s}}}}} \right|}^2}{\varPhi _{{\rm{pump}}}} + {\varPhi _{1{\rm{s}}}}{\varPhi _{\rm{a}}}\varPhi _{{\rm{probe}}}^*{{\rm{e}}^{ - {\rm{i}}\Delta K\varsigma }}} \right), \end{split}\tag{5b}$

      $\begin{split} &\left( {\frac{\partial }{{\partial \zeta }} + \frac{\partial }{{\partial \tau }}} \right){\varPhi _{{\rm{probe}}}} \\ =\;& G\left( {{{\left| {{\varPhi _{\rm{a}}}} \right|}^2}{\varPhi _{{\rm{probe}}}} + {\varPhi _{1{\rm{s}}}}{\varPhi _{\rm{a}}}\varPhi _{{\rm{pump}}}^*{{\rm{e}}^{ - {\rm{i}}\Delta K\varsigma }}} \right),\quad \end{split} \tag{5c} $

      $\begin{split} \left( {\frac{\partial }{{\partial \zeta }} + \frac{\partial }{{\partial \tau }}} \right){\varPhi _{{\rm{1s}}}} =\;& \frac{{{\nu _{1{\rm{s}}}}}}{{{\nu _{\rm{p}}}}}G\Big[ \left( {{{\left| {{\varPhi _{{\rm{pump}}}}} \right|}^2} - {{\left| {{\varPhi _{{\rm{2s}}}}} \right|}^2}} \right){\varPhi _{{\rm{1s}}}}\\ &+ {\varPhi _{{\rm{probe}}}}{\varPhi _{{\rm{pump}}}}\varPhi _{\rm{a}}^{\rm{*}}{{\rm{e}}^{{\rm{i}}\Delta K\varsigma }} \Big],\\[-13pt] \end{split}\tag{5d}$

      $\left( {\frac{\partial }{{\partial \zeta }} + \frac{\partial }{{\partial \tau }}} \right){\varPhi _{2{\rm{s}}}} = \frac{{{\nu _{2{\rm{s}}}}}}{{{\nu _{\rm{p}}}}}G{\left| {{\varPhi _{{\rm{1s}}}}} \right|^2}{\varPhi _{2{\rm{s}}}}. \tag{5e}$

      归一化初始条件为

      $\begin{split} &{\varPhi _{{\rm{probe}}}}\left( {0,\tau } \right) = \sqrt {{r_{{\rm{probe}}}}} {\varPhi _{\rm{p}}}\left( \tau \right),\;\\ &{\varPhi _{{\rm{pump}}}}\left( {0,\tau } \right) = \sqrt {1 - {r_{{\rm{probe}}}}} {\varPhi _{\rm{p}}}\left( \tau \right). \end{split}$

      运用初始条件(6)式对(5)式进行数值求解, 可以得到出射分量j的归一化振幅Φj(1, τ), 则各分量出射光的单脉冲能量为

      $\begin{split} {e_j} =\;& {A_j}{\int {\left| {{E_j}\left( {{l_{\rm{R}}},t} \right)} \right|} ^2}{\rm{d}}t \\ =\;& \frac{{{A_j}{l_{\rm{R}}}{n_j}}}{c}{\left| {{E_{{\rm{p}}\max }}} \right|^2}\int {{{\left| {{\varPhi _j}\left( {1,\tau } \right)} \right|}^2}} {\rm{d}}\tau ,\end{split} $

      式中, Aj为分量j的光束面积.

      各拉曼分量的转化效率定义为输出脉冲能量与入射基频光脉冲能量之比

      $\frac{{{e_j}}}{{{e_{\rm{p}}}}} = \frac{{{A_j}\displaystyle\int {{{\left| {{\varPhi _j}\left( {1,\tau } \right)} \right|}^2}{\rm{d}}\tau } }}{{{A_{\rm{p}}}\displaystyle\int {{{\left| {{\varPhi _{\rm{p}}}\left( \tau \right)} \right|}^2}{\rm{d}}\tau } }}, $

      式中, Ap为基频光的光束面积. 拉曼分量j的归一化转化效率定义为

      ${\eta _j} = \frac{{\displaystyle\int {{{\left| {{\varPhi _j}\left( {1,\tau } \right)} \right|}^2}{\rm{d}}\tau } }}{{\displaystyle\int {{{\left| {{\varPhi _{\rm{p}}}\left( \tau \right)} \right|}^2}{\rm{d}}\tau } }}. $

    • 相位失配参量是影响四波混频强度的一个重要因素. 图2给出了不同rprobe时反斯托克斯光、一阶斯托克斯光和二阶斯托克斯光的归一化转化效率(ηa, η1sη2s)随归一化相位失配参量ΔK的变化关系, 其中G = 90, Wp = 0.2. 当rprobe较小时(rprobe < 0.35), 抽运光的强度大于二阶斯托克斯光的阈值, 因此二阶斯托克斯光具有较大的转化效率. ηa在相位匹配(ΔK = 0)时有最大值, 且η1sη2s随ΔK基本不变. 随着rprobe的增大, 二阶斯托克斯转化减弱, 反斯托克斯光的产生随一阶斯托克斯转化的增加而增加. 当rprobe约为0.35时, ηaη1s达到最大值, η1sη2s的曲线中心出现凹陷. 抽运光强度低于二阶斯托克斯光阈值 (rprobe > 0.35) 时, 二阶斯托克斯光消失, ηaη1srprobe的增大而减小. 当rprobe > 0.37时, ηa的最大值偏离ΔK = 0处, 且ΔK = 0处的ηaη1srprobe的增加迅速减小.

      图  2  G = 90, Wp = 0.2, rprobe取不同值时(a) ηa, (b) η1s和(c) η2s随ΔK的变化

      Figure 2.  (a) ηa, (b) η1s and (c) η2s versus ΔK for different rprobe with G = 90 and Wp = 0.2.

      图2中曲线的变化规律可以用图3图4给出的抽运光、探测光、一阶斯托克斯光和反斯托克斯光的归一化强度在拉曼晶体中的空间演化进行解释, 其中图3rprobe = 0.3, 图4rprobe = 0.39. 对于每一组曲线, |ΔK| = 0, 4和8, G = 90, Wp = 0.2. 根据图5(a)所示的脉冲形状, 当rprobe = 0.3时, 二阶斯托克斯光的产生消耗了入射抽运光和探测光脉冲峰值附近大部分的一阶斯托克斯光, 这导致了峰值附近的反斯托克斯输出很弱. 因此, 将图3中基频光的归一化强度(|Φp|2)设为0.8, 从而偏离了二阶斯托克斯光产生的区域. 相反, 如图5(b)所示, 对于rprobe = 0.39, 没有二阶斯托克斯光产生, 一阶斯托克斯光和反斯托克斯光均产生于入射脉冲峰值附近. 因此图4中, 令|Φp|2等于基频光归一化强度的最大值, 即1.

      图  3  rprobe = 0.3, G = 90, Wp = 0.2时, 抽运光、探测光、一阶斯托克斯光和反斯托克斯光归一化光强随ζ的空间演化 (a) |ΔK| = 0; (b) |ΔK| = 4; (c) |ΔK| = 8

      Figure 3.  Plots of the spatial evolution of pump, probe, first Stokes, and anti-Stokes normalized intensities with rprobe = 0.3, G = 90 and Wp = 0.2: (a) |ΔK| = 0; (b) |ΔK| = 4; and (c) |ΔK| = 8.

      图  4  rprobe = 0.39, G = 90, Wp = 0.2时, 抽运光、探测光、一阶斯托克斯光和反斯托克斯光归一化光强随ζ的空间演化 (a) |ΔK| = 0; (b) |ΔK| = 4; (c) |ΔK| = 8

      Figure 4.  Plots of the spatial evolution of pump, probe, first Stokes, and anti-Stokes normalized intensities with rprobe = 0.39, G = 90 and Wp = 0.2: (a) |ΔK| = 0; (b) |ΔK| = 4; and (c) |ΔK| = 8.

      图  5  G = 90, |ΔK| = 0, Wp = 0.2时抽运光、探测光、一阶斯托克斯光、二阶斯托克斯光和反斯托克斯光的脉冲形状 (a) rprobe = 0.3; (b) rprobe = 0.39

      Figure 5.  Temporal profiles of the pump, probe, first Stokes, second Stokes and anti-Stokes pulses with G = 90, |ΔK| = 0 and Wp = 0.2: (a) rprobe = 0.3; (b) rprobe = 0.39.

      图3图4中, 拉曼晶体中各分量归一化强度的空间演化分为以ζ1ζ2为分界点的3个阶段. 第1个阶段(ζ < ζ1)是受激拉曼散射的积累阶段, 散射光在ζ1点达到一定强度时使抽运光迅速向一阶斯托克斯光转化. 与此同时, 产生的一阶斯托克斯光参与到CARS中, 反斯托克斯光强度在ζ1点也快速增长. 在ζ1点前, |ΔK|越小, CARS效应越强, 造成的斯托克斯散射损耗也就越大, 则ζ1越大; 并且rprobe越大, 探测光越强, 反斯托克斯衰减越大, ζ1也越大. 在第2个阶段(ζ1 < ζ < ζ2), 四波相互作用迅速增强, 抽运光通过SSRS转化为一阶斯托克斯光, 探测光通过CARS转化为反斯托克斯光, 一阶斯托克斯光和反斯托克斯光强在ζ2处达到最大值. |ΔK|越小, 相互作用越强, 导致ζ2处的一阶斯托克斯光和反斯托克斯光强度越大, 剩余抽运光和探测光强度越小. 在第3个阶段 (ζ > ζ2), 反斯托克斯光的强度大于SSRS的阈值, 通过SSRS向探测光转化. 当|ΔK|较大时, 由于ζ2处作为初始散射光的探测光强度大, 反斯托克斯光强度下降更快. 在图3所示的情况下, ζ2均小于1, 因此, 输出的一阶斯托克斯光和反斯托克斯光的强度随|ΔK|的增大而减小. 然而, 当rprobe较大时, 在相位匹配时ζ2大于1 (图4(a)), 因此输出的一阶斯托克斯光和反斯托克斯光的强度最低.

      由以上分析可知, 在某些情况下, 相位匹配时反斯托克斯光的转化效率并不是最高的, 但只要rprobe在合理的范围内, |ΔK| = 0仍然是获得高反斯托克斯转化效率的必要条件. 因此, 在以下的计算中, 令|ΔK|的取值为零.

      图6给出了不同rprobeηa, η1sη2sG的变化关系, 其中|ΔK| = 0, Wp = 0.2. 从图6(a)图6(b)可看出, ηaη1sG的变化规律几乎相同, 这是由于反斯托克斯光产生的前提是一阶斯托克斯光的存在, 如图1(b)所示. 定义Gath, G1sthG2sth分别为ηa, η1sη2s的增益阈值, GaoptG1sopt分别为ηaη1s最大值对应的最佳增益值, 则从图6可看出, GathG1sth, GaoptG1sopt, G1soptG2sth. 可以这样定性地解释: 开始阶段, G较小, 抽运光在拉曼晶体中受激拉曼散射增益很弱, ζ1 > 1, 没有散射光输出. G > Gath(G1sth)时, 抽运光通过SSRS开始向一阶斯托克斯光转化, 与此同时探测光通过CARS开始向反斯托克斯光转化, ηaη1s迅速增长并分别在GaoptG1sopt处达到最大值. 一阶斯托斯克光达到一定强度后进而向二阶斯托克斯光转化, ηaη1sG的增大迅速减小, η2s则随之单调增大, 如图6(c)所示. rprobe越小, 抽运光越强, 各阶散射分量的增益阈值和最佳增益值则越小.

      图  6  |ΔK| = 0, Wp = 0.2, rprobe取不同值时(a) ηa, (b) η1 s和(c) η2 sG的变化

      Figure 6.  (a) ηa, (b) η1s and (c) η2s versus G for different rprobe with |ΔK| = 0 and Wp = 0.2.

      图2图6可知, 当G一定时, 反斯托克斯光的转化效率在一个最佳的rprobe(ropt)时有一最大值, 反之亦是如此. |ΔK| = 0, Wp = 0.2时, 反斯托克斯光最大归一化转化效率ηamax和对应的roptG的变化见图7. roptG的增大单调增长, 这与图6中得到的结果一致. G大于增益阈值后, ηamaxG的增大先迅速增长, 在G = 110, ropt = 0.373时有最大值0.236, 随后缓慢下降. 这是由于反斯托克斯光是由探测光转化而来, 当G < 110时, 反斯托克斯光在低增益时几乎不向探测光转化, rprobe越大, 探测光越强, 则反斯托克斯光越强, 如图8(a)图8(b)所示. 然而, 如图8(c)所示, 当G > 110时, G增大的同时ropt也增大, G的增大使反斯托克斯光向探测光转化增强, ropt的增大使探测光向反斯托克斯光转化增强, 但前者增量大于后者, 因此ηamaxG的增大有所下降.

      图  7  |ΔK| = 0, Wp = 0.2时, roptηamaxG的变化

      Figure 7.  ropt and ηamax versus G with |ΔK| = 0 and Wp = 0.2

      图  8  |ΔK| = 0, Wp = 0.2时抽运光、探测光、一阶斯托克斯光、二阶斯托克斯光和反斯托克斯光的脉冲形状 (a) ropt = 0.270, G = 60; (b) ropt = 0.373, G = 110; (c) ropt = 0.414, G = 160

      Figure 8.  Temporal profiles of the pump, probe, first Stokes, second Stokes, and anti-Stokes pulses with |ΔK |= 0 and Wp = 0.2: (a) ropt = 0.270, G = 60; (b) ropt = 0.373, G = 110; (c) ropt = 0.414, G = 160.

      rprobe一定时, 反斯托克斯光的转化效率在一个最佳的G(Gopt)时有一最大值. |ΔK| = 0, Wp = 0.2时, 反斯托克斯光的最大归一化转化效率ηamax和对应的Goptrprobe的变化曲线如图9所示. 与图7的分析和结论相同, 由于反斯托克斯光是由探测光转化而来, ηamaxrprobe的增大首先增长, 在Gopt = 110, rprobe = 0.373时达到最大值. 当rprobe > 0.373时, rprobe的增大使反斯托克斯光向探测光的转化多于探测光向反斯托克斯光的转化, 因此ηamaxrprobe的增大而下降.

      图  9  |ΔK| = 0, Wp = 0.2时, Goptηamaxrprobe的变化

      Figure 9.  Gopt and ηamax versus rprobe with |ΔK| = 0 and Wp = 0.2

      下面分析Goptrprobe的变化规律. 当rprobe < 0.05时, Goptrprobe的增大而减小, 这是因为当rprobe很小时, 抽运光很强, 大部分的抽运光通过SSRS效应转化成了二阶斯托克斯光, 参与CARS效应的一阶斯托克斯光很弱. Gopt的减小一方面减弱了抽运通道中的SSRS效应, 使得探测通道中的CARS效应增强, 另一方面减弱了探测通道中反斯托克斯光向探测光的转化. 当rprobe > 0.05时, 随着rprobe的增大, 抽运光强度逐渐减弱, Goptrprobe的增大而增大, 增强了抽运通道的SSRS效应以增大一阶斯托克斯光的强度, 使得CARS效应最强而获得最大的反斯托克斯转换效率.

    • 已知基频光脉冲能量ep、基频光脉冲宽度wp、基频光光束面积Ap、拉曼增益系数g以及拉曼晶体的长度lR, 可以估算出基频光脉冲的峰值振幅为${E_{{\rm{p}}\max }} = \sqrt {{e_{\rm{p}}}/\left( {{A_{\rm{p}}}{w_{\rm{p}}}} \right)} $、归一化拉曼增益系数为$G = $$ \dfrac{1}{2}g{\left| {{E_{{\rm{p}}\max }}} \right|^2}{l_{\rm{R}}}$以及归一化基频光脉冲宽度为Wp = wp /tR. 将归一化参量值和探测光所占的能量比rprobe代入归一化耦合波方程组中进行数值求解, 可以计算反斯托克斯光的转化效率.

      如引言所述, 文献[12]采用抽运-探测法在CaCO3晶体中实现了共线相位匹配的皮秒反斯托克斯拉曼频移器. 下面采用文献[12]中的实验结果对本文理论的正确性进行验证. 表1总结了文献[12]中的部分实验参数. 表2为反斯托克斯转化效率的理论值与实验数据的对比结果, 理论结果与实验结果基本一致. 误差产生的原因有: 1) 本文的耦合波理论基于平面波近似, 实际入射的抽运光和探测光的光强在横截面上近似为高斯分布, 且由于透镜的聚焦, 光束有一定的发散角, 在晶体不同位置处光束半径不同; 2) 理论上假设入射抽运光和探测光强度在时间上为高斯分布, 与实际光源产生的脉冲形状有所差异.

      参数参数
      νa/νp1.06Ap/cm22.83 × 10–4
      ν1s/νp0.94lR/cm3.2
      ν2s/νp0.88g/(cm·GW–1)13
      wp/ps20Δk0

      表 1  参考文献[12]中的参数

      Table 1.  Parameters in Ref. [12].

      抽运脉冲
      能量/μJ
      探测脉冲
      能量/μJ
      Grprobe实验
      值/%
      理论
      值/%
      30252060.452.802.54
      3021180.0631.881.86
      26121400.323.504.13

      表 2  不同情况下反斯托克斯转化效率的理论值与实验数据的对比结果

      Table 2.  Comparisons of theoretical and experimental results of anti-Stokes conversion efficiency under different conditions.

      为进一步验证本文理论的正确性, 下面考虑一种更接近实际的情况. 基频光光强在横截面上近似为高斯分布, 且在传播方向上各处的光束半径R均相等. 因此, 基频光的缓变振幅$ E _{\rm{p}}^{\rm{g}}$是径向坐标rt的函数,

      $E_{\rm{p}}^{\rm{g}}(r,t) = E_{\rm{p}}^{\rm{g}}(0,t)\exp\left( { - \frac{{{r^2}}}{{{R^2}}}} \right).$

      为了与平面波近似时得到的结论有可比性, 即平面波近似和高斯近似时基频光具有相同的脉冲能量和光束面积, 高斯近似时基频光在光轴(r = 0)处的峰值振幅$ E _{\rm{pmax}}^{\rm{g}}$(0, t)为平面近似时基频光最大振幅Epmax$\sqrt 2 $[8]. 因此, 高斯近似时基频光的归一化振幅$ \varPhi _{\rm{p}}^{\rm{g}}$(r, τ)表示为

      $\begin{split} \varPhi _{\rm{p}}^{\rm{g}}\left( {r,\tau } \right) =\;& \sqrt 2 \exp\left( { - \frac{{{r^2}}}{{{R^2}}}} \right)\\ &\times\exp \left( { - \frac{{\tau - {\tau _{{\rm{pm}}}}}}{{{W_{\rm{p}}}/2\sqrt {\ln (2)} }}} \right)\exp \left( {{\rm{i}}{\varphi _{\rm{p}}}} \right).\end{split} $

      当考虑光强的横向分布时, 各拉曼分量的归一化振幅$ \varPhi _j^{\rm{g}}$r, ζτ的函数, 则(5)式的归一化初始条件为

      $\begin{split} &\varPhi _{{\rm{probe}}}^{\rm{g}}\left( {r,0,\tau } \right) = \sqrt {{r_{{\rm{probe}}}}} \varPhi _{\rm{p}}^{\rm{g}}\left( {r,\tau } \right),\;\\ &\varPhi _{{\rm{pump}}}^{\rm{g}}\left( {r,0,\tau } \right) = \sqrt {1 - {r_{{\rm{probe}}}}} \varPhi _{\rm{p}}^{\rm{g}}\left( {r,\tau } \right). \end{split}$

      由(3)式和(10)式可以看出高斯近似时归一化拉曼增益系数$G^{\rm g}(r) $与平面波近似时归一化拉曼增益系数G的关系为

      $\begin{split} &{G^{\rm{g}}}\left( r \right) = \frac{g}{2}{\left| {E{{_{\rm{p}}^{\rm{g}}}_{\max }}\left( r \right)} \right|^2}{l_{\rm{R}}}\\ =\;& \frac{g}{2}{\left| {E{{_{\rm{p}}^{\rm{g}}}_{\max }}\left( 0 \right)\exp \left( { - \frac{{{r^2}}}{{{R^2}}}} \right)} \right|^2}{l_{\rm{R}}} \\ =\;& 2G\exp \left( { - \frac{{2{r^2}}}{{{R^2}}}} \right).\end{split} $

      在时间上, 实际激光脉冲的上升沿比下降沿稍陡. 为使基频光脉冲形状更接近实际情况同时便于计算, 令$ \varPhi _{\rm{p}}^{\rm{g}}$(r, τ)的上升沿和下降沿均为高斯函数. 上升沿的宽度(由$ \varPhi _{\rm{p}}^{\rm{g}}$(r, τpm)/2到$ \varPhi _{\rm{p}}^{\rm{g}}$(r, τpm)的归一化时间)与下降沿的宽度(由$ \varPhi _{\rm{p}}^{\rm{g}}$(r, τpm)到$ \varPhi _{\rm{p}}^{\rm{g}}$(r, τpm) /2的归一化时间)的比值为0.8.

      运用(12)式的初始条件对归一化耦合波方程组(5)进行数值求解, 可以得到出射分量j的归一化振幅$ \varPhi _j^{\rm{g}}$(r, 1, τ). 由(8)式, 高斯近似时反斯托克斯光的转化效率$ \eta _{\rm{a}}^{\rm{g}}$

      $\eta _{\rm{a}}^{\rm{g}} = \frac{{\displaystyle\iint {{{\left| {\varPhi _{\rm{a}}^{\rm{g}}\left( {r,1,\tau } \right)} \right|}^2}{\rm{2\pi }}r{\rm{d}}r{\rm{d}}\tau }}}{{\displaystyle\iint {{{\left| {\varPhi _{\rm{p}}^{\rm{g}}\left( {r,\tau } \right)} \right|}^2}{\rm{2\pi }}r{\rm{d}}r{\rm{d}}\tau }}}. $

      在上述高斯近似条件下, 采用表1的实验参数进行计算得出的反斯托克斯转化效率的理论值如表3所列, 与表2中的结果相比可以发现, 两种近似条件下的理论值均与实验结果相符合, 同时也证明了平面波近似时耦合波理论可以正确地反映反斯托克斯拉曼频移器的运转特性.

      抽运脉冲能量/μJ探测脉冲能量/μJ$ \eta _{\rm{a}}^{\rm{g}}$/%
      30252.63
      3021.91
      26124.09

      表 3  高斯近似时反斯托克斯转化效率的理论值

      Table 3.  Theoretical values of anti-Stokes conversion efficiency for Gaussian approximation.

    • 本文在理论上研究了基于抽运-探测法的皮秒反斯托克斯拉曼频移器, 建立了平面波近似下的耦合波方程组. 为使方程组具有普适性, 引入4个无量纲综合参量对方程组进行了归一化. 对方程组数值求解显示该拉曼频移器的性能主要依赖于3个参量: 归一化相位失配参量ΔK、归一化拉曼增益系数G及探测光与基频光的能量比rprobe. 在以往的报道中, 相位匹配(|ΔK| = 0)是皮秒反斯托克斯拉曼频移器的搭建原则[12,14]. 本文通过分析归一化相位失配参量对反斯托克斯转化效率的影响发现, 虽然在相位匹配条件下可以获得最大的反斯托克斯输出, 但前提是当G一定时rprobe在合理的范围之内. 在实际中, 较容易改变的参量为基频光的峰值光强(|Epmax|2)和rprobe, 对于确定的拉曼晶体(glR一定), |Epmax|2决定了G的大小. 因此, 共线相位匹配时, 需选择合适的基频光能量和探测光比例才可以获得高效的反斯托克斯光输出. 需要说明的是, 由于理论上抽运光和探测光的脉冲宽度相同且同步单程通过拉曼晶体, 当其他条件一定时, 基频光的脉冲宽度对反斯托克斯光的转化效率几乎没有影响. 本文提出的归一化耦合波理论可以作为分析基于抽运-探测法的共线反斯托克斯拉曼频移器的理论工具, 辅助激光器的设计以实现超短脉冲(皮秒)反斯托克斯光的最大转化效率.

参考文献 (20)

目录

    /

    返回文章
    返回