二次元多周期1维光子晶体带隙的实质

范希智

(武汉理工大学 物理系，武汉 430070)

引 言

1987年,YABLONOVITCH在研究抑制自发辐射时提出了光子晶体的概念[1]。几乎同时，JOHN在讨论光子局域时也独立地提出了这个概念[2]。光子晶体是将不同介电常数(或折射率)的介质材料在空间按一定的周期排列而形成的一种人为的“晶体”结构。

光子晶体的提出展示了一种新的控制光子的机制，促进了光子技术的进一步研究，给通讯技术、光电子技术等领域的发展和应用带来了新的生机和活力，使得用光子取代电子来传输、处理和存储信息的设想更有希望实现[3]。

光子晶体中的折射率(或介电常数)是空间的周期函数，与半导体材料中的电子在周期性势场作用下形成能带结构相类似，在光子晶体中传播的光子能量也会有能带结构[9]，能带与能带之间会出现光子带隙(或称禁带)，因此，又可将光子晶体称为光子带隙材料或光子半导体。频率(波长)落在带隙中的光子不能在光子晶体中传播，这些光不能进入光子半导体即光子晶体，将被全部反射。

表面看来，1维光子晶体的原始构成材料的折射率只在一个方向上呈周期性变化，因而结构简单，实际上它是一种周期性结构的分层薄膜材料，可采用镀膜方法制备[10]。早在光子晶体概念提出之前，人们已经开始对这种1维周期结构薄膜进行研究，主要用于各种光学器件的增透、反射等的镀膜领域，因此1维光子晶体是一种研究较为成熟的光子晶体，可用于可见光及红外波段。

表征材料光学性能的参量是折射率，而光子带隙材料即光子晶体的折射率有什么特性？从光学角度看，光子晶体的带隙是怎样的情况？本文中将利用分析光学薄膜的干涉矩阵法及相关的概念来讨论二次元1维光子晶体的带隙及其对应的折射率之特征。

1 理论分析

构成1维光子晶体的原始材料的折射率只在1维方向上周期性地变化，一个周期可以由几层不同折射率的原始材料组成。二次元1维光子晶体的周期是最简单的，仅由A和B两种原始材料组成，它们的折射率分别为nA和nB，整个光子晶体是由很多个同样的周期重复排列而成，如ABAB…ABAB，记为(AB)L，L为周期数。由于原料的折射率变化是严格周期性的，因此这种光子晶体可称之为1维本征光子晶体或1维无掺杂光子晶体。目前，这种光子晶体一般附着于玻璃等透明衬底上而不能独立存在，其结构如图1所示。

Fig.1 Model of 1-D photonic crystals structure

分析光子晶体的光学特性的方法有许多[3,11-18]，而对于1维光子晶体来说，采用传输特征矩阵法是非常合适的[12,15,19]，假定光波从折射率为n0的入射介质入射到由均匀、各向同性、非磁性的介质组成的1维光子晶体(多层膜系)的表面上，光子晶体制作在折射率为ng的透明衬底上(如图1中所示的玻璃衬底)，则由传输特征矩阵法得到光子晶体表面的能量反射率为:

(1)

式中,η0为入射介质的导纳,Y是光子晶体与衬底的组合体的等效导纳。

与图2相比，产物在2 997，1595和1 513 cm-1处的3个苯环特征吸收峰依然存在，1 595和1 419 cm-1处为—NO2的伸缩振动特征吸收峰，2 234 cm-1处腈基—C N的伸缩振动特征吸收峰消失。图3中，在3198 cm-1处出现了N—H的伸缩振动特征吸收峰，在1595 cm-1处出现了—C═N的伸缩振动特征吸收峰；由于官能团为共振式，在779，705 cm-1处的2个峰为—NH2的伸缩振动特征吸收峰。上述谱图的特征吸收峰证明产物为对硝基苯甲脒。

(2)

式中,B和C构成光子晶体与衬底的组合传输矩阵的矩阵元素，表示为：

(3)

式中,ηj和ηg分别是第j层膜和衬底介质的导纳。在垂直入射的情况下，导纳与折射率相等：ηi=c/vi(i=0,j,g),vi为入射介质中、第j层膜中和衬底介质中的光速，c为真空中光速。当ηg=1时，Y是纯光子晶体的导纳;在可见光、近红外波段，Y为光子晶体的折射率。

(4)

式中，Tj为第j层薄膜的特征传输矩阵,δj为第j层膜的位相厚度，定义为:

(5)

式中，λ为入射光波的波长,dj为第j层膜的厚度,θj为夹角,nj为第j层膜的折射率,nj=c/vj。

对于图1所示的结构，每个周期内有两层膜，则周期的传输矩阵变为:

(6)

矩阵元分别为：

(7)

很明显，关系式m11m22+m12m21=1成立。则对于结构如图1所示的1维光子晶体，(3)式变为：

(8)

式中,L是指数，表示周期，即存在2L层膜。令E表示单位矩阵，即:

(9)

由(8)式递推得：ML=UL-1(χ)M-UL-2(χ)E，从而有:

ML=

(10)

式中,UL-1(χ)、UL-2(χ)是第L-1阶次、第L-2阶次第2类切比雪夫多项式。于是(3)式变为：

(11)

代入(2)式得:

(12)

由上式可见，等效导纳Y是复数。而Y的实部为:

Yreal=

(13)

利用第2类切比雪夫多项式通项表示式和各阶次第2类切比雪夫多项式之间的递推关系[20]，(13)式变为：

Yreal=

(14)

由上面的讨论可知，等效导纳Y及其实部Yreal都是波长λ的函数，即Y=Y(λ)，Yreal=Yreal(λ)。结合第2类切比雪夫多项式的特性，通过(14)式或(13)式可见，1维光子晶体的周期数达到一定数值时，等效导纳的实部Yreal(λ)在某波长范围内的取值将趋近于或等于零。

2 实验讨论

设构成1维光子晶体周期的两种介质A和B分别是氟化镁(MgF2)和硫化锌(ZnS)，它们的折射率分别为nA=1.38(MgF2)，nB=2.35(ZnS)。每个介质层的光学厚度均为某一参考波长λ0的1/4，即nAdA=nBdB=λ0/4，dA和dB是介质层的几何厚度，与波长λ0对应的频率是ω0。晶体由有限个周期组成，设周期数为10，即介质层的总数为20层。将1维光子晶体制作在折射率ng=1.0的衬底上，并设光从真空垂直入射到光子晶体表面，则入射介质的折射率n0=1.0 (ng=1.0即为真空，如前所述，光学薄膜型光子晶体不能独立存在，一般附着于玻璃等衬底上，这样会形成光子晶体与衬底的组合体，这里设ng=1.0是为了讨论纯光子晶体的情况)。图2a中为这种1维光子晶体的Yreal(λ)曲线。该曲线是根据(14)式或(13)式利用数值计算而绘制的，此时取λ0=510nm，由图可见，大约在440nm～630nm的范围内，等效导纳实部Yreal(λ)是等于零(趋近于零)的。

当实部等于零时，等效导纳Y即变成纯虚数，在此情况下，(1)式中的Y就是纯虚数，这样反射率R=1.0，在相应波长范围内的入射光被全部反射，按照光子晶体带隙的定义，这个波长范围就是带隙。图2b中的光子晶体的反射率曲线是根据(1)式～(5)式进行数值计算而绘制的。由图可见，反射率R=1.0,波长范围也是大约在440nm～630nm之间。

这里讨论的是光学薄膜型的1维光子晶体，如前所述，其等效导纳即为其等效折射率，在光垂直入射时即为光子晶体的折射率[10]。在带隙范围内1维光子晶体的(等效)折射率是虚(等效)折射率。正是由于虚(等效)折射率的出现，处于带隙内的入射光被全部反射，实现光子晶体的功用。

图2是1维10周期的MgF2/ZnS光子晶体在光垂直入射下出现带隙与虚折射率(复折射率的实部为零)相对应的情形。同样，根据(14)式或(13)式利用数值计算而绘制这种光子晶体光在斜入射下的Yreal(λ)曲线，取入射角为θ0=30°，如图3a和图3b所示。在相同的入射角下，根据(1)式～(5)式进行数值计算而绘制光子晶体斜入射时的反射率R(λ)曲线，如图3c和图3d所示。

Fig.2 The curve of real part of the effective admittance and the reflectivity of 1-D photonic crystal with MgF2/ZnS

光在斜入射情形下，可分成p光和s光分别进行讨论，图3a和图3b中分别是p光和s光的Yreal(λ)曲线：Yp,real(λ),Ys,real(λ)；图3c和图3d中分别是p光和s光的R(λ)曲线：Rp(λ),Rs(λ)。由图中R(λ)曲线可见，无论s光还是p光，都出现反射率是100%的波长区域，即出现带隙，p光和s光的带隙范围大约在420nm～580nm和410nm～600nm，而在对应的波长区域内，Yp,real(λ)=0，Ys,real(λ)=0，如图3a和图3b中的Yreal(λ)曲线所示，这说明在带隙内，光子晶体表现出具有虚等效折射率的情况，这表明在斜入射下，也会出现带隙与虚折射率(复折射率的实部为零)相对应的情形。

Fig.3 The curve of real part of the effective admittance and the reflectivity of 1-D photonic crystal with MgF2/ZnS at incident angle of θ0=30°

由此可见，光子晶体出现带隙，其复折射率的实部为零，即折射率成为虚折射率。至于光子晶体带隙大小，它与构成光子晶体的原始材料的折射率及光子晶体之虚等效折射率材料的关系，要根据(1)式～(10)式和第2类切比雪夫多项式的定义出发进行详细的讨论，限于篇幅，作者将于另文发表分析和结果，这里从略。

3 结 论

1维光子晶体在能带结构上存在着带隙，频率或波长处于带隙的光波不能通过或进入1维光子晶体而被全部反射。1维光子晶体的等效折射率(等效导纳)是复数，随入射光波波长而变化，对于带隙范围内的光波而言，1维光子晶体的等效折射率的实部等于或趋近于零，此时变成虚等效折射率，在垂直入射下即是虚折射率。

作者以氟化镁(MgF2)和硫化锌(ZnS)构成的多周期的二次元1维光子晶体为例，分析了这种情况，可以知道,对于频率或波长处于带隙的光波而言，1维光子晶体表现为其等效折射率是虚等效折射率，换言之，1维光子晶体是虚等效折射率材料。金属可以近似看作虚折射率材料，全反射发生时，光疏媒质看作是虚等效折射率材料，而这里讨论的1维光子晶体可以既是虚折射率材料，又是虚等效折射率材料，前者是光垂直入射的情形，后者则是斜入射下的情形。无论哪种，都是处于1维光子晶体的带隙范围内，由此看来，光子晶体带隙的光学本质是：在带隙范围内光子晶体变成虚等效折射率材料。