腔體模式下表面等離子激元共振的仿真研究*

曹 軍，張宇寧

(東南大學電子科學與工程學院，南京210096)

腔體模式下表面等離子激元共振的仿真研究*

曹 軍，張宇寧*

(東南大學電子科學與工程學院，南京210096)

為了分析腔體模式下表面等離子激元共振情況，介紹了電介質球腔體本征模式的多重散射解析方法和有限元數值計算方法(comsol軟件)，運用有限元法建立簡單模型進行數值仿真，簡要分析腔體模式與金屬表面等離子激元的耦合，重點分析了小球半徑與表面等離子波長的關系，提出了改變小球半徑來影響表面等離子波長的構想，并提出了近似的公式。

表面等離子激元;腔體模式;能量耦合;有限元方法

1902年，Wood R W在光學實驗中首次發現了表面等離子激元共振現象［1］。此后幾十年許多科學家對表面等離子激元進行了理論分析與實驗探究。1960年，Stren E A和Farrel R A研究了金屬與介質界面處產生共振的條件，并首次提出了表面等離子激元(Surface Plasmon)的概念。

表面等離子激元是光和金屬表面自由電子的相互作用引起的一種電磁波模式，它局限于金屬與介質表面附近，能夠增強近場。在一般情況下，表面等離子體的波矢量大于光波，所以光波不能直接激發表面等離子。為了成功激勵表面等離子體波，我們需要引入一些特殊結構達到波矢匹配。

1 概述

納米技術成熟之后，表面等離子波受到人們極大的關注，成為目前的研究熱點。惰性金屬薄膜，包含可以共振的自由電子，在滿足一定條件的光激發下，金屬與介質的界面附近會產生強烈的表面等離子波。這種加強和限定金屬表面或納米金屬粒子附近光場的能力已被應用于生物化學傳感，光電子集成器件等多個領域，例如免疫傳感器［2］，熒光傳感器［3］，太陽能電池［4］，等離子激光器［5-6］。

本文討論的是電介質腔體模式與表面等離子體耦合的機制，二者之間的相互作用影響了電介質-金屬表面的光學性能。

當一個電介質腔體被嵌入無限大金屬中時，它可存在一系列本征模式。當外界入射電磁波的頻率與這些本征模式的頻率非常接近時，電介質腔體中很容易被激發出強烈的電磁波。

本文接下來將簡要介紹基于多重散射理論［7-9］的電介質腔體本征模式解析方法，它特別適用于球對稱結構散射中心。我們利用有限元的方法(comsol軟件)建立模型觀察電介質球腔體對表面等離子體的激發情況。有限元是一種數值模擬方法，它采用簡單的問題來近似復雜問題，在有限元內取近似解逼近精確解，實現定性和定量的分析。

2 金屬中電介質腔體本征模式

假設一個半徑為S的電介質小球位于坐標原點，介質球內的材料有一個與頻率有關的復介電常數ε1，周圍環境材料的介電常數為ε2。當如下平面電磁波:

3 建模與仿真

3.1 建立模型

我們建立簡單模型:在金屬層中嵌入電介質。電介質小球的嵌入深度與波長同數量級，這個值小于金屬對于可見光的趨膚深度。金屬表層的上方是空氣，而金屬板的厚度取足夠深，考慮到當前的納米工藝水平，我們設金屬厚度600 nm的立方體。這個厚度遠大于可見光區域金屬的趨膚深度，以保證光不能從下方透射出去。該模型中不存在透射電磁波，所有電磁波通過吸收反映出來。由于硬件條件的限制，我們建立了二維模型仿真，二維的圓代表三維的球。

玻璃作為一種透光率比較好的常用材料一直都活躍在各種科研項目中，我們選取介電常數為2的玻璃小球。金和銀是表面等離子研究方面最常用的貴金屬，它們具有較強的化學惰性，穩定性比較好。而多項研究表明，在同樣的條件下，銀膜可以幫助我們獲得更好的共振特性。我們分別采用銀膜和金膜作為環境介質作了仿真，同樣可以看出銀膜具有一定的優越性。

3.2 應用與分析

多重散射理論對上述模型的計算分析表明，當一個電介質腔體被嵌入無限大的金屬中，腔體中可存在一系列本征模式。上述模型的前3個腔體本征頻率所對應的波長為424 nm、332 nm、288 nm。

我們通過建模仿真，改變頻率觀察腔體中的電場變化，發現430 nm處，腔體電場分布達到一個穩定而較大的狀態(如圖1所示)。可見腔體模式的第1本征頻率所對應的波長確在424 nm附近，且第1本征模式具有最強的響應。隨著模式數的增加，本征模式的響應逐漸減弱。

圖1 第1腔體本征模式

本文中，我們著重關注第1腔體本征模式與表面等離子波的耦合情況。由圖2可以看出，腔體中的電磁場可以在銀-介質表面激發表面等離子波。而隨著電介質小球的深度加大，耦合的情況會不同。

圖2 不同深度的腔體激發表面等離子波的情況

雖然電介質小球嵌入深度越深，越接近于理想腔體模式，但在實際中，我們應該充分考慮入射光線在金屬層中的衰減，即趨膚深度的影響，深度過大，入射光線難以進入腔體中建立腔體模式，耦合效果也會減弱。

不僅電介質小球的嵌入深度可以影響表面等離子波與腔體模式的耦合，而且電介質小球的尺寸大小也會影響激發的表面等離子波的頻率。我們通過改變電介質小球的半徑，觀察到不同波長的表面等離子波受到激發，如圖3所示。

圖3 激發的表面等離子波長與腔體半徑關系(銀膜-玻璃)

對上圖的曲線進行線性擬合，可以得到:

線性回歸模型:λ=9.543R-50。

因此，我們可以用電介質小球的尺寸來控制所激發的表面等離子波的波長。650 nm紅、550 nm綠、450 nm藍光所對應的小球半徑分別為122 nm、109 nm、100 nm。

從圖3可以看出，隨著腔體尺寸的增大，共振激勵的表面等離子波長向長波方向移動，這就是所說的紅移現象［13］。影響紅移的因素有很多，除了量子尺寸效應，主要還有表面效應，即當粒子處于納米尺寸時，粒子表面原子數與總原子數之比隨著粒徑的增加而急劇減小，使得電子在粒子表面和體相之間重新分配，這削弱了鍵強度，引起共振峰紅移。

我們將銀膜換成金膜，再次進行仿真，可以得到圖4所示關系。

圖4 激發的表面等離子波長與腔體半徑關系(金膜-玻璃)

用二次多項式擬合上述曲線更加符合仿真情況，與圖6對比發現，銀膜作為環境材料可以更容易實現對表面等離子波長的線性控制，金膜在此方面略顯不足。

我們將電介質小球的材料換成砷化鎵，環境介質保持為銀膜不變，仿真結果如圖5所示。

圖5 激發的表面等離子波長與腔體半徑關系(銀膜-GaAs)

圖5所示的曲線總的趨勢上符合二次多項式擬合，但是對于450 nm以上的光，控制結果接近于線性。

不同屬性的材料對激發波長的影響是不同的，如何探索最優化的匹配情況，這是值得我們研究的。

4 討論與總結

本文簡要介紹了解析計算金屬中電介質小球電磁模式的多重散射方法，利用有限元方法建立簡單的銀層-電介質小球模型，驗證了腔體模式激發表面等離子的可能性，并分析了嵌入深度和小球尺寸對共振激勵的影響。

電介質小球的嵌入深度越深，越符合理想的腔體模式，但深度越深，到達腔體的光線越弱，我們需要尋找一個合適的深度。小球的尺寸也對激發的表面等離子波長有控制作用，這種控制作用受到材料屬性的強烈影響。在本文所涉及的材料中，銀膜-玻璃小球組合可以實現近似線性控制，這為腔體控制表面等離子體頻率的理念提供了一些參考［14］。此外，電介質球的其他屬性，材料對光的吸收與色散也應該被考慮進去，這些都是要進一步研究的問題。

［1］ Polman A，Atwater H A.Plasmonics:Optics at the Nanoscale［J］. Materials Today，2005，8(1):56-56.

［2］ Liedberg B，Nylander C，Lunstr?m I.Sens［J］.Actuators，1983 (4):299-304.

［3］ Kano H，Kawata S.Opt［J］.Lett，1996，21:1848-1850.

［4］ Ferry V E，Sweatlock L A，Pacifici D，et al.Atwater［J］.Nano Lett，2008(8):4391-4394.

［5］ Okamoto T，H’Dhili F，Kawata S，et al.Appl Phys Lett，2004，85: 3968-3970.

［6］ Zheludev N I，Prosvirnin S L，Papasimakis N，et al.Lasing Spaser［J］.Nature Photonics，2008，2:351-354.

［7］ Stefanou N，Karathanos V，Modinos A.Scattering of electromagnetic waves by periodic structures［J］J Phys:Condens Matter，1992(4): 7389-7400.

［8］ Stefanou N，Yannopapas V，Modinos A.Heterostructures of photonic crystals:frequency bands and transmission coefficients［J］.Computer Physics Communications，1998，113(1):49-77.

［9］ Stefanou N，Yannopapas V，Modinos A，et al.MULTEM2:A new version of a program for transmission and band-structure calculations of photonic crystals［J］.Comput Phys Commun，2000，132: 189-196.

［10］鄭俊娟.金屬-電介質體系中的強吸收及強透射機理研究［D］.北京:中國科學院物理研究所，2005.

［11］鄭俊娟，孫剛.嵌入電介質小球的金屬薄片的電磁波透射特性［J］.物理學報，2010，59(6):4008-4013.

［12］顧小衛，蒙林，孫宜琴，微波爐爐腔內電磁場仿真［J］.電子器件，2008，31(3):803-806.

［13］張妹景.金屬納米顆粒的形狀對表面等離子體共振特性的調制研究［D］.電子科技大學，2010.

［14］江秀明，陳志春，楊紹明.光纖表面等離子體共振傳感器研究進展［J］.傳感技術學報，2003，16(1):74-81.

曹 軍 (1991- )男，漢族，江蘇南通，東南大學電子科學與工程學院，本科，主要研究方向為模擬電路設計，213092298@seu.edu.cn;

張宇寧(1981- )，男，漢族，江蘇南通，博士，副教授，現任教于東南大學。主要從事光電科學與技術方面的研究，包括顯示、照明、光伏等領域，zyn@seu. edu.cn。

Surface Plasmon Polaritons Resonance with Dielectric Spherical Cavity Mode*

CAO Jun，ZHANG Yuning*

(School of Electronic Science and Engineering，Southeast University，Nanjing 210096，China)

Multiple scattering method is introduced to calculate dielectric spherical cavity mode.We use finite element method(comsol multiphysics)to create a simple model to analyze the coupling of cavity mode and surface Plasmon polaritons(SPPs).Simulation results show that cavity size can control the wavelength of SPPs.We focus on the relationship between radius of the ball and wavelength of surface Plasmon polaritions and propose the concept that radius of the ball can be used to affect wavelength of SPPs.Then we propose approximation formula.

surface plasmon polaritons;cavity mode;energy coupling efficiency;finite element method

10.3969/j.issn.1005－9490.2014.02.003

O539

A

1005－9490(2014)02－0182－04

項目來源:江蘇省自然科學基金項目(BK20130629);江蘇省科技支撐工業項目(BE2013087)

2013-06-04修改日期:2013-07-01

EEACC:2810