高级检索

ISSN1001-3806CN51-1125/TN 网站地图

留言板

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

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

飞秒激光致双液滴光学击穿和等离子体分布研究

钱天 陆健 唐懋 张冲 张宏超

引用本文:
Citation:

飞秒激光致双液滴光学击穿和等离子体分布研究

    作者简介: 钱天(1998-),男,硕士研究生,现主要从事激光与液滴相互作用的研究.
    通讯作者: 陆健, lujian@njust.edu.cn
  • 基金项目:

    国家自然科学基金资助项目 11774176

  • 中图分类号: O437.5

Optical breakdown and plasma morphology distribution of double droplets induced by femtosecond laser

    Corresponding author: LU Jian, lujian@njust.edu.cn ;
  • CLC number: O437.5

  • 摘要: 为了研究超短激光脉冲与双液滴相互作用过程中的光学击穿和等离子体分布,基于麦克斯韦方程组和电离速率方程,构建了飞秒激光与双液滴的瞬态耦合模型,使用有限元分析方法,对飞秒激光辐照微米量级双液滴的自由电子密度和光场分布进行了计算,得到了双液滴结构对液滴光学击穿和等离子体变化的影响。结果表明,第2个液滴的击穿阈值约为同等条件下单液滴击穿阈值的35%;等离子体的形态和击穿点的位置随双液滴间距发生变化,且在聚焦区域产生纳米等离子体射流;第2个液滴对激光能量的吸收随着双液滴间距的增加而减少;当分别使用满足击穿阈值的光强入射,双液滴吸收的能量约为单液滴的3%;第2个液滴对激光能量的吸收随光强增大而增大,能量吸收比例最终趋于0.01,仅为单液滴的1.5%。该研究为激光诱导水击穿和激光在大气中的传输提供了一定的参考。
  • 图 1  飞秒激光与双液滴相互作用的模型示意图(xy平面模型)

    Figure 1.  Sketch of femtosecond laser interaction model with double droplets(xy plane model)

    图 2  双液滴与单液滴击穿阈值的比值随半径r的变化

    Figure 2.  Variation of the double-droplets to singledroplet breakdown threshold ratio with droplet radius r

    图 3  不同m值下双液滴等离子体形态的分布变化

    a—m=-0.5  b—m=-0.125  c—m=0.25  d~f—对应m值的激光光场分布

    Figure 3.  Distribution of plasma morphology in double droplets under different m values

    a—m=-0.5  b—m=-0.125  c—m=0.25  d~f—laser field distribution corresponding to m value

    图 4  双液滴最大自由电子密度和总电子数随m值的变化

    Figure 4.  Variations of maximum free electron density and total electron number of double droplets with m value

    图 5  双液滴击穿深度随m值的变化

    Figure 5.  Breakdown depth of double droplets varies with m values

    图 6  双液滴的能量损耗随m值的变化

    Figure 6.  Variations of energy losses of double droplets with m values

    图 7  不同m值下双液滴能量损耗Q和能量损耗比例Q/Qt随激光光强的变化

    Figure 7.  Variation of double droplets energy losses Q and proportion of energy losses Q/Qt with laser intensity at different m values

  • [1] 李蒙, 刘文荣. 雾天多次散射对激光透射仪能见度测量的影响[J]. 激光技术, 2020, 44(4): 503-508.

    LI M, LIU W R. Multiple-scattering effects on the visibility measurement of laser transmissometers in fog[J]. Laser Technoloy, 2020, 44(4): 503-508 (in Chinese). 
    [2]

    RUDENK A, ROSENOW P, HASSON V, et al. Plasma-free water droplet shattering by long-wave infrared ultrashort pulses for efficient fog clearing[J]. Optica, 2020, 7(2): 115-122.
    [3] 陈子琪, 王新兵, 左都罗. CO2激光诱导液滴射流等离子体的实验研究[J]. 激光技术, 2016, 40(6): 888-891.

    CHEN Z Q, WANG X B, ZUO D L. Experimental research of CO2 laser-induced liquid droplet jet flow plasma[J]. Laser Technology, 2016, 40(6): 888-891 (in Chinese). 
    [4]

    KLEIN A L, BOUWHUIS W, VISSER C W, et al. Drop shaping by laser-pulse impact[J]. Physical Review Applied, 2015, 3(4): 044018. doi: 10.1103/PhysRevApplied.3.044018
    [5]

    GELDERBLOM H, LHUISSIER H, KLEIN A L, et al. Drop defor-mation by laser-pulse impact[J]. Journal of Fluid Mechanics, 2016, 794: 676-699. doi: 10.1017/jfm.2016.182
    [6]

    HENIN S, PETIT Y, ROHWETTER P, et al. Field measurements suggest the mechanism of laser-assisted water condensation[J]. Nature Communications, 2011, 2(1): 1-7.
    [7]

    JU J J, LIU J Sh, WANG Ch, et al. Laser-filamentation-induced condensation and snow formation in a cloud chamber[J]. Optics Le-tters, 2012, 37(7): 1214-1216. doi: 10.1364/OL.37.001214
    [8]

    GEINTS Y E, KABANOV A M, MATVIENKO G G, et al. Broadband emission spectrum dynamics of large water droplets exposed to intense ultrashort laser radiation[J]. Optics Letters, 2010, 35(16): 2717-2719. doi: 10.1364/OL.35.002717
    [9] 崔祖文, 向昱霖, 张宇璇, 等. 铀酰溶液中铀含量的激光诱导击穿光谱测量[J]. 激光技术, 2021, 45(3): 331-335.

    CUI Z W, XIANG Y L, ZHANG Y X, et al. Quantitative study on uranium in uranyl solution by laser-induced breakdown spectroscopy[J]. Laser Technology, 2021, 45(3): 331-335 (in Chinese). 
    [10]

    LATZ C, ASSHAUER T, RATHJEN C, et al. Femtosecond-laser assisted surgery of the eye: Overview and impact of the low-energy concept[J]. Micromachines, 2021, 12(2): 122. doi: 10.3390/mi12020122
    [11]

    INHESTER L, HANASAKI K, HAO Y J, et al. X-ray multiphoton ionization dynamics of a water molecule irradiated by an X-ray free-electron laser pulse[J]. Physical Review, 2016, A94(2): 023422.
    [12]

    WOODBURY D, GOFFIN A, SCHWARTZ R M, et al. Self-guiding of long-wave infrared laser pulses mediated by avalanche ionization[J]. Physical Review Letters, 2020, 125(13): 133201. doi: 10.1103/PhysRevLett.125.133201
    [13]

    GEINTS Y E, ZEMLYANOV A A. Phase explosion of a water drop by a femtosecond laser pulse: I. Dynamics of optical breakdown[J]. Atmospheric and Oceanic Optics, 2009, 22(6): 581-589. doi: 10.1134/S1024856009060037
    [14]

    EFIMENKO E S, MALKOV Y A, MURZANEV A A, et al. Femtosecond laser pulse-induced breakdown of a single water microdroplet[J]. Journal of the Optical Society of America, 2014, B31(3): 534-541.
    [15]

    RUDENKO A, MOLONEY J V. Tunable near-to far-infrared optical breakdown in nonlinear interactions of ultrashort laser pulses with water microdroplets in ambient air[J]. Advanced Photonics Research, 2020, 1(2): 2000029. doi: 10.1002/adpr.202000029
    [16]

    BOULAIS E, LACHAINE R, MEUNIER M, et al. Plasma mediated off-resonance plasmonicenhanced ultrafast laser-induced nanocavitation[J]. Nano Letters, 2012, 12(9): 4763-4769. doi: 10.1021/nl302200w
    [17]

    NAZARI M, MIN X, ARONSON M, et al. Plasmon-enhanced pan-microbial pathogen inactivation in the cavitation regime: Selectivity without targeting[J]. ACS Applied Nano Materials, 2019, 2(4): 2548-2558. doi: 10.1021/acsanm.9b00552
    [18]

    DELIBASIC H, PETROVIC V, PETROVIC I. Laser breakdown in water induced by λ=532 nm nanosecond pulses: Analytical calculation of the number density of free electrons[J]. Journal of the Phy-sical Society of Japan, 2020, 89(11): 114501. doi: 10.7566/JPSJ.89.114501
    [19]

    BULAT A F, DYRDA V I, LYSYTSYA M I, et al. Numerical simulation of the stress-strain state of thin-layer rubber-metal vibration absorber elements under nonlinear deformation[J]. Strength of Materials, 2018, 50(3): 387-395. doi: 10.1007/s11223-018-9982-9
    [20]

    ZHANG Ch, TANG M, ZHANG H C, et al. Optical breakdown during femtosecond laser propagation in water cloud[J]. Optics Express, 2019, 27(6): 8456-8475. doi: 10.1364/OE.27.008456
    [21]

    ZHANG Ch, LU J, ZHANG H C, et al. Transient coupling model of plasma and laser field in water[J]. IEEE Journal of Quantum Electronics, 2016, 52(8): 1-9.
    [22]

    JARNAC A, TAMOSAUSKAS G, MAJUS D, et al. Whole life cycle of femtosecond ultraviolet filaments in water[J]. Physical Review A, 2014, 89(3): 033809.
    [23]

    NOACK J, VOGEL A. Laser-induced plasma formation in water at nanosecond to femtosecond time scales: Calculation of thresholds, absorption coefficients, and energy density[J]. IEEE Journal of Quantum Electronics, 1999, 35(8): 1156-1167.
    [24]

    VOGEL A, LINZ N, FREIDANK S, et al. Femtosecond-laser-induced nanocavitation in water: Implications for optical breakdown threshold and cell surgery[J]. Physical Review Letters, 2008, 100(3): 038102.
  • [1] 倪晓武陆建贺安之王平秋马孜周九林 . 高功率激光对光学介质薄膜破坏机理的研究进展. 激光技术, 1994, 18(6): 348-352.
    [2] 陆建倪晓武贺安之 . LSD波点燃和初始阶段传播的光学干涉研究. 激光技术, 1996, 20(2): 73-77.
    [3] 赵琦樊红英李轶国蒋泽伟胡绍云赖庚辛黄燕琳耿旭 . 大能量中空光束大气传输的仿真与实验比对研究. 激光技术, 2014, 38(4): 542-545. doi: 10.7510/jgjs.issn.1001-3806.2014.04.021
    [4] 刘兵王巨胜杨泽后李晓锋樊冬任鹏李斌罗雄冯力天 . 基于非视线红外激光大气散射通信技术研究. 激光技术, 2014, 38(6): 854-858. doi: 10.7510/jgjs.issn.1001-3806.2014.06.028
    [5] 杨军王慧章曦 . 大气湍流对部分相干激光瑞利区间影响的研究. 激光技术, 2016, 40(3): 456-460. doi: 10.7510/jgjs.issn.1001-3806.2016.03.033
    [6] 刘玉丽 . 探测边界层大气温度的转动喇曼激光雷达. 激光技术, 2018, 42(4): 541-544. doi: 10.7510/jgjs.issn.1001-3806.2018.04.020
    [7] 姜楠李晓英牛春晖刘鑫 . 大气湍流对激光空间传输特性影响的实验研究. 激光技术, 2022, 46(5): 708-712. doi: 10.7510/jgjs.issn.1001-3806.2022.05.022
    [8] 刘琳琳杨杰黄见苑克娥尹凯欣胡顺星 . 激光雷达观测淮南大气SO2和NO2浓度廓线实例分析. 激光技术, 2019, 43(3): 353-358. doi: 10.7510/jgjs.issn.1001-3806.2019.03.012
    [9] 石玉立杨丰恺曹念文 . MPL反演南京北郊气溶胶光学厚度准确度的研究. 激光技术, 2015, 39(3): 372-376. doi: 10.7510/jgjs.issn.1001-3806.2015.03.020
    [10] 李志刚王声波郭大浩吴鸿兴戴宇生 . 强激光引导和诱发大气高压放电研究. 激光技术, 1999, 23(6): 332-334.
    [11] 程元丽李思宁王骐 . 激光等离子体和气体放电EUV光刻光源. 激光技术, 2004, 28(6): 561-564.
    [12] 罗锦锋宋世军王平秋刘全喜 . 激光等离子体对硅表面微纳粒子除去机理研究. 激光技术, 2018, 42(4): 567-571. doi: 10.7510/jgjs.issn.1001-3806.2018.04.025
    [13] 陈德章高建波叶景峰路英宾王晟卿光弼王立君陈亦庆 . 纳秒紫外激光空气电离及等离子体导电实验研究. 激光技术, 2008, 32(3): 262-264.
    [14] 朱新旺王新兵傅焰峰卢彦兆石玉华 . CO2激光等离子体极紫外光源收集镜研究. 激光技术, 2010, 34(6): 725-728. doi: 10.3969/j.issn.1001-3806.2010.06.002
    [15] 刘笑兰李晓卿 . 超强激光等离子体中的相对论性朗缪尔孤子. 激光技术, 2013, 37(5): 627-630. doi: 10.7510/jgjs.issn.1001-3806.2013.05.014
    [16] 陈建平倪晓武陆建卞保民王亚伟 . 用于探测激光等离子体冲击波的光纤传感器. 激光技术, 2001, 25(2): 85-90.
    [17] 邹彪陈建平倪晓武 . 强激光等离子体有关物理量的声学测量. 激光技术, 2000, 24(6): 366-369.
    [18] 刘锴范卫星王平秋富星桥杨超韩敬华 . 激光等离子体法清洗微纳颗粒的物态变化研究. 激光技术, 2021, 45(4): 405-410. doi: 10.7510/jgjs.issn.1001-3806.2021.04.001
    [19] 汤明玥李宾中 . 双曲余弦高斯光束在非Kolmogorov湍流中的湍流距离. 激光技术, 2015, 39(4): 581-584. doi: 10.7510/jgjs.issn.1001-3806.2015.04.034
    [20] 王玉诏 . 一种用于激光雷达信号的自适应分段平滑算法. 激光技术, 2014, 38(5): 599-602. doi: 10.7510/jgjs.issn.1001-3806.2014.05.005
  • 加载中
图(7)
计量
  • 文章访问数:  2539
  • HTML全文浏览量:  1790
  • PDF下载量:  25
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-03-24
  • 录用日期:  2022-05-09
  • 刊出日期:  2023-03-25

飞秒激光致双液滴光学击穿和等离子体分布研究

    通讯作者: 陆健, lujian@njust.edu.cn
    作者简介: 钱天(1998-),男,硕士研究生,现主要从事激光与液滴相互作用的研究
  • 南京理工大学 理学院,南京 210094
基金项目:  国家自然科学基金资助项目 11774176

摘要: 为了研究超短激光脉冲与双液滴相互作用过程中的光学击穿和等离子体分布,基于麦克斯韦方程组和电离速率方程,构建了飞秒激光与双液滴的瞬态耦合模型,使用有限元分析方法,对飞秒激光辐照微米量级双液滴的自由电子密度和光场分布进行了计算,得到了双液滴结构对液滴光学击穿和等离子体变化的影响。结果表明,第2个液滴的击穿阈值约为同等条件下单液滴击穿阈值的35%;等离子体的形态和击穿点的位置随双液滴间距发生变化,且在聚焦区域产生纳米等离子体射流;第2个液滴对激光能量的吸收随着双液滴间距的增加而减少;当分别使用满足击穿阈值的光强入射,双液滴吸收的能量约为单液滴的3%;第2个液滴对激光能量的吸收随光强增大而增大,能量吸收比例最终趋于0.01,仅为单液滴的1.5%。该研究为激光诱导水击穿和激光在大气中的传输提供了一定的参考。

English Abstract

    • 超短激光脉冲在透明电介质尤其是在液滴中的光学击穿得到了广泛重视。液滴广泛存在于各种大气环境和生物系统中,在大气传输和各种生物过程中发挥着重要作用。当激光在大气中传输时,液滴集团通常被看作为光通信的阻碍,如大气中的雨雾[1-2]等往往会阻碍激光的传输。一旦激光的强度超过液滴的击穿阈值,液滴就会发生光学击穿,同时在电离过程中产生激光等离子体[3],然后液滴发生形变[4]、位移[5]等一系列动力学现象。这些等离子体的形成会影响激光的传输。当然,光学击穿也可以触发液滴的合并和重新组合,所以飞秒激光诱导的水凝结可以促进核的生长[6-7],甚至可以触发降雨和降雪的可能性。然而,液滴中的激光等离子体也并非是完全有害的,激光诱导击穿光谱通常被用来检测溶液中的成分[8-9]。在生物医学领域中,由于生物组织中含有大量的水分,飞秒激光与生物组织的作用主要是激光与组织中液滴的作用。如激光眼科手术[10]是激光技术在生物医学上的治疗应用。

      在飞秒激光与液滴的相互作用过程中,通过多光子电离[11]和雪崩电离[12]过程,如此强的激光强度在液滴内部产生激光等离子体[13-14],局部改变了液滴中水介质的光学性质。大量的激光能量将被等离子体吸收。激光能量的汇聚可以显著降低液滴的击穿阈值,同时可以在汇聚区域产生光子纳米射流[15],这些光子纳米射流会对水微滴或其它透明电介质产生影响。液滴内纳米电子等离子体通常可以用于癌症细胞的光疗和通过选择性纳米气泡空化灭活细菌或病毒[16-17]。目前已有关于激光诱导等离子体在单个水微滴中击穿的研究,建立了单个液滴光学击穿不同阶段的模型[18-19],研究了单个液滴的击穿阈值和能量吸收[20],然而关于激光与多个液滴的相互作用的研究还有待开展。当激光辐照双液滴时,一个液滴的激光等离子体演化将会有可能干扰另一个液滴附近的光场分布和光学击穿过程。

      本文中主要研究飞秒激光辐照双液滴的光学击穿和等离子体的分布变化。作者基于水介质的电离速率方程和激光光场的非线性传输方程,构建了飞秒激光与双液滴的瞬态耦合模型,对飞秒激光与双液滴的相互作用过程中的光场和电子密度的空间分布进行了计算,研究了双液滴结构对液滴击穿阈值、等离子体形态分布以及激光能量吸收的影响。

    • 激光与液滴等透明介质的相互作用过程中,当激光的强度达到一定的阈值时,激光诱导等离子体通过一些复杂的电离机制产生,主要是多光子电离、隧穿电离和雪崩电离。在水介质的电离实验中,使用的激光强度通常低于隧穿电离所需要的下限,所以本研究中不考虑隧穿电离的影响。激光诱导等离子体的演化可以用自由电子密度ρ随时间变化的速率方程来表示[21-22]

      $ \begin{array}{c} \frac{{\partial \rho }}{{\partial t}} = \left( {\frac{{{\beta _K}}}{{K\hbar \omega }}{I^K} + \frac{{{\sigma _{\rm{a}}}I}}{{{V_{\rm{p}}}}}\rho } \right)\left( {1 - \frac{\rho }{{{\rho _{\rm{n}}}}}} \right) + \\ \nabla \cdot (D\nabla \rho ) - g{\rho ^2} \end{array} $

      (1)

      式中,βKK阶多光子电离截面,σa为雪崩电离截面,Vp为介质的电离势,ρn为中性分子的密度,g为电子的复合速率,ωI分别为激光的频率和强度,$\hbar $为约化普朗克常数,t为时间。值得注意的是,基于各向同性等离子体的假设,全局方程中扩散项转换为一般扩散系数D。同时,雪崩电离的电离截面可以表示为:

      $ {\sigma _{\rm{a}}} = \frac{{{e^2}{\tau _{\rm{c}}}}}{{\left( {c{n_0}{\varepsilon _0}{m_{\rm{e}}}} \right)\left( {1 + {\omega ^2}\tau _{\rm{c}}^2} \right)}} $

      (2)

      式中,τc为自由电子和分子碰撞所需要的平均时间间隔,c为真空中的光速,eme分别为电子电荷量和电子质量,n0ε0分别是线性折射率和真空中的介电常数。

      在忽略谐波产生的情况下,等离子体中单色激光场的特性主要受衍射、光学克尔效应、散焦效应和非线性吸收的影响。为了计算耦合过程中光波的瞬态演化,基于麦克斯韦电磁场理论的非线性波动方程,电场E满足于:

      $ \nabla \times \nabla \times \mathit{\boldsymbol{E}} + \frac{1}{{{c^2}}}\frac{{{\partial ^2}}}{{\partial {t^2}}}\left( {{\varepsilon _{{\rm{RL}}}}\mathit{\boldsymbol{E}}} \right) + {\mu _0}\left( {\frac{{\partial \mathit{\boldsymbol{J}}}}{{\partial t}} + \frac{{{\partial ^2}{\mathit{\boldsymbol{P}}_{{\rm{NL}}}}}}{{\partial {t^2}}}} \right) = 0 $

      (3)

      式中,εRL为线性相对介电常数,μ0为真空磁导率,J是自由电荷的电流密度,PNL是非线性电极化矢量。事实上,等离子体的存在会使空间电流密度产生变化。通过多光子电离形成的等离子体同时会吸收光子能量,产生电流密度JMPA=cn0ε0βKIK-1E。另外,具有很高自由电子密度的等离子体可以吸收和折射光场,产生电流密度JPL=cn0ε0σcρE,其中系数σc是一个复数表达式σc= /(1-iωτc)。于是,总的电流密度为J = JMPA+ JPL。为了简化计算,将上面的电流密度使用P = P′+iJ/ω关系进行归一化处理,P P′分别为电极化矢量和光克尔效应引起的电极化矢量,得到非线性电极化矢量为:

      $ {\mathit{\boldsymbol{P}}_{{\rm{NL}}}} = 2{n_0}{n_2}{\varepsilon _0}I\mathit{\boldsymbol{E}} + {\rm{i}}c{n_0}{\varepsilon _0}\mathit{\boldsymbol{E}}\left( {{\beta _K}{I^{K - 1}} + {\sigma _c}\rho } \right)/\omega $

      (4)

      式中,n2是3阶光克尔折射率。因此,总的相对介电常数为:

      $ {\varepsilon _{\rm{r}}} = {\varepsilon _{{\rm{RL}}}} + 2{n_0}{n_2}I + {\rm{i}}\frac{{c{n_0}{\sigma _{\rm{c}}}\rho }}{\omega } + {\rm{i}}\frac{{c{n_0}{\beta _K}}}{\omega }{I^{K - 1}} $

      (5)

      同时,由于相对介电常数的变化率相比于电场的变化可以忽略不计,所以可以进行以下的近似处理:t(εrE)= Etεr+εrt EεrtE。最后,可以得到非线性波动方程:

      $ \nabla \times \nabla \times \mathit{\boldsymbol{E}} + \frac{1}{{{c^2}}}\frac{\partial }{{\partial t}}\left( {{\varepsilon _{\rm{r}}}\frac{{\partial \mathit{\boldsymbol{E}}}}{{\partial t}}} \right) = 0 $

      (6)

      上述(1)式、(5)式和(6)式组成了激光等离子体的瞬态耦合方程组,用来计算光场的传输及其生成的等离子体瞬态演化过程。

      作者将第2个液滴受到第1个液滴影响的两个液滴视为双液滴结构。如果两个液滴是并联的(两液滴球心的连线沿y轴),则各液滴的光致击穿可以近似认为是相互独立的。否则,如果两个液滴是串联的(两液滴球心的连线沿x轴,或连线在x轴上的投影是非零长的),则第2个液滴的光学击穿过程将会受到第1个液滴的影响。所以,作者主要研究的是串联的两个液滴,飞秒激光与双液滴相互作用的模型如图 1所示。飞秒激光沿x轴方向传播,激光场的偏振方向垂直于xy平面,图 1中的焦点为激光传输过程中全区域全时域上的光场功率密度最大点。为了更好地描述双液滴间距变化对等离子体形态和光场分布产生的影响,作者定义了一个无量纲参数m=d/r,其中d为第2个液滴距第1个液滴焦点的距离(第2个液滴在第1个焦点前为负数,在焦点后为正数),r为双液滴半径。对于双液滴不同的间距或半径变化时,只要m值一定,双液滴附近的光场和等离子体分布都是相似的[20]。双液滴的半径和间距是影响双液滴等离子体分布和激光能量吸收的重要因素。第2个液滴半径越小于第1个液滴时,第2个液滴越容易被击穿;而大于第1个液滴半径或相等时,基本上不发生变化。由于双液滴相对半径的改变对双液滴光学击穿的影响并不明显,所以作者主要研究双液滴间距的影响。

      图  1  飞秒激光与双液滴相互作用的模型示意图(xy平面模型)

      Figure 1.  Sketch of femtosecond laser interaction model with double droplets(xy plane model)

      本文中通过构建飞秒激光与双液滴的耦合模型,求解激光等离子体瞬态耦合方程组,使用COMSOL的有限元方法得到双液滴2维的激光光场和等离子体的瞬态演化过程。在激光瞬态耦合模型中,入射光场符合高斯光束的形式,表达式为:

      $ \begin{array}{c} \mathit{\boldsymbol{E}} = {E_0}\sqrt {\frac{{{w_0}}}{{w(x)}}} \exp \left[ { - 2(\ln 2)\frac{{{{(t - \tau )}^2}}}{{{\tau ^2}}} - {c_0}\frac{{{y^2}}}{{{w^2}(x)}}} \right] \times \\ \exp \left\{ {{\rm{i}}\left[ {\omega t - kx - k\frac{{{y^2}}}{{2{R_w}(x)}} + \mathit{\Phi }(x)} \right]} \right\} \cdot {\mathit{\boldsymbol{e}}_z} \end{array} $

      (7)

      式中,E0为高斯光束中心处电场的振幅,波数k=2πn0/λ,入射激光的波长λ=800 nm,脉冲宽度τ=100 fs,c0为光束的结构参数,这里默认为1,w(x)为光斑半径(w0是光束的束腰半径,取w0=100λ),Rw(x)为高斯光束波面的曲率半径,Φ(x)为高斯光束基模产生的相位因子,ezz方向空间矢量。由于光斑半径远大于液滴的大小,所以入射液滴的激光场可以被认为是平面波。另外,模型中使用的一些介质参数详见参考文献[20]。

    • 在飞秒激光与液滴的光学击穿过程中,飞秒激光产生的热和热弹性应力的结合在液滴中产生相变形成气泡,作者将液滴中气泡的形成作为液滴是否被击穿的判据。气泡的形成是由聚焦中心产生相变的温度决定的,温度的上升与等离子体的自由电子密度有关。以前的研究将液滴局部自由电子密度达到1018 cm-3~1020 cm-3[23],认为液滴发生光学击穿。而对于液滴达到气泡形成的温度来说,此时的临界自由电子密度在1020 cm-3量级[24]。在本文的计算中,将等离子体自由电子密度达到1020 cm-3作为击穿阈值的判据。

      为了探究双液滴对击穿阈值的影响,比较同等条件下单液滴和双液滴的击穿阈值。将第2个液滴放在第1个液滴的焦点处,使得第2个液滴发生击穿,得到双液滴的击穿阈值。对于粒径小于激光波长的液滴,光场在此类液滴附近会发生明显的衍射效应,即光场可以衍射穿过液滴,而对于较大型液滴,折射效应将更加明显,大量的光场能量将更容易被汇聚至液滴内部,导致击穿阈值的降低。如图 2所示,作者使用微米量级粒径的液滴,双液滴与单液滴击穿阈值的比值随着液滴半径的增加而减小,在液滴半径为8 μm后,击穿阈值之比基本趋近于35%。因此,同等条件下双液滴的击穿阈值约为单液滴的35%。

      图  2  双液滴与单液滴击穿阈值的比值随半径r的变化

      Figure 2.  Variation of the double-droplets to singledroplet breakdown threshold ratio with droplet radius r

    • 不同m值下双液滴附近等离子体和光场的空间分布如图 3所示。真实计算的激光入射位置在液滴的左侧7 μm处,图 3a~图 3c中具体时刻为256 fs,272 fs和309 fs,其对应激光诱导等离子体达到整个电离过程中自由电子密度最大值的时刻,同时为了使等离子体的形态分布更加鲜明,对图中自由电子密度ρ(单位为m-3)进行了对数化处理。其中,光场中的最大激光强度Ιmax分别为: 31.1 TW/cm2(见图 3d)、33.5 TW/cm2(见图 3e)和47.3 TW/cm2(见图 3f)。激光脉冲从左向右传播,入射激光的光强为1.5 TW/cm2r=8 μm,此时尺寸参数q=2πr/λ≈60,衍射效应的影响较小。

      图  3  不同m值下双液滴等离子体形态的分布变化

      Figure 3.  Distribution of plasma morphology in double droplets under different m values

      图 3a图 3b中,随着m值的不断增加,第2个液滴的自由电子密度分布从一个横向的“Y”形逐渐变成“I”形,等离子体分布更加的集中,在聚焦区域产生纳米等离子体射流。这种结果反映了光场的分布,如图 3d所示。由于第1个液滴可以被视作为透镜,将激光汇聚在第2个液滴内部发生击穿,且击穿点将会随着m值的增大而不断前移,直至在液滴表面发生击穿。如图 3b图 3c所示,此时第2个液滴在表面发生击穿,同时第2个液滴的自由电子密度区域虽然有所减小,但其形态分布并没有发生太多变化。由图 3f可以发现,此时第2个液滴在第1个液滴焦点后面,第2个液滴受到的光强将会随着m值的再次增大而逐渐减小,直至第2个液滴无法发生击穿。

      不同m值下双液滴最大自由电子密度和总自由电子数变化如图 4所示。在图 4中,总的电子数随着m值的增大而不断地减小。当m值为-0.625~-0.2时,最大自由电子密度基本不变。这是由于液滴的汇聚作用和第2个液滴的折射、吸收等效应,使得第2个液滴击穿点上的光强基本不变;当m值为-0.2~0时,此时击穿点的光强在不断增加,最大自由电子密度随m值增大而增大,并在m=0时达到最大值,可以说明第2个液滴放在第1个液滴的焦点处将更容易被击穿。如图 3e所示,此时m=-0.125,击穿点的光强比图 3d中的光强大;当m>0时,第2个液滴逐渐远离第一个液滴的焦点,此时第2个液滴受到的光强在不断地减弱,使得最大自由电子密度不断减小,直至无法达到击穿条件;当m值达到0.45时,此时双液滴间距为9 μm,第2个液滴达到临界击穿(第2个液滴被击穿时,双液滴间距达到最大)。

      图  4  双液滴最大自由电子密度和总电子数随m值的变化

      Figure 4.  Variations of maximum free electron density and total electron number of double droplets with m value

      对比图 3a图 3b可以发现,自由电子数的分布在不断聚集。此时总的自由电子数在不断地减少,虽然最大自由电子密度在m值为-0.625~0时先基本不变,然后不断增加,但依然可以说明自由电子数存在着不断聚集的现象,只是开始时聚集的程度不够剧烈,没有使最大自由电子密度发生变化。而随着自由电子数的聚集,最大自由电子密度在不断增加,在聚焦区域内形成纳米等离子体射流现象。

      下面讨论第2个液滴击穿点位置的变化。第2个液滴的击穿深度(第2个液滴前表面和击穿点的距离)随m值的变化,如图 5所示。随着m值的增大,击穿深度在不断地变小并且基本呈线性变化。直至当m值达到-0.2时,击穿深度基本在0.1 μm左右保持不变,可以认为此时第2个液滴在表面发生击穿。当m值最小时,击穿深度最大。第2个液滴击穿深度的变化主要是由第1个液滴类透镜的汇聚效果造成的。通过图 3a图 3b中等离子体形态分布的变化,可以更加明显地观察到击穿点的前移。相比于图 4,随着m值的增大,第2个液滴在内部发生击穿并且击穿深度在不断减小时,最大自由电子密度基本不变;当击穿点在表面时,最大自由电子密度在急剧增加,直到m=0时达到最大,然后开始不断地减少。

      图  5  双液滴击穿深度随m值的变化

      Figure 5.  Breakdown depth of double droplets varies with m values

    • 当激光的强度超过液滴的击穿阈值,液滴将会发生光学击穿。同时,大量的激光能量将被产生的等离子体通过逆韧致过程和多光子过程吸收。对于飞秒激光诱导等离子体,由于多光子吸收系数是光强的高阶函数,多光子吸收过程将在光学击穿中占重要地位。这些吸收的能量在电离过程被消耗,然后会转化为液滴的内能和液滴形变的动能。对穿过液滴的光束做圆柱体假设,基于电磁场理论,单个液滴及附近介质的能量损耗Q为:

      $ Q = \int {\int {\int_V {\rm{d}} } } V\int {\rm{d}} t(\mathit{\boldsymbol{J}} \cdot \mathit{\boldsymbol{E}}) = \pi \int {\int {\rm{d}} } y{\rm{d}}x \cdot y\int {\rm{d}} t(\mathit{\boldsymbol{J}} \cdot \mathit{\boldsymbol{E}}) $

      (8)

      飞秒激光穿过液滴的总能量Qt为:

      $ {Q_{\rm{t}}} = \int I (t) \cdot \pi {r^2}{\rm{d}}t = \frac{1}{{2\sqrt {\ln 2} }}{\pi ^{\frac{3}{2}}}{r^2}\tau {I_0} $

      (9)

      式中,I(t)表示随时间变化的光强,I0为入射激光的峰值光强。

      两个液滴的能量损耗Qm值的变化如图 6所示。其中,Q1Q2分别为第1个液滴和第2个液滴的能量损耗。随着m值的增大,第1个液滴的能量损耗Q1基本保持不变,第2个液滴的能量损耗Q2逐渐减少。当m值为-0.625~-0.2时,由于液滴击穿点不断前移,等离子体的分布区域在不断减小,所以第2个液滴的能量损耗在快速地减少,而且下降率大于其它情况。当m值达到0.25时,第2个液滴的能量损耗将小于第1个液滴。因为当第2个液滴逐渐远离第1个液滴的焦点,受到的激光光强在不断地减小,产生的自由电子数在减少,所以第2个液滴对激光能量的吸收也将逐渐减小,直至小于第1个液滴的能量吸收。根据(9)式计算出飞秒激光入射液滴的总能量为321 nJ,则此光强下第2个液滴的最大能量损耗率约为3.3%,而第2个液滴发生临界击穿时的能量损耗仅约为3/10000。对比于单液滴,第2个液滴可以吸收更少的能量而达到光学击穿,约为单液滴能量吸收的1/10。

      图  6  双液滴的能量损耗随m值的变化

      Figure 6.  Variations of energy losses of double droplets with m values

      为了探究单液滴和双液滴分别达到光学击穿时吸收能量的区别,作者分别使用满足击穿阈值的光强入射,计算单液滴和双液滴发生光学击穿时吸收的能量。根据(8)式计算单液滴发生光学击穿时吸收的能量为1.04 nJ。为了使第2个液滴更容易被击穿,将第2个液滴放在m=0的位置上,此时第1个液滴吸收的能量为1.61 pJ,第2个液滴吸收的能量为27.7 pJ,总吸收的能量为29.3 pJ。对于分别使用满足击穿阈值的光强入射,双液滴总吸收的能量约为单液滴的3%。同时,由于双液滴的击穿阈值为单液滴的35%,可以使用更低的激光光强入射双液滴,其将比单液滴吸收更少的能量。

      不同m值下双液滴的能量损耗随光强的变化如图 7a所示。图中Q1为第1个液滴对激光能量的吸收,Q2分别为m=-0.5和m=0.25时第2个液滴对激光能量的吸收。对于一个特定的m值,第2个液滴的能量损耗随着激光光强的增加而增加,但是增长率持续地降低。当入射激光的光强大于3 TW/cm2时,第2个液滴对激光能量的吸收存在明显的饱和作用。在m=-0.5时,第2个液滴的能量损耗趋近于40 nJ;在m=0.25时,第2个液滴的能量损耗趋近于5 nJ。不同m值下双液滴能量损耗比例Q/Qt随光强的变化如图 7b所示。不同m值下第2个液滴的能量吸收随光强变化的趋势是相似的。当入射光强小于液滴的击穿阈值时,液滴对激光能量的吸收可以忽略。当入射光强大于第2个液滴的击穿阈值而小于第1个液滴时,第2个液滴将发生光学击穿吸收大量的激光能量,而第1个液滴没有发生光学击穿或者等离子体自由电子密度不够大,对激光能量的吸收比较微弱。所以此时第2个液滴对激光能量的吸收将会大于第1个液滴的能量吸收。当入射光强大于两个液滴的击穿阈值时,第1个液滴对激光能量的吸收比例随激光强度的增加而不断增加,在达到饱和作用后吸收比例可以达到0.7。第2个液滴对激光能量的吸收比例开始时随激光光强的增加而增加,然而当第1个液滴发生等离子体屏蔽时,第2个液滴的能量吸收比例不断下降,直至趋近于0.01,这仅仅约为单液滴吸收能量比例的1.5%。

      图  7  不同m值下双液滴能量损耗Q和能量损耗比例Q/Qt随激光光强的变化

      Figure 7.  Variation of double droplets energy losses Q and proportion of energy losses Q/Qt with laser intensity at different m values

    • 本文中研究得到波长800 nm、脉宽100 fs的激光脉冲与微米级双液滴相互作用过程中的光学击穿和等离子体的空间分布。计算结果表明,激光辐射串联形态下的双液滴,第2个液滴的击穿阈值约为同等条件下单液滴击穿阈值的35%。当m值为-0.625~0时,双液滴的等离子体形态由“Y”型逐渐变成“I”型,等离子体分布区域在不断的减小,产生纳米等离子体射流。同时,第2个液滴击穿深度也随m值的增大而减小且成线性关系,直至当m=-0.2时,第2个液滴在表面发生光学击穿。另外,第2个液滴对激光能量的吸收随着m值的增大而不断减小,直至低于第1个液滴的能量吸收。在满足液滴击穿阈值时,双液滴将比单液滴吸收更少的能量发生光学击穿。第2个液滴对激光能量的吸收随光强的增大而增大,在入射光强较大时达到饱和,吸收系数将趋于0.01,约为单液滴吸收系数的1.5%。期望作者的工作可以对激光诱导水击穿中激光阈值的降低以及激光在大气中传输的能量吸收提供参考。

参考文献 (24)

目录

    /

    返回文章
    返回