激光与光电子学进展, 2019, 56 (19): 192302, 网络出版: 2019-10-14   

余弦-高斯光束在太赫兹平行平板波导中的传输研究 下载: 1064次

Cos-Gaussian Beams Propagating Inside Terahertz Parallel-Plate Waveguides
作者单位
1 上海大学物理系,上海 200444
2 郑州师范学院物理与电子工程学院, 河南 郑州 450044
摘要
基于平行平板波导的TE1模式理论和电场是近似线偏振的假设,推导了余弦-高斯光束在平行平板波导中的电场传输表达式。在傍轴近似的条件下,余弦-高斯光束可以作为平行平板波导的一个基础解。当频率为0.35 THz时,对余弦-高斯光束在金属铜平行平板波导中的传输进行了数值模拟,结果表明模拟的幅值和相位与根据电场表达式所得结果符合得很好,同时也讨论了z轴方向上模拟的相位与理论结果的偏差和电场散度方程的完备性。
Abstract
Herein, based on the TE1 mode theory of the parallel-plate waveguide and assumption of approximately linearly polarized electric fields, we derive the expression of electric-field propagation of a cos-Gaussian beam in the parallel-plate waveguide. Under the paraxial approximation, the cos-Gaussian beam can be used as a fundamental solution for parallel-plate waveguides. The derived beam is tested by numerical simulation in copper parallel-plate waveguides at a frequency of 0.35 THz. Results show that the simulated amplitude and phase are in good agreement with those from the analytical expression of electric field. The deviation between the simulated phase and theoretical results at the z-axis and completeness of the electric-field divergence equation are also discussed.

1 引言

近年来,频率在0.1~10 THz之间的太赫兹(THz)波受到了广泛关注[1-2],因为其在成像技术[3]、红外传感[4]、通信领域[5-6]和纳聚焦[7-11]方面均有大量的应用,而太赫兹波导[12-16]作为传输太赫兹波的重要器件,在过去的十多年间,也已在理论和实验方面得到了广泛研究。常见的太赫兹波导包括金属腔波导[17-18]、金属平行平板波导(PPWG)[16,19-20]、表面等离激元波导[21-23]和电介质波导[24-26]等。其中PPWG具有低损耗和单TEM(transverse electromagnetic )模传输时几乎没有群速度色散的优点[27],已经被大量应用。例如PPWG可以作为Bragg谐振腔[28],也可以被应用到光谱学[29]的研究中。此外,还有团队利用开槽的PPWG作为漏波天线对太赫兹波进行调制[6]。研究还表明,随着频率增大,PPWG的最低阶横电模式场的衰减逐渐减小[19]

众所周知,来自激光源的普通太赫兹波是呈高斯分布的,并且几乎是线性偏振的。然而,各种太赫兹波导的模态场分布不是准高斯分布,激光源产生的太赫兹波与传统的横电模、横磁模不匹配[18]。因此,激光源和波导模式之间的耦合效率不是很理想。

近年来,余弦-高斯光束[30-34]逐渐进入人们的视野,其在扰动大气中的传输[30-31]以及在介质中的传输[32-33]均有研究。研究表明,经过自然衍射后的余弦-高斯光束最终转变为双曲余弦-高斯光束[30,34],这意味着余弦-高斯光束和双曲余弦-高斯光束在传播过程中是具有互易性的。有人对余弦-高斯光束在Kerr介质[35]和光折变晶体[36]中的传输进行了研究。在之前的工作中,多切片介质填充的PPWG[18]中的余弦-高斯模式已在理论上得到了验证,而余弦-高斯光束在受限金属波导结构中的传输特性尚未得到深入研究。基于此,本文对余弦-高斯光束在波导结构中的传播进行了研究。

2 余弦-高斯光束的分析

图1所示为PPWG的空间结构,坐标轴位于几何中心,此波导在x轴和传输方向z轴上是无限延伸的,两个平板之间的间距为b。此时考虑将边界设置为完美电导体,并且设置两板间填充物为空气。

图 1. 无限长PPWG结构示意图

Fig. 1. Schematic of infinite-length PPWG

下载图片 查看所有图片

首先假设其波导模式只存在一个x方向的电场分量,有如下形式

E=Exex=ψ(x,y,z)exp(-ikzzex,(1)

式中:E为总电场;Ex为假设的唯一电场分量;ψ为慢变振幅复函数;ex为x方向上的单位矢量;kz为传播常数,满足

kz=k02-πb2,(2)

其中k0=2π/λ0,为自由空间中的传播常数,λ0为自由空间中的波长。在PPWG中的波导模式的电场分量应该满足标量亥姆霍兹方程[34]:

2Ex+k02Ex=0(3)

将(1)式所描述的电场方程代入(3)式,可以得到

2ψx2+2ψy2-2ikzψz+πb2ψ=0(4)

由于假设ψ是沿z轴的慢变振幅复函数,所以二阶偏微分项∂2ψ/z2可以被忽略。

根据经典的波导理论[37]可知,PPWG的TE1模的电场分量Exy方向(即受限方向)上呈余弦分布。而在傍轴近似的条件下,在自由空间下标量亥姆霍兹方程的基模解是二维的高斯函数[38-39]。因此,假设标量函数ψy方向受到限制时,(4)式的解在y方向上呈余弦分布,而PPWG在x方向是无限延伸的,标量函数ψx方向上呈一维高斯分布。然后,可以得到试探解,即

ψ=exp-ip(z)-ikzx22q(z)cosπby,(5)

可以称之为余弦-高斯光束,式中:p(z)为与传输距离有关的复数,q(z)为余弦-高斯光束的复光束参数。接下来,将(5)式代入傍轴(4)式,通过对比系数可以得到

-ikzq(z)-2kzp'(z)×x0=0,(6)-kz2q2z-kz21q(z)'×x2=0,(7)

式中:p'(z)为p(z)函数关于坐标变量z的一阶偏导;x0为常系数。因为在xoy平面上的任意一点都要满足(6)式和(7)式,所以它们的系数必须同时为0,则可以得到

q(z)=q0+z,(8)p(z)=-i2lnzq0+1,(9)

式中:q0=iz0;z0w02/λ0,被称为瑞利距离;w0x方向上的光腰宽度。

通过对比自由空间的高斯光束,可以进一步定义光斑半径w(z)、曲率半径R(z)和相位因子φ(z),得到余弦-高斯光束的表达式。值得一提的是,余弦-高斯光束只在x方向上满足高斯分布。接着将(5)式、(8)式和(9)式代入(2)式,可以得到余弦-高斯光束的横电场分量Ex:

Ex(x,y,z)=cosπbyw(z)w0exp·[-ikzz+i2φ(z)-ikzx22R(z)-kzx2k0w2z],(10)

其中,

w(z)=w01+z2z02,(11)R(z)=z1+z02z2,(12)φ(z)=arctanzz0(13)

正如所预测的,余弦-高斯光束表达式和自由空间的高斯光束有类似的数学形式,但也有所不同,例如余弦-高斯光束在PPWG中y方向受到限制而呈现余弦形式,其振幅和相位表达式也与自由空间的高斯光束不同。特别地,(10)式中φ(z)的系数不同于普通的高斯光束,可能是因为余弦-高斯光束在x方向上只是呈一维的高斯分布而不是呈二维的高斯分布。

3 模拟分析

为了验证上述理论的可行性,利用商业有限元软件COMSOL进行数值模拟分析。在这里,对余弦-高斯光束在PPWG中的传输与 (10) 式的结果进行对比分析。

在模拟过程中,将光束频率设为0.35 THz,其在真空中的波长约为0.86 mm,光腰宽度w0=3λ0,平板间间距b为0.85 mm。模拟中,采用金属铜的阻抗边界条件,在远红外和太赫兹波段,金属电介质常数εm可以通过Drude模型[40]得到,金属铜在0.35 THz下的介电常数为-6.4649×105+4.0535 i×106。虽然PPWG在xz方向上是无限延伸的,但是考虑到模拟软件的性能,在x方向上截取了24λ0作为计算窗口,z方向上截取了50λ0作为传输距离,有限长PPWG的三维结构如图2所示。

图 2. 有限长PPWG的三维结构

Fig. 2. Three-dimensional structure of finite-length PPWG

下载图片 查看所有图片

三维数值模拟的场分布如图3(a)所示,可以看到场分布是对称的。另外,也可以得出这样的结论: 当传输距离小于30λ0时,余弦-高斯光束的能量集中在z轴附近。从图3(a)和图3(b)可以看出,随着传输距离的增大,光斑半径在x方向上逐渐增大,这意味着光束在传输过程中能量逐渐分散。同时,电场受y方向边界的严格限制,场分布呈余弦形式。

图 3. 三维数值模拟的场分布图。(a) 0.35 THz下余弦-高斯光束在有限长PPWG中传输的场分布;(b)不同传输距离处横截面的电场强度分布图

Fig. 3. Field distribution in three-dimensional numerical simulation. (a) Field distribution of cos-Gaussian beam propagating in finite-length PPWG at frequency of 0.35 THz; (b) electric-field intensity distribution of cross-sections at different propagating distances

下载图片 查看所有图片

接下来,评估 (10) 式和数值模拟结果在传输过程中幅值和相位的一致性程度。获取的不同截面(xoy平面)的中心电场数据(y=0和z=10λ0, 20λ0, 30λ0, 40λ0, 50λ0),如图4(a)所示。可以看出,模拟电场|Ex_s| (不同的标记)与(10)式计算得到的电场|Ex|(实线)吻合得较好,这两组数据几乎重叠,无法区分。这种极小的偏差是因为(10)式只是一个近似解表达式,而不是严格的解析解。为了描述(10)式的准确性,采用方均根误差(RMSE)来表示模拟数据与公式曲线的符合程度,如图4(b)所示,并将其定义为

ERMSE=i=m(Ex_s,i-Ex,i)2m,(14)

式中:|Ex_s,i|表示图4(a)中不同标记符号对应的数据点;|Ex,i|表示相应实线的数据点;m表示总数据点个数。RMSE的值越小,表示两条曲线符合得越好。从图4(b)中可以看出(10)式的精度是令人满意的。

图 4. 振幅和相位的分析结果。 (a)不同z值下,当y=0 时x方向上的归一化电场强度分布;(b)理论与模拟值之间的RMSE;(c)不同z值下,当y=0 时x方向上的相位分布;(d) z轴处理论和模拟的相位结果对比

Fig. 4. Analysis results of amplitude and phase. (a) Normalized electric-field intensity distribution in x direction for y=0 at different z values; (b) RMSE between theoretical and simulated results; (c) phase distribution of electric field in x direction for y=0 at different z values; (d) comparison between theoretical and simulated results of phase distribution in z-axis

下载图片 查看所有图片

对于相位信息,利用模拟软件自带的arg( )函数来提取相位。当传输距离小于30λ0时,计算得到的窗口的边缘相位曲线会呈现不规则的趋势。这是因为此时光斑半径比计算窗口小很多,窗口边缘的场值几乎为零,这些相位数据便没有意义,可以被忽略。传输到一定距离后,模拟相位曲线与(10)式中的相位吻合得较好。从图4(c)也可看出,等相位面不是一个平面,而是一个曲面。另外在图4(d)中给出了z轴的模拟相位数据(x=y=0),并与(10)式的计算结果进行对比,实线代表(10)式的计算结果结果,圆点代表仿真模拟结果,数据和曲线吻合得很好。

4 讨论

模拟了1.0 THz和3.5 THz下余弦-高斯光束在铜制PPWG中的传输,同时将波导结构缩小相应的尺寸,即改变板间间距b,以保证余弦-高斯光束是单模传输。在1.0 THz和3.5 THz下,金属铜的介电常数εm分别为-5.4894×105+1.2046i×106和-1.8707×105+1.1729 i×105。计算得出z轴上不同的传输距离处的相位数据与 (10) 式计算结果的差值,结果如表1所示。

表1可知z轴模拟的相位和(10) 式的理论值之间的误差与传播距离呈近似线性关系,虽然限于计算机的运算能力不能模拟50λ0之外的传输情况,但是通过估算可以得到0.35 THz下,传输至500λ0之后z轴处的相位误差也仅约为1.2 rad,足以可见(10) 式的准确性。

表 1. 不同频率下z轴不同传输距离处的相位模拟数据与理论的差值

Table 1. Difference in phase data between simulated and theoretical results at different propagating distances of z-axis at different frequencies

zDifference of phase /rad
f=0.35 THzb=0.85 mmf=1 THzb=0.30 mmf=3.5 THzb=0.085 mm
10λ00.0260.0250.016
20λ00.0520.0510.032
30λ00.0750.0710.045
40λ00.0960.0960.056
50λ00.1210.1200.072

查看所有表

另外,需要说明的是,通过填充多层介质的PPWG可以实现余弦-高斯模式,并且该方法已经得到了验证。优化后的高斯光束与填充多层介质的有限PPWG耦合效率约为89.4%[18]。因此,余弦-高斯光束或将可以作为传统高斯光源与PPWG之间一种新的耦合模式。

最后需要阐明的是,在理论推导部分,假定只有一个电场分量Ex,将不满足散度方程Ñ•E=0。在实际的操作模拟中,假定电场主要为x方向上的线偏振,Ey的场值远小于Ex。在假定的Ey=0的情况下,Ez为电场的纵向分量的标量形式,可以表示为

Ez=Ezez=U(x,y,z)exp(-ikzzez,(15)

式中: Ez为电场的纵向分量的矢量形式;U为假设的标量函数;kz为传播常数;ez为z方向上的单位矢量。

在傍轴近似的条件下,∂U/z可以被忽略,可以得到

Ezz=-ikzEz(16)

将(1)式、(15)式和(16)式代入散度方程,Ez有如下形式:

Ez=1ikzExx(17)

通过加入Ez分量,模拟的收敛速度将会加快。

5 结论

基于平行平板波导的TE1模态,在近似线性偏振电场假设的基础上,提出并推导了一种适用于PPWG的余弦-高斯基本光束的电场表达式。对光束在金属铜PPWG中的传输进行数值仿真,验证了表达式具有较高的精度。与传统的高斯光束相比,这些结果为光束在PPWG内部的传输提供了一个全新的视角。余弦-高斯光束在太赫兹平行平板波导耦合和其他太赫兹波导结构中具有广阔的应用前景。

参考文献

[1] Mittleman D M. Perspective:terahertz science and technology[J]. Journal of Applied Physics, 2017, 122(23): 230901.

[2] Zhang L, Zhang M, Liang H W. Realization of full control of a terahertz wave using flexible metasurfaces[J]. Advanced Optical Materials, 2017, 5(24): 1700486.

[3] Hu B B, Nuss M C. Imaging with terahertz waves[J]. Optics Letters, 1995, 20(16): 1716-1718.

[4] MittlemanD.[\s]{1}Sensing[\s]{1}with[\s]{1}terahertz[\s]{1}radiation[M].[\s]{1}Berlin,[\s]{1}Heidelberg:[\s]{1}Springer,[\s]{1}2003.[\s]{1}

[5] Koenig S, Lopez-Diaz D, Antes J, et al. Wireless sub-THz communication system with high data rate[J]. Nature Photonics, 2013, 7(12): 977-981.

[6] Karl N J. McKinney R W, Monnai Y, et al. Frequency-division multiplexing in the terahertz range using a leaky-wave antenna[J]. Nature Photonics, 2015, 9(11): 717-720.

[7] Rusina A, Durach M, Nelson K A, et al. Nanoconcentration of terahertz radiation in plasmonic waveguides[J]. Optics Express, 2008, 16(23): 18576-18589.

[8] Liang H W, Ruan S C, Zhang M, et al. Nanofocusing of terahertz wave on conical metal wire waveguides[J]. Optics Communications, 2010, 283(2): 262-264.

[9] Gao H, Cao Q, Zhu M N, et al. Nanofocusing of terahertz wave in a tapered hyperbolic metal waveguide[J]. Optics Express, 2014, 22(26): 32071-32081.

[10] Seo M A, Park H R, Koo S M, et al. Terahertz field enhancement by a metallic nano slit operating beyond the skin-depth limit[J]. Nature Photonics, 2009, 3(3): 152-156.

[11] Zhan H, Mendis R, Mittleman D M. Superfocusing terahertz waves below λ/250 using plasmonic parallel-plate waveguides[J]. Optics Express, 2010, 18(9): 9643-9650.

[12] Wang K L, Mittleman D M. Metal wires for terahertz wave guiding[J]. Nature, 2004, 432(7015): 376-379.

[13] Awad M M, Cheville R A. Transmission terahertz waveguide-based imaging below the diffraction limit[J]. Applied Physics Letters, 2005, 86(22): 221107.

[14] Wächter M, Nagel M, Kurz H. Metallic slit waveguide for dispersion-free low-loss terahertz signal transmission[J]. Applied Physics Letters, 2007, 90(6): 061111.

[15] Chen D R, Chen H B. A novel low-loss terahertz waveguide: polymer tube[J]. Optics Express, 2010, 18(4): 3762-3767.

[16] Liang H W, Ruan S C, Zhang M. Terahertz surface wave propagation and focusing on conical metal wires[J]. Optics Express, 2008, 16(22): 18241-18248.

[17] Gallot G, Jamison S P. McGowan R W, et al. Terahertz waveguides[J]. Journal of the Optical Society of America B, 2000, 17(5): 851-863.

[18] Wang K, Cao Q, Zhang H F, et al. Cos-Gaussian modal field of a terahertz rectangular metal waveguide filled with multiple slices of dielectric[J]. Optics Communications, 2018, 417: 57-61.

[19] Mendis R, Mittleman D M. An investigation of the lowest-order transverse-electric (TE1) mode of the parallel-plate waveguide for THz pulse propagation[J]. Journal of the Optical Society of America B, 2009, 26(9): A6-A13.

[20] Kim S H, Lee E S, Ji Y B, et al. Improvement of THz coupling using a tapered parallel-plate waveguide[J]. Optics Express, 2010, 18(2): 1289-1295.

[21] Markov A, Guerboukha H, Skorobogatiy M. Hybrid metal wire-dielectric terahertz waveguides: challenges and opportunities [Invited][J]. Journal of the Optical Society of America B, 2014, 31(11): 2587-2600.

[22] Cao Q, Jahns J. Azimuthally polarized surface plasmons as effective terahertz waveguides[J]. Optics Express, 2005, 13(2): 511-518.

[23] He X Y, Wang Q J, Yu S F. Numerical study of gain-assisted terahertz hybrid plasmonic waveguide[J]. Plasmonics, 2012, 7(3): 571-577.

[24] Li H S, Atakaramians S, Lwin R, et al. Flexible single-mode hollow-core terahertz fiber with metamaterial cladding[J]. Optica, 2016, 3(9): 941-947.

[25] Atakaramians S, Afshar V S, Monro T M, et al. Terahertz dielectric waveguides[J]. Advances in Optics and Photonics, 2013, 5(2): 169-215.

[26] 陈颖, 曹景刚, 许扬眉, 等. 双金属挡板金属-电介质-金属波导耦合方形腔的Fano共振传感特性[J]. 中国激光, 2019, 46(2): 0213001.

    Chen Y, Cao J G, Xu Y M, et al. Fano resonance sensing characteristics of metal-dielectric-metal waveguide coupling square cavity with bimetallic baffle[J]. Chinese Journal of Lasers, 2019, 46(2): 0213001.

[27] Mendis R, Grischkowsky D. Undistorted guided-wave propagation of subpicosecond terahertz pulses[J]. Optics Letters, 2001, 26(11): 846-848.

[28] Harsha S S, Laman N, Grischkowsky D. High-Q terahertz Bragg resonances within a metal parallel plate waveguide[J]. Applied Physics Letters, 2009, 94(9): 091118.

[29] Laman N, Harsha S S, Grischkowsky D, et al. High-resolution waveguide THz spectroscopy of biological molecules[J]. Biophysical Journal, 2008, 94(3): 1010-1020.

[30] Eyyubo lu H T, Baykal Y. Analysis of reciprocity of cos-Gaussian and cosh-Gaussian laser beams in a turbulent atmosphere[J]. Optics Express, 2004, 12(20): 4659-4674.

[31] Zhou G Q, Chu X X. Propagation of a partially coherent cosine-Gaussian beam through an ABCD optical system in turbulent atmosphere[J]. Optics Express, 2009, 17(13): 10529-10534.

[32] 黄永超, 蔡达锋, 张廷蓉. 双曲余弦高斯光束在梯度折射率介质中的束腰宽度及其位置[J]. 激光与光电子学进展, 2015, 52(2): 022601.

    Huang Y C, Cai D F, Zhang T R. Beam waist width and waist location of cosh-Gaussian beams in a gradient-index medium[J]. Laser & Optoelectronics Progress, 2015, 52(2): 022601.

[33] 杨玉婷, 张廷蓉, 龚霞. 余弦高斯光束通过含球差分数傅里叶变换系统的传输[J]. 光学学报, 2017, 37(2): 0207001.

    Yang Y T, Zhang T R, Gong X. Propagation of cosine-Gaussian beams in spherically aberrated fractional Fourier transform system[J]. Acta Optica Sinica, 2017, 37(2): 0207001.

[34] Zhu S J, Chen Y H, Wang J, et al. Generation and propagation of a vector cosine-Gaussian correlated beam with radial polarization[J]. Optics Express, 2015, 23(26): 33099-33115.

[35] Chen R P, Ni Y Z, Chu X X. Propagation of a cos-Gaussian beam in a Kerr medium[J]. Optics & Laser Technology, 2011, 43(3): 483-487.

[36] Keshavarz A. Propagation of cos-Gaussian beam in photorefractive crystal[J]. World Academy of Science, Engineering and Technology, 2015, 9(12): 742-745.

[37] Pozar[\s]{1}DM.[\s]{1}Microwave[\s]{1}engineering[M].[\s]{1}Hoboken:[\s]{1}John[\s]{1}Wiley[\s]{1}&[\s]{1}Sons,[\s]{1}2009.[\s]{1}

[38] Kogelnik H, Li T. Laser beams and resonators[J]. Applied Optics, 1966, 5(10): 1550-1567.

[39] Cao Q, Deng X M. Corrections to the paraxial approximation of an arbitrary free-propagation beam[J]. Journal of the Optical Society of America A, 1998, 15(5): 1144-1148.

[40] Ordal M A, Bell R J, Alexander R W, et al. Optical properties of fourteen metals in the infrared and far infrared: Al, Co, Cu, Au, Fe, Pb, Mo, Ni, Pd, Pt, Ag, Ti, V, and W[J]. Applied Optics, 1985, 24(24): 4493-4499.

吴登明, 滕达, 曹清, 白丽华, 李哲. 余弦-高斯光束在太赫兹平行平板波导中的传输研究[J]. 激光与光电子学进展, 2019, 56(19): 192302. Dengming Wu, Da Teng, Qing Cao, Lihua Bai, Zhe Li. Cos-Gaussian Beams Propagating Inside Terahertz Parallel-Plate Waveguides[J]. Laser & Optoelectronics Progress, 2019, 56(19): 192302.

本文已被 1 篇论文引用
被引统计数据来源于中国光学期刊网
引用该论文: TXT   |   EndNote

相关论文

加载中...

关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。
全站搜索
您最值得信赖的光电行业旗舰网络服务平台!