Loading [MathJax]/jax/output/HTML-CSS/jax.js

搜索

x

留言板

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

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

时间延迟双色飞秒激光中H+2的解离动力学研究

王景哲 董福龙 刘杰

黄昊, 牛奔, 陶婷婷, 罗世平, 王颖, 赵晓辉, 王凯, 李志强, 党伟. Sb2Se3薄膜表面和界面超快载流子动力学的瞬态反射光谱分析. 物理学报, 2022, 71(6): 066402. doi: 10.7498/aps.71.20211714
引用本文: 黄昊, 牛奔, 陶婷婷, 罗世平, 王颖, 赵晓辉, 王凯, 李志强, 党伟. Sb2Se3薄膜表面和界面超快载流子动力学的瞬态反射光谱分析. 物理学报, 2022, 71(6): 066402. doi: 10.7498/aps.71.20211714
Huang Hao, Niu Ben, Tao Ting-Ting, Luo Shi-Ping, Wang Ying, Zhao Xiao-Hui, Wang Kai, Li Zhi-Qiang, Dang Wei. Ultrafast carrier kinetics at surface and interface of Sb2Se3 film by transient reflectance. Acta Phys. Sin., 2022, 71(6): 066402. doi: 10.7498/aps.71.20211714
Citation: Huang Hao, Niu Ben, Tao Ting-Ting, Luo Shi-Ping, Wang Ying, Zhao Xiao-Hui, Wang Kai, Li Zhi-Qiang, Dang Wei. Ultrafast carrier kinetics at surface and interface of Sb2Se3 film by transient reflectance. Acta Phys. Sin., 2022, 71(6): 066402. doi: 10.7498/aps.71.20211714

时间延迟双色飞秒激光中H+2的解离动力学研究

王景哲, 董福龙, 刘杰
cstr: 32037.14.aps.73.20241283

Dissociation dynamic study of H+2 in time-delayed two-color femtosecond lasers

Wang Jing-Zhe, Dong Fu-Long, Liu Jie
cstr: 32037.14.aps.73.20241283
Article Text (iFLYTEK Translation)
PDF
HTML
导出引用
  • 通过数值求解薛定谔方程, 计算了具有时间延迟的泵浦及探测飞秒激光联合作用下氢分子离子解离的时间演化动力学. 研究发现, 通过调节探测光的脉宽长度可以对解离过程进行有效的操控; 同时, 结合延迟时间依赖的离子解离动能谱, 可以反演出解离过程中的电子与核的微观动力学行为. 另外, 基于能动量守恒发展了一个描述解离动力学的经典模型, 该模型能够定性地预言延迟时间依赖的解离动能谱. 利用离子动能谱对探测光频率的依赖关系, 提出了一个重构离子核间距的含时演化的方案.
    In recent years, the rapid development of ultrashort pulse laser technology has made it possible to regulate the ionization and dissociation dynamics of atoms and molecules. Among them, the microscopic dynamics of molecular dissociation have always been a hot topic. The phenomenon of molecular dissociation, which is caused by the interaction between femtosecond intense laser fields and H+2 molecules, has attracted widespread attention. Previous theoretical studies on the dissociation of H+2 molecules mainly focused on studying its dissociation dynamics through numerical calculations, with relatively few theoretical models. This paper aims to establish a simple classical model to describe the dissociation dynamics. Firstly, this paper calculates the joint distribution of nuclear energy and electronic energy in the dissociation process of H+2 molecules under the action of pump lasers by numerically solving the Schrödinger equation. The results prove that H+2 molecules initially in the ground state are dissociated into H++H after absorbing a pump photon in the pump light field. Next, this paper studies the dissociation dynamics of H+2 molecules in time-delayed two-color femtosecond lasers. We find that it greatly depends on the specific forms of the pump light and the probe light. By utilizing the dependence of the dissociation kinetic energy release (KER) spectrum on the time delay of the two-color femtosecond lasers, we retrieve the sub-attosecond microscopic dynamic behaviors of electrons and atomic nuclei in the dissociation process. Furthermore, we establish a classical model based on the conservation of energy and momentum to describe the dissociation dynamics. This model can qualitatively predict the ion dissociation KER spectrum depending on the time delay of the two-color femtosecond lasers. The electronic resonant transition between the molecular ground state and the first excited state caused by the probe light will affect the ion kinetic energy spectrum in the dissociation process. Namely, the ion kinetic energy spectrum is dependent on the frequency of the probe laser. By taking advantage of this characteristic, we propose a scheme to reconstruct the evolution of the internuclear distance with time. Our reconstruction results can qualitatively predict the trend of the numerical simulation results, and this scheme may provide some theoretical guidance for experiments.
      PACS:
      64.70.kg(Semiconductors)
      72.20.Jv(Charge carriers: generation, recombination, lifetime, and trapping)
      78.47.jg(Time resolved reflection spectroscopy)
      78.47.da(Excited states)
      通信作者: 董福龙, fldonghb@126.com ; 刘杰, jliu@gscaep.ac.cn
    • 基金项目: 国家自然科学基金委员会-中国工程物理研究院联合基金(批准号: U1930403)和国家自然科学基金(批准号: 12022513, 12404394) 资助的课题.
      Corresponding author: Dong Fu-Long, fldonghb@126.com ; Liu Jie, jliu@gscaep.ac.cn
    • Funds: Project supported by the Joint Fund of the National Natural Science Foundation of China and the China Academy of Engineering Physics (Grant No. U1930403) and the National Natural Science Foundation of China (Grant Nos. 12022513, 12404394).

    Sb2Se3薄膜在可见光区域的光吸收系数达到105 cm–1以上[1], 禁带宽度约为1.1—1.3 eV, 非常接近Si带隙(1.12 eV)[2]. 同时, Sb2Se3具有较高的载流子迁移率, 电子迁移率为15 cm2/(V·s), 空穴迁移率为42 cm2/(V·s)[3]. 因此Sb2Se3是一种非常理想的光伏材料. 近年来, Sb2Se3基太阳能电池发展十分迅速, 其最高光电转换效率达到了10%[4-6]. 与此同时, Sb2Se3在光热、光电探测器、光催化领域的应用也逐渐增多[7-9]. Sb2Se3在上述领域的应用与其光生载流子的复合动力学密切相关. 加深理解Sb2Se3的载流子复合特征, 以及载流子复合与形貌之间的关系对于开发高性能的光伏、光电、光热器件有重要的意义.

    目前关于Sb2Se3薄膜载流子复合动力学的报道包含以下几个方面. Yang等[10]和Wang等[11]分别应用太赫兹泵浦探测技术确定了Sb2Se3自由载流子的寿命介于10—30 ps. 同时Wang等[11]还根据Sb2Se3的纳秒瞬态吸收光谱, 确定其束缚载流子寿命在10 ns量级. 根据Sb2Se3晶体的时间和角分辨光电子谱测量结果, Grad等[12]报道电子由导带高能级向导带底的弛豫时间约为0.4 ps. Ghosh等[13]发现在非晶态Sb2Se3薄膜中浅束缚载流子向深束缚态的弛豫时间约为2.7 ps. 在界面电荷转移方面, Zhang等[14]比较了Sb2Se3/CdS, Sb2Se3/SnO2两种界面处电子转移过程, 证实Sb2Se3/CdS界面的电子转移与电子束缚过程存在着相互竞争, 而Sb2Se3/SnO2界面处良好能级排列可以实现高效的热电子转移. 不同于传统的GaAs和CdTe材料, Sb2Se3虽然是二元化合物半导体, 但是其特殊的低对称一维链状晶体机构, 使其本征缺陷非常复杂 [15]. 因此掌握Sb2Se3载流子复合动力学, 特别是其表面、功能界面的载流子复合动力学信息对开发Sb2Se3基功能器件有重要的参考意义.

    本文选取[hk1]择优取向生长的Sb2Se3薄膜(即一维链状晶体垂直于衬底), 利用表面瞬态反射谱分析了Sb2Se3表面以及Sb2Se3/CdS界面的超快载流子动力学过程. 表面反射测量可以避免体相载流子复合的干扰, 更加灵敏地获得Sb2Se3表面以及Sb2Se3/CdS界面的载流子动力学过程信息. 通过对Sb2Se3表面的相对反射率变化ΔR/R演化过程分析, 得到Sb2Se3表面载流子热化与带隙收缩时间约为0.2—0.5 ps. 实验结果表明载流子热化、带隙收缩过程与载流子浓度和激发光子能量有关. 本文还实验证实, 在Sb2Se3/CdS界面存在自由电子转移和浅束缚电子转移两种电子转移过程.

    将FTO基片置于线列式载气输运设备后, 对腔体进行抽真空至10–2 Pa以下, 设置FTO基片衬底温度为390 ℃, Sb2Se3源温度为480 ℃, 沉积镀膜20 min获得厚度约290 nm的Sb2Se3薄膜. 随后, 在Sb2Se3薄膜表面化学水浴沉积厚度约为50 nm的CdS薄膜(见补充材料图S1 (online)).

    Sb2Se3薄膜吸收光谱表征通过UV-4100 (Hitachi)完成. X射线衍射(X-ray diffraction, XRD)图谱通过X射线衍射谱仪(德国, D8advance, Bruke)测量得到. 通过场发射扫描电子显微镜(Nova Nano SEM450)测量Sb2Se3薄膜表面和断口形貌. 飞秒时间分辨表面瞬态反射谱利用掺钛蓝宝石飞秒激光器(Spectra Physics, Spfire ACE-100 F, 100 fs, 4.5 mJ)泵浦的瞬态吸收光谱仪(Light conversion, HARPIA)测量. 其中探测光的偏振方向为s偏振, 探测光与泵浦光的偏振方向夹角设为54.7°.

    实验采用“双光路法”测量了Sb2Se3薄膜温度改变引起的表面光强反射率的相对变化ΔR/R(T). 808 nm激光经过分光片(T/R = 1∶1)后分为透射光和反射光两部分. 其中反射光光功率P1由功率计1 (北京物科光电技术有限公司, LPE-1A)直接测量. 透射光照在Sb2Se3薄膜表面, 其反射光功率P2由功率计2 (Newport Corporation, 918D-ST-UV)探测. Sb2Se3薄膜温度通过低温恒温器(Oxford Instrument, MicrostatN)改变. 两束光功率比值P2/P1正比于薄膜表面的光强反射率. 相对于室温300 K条件下薄膜光强反射率的相对变化ΔR/R(T)的计算公式为

    ΔR/R(T)=P2P1(T)P2P1(300)P2P1(300), (1)

    其中P2P1(T)P2P1(300)分别为热力学温度T和300 K条件下两束光功率的比值.

    在表面瞬态反射光谱测量中, 探测光的测量深度约为λ/(4πn), 其中n为Sb2Se3复折射率的实部. 本文实验中, 900 nm的探测光具有最大测量深度, 约为16 nm, 对应复折射率的实部n约为4.46[1]. 而波长500 nm的泵浦光在薄膜内的吸收系数α约为4.5×105 cm–1, 产生的最小载流子布居深度约22 nm. 因此本文近似认为Sb2Se3的复折射率在测量深度范围内均匀变化. 为了排除前后表面反射光干涉对测量的影响, 实验还比较了两种不同厚度Sb2Se3薄膜相对反射率变化ΔR/R的动力学曲线, 如图1(a)所示. 可以看出, 在760和860 nm 探测波长处两种Sb2Se3薄膜的ΔR/R动力学曲线基本一致. 这表明, 本文所得Sb2Se3薄膜表面瞬态反射谱不需要考虑样品后表面反射光的干扰.

    图 1 (a)两种厚度Sb2Se3薄膜在760和860 nm处$ {{\Delta }{R}}/{{R}} $归一化曲线比较; (b) Sb2Se3/空气界面$ \dfrac{\partial{R}}{{\partial n}} $和$ \dfrac{\partial{R}}{{\partial k}} $曲线; (c) Sb2Se3/空气界面$ {{\Delta }{R}}/{{R}} $随温度的变化, 其中R为300 K时的界面反射率, $ {\Delta }{R} $为相对300 K条件下界面反射率的变化\r\nFig. 1. (a) Comparisons of $ {{\Delta }{R}}/{{R}} $ curves at probe wavelengths of 760 and 860 nm from Sb2Se3 film of two thicknesses, where the $ {{\Delta }{R}}/{{R}} $ curves have been normalized; (b) $ \dfrac{\partial{R}}{{\partial n}} $ and $ \dfrac{\partial{R}}{{\partial k}} $ curves of Sb2Se3/air interface; (c) dependence of $ {{\Delta }{R}}/{{R}} $ at Sb2Se3/air interface on temperature, where R is the reflectivity of the interface at 300 K, and $ {\Delta }{R} $ is the reflectivity difference with respect to R of the interface at 300 K.
    图 1  (a)两种厚度Sb2Se3薄膜在760和860 nm处ΔR/R归一化曲线比较; (b) Sb2Se3/空气界面RnRk曲线; (c) Sb2Se3/空气界面ΔR/R随温度的变化, 其中R为300 K时的界面反射率, ΔR为相对300 K条件下界面反射率的变化
    Fig. 1.  (a) Comparisons of ΔR/R curves at probe wavelengths of 760 and 860 nm from Sb2Se3 film of two thicknesses, where the ΔR/R curves have been normalized; (b) Rn and Rk curves of Sb2Se3/air interface; (c) dependence of ΔR/R at Sb2Se3/air interface on temperature, where R is the reflectivity of the interface at 300 K, and ΔR is the reflectivity difference with respect to R of the interface at 300 K.

    在空气/Sb2Se3界面, 探测光的光强反射率可以表示为

    R=(n1)2+k2(n+1)2+k2, (2)

    其中k为Sb2Se3薄膜复折射率的虚部. 光激发载流子的多种弛豫和复合过程会影响薄膜复折射率的实部n和虚部k, 并引起空气/Sb2Se3界面的光强反射率R变化. 根据(2)式, 计算反射率R对薄膜复折射率实部n或虚部k的依赖关系, 即RnRk分别为

    Rn=2(n1)(n+1)2+k22[(n1)2+k2](n+1)[(n+1)2+k2]2, (3)
    Rk=2k(n+1)2+k22[(n1)2+k2]k[(n+1)2+k2]2. (4)

    利用文献[1]中Sb2Se3薄膜的复折射率数据, 可以得到RnRk在不同探测波长处的分布, 如图1(b)所示. 随着探测波长的增加, Rk逐渐减小, 而Rn逐渐增大. 当探测波长λ>750 nm时, Rn高达Rk的5倍. 因此在大于750 nm探测光谱范围内, 空气/ Sb2Se3界面的反射率R变化主要源于实部n的改变. 并且在λ>750 nm光谱范围内, Sb2Se3薄膜的复折射率n远大于k, 因此空气/ Sb2Se3界面的相对反射率变化ΔR/R可以近似为

    ΔRR=4Δnn021, (5)

    其中n0为无光照条件下Sb2Se3薄膜复折射率的实部, Δn为光激发引起的薄膜复折射率实部的改变量. (5)式表明, 空气/Sb2Se3界面ΔR/R的时间演化过程能够反映Δn(t)的变化. 光注入载流子后, 引起Sb2Se3薄膜复折射率实部变化Δn(t)的因素包括Drude效应、带隙收缩、能带填充以及晶格温度升高.

    Drude效应是指自由载流子引起复折射率的实部和虚部变化. 在可见-近红外波段, 复折射率虚部的改变一般可以忽略. 自由载流子引起复折射率实部减小, 且变化量与载流子浓度N有关[16]:

    ΔnDrude=e28π2ω2ε0nNμ, (6)

    其中µ为载流子的有效质量.

    光注入载流子后, 载流子与声子耦合会引起Sb2Se3薄膜晶格温度上升, 并引起其表面反射率变化ΔR. 晶格温度上升对ΔR/R的贡献, 可以近似通过热力学温度T时薄膜相对室温(300 K)的反射率变化衡量(详见(1)式). 在不同热力学温度条件下, 测量Sb2Se3薄膜在800 nm附近相对反射率变化ΔR/R, 如图1(c)所示. 以实验所用550 nm激发波长为例, 估算载流子冷却结束后Sb2Se3薄膜晶格上升的最高温度ΔT. 假定载流子冷却后, 高于带隙的能量全部被晶格吸收(暂不考虑自由载流子贡献), 晶格吸收能量为N×(12405501.2)×1.6×1019 J/cm–3, 其中N为载流子浓度, 1.2 eV为文献报道Sb2Se3薄膜的带隙宽度. 根据Sb2Se3薄膜密度ρ = 5.84 g/cm–3, 热容Cp = 124 J/(mol·K)[17], 计算载流子浓度N = 2 × 1020 cm–3时, Sb2Se3薄膜的最大温升ΔT约为20 K. 由图1(c)可知, 由于晶格温度上升引起的ΔR/R最高增幅约0.8%.

    当光注入载流子浓度较高时, 载流子之间的交换和相互作用会引起Sb2Se3薄膜带隙减小, 即带隙收缩效应[18]. 根据Wolff[19]提出的模型, 带隙收缩减小量ΔEg与载流子浓度N13相关:

    ΔEg=e2πε0εr(3Nπ)13, (7)

    其中ε0为真空介电常数, εr为光泵浦前薄膜相对介电常数.

    除了带隙收缩效应, 光生载流子还会引起带填充效应. 带隙收缩与带填充效应均与载流子的浓度和分布有关, 且在明显高于带隙的探测波段, 两者使得折射率实部产生相反的变化[20,21], 因此需要同时考虑这两种效应. 设非光照条件下Sb2Se3薄膜在能量E处的吸收系数为α(E,E0g), 其中E0g为无光注入载流子条件薄膜带隙. 光注入载流子后, 同时考虑带的填充效应、带隙收缩效应时, 薄膜在能量E处吸收系数变为αBF+BGR(E,Eg), 其中新带隙Eg<E0g. 令α0=α(E,E0g),α=αBF+BGR(E,Eg), 根据Kramers-Kronig关系, 在探测光ωp处薄膜复折射率实部变化量Δn可以通过薄膜吸收系数改变量计算:

    ΔnBF+BGR(ωp)=cπP[EgαE2(ωp)2dEE0gα0E2(ωp)2dE], (8)

    其中PEg为柯西积分主值, E为积分变量. 当探测光能量ωp大于带隙时, 满足ωp>E0g. 令ν代替E2(ωp)2的绝对值, 则(8)式可以改写为

    ΔnBF+BGR(ωp)=cπP[E0gEgα[ν]dE+ωpE0gαα0[ν]dE+ωpαα0[+ν]dE]. (9)

    带隙收缩与带的填充((8)式)对折射率变化量Δn的影响显得更加复杂. 根据文献报道, 半导体材料的带隙收缩产生时间常数在亚皮秒时间尺度[22,23], 因此需要同时考虑载流子热化与带隙收缩对材料吸收光谱的影响. 借助图2所示的模型, 本文详细讨论ΔnBF+BGR(ωp)的变化. 首先考虑载流子由非热分布到热载流子分布阶段, 如图2(a)图2(b)所示. (9)式第一项中的α大于0, 且随着载流子热化与带隙收缩, α逐渐增加, 即(9)式第一项的积分结果使得ΔnBF+BGR(ωp)小于0, 且随时间延迟变得更小. 泡利排斥作用使(9)式第三项中的αα0小于0, 即第三项的积分结果使得ΔnBF+BGR(ωp)小于0, 伴随着热流子热化Eωp逐渐减小, 引起ΔnBF+BGR(ωp)更进一步减小. 伴随载流子的热化过程, (9)式第二项中的αα0将会小于0, 即第二项的积分结果使ΔnBF+BGR(ωp)大于0. 但在载流子弛豫到导带底或价带顶前, (9)式中第一项和第三项起主要作用. 综合考虑(9)式中三项的贡献以及Drude效应, 当ωp>Eg时, 总体上Δn<0, 光激发后将首先观测到相对反射率变化ΔR/R<0, 且伴随载流子热化与带隙收缩ΔR/R进一步减小.

    图 2 Sb2Se3薄膜载流子弛豫过程示意图\r\nFig. 2. Schematic diagram of carrier relaxation in Sb2Se3 film.
    图 2  Sb2Se3薄膜载流子弛豫过程示意图
    Fig. 2.  Schematic diagram of carrier relaxation in Sb2Se3 film.

    在第2阶段, 载流子的冷却过程, 如图2(b)图2(c)所示, 折射率实部改变量Δn的正负性与载流子浓度和探测波长能量ωp有关. 在载流子浓度、分布一定的条件下, 存在一个能量过渡点Etransition. 当ωp<Etransition时, 带的填充效应与带隙收缩均引起折射率实部Δn小于0; 当ωp>Etransition带隙收缩使Δn小于0, 而带的填充效应使Δn大于0[21,24].

    根据上述分析可知, 影响Sb2Se3薄膜相对反射率变化ΔR/R的因素包括Drude效应、晶格温度上升效应、带隙收缩效应和带隙填充效应. 这4种过程具有时间先后顺序, 如Drude效应在光注入载流子后其影响一直存在, 带隙收缩和填充效应的影响在载流子热化后更加显著, 而晶格温度上升过程则在载流子冷却过程同步发生. 当探测波长ωp<Etransition且接近带隙时, 带的填充效应与带隙收缩对折射率变化量Δn的影响一致, 在此范围内测量ΔR/R可以获得载流子热化过程信息. 当探测波长ωp>Etransition时, 带的填充效应(Δn>0)与带隙收缩(Δn<0)对折射率实部n的作用相反; 本文实验中载流子浓度在10191020cm–3范围内, 带的填充效应影响更加显著[20,25,26], 在探测波长ωp处可以观测ΔR/R由负值转为正值过程. 相应地对ΔR/R的符号改变过程分析, 可以得到材料载流子冷却过程信息.

    图3(a)为Sb2Se3薄膜和Sb2Se3/CdS薄膜的紫外可见吸收光谱. 由于CdS层非常薄, Sb2Se3薄膜和Sb2Se3/CdS薄膜的吸收光谱基本重合. 图3(b)为Sb2Se3薄膜与Sb2Se3/CdS异质结的XRD图谱, 符合文献报道的Sb2Se3物相结构[1,27]. 由图3(c)图3(d)可看出, FTO上沉积的Sb2Se3薄膜呈现为高密度的纳米棒, 其中棒的直径大多介于100—200 nm之间, 棒长约为290 nm.

    图 3 (a) Sb2Se3与Sb2Se3/CdS薄膜的吸收光谱; (b) Sb2Se3与Sb2Se3/CdS薄膜的XRD图谱; (c) Sb2Se3薄膜的表面形貌; (d) Sb2Se3薄膜的截面形貌\r\nFig. 3. (a) Absorbance spectra of Sb2Se3 and Sb2Se3/CdS film; (b) XRD diffraction spectra of Sb2Se3 and Sb2Se3/CdS film; (c) surface morphologies of Sb2Se3 film; (d) cross-sectional morphology of Sb2Se3 film.
    图 3  (a) Sb2Se3与Sb2Se3/CdS薄膜的吸收光谱; (b) Sb2Se3与Sb2Se3/CdS薄膜的XRD图谱; (c) Sb2Se3薄膜的表面形貌; (d) Sb2Se3薄膜的截面形貌
    Fig. 3.  (a) Absorbance spectra of Sb2Se3 and Sb2Se3/CdS film; (b) XRD diffraction spectra of Sb2Se3 and Sb2Se3/CdS film; (c) surface morphologies of Sb2Se3 film; (d) cross-sectional morphology of Sb2Se3 film.

    图4(a)为550 nm波长激发条件下, Sb2Se3薄膜的相对反射率变化ΔR/R的二维图, 图中虚线标识的典型时间延迟的瞬态反射谱和典型探测波长的ΔR/R动力学曲线如图4(b)图4(c)所示. 可以看出, 在500—900 nm的探测光谱范围内, 光注入载流子后ΔR/R首先均减小, 且峰值出现在650 nm附近、860 nm附近. 随着时间延迟增加, 700—900 nm探测光谱范围内ΔR/R由负转为正, 峰值出现在760 nm附近. 而波长小于700 nm的光谱范围内ΔR/R始终保持为负值. 将本文Sb2Se3薄膜的瞬态反射光谱与GaAs的表面瞬态反射谱对比[26,28], 两者具有明显相似之处. 第一, 光注入载流子后, 在能量大于吸收边的光谱范围内首先出现ΔR/R减小(ΔR/R<0). 第二, 随着时间延迟的增加, 能量大于吸收边部分光谱存在ΔR/R由负转为正的动力学过程. 由图4(b)可以看出, 在时间延迟0.8 ps时760 nm附近信号强度明显弱于860 nm附近信号. 根据2.2节表面瞬态反射模型分析, 可以确定探测波长760 nm的光子能量大于过渡能量位置Etransition, 探测波长860 nm 的光子能量小于过渡能量位置Etransition. 在探测波长760 nm附近, 带的填充效应和带隙收缩的作用相互抵消, 而860 nm附近两者作用相互叠加, 因此在时间延迟小于几皮秒时间范围内探测波长860 nm的信号明显强于760 nm处的信号幅度. 时间延迟大于10 ps后, 探测波长760 nm处ΔR/R的过零行为更显著(图4(c)). 这也进一步表明760 nm附近带隙收缩和带的填充效应的作用相反. 当热载流子冷却后带隙填充效应占主要作用, 从而该探测波长附近ΔR/R出现更加明显的过零值行为. 而探测波长860 nm附近相对较弱的ΔR/R过零值行为, 应主要来自晶格温度上升的影响. 根据2.2节中的理论分析, 利用探测波长860 nm处ΔR/R的上升沿(图4(c), 动力学过程I)分析载流子的热化与带隙收缩过程. 利用探测波长760 nm处ΔR/R的恢复过程(图4(c), 动力学过程II)分析载流子的冷却过程.

    图 4 (a) Sb2Se3薄膜的相对反射率变化$ {{\Delta }{R}}/{{R}} $的二维图像(激发波长550 nm, 载流子浓度1.55 × 1020cm–3); (b)典型时间延迟对应的$ {{\Delta }{R}}/{{R}} $谱; (c)典型探测波长对应的$ {{\Delta }{R}}/{{R}} $动力学曲线\r\nFig. 4. (a) Two-dimensional image of relative reflectivity change $ {{\Delta }{R}}/{{R}} $ of Sb2Se3 film excited by 550 nm laser (300 nJ); (b) $ {{\Delta }{R}}/{{R}} $ with various probe wavelengths at three typical time delays; (c) evolutions of $ {{\Delta }{R}}/{{R}} $ of four typical probe wavelengths
    图 4  (a) Sb2Se3薄膜的相对反射率变化ΔR/R的二维图像(激发波长550 nm, 载流子浓度1.55 × 1020cm–3); (b)典型时间延迟对应的ΔR/R谱; (c)典型探测波长对应的ΔR/R动力学曲线
    Fig. 4.  (a) Two-dimensional image of relative reflectivity change ΔR/R of Sb2Se3 film excited by 550 nm laser (300 nJ); (b) ΔR/R with various probe wavelengths at three typical time delays; (c) evolutions of ΔR/R of four typical probe wavelengths

    比较探测波长500, 630, 760 和860 nm处的ΔR/R上升沿(图4(c), 零值负最大值)可以看出, 上升时间随波长增加明显延长. 从时间尺度上看, 与该现象对应的载流子过程发生在光激发后1 ps之内, 应该与载流子的热化和带隙收缩过程有关[23]. 图5(a)展示了探测波长860 nm处ΔR/R的上升沿对激发光子能量的依赖关系. 不同激发波长对应的光注入载流子浓度约1.0 × 1020 cm–3, 以降低俄歇复合等带来的差异. 由图5(a)可以看出, 随着激发光子能量增加, 860 nm处ΔR/R信号上升时间逐渐延长. 利用高斯型仪器相应函数和e指数函数的卷积, 对探测波长860 nm处ΔR/R信号上升沿进行拟合. 当激发波长为700和760 nm时, 信号上升沿只需要90 fs (半高全宽)的仪器响应函数即可良好拟合. 激发波长为650, 600, 550及500 nm时, 860 nm处ΔR/R的上升沿需要加入额外的e指数过程才能良好拟合, 对应的拟合时间常数则表示载流子的热化与带隙收缩时间. 当激发光子介于550—650 nm之间时, 载流子热化过程和带隙收缩时间约0.20 ps.

    图 5 (a)不同激发光子能量条件下860 nm探测波长处$ {{\Delta }{R}}/{{R}} $动力学曲线比较; (b)不同激发光子能量条件下载流子热化和带隙收缩时间; (c)不同载流子浓度条件下860 nm探测波长处$ {{\Delta }{R}}/{{R}} $动力学曲线比较(激发波长550 nm); (d)不同载流子浓度条件下载流子热化和带隙收缩时间\r\nFig. 5. (a) Comparisons of kinetic curves of $ {{\Delta }{R}}/{{R}} $ at probe wavelength 860 nm with different excitation photon energies; (b) carrier thermalization and band gap renormalization times with different excitation photon energies; (c) comparisons of kinetic curves of $ {{\Delta }{R}}/{{R}} $ with different carrier concentrations (excitation wavelength of 550 nm); (d) carrier thermalization and band gap renormalization time with different carrier concentrations.
    图 5  (a)不同激发光子能量条件下860 nm探测波长处ΔR/R动力学曲线比较; (b)不同激发光子能量条件下载流子热化和带隙收缩时间; (c)不同载流子浓度条件下860 nm探测波长处ΔR/R动力学曲线比较(激发波长550 nm); (d)不同载流子浓度条件下载流子热化和带隙收缩时间
    Fig. 5.  (a) Comparisons of kinetic curves of ΔR/R at probe wavelength 860 nm with different excitation photon energies; (b) carrier thermalization and band gap renormalization times with different excitation photon energies; (c) comparisons of kinetic curves of ΔR/R with different carrier concentrations (excitation wavelength of 550 nm); (d) carrier thermalization and band gap renormalization time with different carrier concentrations.

    特别是利用500 nm光子激发时, 载流子的热化时间、带隙收缩时间增加至0.52 ps. 这可能与载流子被激发到不同的能级位置有关. 文献[29]报道, 自旋-轨道耦合和晶体场的作用使Sb2Se3价带中存在4个子带. Grad等[12]利用时间与角分辨光电子发射谱技术进行测量, 3 eV光子激发条件下Sb2Se3薄膜的载流子热化时间约为0.4 ps, 略小于本文所测结果. 这归因于光电子发射谱只反映自由电子的热化过程, 而表面瞬态反射谱则与自由电子和空穴的分布均相关. 本文还测量了不同载流子浓度对载流子热化、带隙收缩时间的影响, 如图5(c)图5(d)所示. 从图5(d)可以看出, 随着光载流子浓度增加, 载流子热化、带隙收缩时间由0.13 ps逐渐增加至0.42 ps (拟合用高斯型仪器响应函数半高全宽90 fs). 这可能与高浓度载流子对电子与LO声子作用的屏蔽效应有关[30].

    完成热化过程的载流子在导带底和价带顶部形成热载流子(载流子温度大于晶格温度). 半导体材料的热载流子冷却过程是描述载流子复合的重要参数[31], 同时对开发基于Sb2Se3薄膜的热载流子光伏器件有重要的意义[32]. 根据前文分析可知, 载流子的冷却过程会加强带的填充效应, 引起探测波长760 nm处ΔR/R由负最大值-恢复-转为正值的过程(动力学过程II). 图6(a)为光注入载流子浓度近似相等条件下, 不同激发波长对应的760 nm处ΔR/R的恢复过程(动力学过程II). ΔR/R由负最大值-恢复-转为正值过程持续时间约为30 ps. Sb2Se3薄膜的自由载流子寿命约为10—35 ps[11,12], 因此动力学过程II涉及到了载流子的复合过程、带隙收缩的恢复过程、载流子冷却过程以及晶格温度上升过程. 其中载流子的冷却与晶格温度上升过程直接关联, 且两者对ΔR/R的影响一致(图1(c)), 因此可以作为一个过程考虑. 载流子浓度近似相等条件下, 载流子的复合过程、带隙收缩恢复过程规律相同, 因此比较不同激发波长时ΔR/R的动力学过程II, 就可以得到载流子冷却过程与激发波长的关系. 对图6(a)中550 nm (载流子浓度10.32 × 1019 cm–3), 600 nm (载流子浓度10.34 × 1019 cm–3), 650 nm (载流子浓度9.59 × 1019 cm–3), 760 nm (载流子浓度9.91 × 1019 cm–3)对应的ΔR/R的动力学过程II进行单e指数拟合, 得到衰减寿命τ(图6(b)). 760 nm激发光子超过带隙能量较小, 对应的衰减寿命τ最小, 约3.0 ps. 而激发光子550, 600, 650 nm存在较大的超带隙能量, 对应的衰减寿命τ介于3.4—3.8 ps. 由于动力学过程II的衰减寿命τ包含了载流子冷却、载流子复合、带隙收缩恢复过程的影响, 因此, 本文只能估测Sb2Se3薄膜的载流子冷却时间不超过3.0—4.0 ps.

    图 6 (a)不同激发光子能量条件下探测波长760 nm处${{\Delta }{R}}/{{R}}$的动力学过程II比较; (b)不同激发光子能量条件下${{\Delta }{R}}/{{R}}$的动力学过程II的衰减寿命\r\nFig. 6. (a) Comparisons of kinetics II of ${{\Delta }{R}}/{{R}}$ at 760 nm with different excitation photon energies; (b) decay lifetime of kinetics II of ${{\Delta }{R}}/{{R}}$ with different excitation energies.
    图 6  (a)不同激发光子能量条件下探测波长760 nm处ΔR/R的动力学过程II比较; (b)不同激发光子能量条件下ΔR/R的动力学过程II的衰减寿命
    Fig. 6.  (a) Comparisons of kinetics II of ΔR/R at 760 nm with different excitation photon energies; (b) decay lifetime of kinetics II of ΔR/R with different excitation energies.

    最后, 讨论Sb2Se3与CdS界面之间的电子转移过程. CdS是目前Sb2Se3太阳能电池最常用的电子传输层材料[33,34]. 因此理解Sb2Se3/CdS界面电子转移的特征对于器件、材料方面的研究有重要意义. 图7(a)为3个典型时间延迟条件下Sb2Se3/CdS界面的ΔR/R谱. 比较Sb2Se3/CdS与Sb2Se3薄膜(图4(b))的ΔR/R谱可以发现, Sb2Se3/CdS界面ΔR/R谱在495 nm附近出现一个新的谱峰. 为了验证495 nm附近的新谱峰是否与Sb2Se3/CdS界面电子转移过程有关, 本文测量了玻璃衬底沉积的CdS的表面瞬态反射光谱(图S3 (online)). 可以看出, CdS层ΔR/R谱确实在495 nm附近存在负峰, 且峰位与其吸收边接近(图S2 (online)). 这表明, 在Sb2Se3/CdS界面存在明显的电子转移过程. 在载流子浓度1.44 × 1020 cm–3条件下, Sb2Se3与Sb2Se3/CdS 在495 nm处ΔR/R动力学曲线对比如图7(b)所示. 在时间延迟0.95 ps时, Sb2Se3ΔR/R信号已经开始恢复, 而由于电子转移过程Sb2Se3/CdS的ΔR/R信号继续减小(图7(b)红圈所示), 至1.6 ps时达到最小值. 这表明光激发后1.6 ps 时Sb2Se3/CdS界面处已经存在明显电子转移过程. 另外, Sb2Se3与CdS之间的电子转移过程降低了Sb2Se3的自由载流子寿命, 也延缓了760 nm处ΔR/R的动力学过程II的过零时间(由负值转为正值, 图S4 (online)). 这里需要说明, CdS的表面瞬态反射谱(图S3 (online))在760 nm处无明显信号, 因此动力学过程II过零时间的延缓源于Sb2Se3自由载流子寿命减小(详细分析见补充材料(online)).

    图 7 (a)典型时间延迟条件下Sb2Se3/CdS的相对反射率变化$ {{\Delta }{R}}/{{R}} $谱(激发波长650 nm); (b) Sb2Se3和Sb2Se3/CdS在495 nm处${{\Delta }{R}}/{{R}}$动力学曲线对比\r\nFig. 7. (a) Relative reflectance change $ {{\Delta }{R}}/{{R}} $ spectra of Sb2Se3/CdS at three time delays (excitation wavelength of 650 nm); (b) comparison of kinetic curves of $ {{\Delta }{R}}/{{R}} $ at 495 nm for Sb2Se3 and Sb2Se3/CdS.
    图 7  (a)典型时间延迟条件下Sb2Se3/CdS的相对反射率变化ΔR/R谱(激发波长650 nm); (b) Sb2Se3和Sb2Se3/CdS在495 nm处ΔR/R动力学曲线对比
    Fig. 7.  (a) Relative reflectance change ΔR/R spectra of Sb2Se3/CdS at three time delays (excitation wavelength of 650 nm); (b) comparison of kinetic curves of ΔR/R at 495 nm for Sb2Se3 and Sb2Se3/CdS.

    图4(a)可知, Sb2Se3薄膜的相对反射率变化ΔR/R持续到纳秒量级. 考虑到Sb2Se3薄膜的自由载流子寿命约30—50 ps, 长时间尺度的ΔR/R动力学过程应与束缚载流子有关. 根据文献[14]报道, Sb2Se3薄膜的自由载流子寿命τ = 31.7 ps, Sb2Se3/dS界面的电子转移时间约为τET = 49.3 ps, 可以估计Sb2Se3/CdS界面处的电子转移效率φ约为39.1%:

    φ=11+τET/τ. (10)

    由此可以推断, 在Sb2Se3/CdS界面处约60.9%的自由电子转化为束缚电子. 而CdS是目前高效(光电转换效率大于7%)Sb2Se3太阳能电池普遍使用的电子传输层材料. 除了自由电子转移过程, 在Sb2Se3/CdS界面还很可能存在束缚电子转移过程.

    图8为Sb2Se3/CdS与Sb2Se3薄膜在760 nm处的ΔR/R在500 ps时间尺度内的对比. 首先对比4种载流子浓度条件下Sb2Se3ΔR/R动力学. 在低载流子浓度(图8(a), 4.79×1019 cm–3)条件下, ΔR/R动力学曲线在20—200 ps时间延迟范围内出现了缓慢上升. 当载流子浓度升高至9.59× 1019, 1.44 × 1020 cm–3时(图8(b)图8(c)), 在30 ps时间延迟后ΔR/R动力学曲线先下降, 后转为缓慢上升, 约持续至300 ps. 当载流子浓度进一步升高至1.92 × 1020 cm–3时, ΔR/R动力学曲线的衰减过程更加明显, 且提前至9 ps时间延迟时发生. 提高载流子浓度, 引入的高阶复合过程降低了Sb2Se3自由载流子寿命. 由于带隙收缩效应与带的填充效应相互竞争, 其中带隙收缩在760 nm产生负的ΔR/R, 而带的填充效应在760 nm处产生正的ΔR/R. 带的填充效应与载流子浓度N相关, 而带的收缩效应与载流子浓度N13相关, 载流子浓度减小对带的填充效应影响更显著, 因此减小载流子寿命可以导致ΔR/R衰减加快. 而ΔR/R曲线的缓慢上升过程开始于自由载流子复合的后期. 同时考虑到Sb2Se3自由载流子寿命在10 ps量级, 这种缓慢的上升过程应与束缚载流子的弛豫过程有关. 在1.92 × 1020 cm–3载流子浓度条件下, ΔR/R曲线没有出现缓慢上升过程, 这可能是源于Sb2Se3缺陷态的饱和效应.

    图 8 4种载流子浓度条件下Sb2Se3与Sb2Se3/CdS在760 nm处的$ {{\Delta }{R}}/{{R}} $的动力学曲线对比 (a) $ 4.79\times {10}^{19} $ cm–3; (b) 9.59 × 1019 cm–3; (c) 1.44 × 1020 cm–3; (d) 1.92 × 1020 cm–3\r\nFig. 8. Comparisons of kinetics of $ {{\Delta }{R}}/{{R}} $ at 760 nm for Sb2Se3 and Sb2Se3/CdS under four different carrier concentrations of (a) $ 4.79\times {10}^{19} $ cm–3, (b) 9.59 × 1019 cm–3, (c) 1.44 × 1020 cm–3, (d) 1.92 × 1020 cm–3
    图 8  4种载流子浓度条件下Sb2Se3与Sb2Se3/CdS在760 nm处的ΔR/R的动力学曲线对比 (a) 4.79×1019 cm–3; (b) 9.59 × 1019 cm–3; (c) 1.44 × 1020 cm–3; (d) 1.92 × 1020 cm–3
    Fig. 8.  Comparisons of kinetics of ΔR/R at 760 nm for Sb2Se3 and Sb2Se3/CdS under four different carrier concentrations of (a) 4.79×1019 cm–3, (b) 9.59 × 1019 cm–3, (c) 1.44 × 1020 cm–3, (d) 1.92 × 1020 cm–3

    进一步对比Sb2Se3与Sb2Se3/CdS在760 nm处的ΔR/R的动力学曲线(图8(a)(c)), Sb2Se3/CdS信号中没有束缚载流子弛豫过程引发的ΔR/R缓慢上升过程, 而是信号幅度的逐渐降低. 即使是在1.92 × 1020 cm–3载流子浓度条件下, Sb2Se3/CdS信号也显示了相对较大的衰减速率(图8(d)). 本文推测Sb2Se3与CdS之间束缚电子转移是引起上述动力学过程差异的主要原因. 考虑到浅束缚电子在弛豫过程可以伴随空间扩散, 本文认为 Sb2Se3与CdS之间可以产生浅束缚电子转移过程.

    根据表面瞬态反射谱测量, 光激发后Sb2Se3/CdS界面处存在载流子过程, 如图9所示. 电子和空穴经过热化后, 形成热电子和热空穴分布. 热电子和热空穴在冷却过程中, 同时发生电子、空穴的束缚过程, 以及界面处的电子转移过程. 其中电子可以由导带直接发生转移, 或在浅束缚态上发生转移. Sb2Se3具有短的自由载流子寿命和相对较长的束缚载流子寿命, 因此界面处束缚载流子的收集对提升Sb2Se3太阳能电池的光电转换效率具有重要意义. 与束缚载流子转移相关的输运过程、界面转移机制分析需要进一步深入讨论.

    图 9 Sb2Se3/CdS界面处的主要载流子过程(CB, 导带; VB, 价带)\r\nFig. 9. Main carrier processes at the interface of Sb2Se3/CdS (CB, conduction band; VB, valence band).
    图 9  Sb2Se3/CdS界面处的主要载流子过程(CB, 导带; VB, 价带)
    Fig. 9.  Main carrier processes at the interface of Sb2Se3/CdS (CB, conduction band; VB, valence band).

    本文测量了Sb2Se3薄膜以及Sb2Se3/CdS界面的飞秒时间分辨表面瞬态反射谱. 基于Drude效应、带隙收缩、带的填充效应以及晶格温度上升效应, 提出了相对反射率变化ΔR/R(光谱范围750—900 nm)的动力学行为解释方法. 根据ΔR/R的时间演化特征, 确定Sb2Se3载流子热化、带隙收缩时间介于0.2—0.5 ps之间, 且随载流子浓度升高、激发光子能量增加而延长. 还估计了Sb2Se3的热载流子冷却时间约为3—4 ps, 与激发光子能量有关. 利用表面瞬态反射谱, 在Sb2Se3/CdS界面观测到了自由电子转移光谱特征, 证实在光激发后1.6 ps内界面处已经产生明显的自由电子转移. 同时, 实验证实了Sb2Se3/CdS界面存在浅束缚电子转移过程. 浅束缚电子收集对于提升Sb2Se3太阳能电池的光电转换效率有重要意义. 浅束缚电子界面转移涉及的输运过程以及界面转移机制将在下一步工作中讨论.

    [1]

    Alnaser A S, Tong X M, Osipov T, et al. 2004 Phys. Rev. A 93 183202Google Scholar

    [2]

    Manschwetus B, Nubbemeyer T, Gorling K, Steinmeyer G, Eichmann U, Rottke H, Sandner W 2009 Phys. Rev. Lett. 102 113002Google Scholar

    [3]

    Mi Y H, Peng P, Camus N, et al. 2020 Phys. Rev. Lett. 125 173201Google Scholar

    [4]

    Pan S Z, Zhang W B, Li H, et al. 2021 Phys. Rev. Lett. 126 063201Google Scholar

    [5]

    Guo Z N, Zhang Z H, Deng Y K, Wang J G, Ye D F, Liu J, Liu Y Q 2024 Phys. Rev. Lett. 132 143201Google Scholar

    [6]

    张颖, 王兴, 徐忠锋, 任洁茹, 张艳宁, 周贤明, 梁昌慧, 张小安 2024 物理学报 73 023101Google Scholar

    Zhang Y, Wang X, Xu Z F, Ren J R, Zhang Y N, Zhou X M, Liang C H, Zhang X A 2024 Acta Phys. Sin. 73 023101Google Scholar

    [7]

    骆炎, 余璇, 雷建廷, 陶琛玉, 张少锋, 朱小龙, 马新文, 闫顺成, 赵晓辉 2024 物理学报 73 044101Google Scholar

    Luo Y, Yu X, Lei J T, Tao C Y, Zhang S F, Zhu X L, Ma X W, Yan S C, Zhao X H 2024 Acta Phys. Sin. 73 044101Google Scholar

    [8]

    Jin W W, Wang C C, Zhao X G, et al. 2024 Chin. Phys. Lett. 41 053101Google Scholar

    [9]

    Bucksbaum P H, Zavriyev A, Muller H G, Schumacher D W 2019 Phys. Rev. Lett. 64 1883Google Scholar

    [10]

    Frasinski L J, Posthumus J H, Plumridge J, Codling K, Taday P F, Langley A J 1999 Phys. Rev. Lett. 83 3625Google Scholar

    [11]

    Jolicard G, Atabek O 1992 Phys. Rev. A 46 5845Google Scholar

    [12]

    Posthumus J H, Plumridge J, Frasinski L J, et al. 2000 J. Phys. B: At. Mol. Opt. Phys. 33 L563Google Scholar

    [13]

    Niikura H, Légaré F, Hasbani R, Ivanov M Y, Villeneuve D M, Corkum P B 2003 Nature 421 826Google Scholar

    [14]

    Staudte A, Pavičić D, Chelkowski S, et al. 2007 Phys. Rev. Lett. 98 073003Google Scholar

    [15]

    Xu H, Li Zhi C, He F, Wang X, Atia T N A, Kielpinski D, Sang R T, Litvinyuk I V 2017 Nat. Commun. 8 15849Google Scholar

    [16]

    Hanus V, Kangaparambil S, Larimian S, et al. 2019 Phys. Rev. Lett. 123 263201Google Scholar

    [17]

    Li X K, Yu X T, Ma P, Zhao X N, Wang C C, Luo S Z, Ding D J 2022 Chin. Phys. B 31 103304Google Scholar

    [18]

    Leth H A, Madsen L B, Mølmer K 2010 Phys. Rev. A 81 053409Google Scholar

    [19]

    Leth H A, Madsen L B, Mølmer K 2010 Phys. Rev. A 81 053410Google Scholar

    [20]

    Liu K L, Barth I 2021 Phys. Rev. A 103 013103Google Scholar

    [21]

    Sami F, Vafaee M, Shokri B 2016 J. Phys. B: At. Mol. Opt. Phys. 44 165601Google Scholar

    [22]

    Zhao M M, Li L H, Si B W, Wang B B, Fu B N, Han Y C 2022 Chin. Phys. Lett. 39 083401Google Scholar

    [23]

    Hu T C, Zhu S K, Zhao Y N, et al. 2022 Chin. Phys. B 31 047901Google Scholar

    [24]

    Pavicic D, Kiess A, Hansch T W, Figger H 2005 Phys. Rev. Lett. 94 163002Google Scholar

    [25]

    Magrakvelidze M, He F, Niederhausen T, Litvinyuk I V, Thumm U 2009 Phys. Rev. A 79 033410Google Scholar

    [26]

    Kling M F, Siedschlag C, Verhoef A J, et al. 2006 Science 312 246Google Scholar

    [27]

    Esry B D, Sayler A M, Wang P Q, Carnes K D, BenItzhak I 2006 Phys. Rev. Lett. 97 013003Google Scholar

    [28]

    Guo W, Lu X Q, Zhao D, Wang X L 2014 Phys. Scr. 89 025401Google Scholar

    [29]

    Krausz F, Ivanov M 2009 Rev. Mod. Phys. 81 163Google Scholar

    [30]

    Feng L Q 2015 Phys. Rev. A 92 053832Google Scholar

    [31]

    Roudnev V, Esry B D, Itzhak I B 2004 Phys. Rev. Lett. 93 163601Google Scholar

    [32]

    Gibson G N, Li M, Guo C, Neira J 1997 Phys. Rev. Lett. 79 2022Google Scholar

    [33]

    Alnaser A S, Ulrich B, Tong X M, et al. 2005 Phys. Rev. A 72 030702Google Scholar

    [34]

    Hua J J, Esry B D 2009 Phys. Rev. A 80 013413Google Scholar

    [35]

    Benis E P, Bakarezos M, Papadogiannis N A, et al. 2012 Phys. Rev. A 86 043428Google Scholar

    [36]

    Hu H T, Xu H, Bai Y, Sang R T, Litvinyuk I V, Liu P, Li R X 2016 Phys. Rev. A 94 053415Google Scholar

    [37]

    Fischer B, Kremer M, Pfeifer T, et al. 2010 Phys. Rev. Lett. 105 223001Google Scholar

    [38]

    Jia Z M, Zeng Z N, Li R X, Xu Z Z, Deng Y P 2014 Phys. Rev. A 89 023419Google Scholar

    [39]

    Zhang J, Pan X F, Du H, Xu T T, Guo J, Liu X S 2017 Opt. Commun. 382 495Google Scholar

    [40]

    Liu K L, Zhang Q B, Lu P X 2012 Phys. Rev. A 86 033410Google Scholar

    [41]

    Wanie V 2016 J. Phys. B: At. Mol. Opt. Phys. 49 025601Google Scholar

    [42]

    Balint K G G 2015 Theory of Molecular Collisions (Cambridge: Royal Society of Chemistry

    [43]

    Lu R F, Zhang P Y, Han K L 2008 Phys. Rev. E 77 066701Google Scholar

    [44]

    Lehtovaara L, Toivanen J, Eloranta J 2007 J. Comput. Phys. 221 148Google Scholar

    [45]

    Feit M D, Fleck Jr J A, Steiger A 1982 J. Comput. Phys. 47 412Google Scholar

    期刊类型引用(0)

    其他类型引用(1)

  • 图 1  不同电子本征态对应的核间距依赖的分子势能面

    Fig. 1.  Molecular potential energy surfaces as a function of the internuclear distance for different electronic eigenstates.

    图 2  单色泵浦光作用下解离动力学稳定时的核动能与电子能量的联合分布

    Fig. 2.  Joint distribution of the nuclear and electron energy after the dissociation caused by the monochromatic pump laser.

    图 3  序列双色激光中H+2的延迟时间依赖的离子解离动能谱 (a)经(10)式计算所得的H+2的离子解离动能谱, 其中探测光τ2=2T2, λ2=580 nm; (b)—(d) 与(a)相同, 区别仅在于探测光脉冲时间为τ2=4T2,6T28T2. 图中实线是经典模型的计算结果

    Fig. 3.  Time-dependent dissociation kinetic energy spectra of H+2 in sequential two-color femtosecond lasers: (a) The dissociation kinetic energy spectra of H+2 calculated by Eq. (10), in which τ2=2T2 and λ2=580 nm; (b)−(d) the same as panel (a), but τ2=4T2,6T28T2, respectively. The solid lines are the results calculated by the classical model.

    图 4  (10)式计算得到的不同探测光波长下H+2的延迟时间依赖的离子解离动能谱 (a)探测光脉冲时间为τ2=6T2, 波长为λ2=180 nm; (b)—(d)与(a)相同, 区别仅在于探测光波长为λ2=288,410,580nm

    Fig. 4.  Dissociation KER spectra calculated by Eq. (10) as a function of td: (a) τ2=6T2 and λ2=180 nm; (b)−(d) the same as (a), but λ2=288,410,580nm, respectively.

    图 5  利用动能谱重构的核间距的时间演化与数值结果对比

    Fig. 5.  Comparison between the reconstructed time evolution of the internuclear distance and the numerical simulation results.

    表 1  利用波长依赖的动能谱重构出的H+2解离过程中核间距的时间演化

    Table 1.  Reconstructed time evolution of the internuclear distance in the dissociation process of H+2 utilizing the wavelength-dependence KER spectra.

    λ2/nm ω2/a.u. R/a.u. td/a.u. t=td+T2 t Δt R ΔR
    180 [0.21093, 0.2953] [2.5,3.05] [0,20] [25,45] 35 10 2.75 0.25
    288 [0.13183, 0.18456] [3.3,3.85] [0,47] [40,87] 63.5 23.5 3.6 0.3
    410 [0.0926, 0.12964] [3.9,4.45] [0,60] [57,117] 87 30 4.2 0.3
    580 [0.06546, 0.09164] [4.45,5.0] [20,70] [100,150] 125 25 4.7 0.3
    下载: 导出CSV
  • [1]

    Alnaser A S, Tong X M, Osipov T, et al. 2004 Phys. Rev. A 93 183202Google Scholar

    [2]

    Manschwetus B, Nubbemeyer T, Gorling K, Steinmeyer G, Eichmann U, Rottke H, Sandner W 2009 Phys. Rev. Lett. 102 113002Google Scholar

    [3]

    Mi Y H, Peng P, Camus N, et al. 2020 Phys. Rev. Lett. 125 173201Google Scholar

    [4]

    Pan S Z, Zhang W B, Li H, et al. 2021 Phys. Rev. Lett. 126 063201Google Scholar

    [5]

    Guo Z N, Zhang Z H, Deng Y K, Wang J G, Ye D F, Liu J, Liu Y Q 2024 Phys. Rev. Lett. 132 143201Google Scholar

    [6]

    张颖, 王兴, 徐忠锋, 任洁茹, 张艳宁, 周贤明, 梁昌慧, 张小安 2024 物理学报 73 023101Google Scholar

    Zhang Y, Wang X, Xu Z F, Ren J R, Zhang Y N, Zhou X M, Liang C H, Zhang X A 2024 Acta Phys. Sin. 73 023101Google Scholar

    [7]

    骆炎, 余璇, 雷建廷, 陶琛玉, 张少锋, 朱小龙, 马新文, 闫顺成, 赵晓辉 2024 物理学报 73 044101Google Scholar

    Luo Y, Yu X, Lei J T, Tao C Y, Zhang S F, Zhu X L, Ma X W, Yan S C, Zhao X H 2024 Acta Phys. Sin. 73 044101Google Scholar

    [8]

    Jin W W, Wang C C, Zhao X G, et al. 2024 Chin. Phys. Lett. 41 053101Google Scholar

    [9]

    Bucksbaum P H, Zavriyev A, Muller H G, Schumacher D W 2019 Phys. Rev. Lett. 64 1883Google Scholar

    [10]

    Frasinski L J, Posthumus J H, Plumridge J, Codling K, Taday P F, Langley A J 1999 Phys. Rev. Lett. 83 3625Google Scholar

    [11]

    Jolicard G, Atabek O 1992 Phys. Rev. A 46 5845Google Scholar

    [12]

    Posthumus J H, Plumridge J, Frasinski L J, et al. 2000 J. Phys. B: At. Mol. Opt. Phys. 33 L563Google Scholar

    [13]

    Niikura H, Légaré F, Hasbani R, Ivanov M Y, Villeneuve D M, Corkum P B 2003 Nature 421 826Google Scholar

    [14]

    Staudte A, Pavičić D, Chelkowski S, et al. 2007 Phys. Rev. Lett. 98 073003Google Scholar

    [15]

    Xu H, Li Zhi C, He F, Wang X, Atia T N A, Kielpinski D, Sang R T, Litvinyuk I V 2017 Nat. Commun. 8 15849Google Scholar

    [16]

    Hanus V, Kangaparambil S, Larimian S, et al. 2019 Phys. Rev. Lett. 123 263201Google Scholar

    [17]

    Li X K, Yu X T, Ma P, Zhao X N, Wang C C, Luo S Z, Ding D J 2022 Chin. Phys. B 31 103304Google Scholar

    [18]

    Leth H A, Madsen L B, Mølmer K 2010 Phys. Rev. A 81 053409Google Scholar

    [19]

    Leth H A, Madsen L B, Mølmer K 2010 Phys. Rev. A 81 053410Google Scholar

    [20]

    Liu K L, Barth I 2021 Phys. Rev. A 103 013103Google Scholar

    [21]

    Sami F, Vafaee M, Shokri B 2016 J. Phys. B: At. Mol. Opt. Phys. 44 165601Google Scholar

    [22]

    Zhao M M, Li L H, Si B W, Wang B B, Fu B N, Han Y C 2022 Chin. Phys. Lett. 39 083401Google Scholar

    [23]

    Hu T C, Zhu S K, Zhao Y N, et al. 2022 Chin. Phys. B 31 047901Google Scholar

    [24]

    Pavicic D, Kiess A, Hansch T W, Figger H 2005 Phys. Rev. Lett. 94 163002Google Scholar

    [25]

    Magrakvelidze M, He F, Niederhausen T, Litvinyuk I V, Thumm U 2009 Phys. Rev. A 79 033410Google Scholar

    [26]

    Kling M F, Siedschlag C, Verhoef A J, et al. 2006 Science 312 246Google Scholar

    [27]

    Esry B D, Sayler A M, Wang P Q, Carnes K D, BenItzhak I 2006 Phys. Rev. Lett. 97 013003Google Scholar

    [28]

    Guo W, Lu X Q, Zhao D, Wang X L 2014 Phys. Scr. 89 025401Google Scholar

    [29]

    Krausz F, Ivanov M 2009 Rev. Mod. Phys. 81 163Google Scholar

    [30]

    Feng L Q 2015 Phys. Rev. A 92 053832Google Scholar

    [31]

    Roudnev V, Esry B D, Itzhak I B 2004 Phys. Rev. Lett. 93 163601Google Scholar

    [32]

    Gibson G N, Li M, Guo C, Neira J 1997 Phys. Rev. Lett. 79 2022Google Scholar

    [33]

    Alnaser A S, Ulrich B, Tong X M, et al. 2005 Phys. Rev. A 72 030702Google Scholar

    [34]

    Hua J J, Esry B D 2009 Phys. Rev. A 80 013413Google Scholar

    [35]

    Benis E P, Bakarezos M, Papadogiannis N A, et al. 2012 Phys. Rev. A 86 043428Google Scholar

    [36]

    Hu H T, Xu H, Bai Y, Sang R T, Litvinyuk I V, Liu P, Li R X 2016 Phys. Rev. A 94 053415Google Scholar

    [37]

    Fischer B, Kremer M, Pfeifer T, et al. 2010 Phys. Rev. Lett. 105 223001Google Scholar

    [38]

    Jia Z M, Zeng Z N, Li R X, Xu Z Z, Deng Y P 2014 Phys. Rev. A 89 023419Google Scholar

    [39]

    Zhang J, Pan X F, Du H, Xu T T, Guo J, Liu X S 2017 Opt. Commun. 382 495Google Scholar

    [40]

    Liu K L, Zhang Q B, Lu P X 2012 Phys. Rev. A 86 033410Google Scholar

    [41]

    Wanie V 2016 J. Phys. B: At. Mol. Opt. Phys. 49 025601Google Scholar

    [42]

    Balint K G G 2015 Theory of Molecular Collisions (Cambridge: Royal Society of Chemistry

    [43]

    Lu R F, Zhang P Y, Han K L 2008 Phys. Rev. E 77 066701Google Scholar

    [44]

    Lehtovaara L, Toivanen J, Eloranta J 2007 J. Comput. Phys. 221 148Google Scholar

    [45]

    Feit M D, Fleck Jr J A, Steiger A 1982 J. Comput. Phys. 47 412Google Scholar

  • [1] 钟振祥. 氢分子离子超精细结构理论综述. 物理学报, 2024, 73(20): 203104. doi: 10.7498/aps.73.20241101
    [2] 郑悦, 张宇璇, 孙少华, 丁鹏基, 胡碧涛, 刘作业. 飞秒激光脉冲对N2分子非绝热准直的调控. 物理学报, 2023, 72(6): 064203. doi: 10.7498/aps.72.20222112
    [3] 俞祖卿, 杨魏吉, 何峰. H2+在强激光脉冲作用下的电离率和原子核间距的关系. 物理学报, 2016, 65(20): 204202. doi: 10.7498/aps.65.204202
    [4] 姚云华, 卢晨晖, 徐淑武, 丁晶新, 贾天卿, 张诗按, 孙真荣. 飞秒激光脉冲整形技术及其应用. 物理学报, 2014, 63(18): 184201. doi: 10.7498/aps.63.184201
    [5] 杨青, 杜广庆, 陈烽, 吴艳敏, 欧燕, 陆宇, 侯洵. 时间整形飞秒激光诱导熔融硅表面纳米周期条纹的电子动力学研究. 物理学报, 2014, 63(4): 047901. doi: 10.7498/aps.63.047901
    [6] 姚洪斌, 张季, 彭敏, 李文亮. H2+在强激光场中的解离及其量子调控的理论研究. 物理学报, 2014, 63(19): 198202. doi: 10.7498/aps.63.198202
    [7] 王文亭, 张楠, 王明伟, 何远航, 杨建军, 朱晓农. 飞秒激光烧蚀金属靶的冲击温度. 物理学报, 2013, 62(21): 210601. doi: 10.7498/aps.62.210601
    [8] 王文亭, 张楠, 王明伟, 何远航, 杨建军, 朱晓农. 飞秒激光烧蚀固体靶的冲击压强. 物理学报, 2013, 62(17): 170601. doi: 10.7498/aps.62.170601
    [9] 刘小亮, 孙少华, 曹瑜, 孙铭泽, 刘情操, 胡碧涛. 飞秒激光低压N2等离子体特性的实验研究. 物理学报, 2013, 62(4): 045201. doi: 10.7498/aps.62.045201
    [10] 郭凯敏, 高 勋, 郝作强, 鲁毅, 孙长凯, 林景全. 空气中飞秒激光等离子体荧光辐射光谱研究. 物理学报, 2012, 61(7): 075212. doi: 10.7498/aps.61.075212
    [11] 朱竹青, 王晓雷. 飞秒激光空气等离子体发射光谱的实验研究. 物理学报, 2011, 60(8): 085205. doi: 10.7498/aps.60.085205
    [12] 高勋, 宋晓伟, 郭凯敏, 陶海岩, 林景全. 飞秒激光烧蚀硅表面产生等离子体的发射光谱研究. 物理学报, 2011, 60(2): 025203. doi: 10.7498/aps.60.025203
    [13] 郭凯敏, 高勋, 薛念亮, 赵振明, 李海军, 鲁毅, 林景全. 飞秒激光等离子体单丝导电性能的空间分辨研究. 物理学报, 2011, 60(10): 105203. doi: 10.7498/aps.60.105203
    [14] 唐小锋, 牛铭理, 周晓国, 刘世林. 基于阈值光电子-光离子符合技术的分子离子光谱和解离动力学研究. 物理学报, 2010, 59(10): 6940-6947. doi: 10.7498/aps.59.6940
    [15] 赵红敏, 王鹿霞. 异质结中桥分子电子转移的飞秒激光控制研究. 物理学报, 2009, 58(2): 1332-1337. doi: 10.7498/aps.58.1332
    [16] 王晓雷, 张 楠, 赵友博, 李智磊, 翟宏琛, 朱晓农. 飞秒激光激发空气电离的阈值研究. 物理学报, 2008, 57(1): 354-357. doi: 10.7498/aps.57.354
    [17] 仲佳勇, 李玉同, 鲁 欣, 张 翼, Bernhardt Jens, 刘 峰, 郝作强, 张 喆, 于全芝, 陈 民, 远晓辉, 梁文锡, 赵 刚, 张 杰. 空气中单个激光等离子体通道的形成条件. 物理学报, 2007, 56(12): 7114-7119. doi: 10.7498/aps.56.7114
    [18] 李成斌, 贾天卿, 孙海轶, 李晓溪, 徐世珍, 冯东海, 王晓峰, 葛晓春, 徐至展. 飞秒激光对氟化镁烧蚀机理研究. 物理学报, 2006, 55(1): 217-220. doi: 10.7498/aps.55.217
    [19] 朱频频, 刘建胜, 徐至展. Ar原子团簇与飞秒强激光相互作用产生的高能离子计算. 物理学报, 2004, 53(3): 803-807. doi: 10.7498/aps.53.803
    [20] 何 峰, 余 玮, 陆培祥. 飞秒强激光作用下线性等离子体层中光场和电子密度的自洽分布. 物理学报, 2003, 52(8): 1965-1969. doi: 10.7498/aps.52.1965
  • 期刊类型引用(0)

    其他类型引用(1)

计量
  • 文章访问数:  1375
  • PDF下载量:  57
  • 被引次数: 1
出版历程
  • 收稿日期:  2024-09-12
  • 修回日期:  2024-10-18
  • 上网日期:  2024-11-12
  • 刊出日期:  2024-12-20

/

返回文章
返回