金属光栅透射光谱的时域有限差分法模拟.docxVIP

金属光栅透射光谱的时域有限差分法模拟 1. 维金属光栅的透射特性 1998年后，研究发现，当一组光垂直排列在金属膜上时，一些特定波长的传输率有很大提高。根据经典的粒度理论，当一组光通过直径为d（d）的小波时，传输率应约为4%。也就是说，当波长约为900nm的波长通过直径为150nm的小波时，传输率应约为0.001，且试验结果约为4%（源测量值）。然后，schreter等人使用散射矩阵法提供了光透光的一维金属光束的数值解。这表明，即使是电子氧化物也具有这种增强效果。我们称它为一种齿裂现象。由于光场局域、微腔测量电动力学、高密度数据存储、近场光学等领域的巨大应用潜力，导致了对该杂志性能特性的研究。随着这一改进的物理起源的研究，我们可以提出这样一种新的理论。由此提出的光谱和北海贝森等人的二维网格矩阵的结构不同，光在狭窄缝中的模型也不同。然而，由于透射光的特性，如伯森等人的实验结果一样，三维网格的成像可以扩展到分析伯森等人实验模型的基础，这已经成为许多人选择的研究对象。目前，形成了两种代表性观点。首先，金属膜表面离子激元的共振激励和耦合产生的这种增强效果被称为spw（表面plagoonw）理论。然而，该理论很难根据金属膜厚度的增加来解释反射增强的波长。另一种观点认为，f-p腔的作用是导致这种现象的原因。这被称为f-p腔理论。然而，该理论忽略了f-p腔的随机性对透射增强波长的影响。 本文用时域有限差分(FDTD)法模拟平面光(400—2000nm)正入射通过一维金属光栅的透射特性. 通过详细讨论金属厚度等对透射光场的影响及光栅狭缝内的光场分布特性,以期对Ebbesen现象给出物理解释. 2. nix,iy,nt的干计算n 计算所用模型如图1所示.图中金属膜附在基底上,光栅上方为空气.如果不作特殊说明,基底的相对介电常量为2.银作为一种典型的贵重金属,普遍作为研究金属及其结构光学特性的代表,Ebbesen等人在实验中同样也首选金属银膜,主要原因在于它可使增强效应更为显著.因此,本文也以金属银作为研究对象,其介电常量根据已知的分立介电常量拟合而得.光栅的周期为a=900nm,光栅狭缝(即空气缝)宽度d=60nm,狭缝处金属厚度为零,非狭缝处银膜厚度用h表示.假设光栅平面无限大,入射波为TM(transverse-magnetic)波,磁场方向沿狭缝方向. 计算中虚框上下边界采用二次吸收边界条件,左右边界采用周期性边界条件.将计算空间沿x,y方向分成135×45个网格单元,Δx,Δy分别为网格单元沿x,y方向上的长度,Δx=Δy=20nm,时间增量Δt=Δx/(2c),c为真空中的光速,总计算时间为16个入射光周期,(i,j)为格点坐标.任意一个时间和空间的函数表示为 Fn(i,j)=F(iΔx,jΔy,nΔt).(1)Fn(i,j)=F(iΔx,jΔy,nΔt).(1) 根据Maxwell方程组,在无磁损耗介质中,空间各点的电场、磁场的各分量可写成 En+1x(i+1/2,j)=C(i+1/2,j){-D(i+1/2,j)×Enx(i+1/2,j)+[Ηn+1/2z(i+1/2,j+1/2)-Ηn+1/2z(i+1/2,j-1/2)]/Δy},(2)En+1y(i,j+1/2)=C(i,j+1/2){-D(i,j+1/2)×Eny(i,j+1/2)-[Ηn+1/2z(i+1/2,j+1/2)-Ηn+1/2z(i-1/2,j+1/2)]/Δx},(3)Ηn+1/2z(i+1/2,j+1/2)=C′(i+1/2,j+1/2)×{-D′(i+1/2,j+1/2)Ηn-1/2z(i+1/2,j+1/2)+[Enx(i+1/2,j+1)-Enx(i+1/2,j)]/Δy-[Eny(i+1,j+1/2)-Eny(i,j+1/2)]/Δx},(4)En+1x(i+1/2,j)=C(i+1/2,j){?D(i+1/2,j)×Enx(i+1/2,j)+[Hn+1/2z(i+1/2,j+1/2)?Hn+1/2z(i+1/2,j?1/2)]/Δy},(2)En+1y(i,j+1/2)=C(i,j+1/2){?D(i,j+1/2)×Eny(i,j+1/2)?[Hn+1/2z(i+1/2,j+1/2)?Hn+1/2z(i?1/2,j+1/2)]/Δx},(3)Hn+1/2z(i+1/2,j+1/2)=C′(i+1/2,j+1/2)×{?D′(i+1/2,j+1/2)Hn?1/2z(i+1/2,j+1/2)+[Enx(i+1/2,j+1)?Enx(i+1/2,j)]/Δy?[Eny(i+1,j+1/2)?Eny(i,j+1/2)]/Δx},(4) 式中 C(i,j)={√μ0ε0[σ(i,j)2+ε(i,j)Δt]}-1,(5)D(i,j)=√μ0ε0[σ(i,j)