快速傅里叶变换在阿秒束线光路稳定控制中的应用

江昱佼; 高亦谈; 黄沛; 赵昆; 许思源; 朱江峰; 方少波; 滕浩; 侯洵; 魏志义

doi:10.7498/aps.68.20191164

摘要

Abstract

作者及机构信息

Authors and contacts

文章全文

1. 引　言
2. 算法介绍
2.1 快速傅里叶变换
2.2 傅里叶变换在稳定反馈中的应用
3. 实验应用及锁定结果
3.1 使用单频连续激光进行抽运探测光路的锁定
3.2 倍频光谱干涉锁定载波包络相位
4. 结　论

参考文献

摘要

本文报道了将快速傅里叶变换算法应用于抽运探测系统和倍频光谱干涉系统(f-2f interferometry), 对光路进行反馈控制的原理和结果, 分别得到相对臂长抖动均方根1.24 nm (对应时间为4.1 as)的抽运-探测光路锁定和积分时间3 ms下相对相移均方根227 mrad的慢环载波包络相位(carrier envelop phase, CEP)锁定. 这样的锁定精度可以保证产生阿秒脉冲的飞秒激光脉冲拥有稳定的CEP, 并且为后续阿秒抽运探测提供了稳定的实验条件.

关键词:

Abstract

With the unveiling of molecular and atomic dynamics, scientists crave finer and faster tools to communicate with the microworld. Attosecond pump-probe enjoys its reputation as the fastest camera, hinting ultrafast movements in the delay graph. To employ this camera, the stability and delay control should have very great accuracy comparable to the camera resolution. It is also of significant importance for stabilizing the carrier envelope phase (CEP) in few-cycle laser field. When dealing with a huge quantity of data, conventional Fourier transform algorism is challenging in high-speed control. Here we put forward the efficient calculation method, fast Fourier transform (FFT) algorism in Mach-Zehnder interferometer for arm length locking and f-2f for CEP locking. In the interferometer locking, 532 nm continuous wave laser is used in the Mach-Zehnder interferometer, and the phase of the FFT term corresponding to the delay between the two arms of the interferometer serves as a feedback signal on piezo transducer (PZT) in the delay line to reduce the change of the arm length. In the CEP control experiment, data to be analyzed are the f-2f spectrum interference fringes recorded by the spectrometer. The CEP values are obtained from the first order of FFT module output of the integrated spectrum interference fringes, and a labview program examines the relative phase drift and sends a feedback voltage signal to the PZT through the proportion integration differentiation module to compensate slow CEP drift after the chirped pulse amplification system by changing the insert length of a prism pair. The results show that the root mean square (RMS) of the arm length difference is 1.24 nm (4.1 attosecond for light to travel) per meter in the interferometer locking over 12 h, and the RMS of CEP is 227 mrad under 3 ms integration time in the CEP locking over 20 min. These results are able to meet the requirement of the accuracy for attosecond pulse generation and attosecond pump-probe experiments. We also use FFT to stabilize the CEP and relative time simultaneously in the waveform synthesis for 8 h (Huang P, Fang S, Gao Y, Zhao K, Hou X, Wei Z 2019 Appl. Phys. Lett. 115 031102), the phase-locking system results in a CEP stability of 280 mrad and a relative time stability of 110 as at a repetition rate of 1 kHz. These results imply that the FFT is versatile and reliable in ultrafast control.

Keywords:

PACS:

42.25.Hz	(Interference)
42.55.-f	(Lasers)
02.30.Nw	(Fourier analysis)
07.05.Hd	(Data acquisition: hardware and software)

作者及机构信息

1.
西安电子科技大学物理与光电工程学院, 西安　710071

2.
中国科学院物理研究所, 北京凝聚态物理国家研究中心, 北京　100190

3.
中国科学院大学, 北京　100049

4.
中国科学院西安光学精密机械研究所, 西安　710119

通信作者: 赵昆, zhaokun@iphy.ac.cn ; 魏志义, zywei@iphy.ac.cn

基金项目: 国家重点研发计划(批准号: 2017 YFB0405202, 2017 YFC0110301)、国家自然科学基金重大项目(批准号: 61690221)、国家自然科学基金重点项目(批准号: 11434016, 91850209)、国家自然科学基金(批准号:11574384, 61575219, 11774277)、中国科学院仪器研制项目(批准号: YZ201658)、中国科学院前沿科学重点研究计划(批准号: QYZDJ-SSW-JSC006)和中国科学院青年创新促进会(批准号: 2018007)资助的课题

Authors and contacts

1.
School of Physics and Optoelectronic Engineering, Xidian University, Xi’an 710071, China

2.
Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China

3.
University of Chinese Academy of Sciences, Beijing 100049, China

4.
Xi’an Institute of Optics and Precision Mechanics, Chinese Academy of Sciences, Xi’an 710119, China

Corresponding author: Zhao Kun, zhaokun@iphy.ac.cn ; Wei Zhi-Yi, zywei@iphy.ac.cn

Funds: Project supported by the National Key Research and Development Program of China (Grant Nos. 2017YFB0405202, 2017YFC0110301), the Major Program of the National Natural Science Foundation of China (Grant No. 61690221), the Key Program of the National Natural Science Foundation of China (Grant Nos. 11434016, 91850209), the National Natural Science Foundation of China (Grant Nos. 11574384, 61575219, 11774277), the Instrument Developing Project of the Chinese Academy of Sciences (Grant No. YZ201658), the Frontier Science Key Research Project of the Chinese Academy of Sciences (Grant No. QYZDJ-SSW-JSC006), and the Youth Innovation Promotion Association, Chinese Academy of Sciences, China (Grant No. 2018007)

文章全文

1. 引　言

20世纪80年代, 皮秒和飞秒超快激光飞速发展, 为电子学和光学构筑了桥梁^[1], 带来了大量的技术革新, 例如位于光频率和电子频率之间的太赫兹波段的产生^[2]. 对于更快的电子运动, 例如忽略量子效应时周期约为150 as的电子在玻尔半径上绕氢原子核的运动, 飞秒激光无法探测. 2001年, 阿秒脉冲首次在实验上得到证实^[3], 经过十几年的参数提升^[4-6], 阿秒激光不断将人类带入全新的超快世界, 阿秒条纹相机(attosecond streak camera)^[7]、离子速度成像^[8]、冷靶反冲离子动量谱^[9]实验让原子尺度下电子动力学过程的时域演化研究得以实现. 近年来固体高次谐波^[10]、凝聚态材料条纹相机(例如使用时间飞行谱仪和阿秒瞬态吸收光谱方法^[11])的发展将阿秒研究推向下一代进程, 丰富的电子带内运动、带间跃迁、干涉以及和晶格作用(例如描述晶格整体运动的慢声子和晶格独立振动的快声子)被研究和验证. 所有的这些高分辨率的实验都要求使用稳定的阿秒脉冲或者脉冲串, 用于抽运的激光必须和探测光保持极高精度的可调延时, 无论是载波包络相位(carrier envelop phase, CEP)锁定^[12]、平衡光学互相关测量技术^[13], 还是后期抽运探测实验^[14]的实施都要对光学器件的振动和实验室环境导致的相位慢漂进行控制. 2016年, Kaldun等^[15]使用抽运探测的方法在时间上重建了自电离共振的演化, 他们的臂长差抖动在时间上约10 as. 2019年, Ishii等^[16]在薄片组展宽后锁定CEP, 在2000 s内单发(采集间隔1 s)测量的CEP抖动的均方根为101 mrad. 本文介绍一种通用的锁定机制, 包括原理和在实验室中的实现.

2. 算法介绍

2.1 快速傅里叶变换

在分析时间或空间域的连续和离散信号时, 为了透过杂乱的信号看到更直观的本质信息, 并且找出想要研究的频率成分进行解析或者加密, 通常要对时间域的信号进行傅里叶变换. 设有限离散函数x(n), n = 0, 1, 2, ···, N－1, 则其傅里叶变换函数为

$X\left( k \right) = \mathop \sum\limits_{n = 0}^{N - 1} x\left( n \right){{\rm{e}}^{ - {\rm{i}}2{\text{π}}kn/N}},\;k = 0,1,2, \cdots ,N - 1,$

(1)

X(k)的傅里叶逆变换为x(n),

$x\left( n \right) = \frac{1}{N}\mathop \sum\limits_{k = 0}^{N - 1} X\left( k \right){{\rm{e}}^{{\rm{i}}2{\text{π}}kn/N}},n = 0,1,2, \cdots ,N - 1,$

(2)

其中虚数单位 ${\rm{i}} = \sqrt { - 1}$ .

这种离散傅里叶变换和逆变换的时间复杂度是N². 实际应用中采样点经常很多, 庞大的计算量导致计算效率低下, 可能导致数据堆积. 为了简化计算, 出现了若干种分支计算方式, 例如离散余弦变换和离散正弦变换, 在这里当N很大时计算量依然巨大, 基本的傅里叶变换不现实.

1965年, Cooley和Tukey^[17]充分利用了傅里叶变换的周期性和对称性以及蝴蝶操作对展开项合并, 大幅度提高了计算效率, 将时间复杂度O(N²)的运算量降低为O(2Nlog₂N)量级. 应用广泛的计算软件例如matlab和labview均采用了快速傅里叶变换算法.

2.2 傅里叶变换在稳定反馈中的应用

对于数据处理软件来说, 较为稳定的干涉条纹在某一时间t₀下为一个一维数组, 由f(x, φ_x, t₀)表示, 可以写成不同频率成分g(ω, x, φ_{ω, x}, t)在t₀时刻对ω的积分,

$f\left( {x,{\varphi _x},{t_0}} \right) = \int_0^{{\omega _{\max }}} {g\left( {\omega ,x,{\varphi _{\omega ,x}},{t_0}} \right){\rm{d}}\omega } .$

(3)

为了稳定干涉条纹, 需要对条纹中主要的信息频率成分进行研究. 利用f(x, φ_x, t₀)的傅里叶变换F(ω, ϕ_ω, t₀), 可以从一个随时间变化的干涉条纹得到随时间变化的功率谱(实部)和相位谱(虚部). 功率谱中最强的零级成分(直流成分)和较弱的成分均可看为背景噪声, 选择高信噪比的一级频率成分, 对此段功率谱峰值位置对应的频率进行锁定, 可以保证干涉条纹的精确锁定. 此时一级频率成分的傅里叶逆变换为 $f_1(x, \varphi _{x,{\omega _1}}, t_0)$ , ω₁为一级中心频率. 若干涉条纹精确锁定, 不同时间t下 $\varphi _{{\omega _1}}$ 均接近, 平均值 ${\bar \varphi _{{\omega _1}}}$ 将是一个稳定值, 反之亦然.

3. 实验应用及锁定结果

3.1 使用单频连续激光进行抽运探测光路的锁定

极紫外(extreme ultra-violet, XUV)抽运探测在研究分子原子能级结构和动力学有着广泛的应用^[18], 在实验中要求抽运光和探测光的延时可调, 步长精确至阿秒量级, 这就要求两臂相对位置精确锁定. 为此, 我们在光路中加入较弱的532 nm 连续激光.

当频率为ω, 传播常数为k的平面波

${{E}}\left( t \right) = {{E}}\cos \left( {\omega t - {kr} } \right),$

(4)

通过马赫-曾德尔干涉仪后, 干涉平面的电场强度可表示为

${{E}}\left( t \right) = {{{E}}_1}\cos \left( {\omega t - {{{k}}_1}{{{r}}_1}} \right) + {{{E}}_2}\cos \left[ {\omega \!\left( {t + \tau } \right) -{{{k}}_2}{{{r}}_2} } \right],$

(5)

这里E₁和E₂分别为两路光电场矢量, τ是两臂光程差. 合成场的光强为

$I = {I_1} + {I_2} + 2\sqrt {{I_1}{I_2}} \cos \theta \cos \varphi ,$

(6)

$\varphi = {{{k}}_2}{{{r}}_2} - {{{k}}_1}{{{r}}_1},$

(7)

其中, θ是E₁和E₂的夹角, φ是θ和τ的函数. 所以以上方程可写为

$I = {I_1} + {I_2} + f\left( {\theta ,\tau } \right),$

(8)

可见光强包含非相干叠加项I₁和I₂, 而干涉项f(θ, τ)是θ和τ的函数. 当干涉仪处于较为稳定的状态时, θ可看作保持不变, 干涉条纹的变化是由于τ的改变, 当τ锁定, 干涉条纹保持稳定. 也就是说, 可以通过干涉条纹的稳定性来评估光路锁定的效果.

实验装置如图1, 分束片分别将近红外(near infrared ray, NIR)脉冲激光和扩束后的532 nm 连续激光分成两路, 一路较强的NIR用来产生XUV光, 金属膜滤除残留的NIR, XUV经过轮胎镜聚焦到时间飞行谱仪(time of flight, TOF)入口附近; 另一路较弱的NIR作为抽运信号通过用于调节本路光臂长的精密平移台, 经过打孔镜与XUV复合, 进入时间飞行谱仪. 与此同时, 一路绿光从金属膜边缘通过, 在打孔镜表面反射, 另一路伴随着作为抽运光的NIR穿过打孔镜, 这两束绿光经过透镜放大后光斑在电荷耦合器件(CCD)上重合. 这个系统中, CCD是软件控制系统的信号输入硬件, 调节光程的压电陶瓷(piezo transducer, PZT)精密平移台根据软件反馈对光路相对长度精确锁定并调节延时.

$图 1 抽运探测实验示意图 (m1−m15, 平面高反镜; f1, 平凹透镜; f2, f3, 凸透镜)\r\nFig. 1. Sketch of the pump-probe setup (m1−m15, flat mirror; f1, flat-concave lens; f2, f3, convex lens).$

图 1 抽运探测实验示意图 (m1−m15, 平面高反镜; f1, 平凹透镜; f2, f3, 凸透镜)

Fig. 1. Sketch of the pump-probe setup (m1−m15, flat mirror; f1, flat-concave lens; f2, f3, convex lens).

下载: 全尺寸图片幻灯片

为了对光路锁定系统的可行性进行验证, 搭建了如图2所示的测试光路. 测试光路与抽运探测实验光路原理相同, 一束532 nm绿光扩束后进入马赫增德尔干涉仪, M3和M4放在平移台上调节臂长, 偏振片P和半波片配合调节光强. 另一路M7放在PZT平移台上, 两路光通过打孔镜复合, 聚焦在CCD上.

$图 2 用532 nm连续激光干涉锁定光路示意图(M1, 凹透镜; M2, 凸透镜; BS, 分束片; M3−M10, 平面银镜; P, 偏振片; ${\rm{\lambda }}/2$, 半波片; FL, 凸透镜; CW, 连续)\r\nFig. 2. Arm length stabilization of a Mach-Zehnder interferometer (M1, concave lens; M2, convex lens; M3−M10, plane mirrors; P, polarizer; ${\rm{\lambda }}/2$, half wave plate; FL, focusing lens; CW, continuous wave).$

图 2 用532 nm连续激光干涉锁定光路示意图(M1, 凹透镜; M2, 凸透镜; BS, 分束片; M3−M10, 平面银镜; P, 偏振片;

${\rm{\lambda }}/2$

, 半波片; FL, 凸透镜; CW, 连续)

Fig. 2. Arm length stabilization of a Mach-Zehnder interferometer (M1, concave lens; M2, convex lens; M3−M10, plane mirrors; P, polarizer;

${\rm{\lambda }}/2$

, half wave plate; FL, focusing lens; CW, continuous wave).

下载: 全尺寸图片幻灯片

编写了labview程序实现光路的反馈控制. 该程序通过采集卡取得CCD上的干涉条纹阵列, 对条纹做纵向积分得到随时间变化的f(x, φ_x, t), 将计算得到的 ${\bar \varphi _{{\omega _1}}}$ 偏移量 $\Delta {\bar \varphi _{{\omega _1}}}$ 输入到比例-积分-微分(proportion integration differentiation, PID)控制模块和数模转换器, 获得输出到PZT平移台的电压值; 通过调整x范围、纵向积分范围、快速傅里叶变换一级频率选择范围、PID参数等降低相位的快速抖动和慢漂导致的偏移量 $\Delta {\bar \varphi _{{\omega _1}}}$ , 从而得到稳定的 ${\bar \varphi _{{\omega _1}}}$ , 并且可以长时间锁定. 图3为12 h的锁定结果, 锁定后两臂(单臂长约1 m)光程差抖动的均方根(root-mean-square, RMS)为1.23 nm (14.6 mrad). 如此精密的锁定, 不仅可以在阿秒产生测量装置中发挥重要的作用^[19], 还能给非共线抽运探测实验提供阿秒量级的稳定性^[20].

$图 3 用532 nm连续激光干涉锁定12 h的稳定性, 其中均方差为14.6 mrad, 插图为未锁定下的相对相位漂移\r\nFig. 3. Interferometer locking stability over 12 hours. Inset is the relative phase drift when armlength is not locked.$

图 3 用532 nm连续激光干涉锁定12 h的稳定性, 其中均方差为14.6 mrad, 插图为未锁定下的相对相位漂移

Fig. 3. Interferometer locking stability over 12 hours. Inset is the relative phase drift when armlength is not locked.

下载: 全尺寸图片幻灯片

3.2 倍频光谱干涉锁定载波包络相位

CEP是脉冲包络最高点附近载波电场的峰值相对于脉冲强度包络最高点处的相位偏移. 对于中心频率ω₀, 脉宽τ的高斯脉冲的电场E(t)可表示为

$E\left( t \right) = {{\rm{e}}^{ -\textstyle \frac{{2\ln 2}}{{{\tau ^2}}}{t^2}}}{\rm{cos}}\left( {{\omega _0}t + {\varphi _{{\rm{CE}}}}} \right),$

(9)

其中φ_CE即为CEP, 该值直接影响周期量级脉冲的实际电场形状, 如图4所示. 因此在周期量级脉冲与物质相互作用的过程^[21], 如高次谐波以及阿秒脉冲的产生、阈上电离^[22]、探究原子分子运动过程等研究中, CEP及其稳定性是重要参数.

$图 4 在同一脉冲包络下不同CEP值所对应的实际电场(CEP对少周期脉冲电场实际形状影响显著)\r\nFig. 4. Actual electric field of a few cycle laser pulse under different CEP values, which affect the electric field significantly.$

图 4 在同一脉冲包络下不同CEP值所对应的实际电场(CEP对少周期脉冲电场实际形状影响显著)

Fig. 4. Actual electric field of a few cycle laser pulse under different CEP values, which affect the electric field significantly.

下载: 全尺寸图片幻灯片

CEP会在激光的传播过程中产生抖动, 抖动源来自激光产生、放大和传播过程中的量子噪声和空气流动、温度变化、机械振动等带来的环境噪声. 若要稳定脉冲的CEP值则需要对高频与低频抖动分别进行锁定. 对于高频抖动, 较为常见是利用脉冲通过周期极化晶体后产生的拍频信号(即载波包络相移CEO)与锁相环电路中的参考信号鉴相, 锁定到脉冲重复频率的某一分频处, 使得CEP呈现周期性锁定^[23], 该过程多用于激光振荡器. 而低频抖动的锁定主要用于激光放大器. 对于脉冲能量在mJ量级, 重复频率kHz量级的放大激光而言, CEP的值可以通过阈上电离法^[24]与光学干涉^[25]的办法获得并反馈控制. 阈上电离法是通过周期量级脉冲与气体分子相互作用, 并对电离出的电子飞行路径进行分析, 能够获取光脉冲内电场分布并进而推导出CEP. 但由于该方案所需系统庞大, 造价昂贵, 因此并未得到广泛应用. 而光学干涉法是利用单个脉冲内部不同频率成分之间产生光谱干涉条纹以提取CEP的相对值变化, 该方法的实验装置要简易许多. 目前最常见的是倍频光谱(f-2f)干涉法, 将光谱覆盖一个倍频程的超连续光的长波部分倍频, 与基频光的相应频率(短波)部分拍频得到光谱干涉信号, 本研究正是基于这种方法.

对高斯脉冲电场E(t), 设脉宽 $\tau = \sqrt {2\ln 2}$ 以简化计算, 并对其进行傅里叶变换得

$E\left( \omega \right) = \sqrt {2{\text{π}}} {{\rm{e}}^{ -\textstyle{{\left( {\frac{{\omega - {\omega _0}}}{2}} \right)}^2}}}{{\rm{e}}^{{\rm{i}}{\varphi _{{\rm{CE}}}}}} = A\left( \omega \right){{\rm{e}}^{{\rm{i}}{\varphi _{{\rm{CE}}}}}}.$

(10)

考虑由激光电场分量E(ω)引入的电偶极矩, 选取一阶和二阶项^[26]

${P^{\left( 1 \right)}}\left( \omega \right) = {\chi ^{\left( 1 \right)}}E\left( \omega \right) = {\chi ^{\left( 1 \right)}}A\left( \omega \right){{\rm{e}}^{{\rm{i}}{\varphi _{{\rm{CE}}}}}},$

(11)

$\begin{split} \, & {P^{\left( 2 \right)}}\left( \omega \right) = {\chi ^{\left( 2 \right)}}\mathop \int \nolimits_{ - \infty }^\infty E\left( {{\omega _1}} \right)E\left( {\omega - {\omega _1}} \right){\rm{d}}{\omega _1}\\ & \quad= {{\rm{e}}^{{\rm{i}}2{\varphi _{{\rm{CE}}}}}}{\chi ^{\left( 2 \right)}}\mathop \int \nolimits_{ - \infty }^\infty A\left( {{\omega _1}} \right)A\left( {\omega - {\omega _1}} \right){\rm{d}}{\omega _1}. \end{split}$

(12)

当短波部分与长波部分的二次谐波光谱干涉, 满足相位匹配时, 由麦克斯韦方程可得二次谐波相位滞后基频波 ${\text{π}}/2$ , 也就是 ${E_{{\rm{SH}}}}\left( \omega \right) \propto {\rm{i}}{P^{\left( 2 \right)}}\left( \omega \right)$ , 基频光通过固定的色散材料时的长波短波之间会产生一定的延迟, 设为τ_d, 因此 ${E_{\rm{F}}}\left( \omega \right) \propto {{\rm{e}}^{ - {\rm{i}}\omega {\tau _{\rm{d}}}}}{P^{\left( 1 \right)}}\left( \omega \right)$ , 可得干涉条纹的光谱强度

$\begin{split} & I\left( \omega \right) \propto {\left| {{{\rm{e}}^{ - {\rm{i}}\omega {\tau _{\rm{d}}}}}{P^{\left( 1 \right)}}\left( \omega \right) + {\rm{i}}{P^{\left( 2 \right)}}\left( \omega \right)} \right|^2}\\ \propto \,&{I_{\rm{F}}}\left( \omega \right) + {I_{{\rm{SH}}}}\left( \omega \right) + 2{\left[ {{I_{\rm{F}}}\left( \omega \right){I_{{\rm{SH}}}}\left( \omega \right)} \right]^{\textstyle\frac{1}{2}}}\\ \times \,& \cos \left[ {\omega {\tau _{\rm{d}}} + {\varphi _{{\rm{CE}}}} + \left( {{\text{π}}/2} \right)} \right], \end{split}$

(13)

其中包含了强度叠加项和干涉项, 相位与ωτ_d + φ_CE有关, 对给定的波长, 且固定的光路下, 可认为延迟时间τ_d不变, 影响干涉条纹位置的变量是φ_CE. 对I(ω)进行傅里叶变换可知

$F\left( {I\left( \omega \right)} \right) \propto {{\rm{e}}^{{{\textstyle\left( {\frac{{k - {\tau _{\rm{d}}}/\left( {2{\text{π}}} \right)}}{2}} \right)}^2}}}{{\rm{e}}^{{\rm{i}}{\varphi _{{\rm{CE}}}}}},$

(14)

${\rm{arg}}\left[ {F\left( {I\left( \omega \right)} \right){|_{k = {\tau _{\rm{d}}}/\left( {2{\text{π}}} \right)}}} \right] \propto {\varphi _{{\rm{CE}}}}.$

(15)

傅里叶变换后, 当 $k = {\tau _{\rm{d}}}/\left( {2{\text{π}}} \right)$ 时该项所对应的虚部值即包含所研究光场的CEP信息.

线偏振高斯脉冲在介质中传播时, 其电场可表示为^[27]

$E\left( {z,t'} \right) = {E_0}{{\rm{e}}^{ - \textstyle\frac{{2{\rm{ln}}2}}{{{\tau ^2}}}{{\left( {t'} \right)}^2}}}{\rm{cos}}\left( {{\omega _0}t' + \Delta \varphi + {\varphi _{{\rm{CE}}}}} \right),$

(16)

式中φ_CE为初始相位; $\Delta \varphi = {\omega _0}\left( {{1}/{{{v_{\rm{g}}}}} - {1}/{{{v_{\rm{p}}}}}} \right)z$ , 其中ω₀是载波中心频率, 群速度v_g是脉冲包络的传播速度, 相速度v_p是载波传播速度, z是脉冲传播距离; ∆φ直接影响载波电场的相位值. 且在脉冲在材料中传播时满足:

$\frac{1}{{{v_{\rm{g}}}}} = \frac1{c}{{\left( {n + {\omega _0}\frac{{\partial n}}{{\partial {\omega _0}}}} \right)}},$

(17)

${1}/{{{v_{\rm{p}}}}} = n/c,$

(18)

综上可以得出 $\Delta \varphi = \dfrac{{\omega _0^2}}{{\rm{c}}}\dfrac{{\partial n}}{{\partial {{\rm{\omega }}_0}}}z$ , 其中n为材料的折射率, c为真空中光速. 因此选用合适的材料, 改变光在材料中的传播距离, 便能对CEP进行调控. 我们的实验系统使用LAK16 玻璃材料的棱镜作为色散介质, 计算可得当棱镜插入量改变为 17.5 μm时, CEP改变 2π 弧度.

我们的锁定系统先是通过自主搭建的f-2f干涉仪, 结构见图5. 将25 fs, 中心波长780 nm的红外光经光谱展宽后, 通过尖劈分少部分能量聚焦在BBO上产生倍频光, 并利用光谱仪采集倍频部分与基频光中相同频率成分形成的干涉信号以获取CEP信息. f-2f干涉仪中加入了选偏元件, 通过转动半波片和偏振棱镜可以保证通过BBO后偏振方向旋转了 ${{\text{π}}}/{2}$ 的倍频光与基频光在同一方向上的分量相近, 从而得到调制最深的干涉条纹, 有利于提高数据处理精度. 随后将光谱仪得到的干涉数据送入labview程序, 截取有效的深调制干涉范围, 对强度-频率数据进行快速傅里叶变换. 某时刻下原始数据无论以时间、频率、或者位置为横坐标, 都看成上文提到的f(x, φ_x, t), 也就是在不同时间t下的测得横坐标用x表示的光谱图, 方便从数学上理解. 我们的目的是使得原始数据稳定, 也就是f(x, φ_x, t)不随时间t改变. 现实情况是每个t都会产生一张光谱图, 不同x处相位不同. 程序在每个时刻下在有效的x范围内处理光谱强度, 提取出含有CEP信息的相位项, 和作为参考值的初始相位对比. 如果能保持这个相位不变, 那么随着t的改变, 有效的x范围内光谱相位就基本保持不变. 具体算法是对时刻t'的干涉条纹谱进行快速傅里叶变换, f(x, φ_x)的变换表示为F(ω, ϕ_ω), 取一阶项, 选择一个有效参考点减少随机性, 例如强度最高处对应的点, 也就是频域中频率成分最多的点, 横坐标为ω₁, 取其关于参考相位ϕ_ω₁(t₀)的相位差 $\Delta {\phi _{{\omega _1}}}\left( {t'} \right)$ , 经过PID给出电压信号反馈到控制色散元件的PZT平移台, 来抑制 $\Delta {\phi _{{\omega _1}}}\left( {t'} \right)$ , 将由于机械元件振动和环境涨落对CEP造成的影响通过改变棱镜对的插入量弥补回来. 对于我们的锁定环路来说, 光谱干涉信号越强, 信噪比也会越强, 有利于增强CEP锁定的稳定性.

$图 5 f-2f装置光路图(采用共线设计以避免延时变化)\r\nFig. 5. f-2f setup diagram, where collinear design is applied to avoid delay variation.$

图 5 f-2f装置光路图(采用共线设计以避免延时变化)

Fig. 5. f-2f setup diagram, where collinear design is applied to avoid delay variation.

下载: 全尺寸图片幻灯片

图6为锁定结果, 光谱仪积分时间3 ms, 在20 min的锁定时间内CEP的抖动仅为227 mrad, 显示反馈系统在CEP的慢环锁定上有良好的效果.

$图 6 f-2f锁定CEP20分钟的效果(均方差为227 mrad)\r\nFig. 6. CEP offset over 20 minutes when CEP is locked using f-2f.$

图 6 f-2f锁定CEP20分钟的效果(均方差为227 mrad)

Fig. 6. CEP offset over 20 minutes when CEP is locked using f-2f.

下载: 全尺寸图片幻灯片

4. 结　论

本文分析了快速傅里叶变换在运算效率上的优势, 基于这种高速运算, 补偿机械振动和环境涨落引入的噪声得以实现. 利用532 nm连续激光锁定抽运探测光路中马赫-曾德尔干涉仪两臂相对光程, 通过对干涉条纹数据进行快速傅里叶变换, 计算出频域一阶项的平均相位, 由labview程序反馈到控制相对臂长的PZT上实现了相对相移的RMS为14.6 mrad的锁定, 对应时间4.1 as. 还将快速傅里叶变换应用于f-2f系统, 通过光谱干涉信号傅里叶变换的一阶项峰值处的相位偏移量, 得出需要的反馈电压, 由labview程序控制棱镜对插入量从而补偿CEP慢漂, 在光谱仪3 ms积分时间下得到了RMS为227 mrad的锁定结果. 以上结果表明快速傅里叶变换在反馈控制方面的应用有着响应速度快, 运算精度高等优势. 最近, 本研究组基于f-2f自参考装置在相干合成光场中同时实现了CEP与相对延迟(relative timing)的8 h长时间锁定, 同样基于快速傅里叶变换算法^[28]. 我们可以提供阿秒量级的相对时间稳定性, 这在当今使用不同手段观测亚原子尺度电子运动, 甚至通过光场调控驱动电子成为热门的大环境下显得尤为重要, 不仅仅是物理学, 在信息存取、化学、生物学、高分子材料方面都有广泛的应用.

参考文献

[1]	Bloembergen N 1999 Rev. Mod. Phys. 71 283 Google Scholar
[2]	Auston D 1988 Top. Appl. Phys. 60 183
[3]	Drescher M, Hentschel M, Kienberger R, Tempea G, Spielmann C, Reider G, Corkum P, Krausz F 2001 Science 291 5510
[4]	Zhao K, Zhang Q, Chini M, Wu Y, Wang X W, Chang Z 2012 Opt. Lett. 37 3891 Google Scholar
[5]	Li J, Ren X, Yin Y, Zhao K, Chew A, Cheng Y, Cunningham E, Wang Y, Hu S, Wu Y, Chini M, Chang Z 2017 Nat. Commun. 8 186
[6]	Gaumnitz T, Jain A, Pertot Y, Huppert M, Jordan I, Ardana-Lamas F, Wörner H 2017 Opt. Express 25 27506 Google Scholar
[7]	Itatani J, Quéré F, Yudin G, Ivanov M, Krausz F, Corkum P 2002 Phys. Rev. Lett. 88 173903
[8]	Eppink A, Parker D 1997 Rev. Sci. Instrum. 68 3477 Google Scholar
[9]	Dörner R, Mergel V, Jagutzki O, Spielberger L, Ullrich J, Moshammer R, Schmidt-Bögcking H 2000 Phys. Rep. 330 95 Google Scholar
[10]	Luu T, Garg M, Kruchinin S, Moulet A, Hassan M, Goulielmakis E 2015 Nature 521 498 Google Scholar
[11]	Goulielmakis E, Loh Z, Wirth A, Santra R, Rohringe N, Yakovlev V, Zherebtsov S, Pfeifer T, Azzeer A, Kling M, Leone S, Krausz F 2010 Nature 466 739 Google Scholar
[12]	Jones D, Diddams S, Ranka J, Stentz A, Windeler S, Hall J, Cundiff S 2000 Science 288 635 Google Scholar
[13]	Schibli T, Kim J, Kuzucu O, Gopinath J, Tandon S, Petrich G, Kolodziejski L, Fujimoto J, Ippen E, Kaertner F 2003 Opt. Lett. 28 947
[14]	Chini M, Wang X, Cheng Y, Wu Y, Zhao K, Zhang Q, Cunningham E, Wang Y, Zang H, Chang Z 2012 25 th IEEE Photonics Conference Burlingame, CA, USA, September 23–27, 2012 p622
[15]	Kaldun A, Blättermann A, Stooß V, Donsa S, Wei H, Pazourek R, Nagele S, Ott C, Lin C, Burgdörfer J, Pfeifer T 2016 Science 354 738
[16]	Ishii N, Xia P, Kanai T, Itatani J 2019 Opt. Express 27 11447 Google Scholar
[17]	Cooley J, Tukey J 1965 Math. Comp. 19 297 Google Scholar
[18]	Seres E, Seres J, Serrat C, Namba S 2016 Phys. Rev. B 94 165125 Google Scholar
[19]	Feng X, Gilbertson S, Mashiko H, Wang H, Khan S, Chini M, Wu Y, Zhao K, Chang Z 2009 Phys. Rev. Lett. 103 183901 Google Scholar
[20]	Pertot Y, Schmidt C, Matthews M, Chauvet A, Huppert M, Svoboda V, Conta A, Tehlar A, Baykusheva D, Wolf J, Wörner H 2017 Science 355 264 Google Scholar
[21]	Krausz F, Ivanov M 2009 Rev. Mod. Phys. 81 163 Google Scholar
[22]	Rathje T, Johnson N, Moller M, Sussmann F, Adolph D, Kubel M, Kienberger R, Kling M, Paulus G, Sayler A 2012 J. Phys. B 45 074003
[23]	Yu Z, Han H, Xie Y, Teng H, Wang Z, Wei Z 2016 Chin. Phys. B 25 044205 Google Scholar
[24]	Wittmann T, Horvath B, Helml W, Schatzel M, Gu X, Cavalieri A, Paulus G, Kienberger R 2009 Nat. Phys. 5 5 Google Scholar
[25]	Xu L, Spielmann C, Poppe A, Brabec T, Krausz F, Hansch T 1996 Opt. Lett. 21 24 Google Scholar
[26]	Kakehata M, Takada H, Kobayashi Y, Torizuka K, Fujihira Y, Homma T, Takahashi H 2001 Opt. Lett. 26 1436
[27]	张伟 2013 博士学位论文 (北京: 中国科学院大学) Zhang W 2013 Ph. D. Dissertation (Beijing: University of Chinese Academy of Sciences) (in Chinese)
[28]	Huang P, Fang S, Gao Y, Zhao K, Hou X, Wei Z 2019 Appl. Phys. Lett. 115 031102 Google Scholar

图 1 抽运探测实验示意图 (m1−m15, 平面高反镜; f1, 平凹透镜; f2, f3, 凸透镜)

Fig. 1. Sketch of the pump-probe setup (m1−m15, flat mirror; f1, flat-concave lens; f2, f3, convex lens).

下载: 全尺寸图片幻灯片

图 2 用532 nm连续激光干涉锁定光路示意图(M1, 凹透镜; M2, 凸透镜; BS, 分束片; M3−M10, 平面银镜; P, 偏振片; ${\rm{\lambda }}/2$ , 半波片; FL, 凸透镜; CW, 连续)

Fig. 2. Arm length stabilization of a Mach-Zehnder interferometer (M1, concave lens; M2, convex lens; M3−M10, plane mirrors; P, polarizer; ${\rm{\lambda }}/2$ , half wave plate; FL, focusing lens; CW, continuous wave).

下载: 全尺寸图片幻灯片

图 3 用532 nm连续激光干涉锁定12 h的稳定性, 其中均方差为14.6 mrad, 插图为未锁定下的相对相位漂移

Fig. 3. Interferometer locking stability over 12 hours. Inset is the relative phase drift when armlength is not locked.

下载: 全尺寸图片幻灯片

图 4 在同一脉冲包络下不同CEP值所对应的实际电场(CEP对少周期脉冲电场实际形状影响显著)

Fig. 4. Actual electric field of a few cycle laser pulse under different CEP values, which affect the electric field significantly.

下载: 全尺寸图片幻灯片

图 5 f-2f装置光路图(采用共线设计以避免延时变化)

Fig. 5. f-2f setup diagram, where collinear design is applied to avoid delay variation.

下载: 全尺寸图片幻灯片

图 6 f-2f锁定CEP20分钟的效果(均方差为227 mrad)

Fig. 6. CEP offset over 20 minutes when CEP is locked using f-2f.

下载: 全尺寸图片幻灯片

[1]	Bloembergen N 1999 Rev. Mod. Phys. 71 283 Google Scholar
[2]	Auston D 1988 Top. Appl. Phys. 60 183
[3]	Drescher M, Hentschel M, Kienberger R, Tempea G, Spielmann C, Reider G, Corkum P, Krausz F 2001 Science 291 5510
[4]	Zhao K, Zhang Q, Chini M, Wu Y, Wang X W, Chang Z 2012 Opt. Lett. 37 3891 Google Scholar
[5]	Li J, Ren X, Yin Y, Zhao K, Chew A, Cheng Y, Cunningham E, Wang Y, Hu S, Wu Y, Chini M, Chang Z 2017 Nat. Commun. 8 186
[6]	Gaumnitz T, Jain A, Pertot Y, Huppert M, Jordan I, Ardana-Lamas F, Wörner H 2017 Opt. Express 25 27506 Google Scholar
[7]	Itatani J, Quéré F, Yudin G, Ivanov M, Krausz F, Corkum P 2002 Phys. Rev. Lett. 88 173903
[8]	Eppink A, Parker D 1997 Rev. Sci. Instrum. 68 3477 Google Scholar
[9]	Dörner R, Mergel V, Jagutzki O, Spielberger L, Ullrich J, Moshammer R, Schmidt-Bögcking H 2000 Phys. Rep. 330 95 Google Scholar
[10]	Luu T, Garg M, Kruchinin S, Moulet A, Hassan M, Goulielmakis E 2015 Nature 521 498 Google Scholar
[11]	Goulielmakis E, Loh Z, Wirth A, Santra R, Rohringe N, Yakovlev V, Zherebtsov S, Pfeifer T, Azzeer A, Kling M, Leone S, Krausz F 2010 Nature 466 739 Google Scholar
[12]	Jones D, Diddams S, Ranka J, Stentz A, Windeler S, Hall J, Cundiff S 2000 Science 288 635 Google Scholar
[13]	Schibli T, Kim J, Kuzucu O, Gopinath J, Tandon S, Petrich G, Kolodziejski L, Fujimoto J, Ippen E, Kaertner F 2003 Opt. Lett. 28 947
[14]	Chini M, Wang X, Cheng Y, Wu Y, Zhao K, Zhang Q, Cunningham E, Wang Y, Zang H, Chang Z 2012 25 th IEEE Photonics Conference Burlingame, CA, USA, September 23–27, 2012 p622
[15]	Kaldun A, Blättermann A, Stooß V, Donsa S, Wei H, Pazourek R, Nagele S, Ott C, Lin C, Burgdörfer J, Pfeifer T 2016 Science 354 738
[16]	Ishii N, Xia P, Kanai T, Itatani J 2019 Opt. Express 27 11447 Google Scholar
[17]	Cooley J, Tukey J 1965 Math. Comp. 19 297 Google Scholar
[18]	Seres E, Seres J, Serrat C, Namba S 2016 Phys. Rev. B 94 165125 Google Scholar
[19]	Feng X, Gilbertson S, Mashiko H, Wang H, Khan S, Chini M, Wu Y, Zhao K, Chang Z 2009 Phys. Rev. Lett. 103 183901 Google Scholar
[20]	Pertot Y, Schmidt C, Matthews M, Chauvet A, Huppert M, Svoboda V, Conta A, Tehlar A, Baykusheva D, Wolf J, Wörner H 2017 Science 355 264 Google Scholar
[21]	Krausz F, Ivanov M 2009 Rev. Mod. Phys. 81 163 Google Scholar
[22]	Rathje T, Johnson N, Moller M, Sussmann F, Adolph D, Kubel M, Kienberger R, Kling M, Paulus G, Sayler A 2012 J. Phys. B 45 074003
[23]	Yu Z, Han H, Xie Y, Teng H, Wang Z, Wei Z 2016 Chin. Phys. B 25 044205 Google Scholar
[24]	Wittmann T, Horvath B, Helml W, Schatzel M, Gu X, Cavalieri A, Paulus G, Kienberger R 2009 Nat. Phys. 5 5 Google Scholar
[25]	Xu L, Spielmann C, Poppe A, Brabec T, Krausz F, Hansch T 1996 Opt. Lett. 21 24 Google Scholar
[26]	Kakehata M, Takada H, Kobayashi Y, Torizuka K, Fujihira Y, Homma T, Takahashi H 2001 Opt. Lett. 26 1436
[27]	张伟 2013 博士学位论文 (北京: 中国科学院大学) Zhang W 2013 Ph. D. Dissertation (Beijing: University of Chinese Academy of Sciences) (in Chinese)
[28]	Huang P, Fang S, Gao Y, Zhao K, Hou X, Wei Z 2019 Appl. Phys. Lett. 115 031102 Google Scholar

[1]	李月, 李竣, 薛正跃, 王晶晶, 王贵师, 高晓明, 谈图. 本振光功率锁定方法应用于激光外差辐射计的研究. 物理学报, 2023, 72(9): 093201. doi: 10.7498/aps.72.20230261
[2]	陈经纬, 罗斌, 曾小明, 母杰, 王逍. 光参量啁啾脉冲放大数值模拟平台中超短脉冲聚焦模拟算法. 物理学报, 2023, 72(9): 094204. doi: 10.7498/aps.72.20222387
[3]	徐昭, 周昕, 白星, 李聪, 陈洁, 倪洋. 基于深度学习的相位截断傅里叶变换非对称加密系统攻击方法. 物理学报, 2021, 70(14): 144202. doi: 10.7498/aps.70.20202075
[4]	徐靖翔, 孔明, 许新科. 基于旋转不变技术信号参数估计的激光扫频干涉测量方法. 物理学报, 2021, 70(3): 034205. doi: 10.7498/aps.70.20201135
[5]	杨增强, 张力达. 红外激光载波包络相位对氦原子的极紫外光(XUV)吸收谱的量子调控研究. 物理学报, 2015, 64(13): 133203. doi: 10.7498/aps.64.133203
[6]	尤良芳, 令维军, 李可, 张明霞, 左银燕, 王屹山. 基于单个BBO晶体载波包络相位稳定的高效率光参量放大器. 物理学报, 2014, 63(21): 214203. doi: 10.7498/aps.63.214203
[7]	张岩, 于旭东, 邸克, 李卫, 张靖. 压缩态光场平衡零拍探测的位相锁定. 物理学报, 2013, 62(8): 084204. doi: 10.7498/aps.62.084204
[8]	田原野, 魏珊珊, 郭福明, 李苏宇, 杨玉军. 共振条件下载波包络相位效应对阈上电离谱的影响. 物理学报, 2013, 62(15): 153202. doi: 10.7498/aps.62.153202
[9]	周庆勇, 姬剑锋, 任红飞. 非等间隔计时数据的X射线脉冲星周期快速搜索算法. 物理学报, 2013, 62(1): 019701. doi: 10.7498/aps.62.019701
[10]	韩庆生, 乔耀军, 李蔚. 基于全光时域分数阶傅里叶变换的光脉冲最小损伤传输新方法. 物理学报, 2011, 60(1): 014219. doi: 10.7498/aps.60.014219
[11]	王建良, 张春梅, 宋立伟, 冷雨欣. 双光路测量红外飞秒激光脉冲的载波包络相位稳定性. 物理学报, 2009, 58(6): 3966-3970. doi: 10.7498/aps.58.3966
[12]	盖琦, 王明伟, 李智磊, 翟宏琛. 基于离散四元数傅里叶变换的双随机相位加密技术. 物理学报, 2008, 57(11): 6955-6961. doi: 10.7498/aps.57.6955
[13]	朱江峰, 杜强, 王向林, 滕浩, 韩海年, 魏志义, 侯洵. 飞秒钛宝石放大激光脉冲的载波包络相位测量与控制. 物理学报, 2008, 57(12): 7753-7757. doi: 10.7498/aps.57.7753
[14]	邓玉强, 曹士英, 于靖, 徐涛, 王清月, 张志刚. 小波变换提取放大超短脉冲载波-包络相位的研究. 物理学报, 2008, 57(11): 7017-7021. doi: 10.7498/aps.57.7017
[15]	雷兵, 冯莹, 刘泽金. 利用全光纤耦合环实现三路光纤激光器的相位锁定. 物理学报, 2008, 57(10): 6419-6424. doi: 10.7498/aps.57.6419
[16]	韩海年, 赵研英, 张炜, 朱江峰, 王鹏, 魏志义, 李师群. PPLN晶体差频测量飞秒激光脉冲的载波包络相移. 物理学报, 2007, 56(5): 2756-2759. doi: 10.7498/aps.56.2756
[17]	毕磊, 包景东. 非线性耗散对亚稳态系统量子衰变速率的影响. 物理学报, 2007, 56(4): 1919-1923. doi: 10.7498/aps.56.1919
[18]	邓玉强, 王清月, 吴祖斌, 张志刚. 载波-包络相位对于基频光与其自身倍频光脉冲合成的影响. 物理学报, 2006, 55(2): 737-742. doi: 10.7498/aps.55.737
[19]	韩海年, 魏志义, 张军, 聂玉昕. 飞秒钛宝石激光脉冲的载波包络相移测量研究. 物理学报, 2005, 54(1): 155-158. doi: 10.7498/aps.54.155
[20]	马建国, 郭应焕, 王蕴玉. 正电子湮没寿命谱分析的快速傅里叶变换解法. 物理学报, 1994, 43(4): 547-554. doi: 10.7498/aps.43.547

计量

文章访问数: 9679
PDF下载量: 90

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

搜索

留言板

快速傅里叶变换在阿秒束线光路稳定控制中的应用