全文HTML
--> --> -->基于微腔的光频梳是利用级联的四波混频效应, 在光学微腔中产生一系列等间距的频率成分, 从而将单一的抽运频率扩展到一个非常超宽光谱范围. 微腔光频梳具有体积小、结构简单、易实现光学集成等优点. 通过理论研究证实, 微腔可工作于多种模式, 例如: 光孤子、图灵环、暗孤子、混沌等模式[10], 并且这些工作模式也在实验中被观察到[11-15]. 目前, 在光学微腔的理论研究中, 大多都是以Lugiato-Lefeve方程 (Lugiato-Lefeve equation, LLE)为理论模型, 通过求解LLE, 得到稳态时腔内的光场分布[16]. 关于光学微腔中二次谐波的产生, 已经在实验中被观察到. Xue等[17]通过四波混频的方法在正色散的光学微腔中同时观察到了Kerr光频梳和二次谐波的产生, 此外还建立了光学微腔中二次谐波产生的理论模型, 理论分析和实验结果相一致. Lin等[18]也在片上集成的光学微腔中实现了高效率的二次谐波产生.
然而, 目前对微腔光频梳的理论研究采用的LLE, 只包含了腔内的损耗、色散、抽运项, 并没有考虑微腔中的非线性效应[19]. 但是由于光学微腔本身体积很小, 腔内光功率密度较高, 因而腔内的非线性效应是无法被忽略的. 当考虑微腔内的二阶和三阶非线性效应时, 腔内会由于非线性效应产生二次谐波[20]. 因此, 本文将在二阶、三阶非线性效应的基础上, 分析光学微腔内二次谐波的产生情况, 并研究各个参数对腔内基频光和倍频光光场产生的影响.

Figure1. Structure of the SiN microresonator.




通过求解(1)式和(2)式, 可以得到微腔内基频和倍频光场的分布情况. 为了能够清楚腔内光场分布随时间的演化过程, 采用分布傅里叶法求解(1)式和(2)式[19].
3.1.光学微腔中二次谐波的产生
以SiN为材料的光学微腔为正色散. SiN材料的各项参数如下:
首先, 忽略抽运频率和微腔本征频率之间的频率失谐, 即失谐参量δ1 = 0, 连续抽运光的强度Ein = 1 W1/2. 在非线性效应的作用下, 腔内有二次谐波产生, 二次谐波随着光在腔内循环次数产生的光场分布如图2(a)所示, 可见初始产生的倍频光场较弱, 随着循环次数的增加, 腔内的倍频光场不断变强. 图2(b)和(c)表示腔内基频光和倍频光功率随循环次数不断增加, 并最终稳定于某一值, 此时腔内光场相对稳定, 稳定后基频光和倍频光光场分布如图2(d)和(e)所示, 基频光在腔内形成了一个较宽的脉冲, 且脉冲顶端较为平坦而倍频光场则类似正弦信号, 其功率和基频光相比相对较弱. 图2(f)和(g)是对应的光谱, 基频光和倍频光都在微腔内形成了梳状的光谱, 但倍频光谱宽要小于基频光的谱宽.

Figure2. Second harmonic waves generates inside the SiN microresonator (δ1 = 0, Ein = 1 W1/2): (a) Evolution of the second harmonic waves; (b) curves of the fundamental frequency waves power and (c) the second frequency waves power vary with the round trip number; (d) stationary distribution of the fundamental frequency waves light field and (e) the second frequency waves light field; (f) spectra of the fundamental frequency waves and (g) the second frequency waves.
以上是失谐参数δ1 = 0时腔内光场的演化情况. 下面讨论失谐参数δ1不为0的情况. 假设δ1 = 0.1, 抽运光强Ein = 1 W1/2. 通过求解(1)式和(2)式, 得到微腔内基频光和倍频光功率变化如图3(a)和(b)所示. 对比图2(b)和(c)可以发现, 当微腔的抽运频率和本征频率发生失谐时, 腔内的基频光和倍频光功率在光循环过程中发生振荡, 且经过多次振荡后腔内最终稳定的光功率和图2(b)和(c)中的光功率相比, 强度明显减弱, 说明频率失谐会对腔内的光功率造成明显的损耗. 当腔内光功率稳定后, 腔内基频光和倍频光光场分布随循环次数的演化如图3(c)和(d)所示, 其中任意时刻光场的分布如图3(e)和(f)所示. 可见, 基频光在腔内最终形成了稳定的暗孤子, 而倍频光在腔内的光场仍然近似正弦信号, 但正弦信号的相位随着循环次数而变化.

Figure3. Light field evolution in the microresonator with the frequency detuning (δ1 = 0.1 and Ein = 1 W1/2): (a) Curves of the fundamental frequency waves power and (b) the second frequency waves power vary with the round trip number; evolution of (c) the fundamental frequency waves (d) the second frequency waves after the light fields are stable; stationary distribution of (e) the fundamental frequency waves light field and (f) the second frequency waves light field.
通过以上分析可知, 微腔的频率失谐参数δ1会引起腔内基频光和倍频光功率的损耗, 并导致腔内光场分布发生变化. 因此, 接下来将分析失谐参数δ1对微腔内二次谐波产生的影响.
2
3.2.失谐参量对微腔内二次谐波的影响
假设抽运光场Ein = 1 W1/2不变, 依次增加失谐参数δ1, 得到微腔内基频光和倍频光功率变化曲线如图4所示. 图4(a)和(b)是δ1取不同值时, 基频光在腔内循环过程中的功率变化曲线, 随着失谐参数δ1的不断增加, 基频光的功率在趋于稳定的过程中振荡越来越明显, 当δ1大于0.2后(图4(b)), 功率的振荡变得非常剧烈. 因此, 抽运光频率和微腔本征频率之间的失谐会导致腔内光场分布的不稳定. 图4(c)和(d)是失谐参数δ1对腔内倍频光的影响. 和基频光的变化趋势一致, 倍频光功率在腔内经历先增大再减小的变化过程, 随着δ1的增加, 曲线下降沿的抖动也越来越剧烈, 腔内最终稳定存在的倍频光功率总体呈下降趋势. 可以看到, δ1 = 1时, 微腔内基频光和倍频光的功率变化曲线都出现了周期性的抖动, 腔内光场趋于稳定后, 光功率仍然保持周期性振荡. 图5(a)和(b)给出了腔内基频和倍频光场分布的演化, 光场分布随腔内光场循环的次数进行周期性的变化. 在光场变化周期中, 每隔100次循环选取一个时间点, 画出这些时间点腔内的光场分布的形状, 结果如图5(c)和(d)所示. 由图5(c)可见, 在某些时刻基频光在腔内可形成单脉冲形式的分布(图5(c)②③), 但也仍有脉冲分裂的情况存在(图5(c)④⑤). 图5(d)为倍频光在各个时刻的光场分布情况, 可以看出, 当失谐参量δ1 = 1时倍频光在腔内仍然近似正弦分布, 正弦波的相位随时间不断变化.
Figure4. Influence of the frequency detuning δ1 on the power change curves: (a) Curves of the power variation for the fundamental frequency waves, 0.02 ≤ δ1 ≤ 0.08; (b) curves of the power variation for the fundamental frequency waves, 0.2 ≤ δ1 ≤ 0.8; (c) curves of the power variation for the second harmonic waves, 0.02 ≤ δ1 ≤ 0.08; (d) curves of the power variation for the second harmonic waves, 0.2 ≤ δ1 ≤ 0.8.

Figure5. Stable evolution of the dual light fields when δ1 = 1: (a) Evolution of the fundamental frequency waves; (b) evolution of the second harmonic waves; (c) intensity profiles of the fundamental frequency waves at six different moments within a period, the waveforms are plotted every hundred times; (d) intensity profiles of the second harmonic waves at six different moments within a period, the waveforms are plotted every hundred times.
研究过程中, 我们进一步增大了失谐参量δ1的取值, 结果发现当δ1的值过大时, 腔内的基频光和倍频光的光场均呈现无规则的变化, 如图6(a)和(b)所示. 经过光场在腔内的长时间演化, 微腔中的光场在时间和空间上都没有形成稳定分布, 此时微腔内处于混沌状态[22]. 任取某一时刻, 分别绘制基频光和倍频光的光场分布, 如图6(c)和(d)所示, 图6(e)和(f)分别为其对应的光谱. 可见, 混沌模式的光场表现为光功率随时间无规律地变化, 其相对应的光谱也是无规则的, 每一个梳齿对应的功率也是时刻发生变化的. 产生混沌的原因为失谐参量δ1的值过大, 导致腔内光场的分布不稳定, 在色散以及高阶非线性效应的作用下, 腔内开始随机地出现毛刺以及功率较高的尖峰[23]. 因此, 在利用光学微腔产生二次谐波光梳的过程中, 要尽量避免失谐参量δ1过大, 即抽运频率要尽量接近微腔的本征谐振频率, 以免在腔内产生混沌状态.

Figure6. Chaos inside the microresonator, when the value of detuning parameter is too large: (a) Chaos of the fundamental frequency waves; (b) chaos of the second harmonic waves; (c) intensity profile of the fundamental frequency waves in a moment; (d) intensity profile of the second harmonic waves in a moment; (e) spectrum of the fundamental frequency waves; (f) spectrum of the second harmonic waves.
失谐参量δ1是抽运光频率和微腔谐振频率之间的失谐, δ1的值越大, 表示抽运与微腔谐振之间的频率失谐越严重. 这种失谐会导致腔内光场不同的模式之间发生相互耦合, 因而使腔内基频和倍频光场的不稳定. 特别是当δ1值过大时, 腔内光场的模式发生混乱, 因而产生不规则的光场分布. 根据以上分析, δ1的值越小, 不仅可以有效避免混沌的现象产生, 同时还可获得较强的倍频光功率以及腔内稳定的光场分布.
2
3.3.抽运功率对微腔内二次谐波的影响
微腔的抽运功率同样也会对腔内基频光和倍频光的功率以及光场分布产生影响. 我们保持失谐参量δ1 = 0.1不变, 在此基础上, 不断增加抽运强度Ein, 得到微腔内基频光和倍频光功率变化曲线分别如图7(a)和(b)所示. 可明显看出, 随着抽运功率的不断增大, 腔内最终稳定后光功率也会越大, 且抽运功率的大小几乎不会影响腔内形成稳定光场所消耗的时间. 因此, 为了获得更高功率的倍频光, 可以增强抽运光的功率.
Figure7. Influence of the pump power on the power change curves: (a) Power variation for the fundamental frequency waves; (b) power variation for the second harmonic waves.
然而抽运功率的增大也会导致腔内光场的不稳定. 将抽运光场强度Ein增大到100 W1/2, 腔内基频光和倍频光功率变化曲线如图8(a)和(b)所示, 可见, 腔内的光功率经过大幅振荡衰减后趋于稳定, 稳定时期的光功率仍然保持周期性的小幅振荡. 当光功率稳定振荡后, 基频光和倍频光的光场演化如图8(c)和(d)所示, 光场分布同样呈周期性变化, 并且可以认为在任意时刻, 基频光和倍频光的光场在腔内是均匀分布的.

Figure8. Evolution of the light field in the microresonator at Ein = 100 W1/2: (a) Curve of power variation for the fundamental frequency waves; (b) curve of power variation for the second harmonic waves; (c) periodic evolution of the fundamental frequency waves; (d) periodic evolution of the second harmonic waves.
在此基础上进一步增大抽运光强度, 抽运光场Ein分别为800, 1000和1200 W1/2时, 基频光和倍频光功率变化曲线如图9所示. 当Ein = 800 W1/2时, 基频光的功率在腔内经过功率振荡后, 最终大约稳定于20 kW, 而倍频光功率在循环过程中增加, 达到峰值后迅速衰减, 然后趋向于稳定, 稳定部分的光功率也同样保持周期性振荡. 这是由于抽运光功率过高, 导致腔内二阶非线性效应明显, 较多的基频光转化为倍频光, 因而腔内累积的倍频光功率越来越大, 因此光在腔内循环过程中, 倍频光的功率有明显的增长. 当倍频光功率达到一定值时, 由频率失谐、色散、损耗引起的倍频光功率降低, 超过了由非线性效应产生的倍频光强度, 因而倍频光在功率增大到一定程度后会出现功率迅速衰减的现象, 所以在图9(b)中倍频光的功率会出现一个脉冲形式的变化曲线. 由于倍频光是由基频光通过非线性效应产生的, 因而在图9(a)中对应了一个向下的脉冲形状. 最终腔内形成的稳定光场, 基频光和倍频光功率都会在一定范围内等幅振荡 (如图9(b)中放大部分的曲线所示). 当Ein = 1000 W1/2和Ein = 1200 W1/2时, 腔内功率变化曲线如图9(c)—(f)所示. 抽运功率增加, 在图9(d)和(f)中, 倍频光功率变化曲线中出现了两个脉冲, 并且抽运光功率越强, 曲线中第二个脉冲的强度更大, 因此对应的基频光功率振荡也更加明显. 从LLE的(1)式和(2)式可以看出, 非线性效应的强弱和腔内光功率的大小是相关的, 腔内光功率越大, 由自相位调制和交叉相位调制引起的倍频效率也越高. 因此, 当抽运功率较弱时, 腔内基频和倍频光强也处于相对较弱的水平, 增益、损耗、色散、非线性之间相互作用可以达到相对平衡的状态, 最终基频和倍频光场均会趋于稳定. 而当抽运功率较强时, 腔内增益和非线性均表现为较强的作用, 此时腔内光场难以达到稳定状态. 因此, 光学微腔抽运功率过大, 会导致腔内光场处于不稳定的状态. 为保证腔内光场和光功率的稳定, 需要避免抽运功率过强的情况.

Figure9. Curves of the power variation for the fundamental frequency waves and the second harmonic waves: Power variation for (a) the fundamental frequency waves and (b) the second harmonic waves at Ein = 800 W1/2; power variation for (c) the fundamental frequency waves and (d) the second harmonic waves at Ein = 1000 W1/2; power variation for (e) the fundamental frequency waves and (f) the second harmonic waves at Ein = 1200 W1/2.
2
3.4.微腔大小对腔内光场的影响
以上都是在微腔FSR = 231.3 GHz (对应的微腔半径为100 μm)的条件下进行的讨论. 但是研究中发现, 微腔半径的大小同样会对腔内光场产生影响. 设微腔FSR = 300 GHz, δ1 = 0.1, Ein = 100 W1/2, 得到腔内基频光和倍频光功率的随时间变换, 最终都收敛于某一值, 如图10(a), (b)所示, 即在围墙内可以得到稳定的基频、倍频的光场分布. 当光功率稳定后, 光场演化情况如图10(c), (d)所示, 由于微腔半径增大, 在腔内基频光与倍频光之间非线性作用的距离变长, 在较强的非线性作用下, 基频光和倍频光最终在腔内演化出16个脉冲, 在对应的光谱图中(图10(e), (f)), 每隔16 FSR, 就会有一个功率较大的梳齿, 此时微腔工作于“图灵环”模式(turing patterns)[17].
Figure10. Turning patterns in the microresonator, when FSR = 300 GHz (Ein = 100 W1/2, δ1 = 0.1): (a) Curves of the power variation for the fundamental frequency waves; (b) curves of the power variation for the second harmonic waves; (c) evolution of the fundamental frequency waves; (d) evolution of the second harmonic waves; (e)spectra of the fundamental harmonic waves; (f) spectra of the second harmonic waves.
研究过程中发现, 可工作于“图灵环”模式的微腔FSR的取值范围较小. 将微腔的FSR减小为298 GHz, 腔内最终稳定的基频光场分布如图11(a)所示, 基频光的光场在腔内仍然为多脉冲的形式, 但脉冲的强度受到了调制, 这种调制在倍频光的光场分布中更为明显 (如图11(b)所示), 多脉冲的光场受到正弦信号的调制. 由于光场中包含了多个频率成分, 因此其光谱图中噪声较多, 在主梳齿附近会出现很多的毛刺. 进一步减小微腔的FSR, 当减小为295 GHz时, “图灵环”工作模式消失, 基频光和倍频光的光场均演化为正弦分布, 且正弦信号的相位随时间变化, 光谱也转化为很窄的梳状光谱(如图12所示).

Figure11. Stable intensity profile and spectra of the fundamental frequency waves for FSR = 298 GHz: (a) Intensity profile of the fundamental frequency waves; (c) spectrum of the fundamental frequency waves; (b) intensity profile of second harmonic waves; (d) spectrum of second harmonic waves.

Figure12. Stable intensity profile and spectra of the fundamental frequency waves for FSR = 295 GH: (a) Evolution of the fundamental frequency waves; (b) evolution of the second frequency waves; (c) intensity profile of the fundamental frequency waves; (d) spectrum of the fundamental frequency waves; (e) intensity profile of the second harmonic waves; (f) spectrum of the second harmonic waves.
研究结果表明, 当抽运频率和微腔谐振频率之间没有发生频率失谐, 即失谐参量δ1 = 0.1时, 腔内的基频光和倍频光功率不断增强, 并最终达到稳定, 稳定后的基频光场为平顶脉冲的形式, 而倍频光场呈正弦分布. 当抽运频率和谐振频率之间发生频率失谐后, 腔内基频和倍频光功率在演化过程中出现振荡, 但最终稳定于某一值. 稳定后的光场也发生变化, 基频光场变为暗孤子的形式, 倍频光场仍然为正弦分布. 失谐参量δ1增加, 将导致腔内基频和倍频光功率在演化过程中振荡变强, 且最终稳定的光功率变弱, 稳定后的光场分布为周期性变化. δ1进一步增加会使得微腔光场处于混沌状态. 因此, 为保证微腔中二次谐波的产生以及腔内光场的稳定, 实验中应尽量减小抽运频率与微腔振荡频率之间的频率失谐.
当微腔抽运功率较小时, 增大抽运功率可有效增强腔内基频和倍频光场的强度; 而当抽运光强较强时, 会导致腔内色散以及非线性效应过强, 最终稳定的光场仍然呈周期性变化, 且抽运光功率越强, 光功率的演化曲线振荡越强. 因此, 为保证腔内光场和光功率的稳定, 需要避免抽运功率过强的情况.
光学微腔的尺寸同样会对腔内光场分布产生影响. 当选取微腔FSR = 300 GHz, 失谐参量δ1 = 0.1, 抽运强度Ein = 100 W1/2, 微腔工作于“图灵环”状态, 腔内基频和倍频光场均演化出稳定的16个脉冲, 相对应光谱每隔16 FSR, 就会有一个功率较大的梳齿. 当FSR减小后, 腔内多脉冲形式的光场受到了正弦信号的调制, 光谱中也出现较多的毛刺. 因此, 选择合适的微腔尺寸, 可以有效选择微腔的工作模式, 得到需要的光谱. 该理论分析结果对于实验中微腔参数和抽运参数的选取有重要意义. 此外, 本文中采用的微腔光场演化模型是二维模型, 忽略了光场在垂直方向上(y方向)的分布. 为了更精确地分析腔光场的分布和演化情况, 需要考虑光场在垂直方向分布情况时, 因而要引入三维的光场演化模型, 此时抽运光场是关于y的表达式, 腔内基频和倍频光场也和y有关. 关于微腔内基频和倍频光场三维模型的建立和分析, 将在后续的工作中进一步展开.