作者简介:刘 强(1983—),男,深圳大学助理教授、博士.研究方向:太赫兹器件.E-mail:qliu@szu.edu.cn
中文责编:方 圆; 英文责编:溯 心
1)深圳大学电子与信息工程学院,广东深圳 518060; 2)深圳大学太赫兹技术研究中心,广东深圳 518060; 3)深圳大学光电子器件与系统教育部/广东省重点实验室,广东深圳 518060
1)College of Electronic and Information Engineering, Shenzhen University, Shenzhen 518060, Guangdong Province, P.R.China2)THz Technical Research Center, Shenzhen University, Shenzhen 518060, Guangdong Province, P.R.China3)Key Laboratory of Optoelectronics Devices and Systems of Ministry of Education and Guangdong Province, Shenzhen University, Shenzhen 518060, Guangdong Province, P.R.China
nonlinear optics; optical frequency-difference effect; whispering gallery mode; surface plasmonic mode; sub-wavelength confinement; terahertz wave generation
DOI: 10.3724/SP.J.1249.2019.02140
提出并研究基于回音壁共振模和表面等离子激元共振模的混合微腔太赫兹产生系统,运用有限元方法模拟实现了室温下11.7 THz的太赫兹波产生,所获得的太赫兹波模体积小,超过衍射极限,亚波长局域使得Purcell因子达到7.2×103(砷化镓中太赫兹波长的负三次方). 研究表明,该结构的太赫兹波转化效率可达到量子极限(~5.5%),太赫兹波的最高输出功率可达166.3 mW. 该系统可在室温下工作,转化效率接近量子极限,太赫兹输出功率较高,且易于光学集成.
By utilizing the nonlinear optical frequency-difference effect, an efficient terahertz(THz)wave generator based on a hybrid micro-cavity of whispering gallery(WG)mode and surface plasmonic mode is proposed and investigated. Simulations based on finite element method are performed. It is found that the THz wave at the frequency 11.7 THz can be produced at room temperature, and the terahertz modal volume is small, which breaks the diffraction limit, the sub-wavelength confinement in turn produces a high Purcell factor up to 7.2×103(one in cubic of the THz wavelength in GaAs). Further studies show that the conversion efficiency can reach the quantum limit(~5.5%), and the maximum output of the generated THz wave is as high as 166.3 mW. The advantages of the design include: capability for operation at room temperature, high conversion efficiency that approaches to the quantum limit, high output power for the THz wave, and convenience for all-optical integration.
理论上,通过两束光的非线性光学差频可在室温下产生太赫兹波[1-2]. 采用三模式共振腔系统可实现光与非线性介质的强相互作用,提高太赫兹波的转化效率. 其中,高品质因子Q谐振腔是一个重要因素,因为Q值越大,光学模式场和非线性介质的相互作用时间越长,从而可以提高太赫兹波的转化效率. 另一方面,非线性混频理论表明,太赫兹波的转化效率与Purcell因子(Purcell因子=Q/Vm, Vm为模体积)成正比[2],因此,也需要高品质因子谐振腔. 但Q值越大,谐振模式的谱线宽度越窄,系统工作特性对制作加工过程中的误差越敏感. 文献[3]提出一种通过调节环形腔中部分金属尖端的高度,实现谐振频率移动的改良办法,然而在微系统集成方面该方法仍具有一定局限性.
基于表面等离激元共振的微腔和波导结构可能为上述问题的解决带来新思路[4-5], 因为这种微腔的模体积非常小,能达到超越衍射极限的深亚波长尺寸. 由此,即使一个Q值较低的表面等离激元微腔,也能实现大的Purcell因子,这样就对制作工艺引起的工作性能误差有较大的空间容忍度. 例如,文献[4]报道一种基于表面等离激元共振的纳米腔系统,尽管其光学Q值仅约330,但由于其深亚波长结构的模体积仅约为5.6×10-3 μm3,所获得的Purcell因子高达60 000(λ/nair)-3. 其中, λ为光波长; nair为空气的折射率.
本研究提出一种基于回音壁共振模与表面等离激元共振模的混合微腔太赫兹产生系统,并利用COMSOL软件[6-7]模拟实现了高效率太赫兹光源.
基于回音壁共振模与表面等离激元共振模的混合微腔太赫兹产生系统示于图1. 其中,绿色部分为本征砷化镓; 灰色部分为空气间隙; 蓝色部分为掺杂砷化镓. 该系统支持三模式共振:红外泵浦光模式(λ1=1 462.7 nm)、红外辅助泵浦光模式(λ2=1 551.0 nm)以及太赫兹波谐振模式,故称之为混合微腔结构. 当满足光学差频效应中的相位匹配条件时,可输出频率为11.7 THz的电磁波(λTHz=25.7 μm). 值得注意的是,所述的红外泵浦光模式和辅助泵浦光模式通过图1的直波导耦合进入环形腔. 为使环形腔在垂直方向上对太赫兹波产生局域作用,选取绿色环形腔的高度H1=5 μm,近似为GaAs介质中太赫兹波的半波长.
模拟过程中的具体结构参数如下:环形腔外环半径Rex=9 484 nm, 内环半径Rin=1 490 nm,掺杂砷化镓环对应的半径Rd=1 388 nm,掺杂区与本征半导体之间的空气间隙g=102 nm. 注意到环形腔外的直波导与环形腔之间也有一个空气缝隙,其作用是耦合输入与输出,此缝隙宽为w=120 nm,直波导的宽度W=1 μm. 此外,混合微腔系统的厚度H1=5 μm,为节约计算时间,衬底(图1中黄色部分)的厚度H2设为200 nm,衬底半径假设为无限大. 另外,图1中直波导的下(上)端口为输入(输出)端口,计算区域的其他部分均为空气; 为避免光在边界处发生反射,除端口外的边界均设置为理想匹配层(perfectly matched layer, PML).
模拟计算表明,太赫兹波以亚波长局域的方式被限制在混合微腔系统中,太赫兹模体积仅为2.71 μm3或4.3×10-3(λTHz/nGaAs)3. 其中, nGaAs=3为11.7 THz频率处本征GaAs的折射率. 注意到本征GaAs的介电常数实部为9,而虚部仅为实部的0.18%,因此在模拟分析中忽略了介电常数的虚部. 计算得到的太赫兹波谐振Q值为31,Purcell因子为Q/Vm= 7.2× 103(λTHz/nGaAs)-3.光学差频的本质是二阶非线性效应,数学上可引入极化向量P(2)来描述
P(2)=(d11 d12 d13 d14 d15 d16
d21 d22 d23 d24 d25 d26
d31 d32 d33 d34 d35 d36)×(E2x
E2y
E2z
EyEz
ExEz
ExEy)(1)
其中, dij(i=1, 2, 3, j=1,2,3,…,6)是极化率张量矩阵元,由非线性材料决定; Ex、 Ey及Ez分别为x, y及z方向的电场强度. GaAs属于4-3m群, d14=d25=d36=d=χ/2=274 pm/V[8],其他张量矩阵元素均为0, χ为二阶极化系数.
在x、 y、 z方向上的极化分量为
{P(THz)x=d·EyEz
P(THz)y=d·ExEz
P(THz)z=d·ExEy(2)
根据泵浦光TE(只有z方向的电场, x和y方向都没有电场分布)和TM(只有z方向的磁场, x和y方向都没有磁场分布)两种模式,方程(2)所给出的极化分量可得到进一步简化.
考虑到需要在系统中局域太赫兹波,本研究选择利用在太赫兹频段具有负介电常数的材料,这里选用掺杂GaAs.
根据Drude模型,掺杂GaAs材料的介电常数[9]为
ε(ω)=ε(∞)+(iσ)/(ωε0)(3)
其中, ε(∞)=10.89, 是砷化镓在频率无穷大时的介电常数[9]; ε0≈8.854×10-12 F/m 是真空的介电常数; σ是复电导率,且
σ=(106 Ne2 τ/m*)/(1-iωTHzτ)(4)
其中, ωTHz为太赫兹波的角频率; N为GaAs中的载流子浓度(单位:cm-3); e=1.6 ×10-19 C,是电子电荷量; m*=1.39×10-32 kg,为载流子有效质量[10]; τ为载流子的弛豫时间,且
τ=m*μm/e(5)
对于n型掺杂的砷化镓,载流子的移动率μm可以近似表示为[10]
μm=(0.94)/(1.0+(10-17 N)1/2)(6)
图2给出两种n型掺杂浓度下砷化镓介电常数的实部与虚部. 考虑到本研究所模拟太赫兹波的频率集中在11.7 THz附近,因此需要找到较为合适的掺杂浓度以及相对应的砷化镓介电常数εT(如较大的实部绝对值,较小的虚部). 以下的场分布模拟实验发现,在11.7 THz附近的优化掺杂浓度为N=8.561×1018 cm-3,对应的εT = -81.3+148.4i.
为对环形微腔进行表征,以往研究分别采用了两个基本整数:角节点数mi和径向节点数pi(i=1,2,3, 模式序数1和2分别代表泵浦光和辅助泵浦光,3代表THz波谐振模式). 对于光学差频(difference frequency generation, DFG),环形微腔中的角节点数mi应该满足相位匹配条件[1-2]
Δm=m1-m2-m3=±2(7)
DFG过程遵守能量守恒定律
ωTHz=ω1-ω2(8)
太赫兹波的转化效率不超过量子极限效率ηlim, 即
ηTHz≤ηlim(9)
量子极限效率ηlim定义为
ηlim=ΓTΓ1ωTHz/ω1(10)
其中, Γ1(ΓT)代表泵浦光(太赫兹波)的耦合效率,均为从直波导到微腔之间的耦合效率.
量子极限效率下的临界泵浦光功率Pcr为[11]
Pcr=(ω1)/(4|β|2Q1Q2QTΓ1)(11)
其中, β为模式交叠因子,定义为
β=1/4×
(∫d3r∑i, j, kε0 χ(2)ijkE*1,i(E2, jET, k+ET, jE2, k))/((∫d3r·ε1|E1|2)1/2(∫d3r·ε2|E2|2)1/2(∫d3r·εT|ET|2)1/2)(12)
其中, |β|2单位是J-1; Q1, Q2 和QT 分别表示混合微腔系统中泵浦光、辅助泵浦光和太赫兹波的谐振Q值. 由于临界泵浦光的功率与模式交叠以及3个模式的Q值之积成反比,因而需要提高模式交叠和3个模式的光学Q值,以获得较小的临界泵浦光输入功率,可在单光子探测方面获得应用[12].
在量子极限转化条件下,即转化效率达到式(10)的定义值时,根据耦合模理论(coupled mode theory)[8],输入辅助泵浦光的振幅s+2与输入泵浦光的振幅s+1满足
|s+2|2=(ω2|1-4/(ω1)|β|2Q1Q2QTΓ1|s+1|2|2)/(16|β|2Q1Q2QTΓ2)(13)
即
(|s+2|2)/(Pcr)=(ω2)/(4ω1)(Γ1)/(Γ2)|1-(|s+1|2)/(Pcr)|2(14)
由式(14)可见,当输入泵浦光的功率|s+1|2等于临界功率Pcr时,输入的辅助泵浦光的功率|s+2|2为0,即泵浦光完全转化为太赫兹波. 此现象称为泵浦光的完全下转化(complete down-conversion). 在接下来有关太赫兹功率输出的讨论中,需协同考虑辅助泵浦光、太赫兹波的上转化(up-conversion)以及泵浦光的下转化(down-conversion)两种情况.
为增强非线性效应,获得高的非线性DFG效果,需要提高Q值和模式交叠因子值β. 首先,需对泵浦光和辅助泵浦光这两个光学模式分别实现高Q值场局域. 其次,通过调整掺杂区和本征半导体之间的空气间隙宽度,使泵浦光、辅助泵浦光和太赫兹波同时满足相位匹配条件式(7)和能量守恒条件式(8),从而产生太赫兹波输出. 其中,THz波在系统中通过表面等离激元共振而获得亚波长局域,从而实现了高模式交叠因子β. 数值计算表明,在量子极限转换效率下,由式(11)所需的泵浦光功率为875.5 mW.
为显示泵浦光、辅助泵浦光和太赫兹波之间的模式交叠,分别研究泵浦光、辅助泵浦光和太赫兹波在混合微腔中的场分布情况. 图3给出泵浦光(TE模式)在耦合波导及环形腔中的电场分布情况,此处λ1=1 462.7 nm,角节点数m1=119, 径向节点数p1= 4. 值得注意的是,这里的输入端位于直波导的下端口,输出端为直波导的上端口. 模拟过程中为避免光在边界处发生反射,除端口外的边界均设置为PML.
图3 TE模式下泵浦光的电场E1,z分布图
Fig.3 Electric field profile(E1,z)of the pump at TE mode
图5给出太赫兹波(TM模式)在耦合波导及环形腔中的磁场分布情况,此处频率fT=11.7 THz,角节点数mT=1, 径向节点数p3=3.
图4 TM模式辅助泵浦光的磁场H2,z分布图
Fig.4 Magnetic field profile(H2,z)of the assisting pump at TM mode
图5 TM模式下太赫兹波磁场HT,z分布图
Fig.5 Magnetic field profile(HT,z)of THz at TM mode
通过式(12)可知, E*1,z(E2,xET,y+ET,xE2,y)在三维空间的积分正比于模式交叠因子β. 图6为β在空间的分布情况,计算显示在相位匹配条件下,|β|2=16 J-1. 其次, β的最大值出现在靠近环形腔的边缘区域,而旁边的直波导中没有分布. 我们还注意到,在本征GaAs与掺杂GaAs之间的狭窄空气间隙区域(g=102 nm), β≈0, 这是因为一般的光学回音壁模式无法存在于远离边缘的区域.
考虑到太赫兹波长远大于波导的宽度,因此,以下关于THz波电场能量密度分布的数值模拟中,忽略波导处的THz波电场分布. 为预测太赫兹波在混合腔中的品质因子QT, 在如图7所示的4个边界处都添加了PML.
图7 THz波电场能量密度的横截面侧视图(切平面过环形腔对称轴,且与z轴平行)
Fig.7 Cross-sectional side view of THz mode profile in terms of electric energy density with the symmetry axis of the ring cavity on the cutting plane and parallel to the z axis
计算得QT=31, 模体积超越衍射极限,仅为4.3×10-3(λTHz/nGaAs)3. 值得注意的是,通过引入掺杂砷化镓以及本征砷化镓与掺杂砷化镓区域之间的狭窄空气间隙(g=102 nm, 见图1中灰色区域),太赫兹波被限制在具有亚波长尺度的微腔中(如图1中灰色区域,或见图7中该狭窄空气间隙中的电场分布).
假设泵浦光(ω1)和辅助泵浦光(ω2)的输入功率相等. 计算显示,泵浦光、辅助泵浦光和太赫兹波的谐振Q值分别为Q1=4.7×105、 Q2=1.6×106和QT=31. 波导与微腔之间,泵浦光和辅助泵浦光的耦合效率分别为Γ1=0.963 7和Γ2=0.849 1. 考虑到太赫兹波的自由空间波长远大于混合腔的几何尺寸,且当QT较小时,太赫兹波更易于逃逸出腔外,使得太赫兹输出功率损耗仅为2.6%,这里假设太赫兹波的输出耦合效率为ΓT=1.
通过耦合模理论[8]分析,由式(11)~式(14)可以计算得到图8和图9,其分别为基于DFG的太赫兹波输出功率与泵浦光输入功率的关系曲线. 注意到图8中太赫兹波输出功率存在一个峰值,之后输出呈现饱和. 由式(10)可计算得到量子极限转化效率ηlim≈5.5%,对应的临界输入功率为Pcr=875.5 mW; 此情况(如图8中虚线)对应的THz输出功率为48.2 mW. THz波的最大输出功率为0.19Pcr(或166.3 mW),对应的泵浦光的输入功率约为6Pcr,太赫兹转化效率下降为3.2%.
图8 输出THz波功率与输入泵浦光的功率之间的关系
Fig.8 The relationship between the normalized output power of THz wave and the input pump power
图9 小信号输入条件(PP<<Pcr)下,输出THz波功率PT与输入泵浦光功率平方之间的关系.
Fig.9 The relationship between the output power of THz wave and the input power of pumping wave when PP<<Pcr
最后,与同样采用环形腔进行三模耦合的文献[14]不同,本研究所探讨的混合腔可通过调节n掺杂浓度改变有效折射率,进而实现混合腔中THz谐振频率的变化,因而在频率调谐方面提供了另一种自由度; 其次,本研究直观显示出三模DFG发生的区域,即模式交叠区域; 最后,本研究对输出THz波功率与输入泵浦光的功率之间的关系,特别是饱和信号和小信号输入条件下的关系作了较深入分析.
本研究提出并研究一种利用DFG效应的基于回音壁共振模和表面等离激元共振模的混合微腔太赫兹产生系统,旨在实现常温下的高量子效率太赫兹光源. 应用相位匹配条件和Drude模型设计泵浦光、辅助泵浦光和太赫兹光三模共振混合微腔结构. 所获得的太赫兹模体积远超衍射极限,使Purcell因子Q/Vm高达约7.2×103(λTHz/nGaAs)-3. 数值结果表明,该结构可将THz波转化效率提高到量子极限(输入功率875.5 mW时, PT/Pp=5.5%, THz波输出功率最高达166.3 mW). 此外,还分别针对饱和信号及小信号输入情况下的输出功率进行讨论. 高Purcell因子、高转换效率以及较高的THz输出功率,结合室温工作方式、易于实现光学集成和制造等优点,使该太赫兹波产生系统具有潜在应用价值.
致 谢:衷心感谢美国诺福克州立大学Sacharia Albin教授的悉心指导!
深圳大学学报理工版
JOURNAL OF SHENZHEN UNIVERSITY SCIENCE AND ENGINEERING
(1984年创刊 双月刊)
主 管 深圳大学
主 办 深圳大学
编辑出版 深圳大学学报理工版编辑部
主 编 李清泉
国内发行 深圳市邮电局
国外发行 中国国际图书贸易集团有限公司(北京399信箱)
地 址 北京东黄城根北街16号
邮 编 100717
电 话 0755-26732266
0755-26538306
Email journal@szu.edu.cn
标准刊号 ISSN 1000-2618
CN 44-1401/N