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

曹 軍，張宇寧

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

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

曹 軍，張宇寧*

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

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

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

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

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

1 概述

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

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

當(dāng)一個(gè)電介質(zhì)腔體被嵌入無限大金屬中時(shí)，它可存在一系列本征模式。當(dāng)外界入射電磁波的頻率與這些本征模式的頻率非常接近時(shí)，電介質(zhì)腔體中很容易被激發(fā)出強(qiáng)烈的電磁波。

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

2 金屬中電介質(zhì)腔體本征模式

假設(shè)一個(gè)半徑為S的電介質(zhì)小球位于坐標(biāo)原點(diǎn)，介質(zhì)球內(nèi)的材料有一個(gè)與頻率有關(guān)的復(fù)介電常數(shù)ε1，周圍環(huán)境材料的介電常數(shù)為ε2。當(dāng)如下平面電磁波:

3 建模與仿真

3.1 建立模型

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

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

3.2 應(yīng)用與分析

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

我們通過建模仿真，改變頻率觀察腔體中的電場(chǎng)變化，發(fā)現(xiàn)430 nm處，腔體電場(chǎng)分布達(dá)到一個(gè)穩(wěn)定而較大的狀態(tài)(如圖1所示)。可見腔體模式的第1本征頻率所對(duì)應(yīng)的波長(zhǎng)確在424 nm附近，且第1本征模式具有最強(qiáng)的響應(yīng)。隨著模式數(shù)的增加，本征模式的響應(yīng)逐漸減弱。

圖1 第1腔體本征模式

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

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

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

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

圖3 激發(fā)的表面等離子波長(zhǎng)與腔體半徑關(guān)系(銀膜-玻璃)

對(duì)上圖的曲線進(jìn)行線性擬合，可以得到:

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

因此，我們可以用電介質(zhì)小球的尺寸來控制所激發(fā)的表面等離子波的波長(zhǎng)。650 nm紅、550 nm綠、450 nm藍(lán)光所對(duì)應(yīng)的小球半徑分別為122 nm、109 nm、100 nm。

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

我們將銀膜換成金膜，再次進(jìn)行仿真，可以得到圖4所示關(guān)系。

圖4 激發(fā)的表面等離子波長(zhǎng)與腔體半徑關(guān)系(金膜-玻璃)

用二次多項(xiàng)式擬合上述曲線更加符合仿真情況，與圖6對(duì)比發(fā)現(xiàn)，銀膜作為環(huán)境材料可以更容易實(shí)現(xiàn)對(duì)表面等離子波長(zhǎng)的線性控制，金膜在此方面略顯不足。

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

圖5 激發(fā)的表面等離子波長(zhǎng)與腔體半徑關(guān)系(銀膜-GaAs)

圖5所示的曲線總的趨勢(shì)上符合二次多項(xiàng)式擬合，但是對(duì)于450 nm以上的光，控制結(jié)果接近于線性。

不同屬性的材料對(duì)激發(fā)波長(zhǎng)的影響是不同的，如何探索最優(yōu)化的匹配情況，這是值得我們研究的。

4 討論與總結(jié)

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

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

［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］鄭俊娟.金屬-電介質(zhì)體系中的強(qiáng)吸收及強(qiáng)透射機(jī)理研究［D］.北京:中國(guó)科學(xué)院物理研究所，2005.

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

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

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

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

曹 軍 (1991- )男，漢族，江蘇南通，東南大學(xué)電子科學(xué)與工程學(xué)院，本科，主要研究方向?yàn)槟M電路設(shè)計(jì)，213092298@seu.edu.cn;

張宇寧(1981- )，男，漢族，江蘇南通，博士，副教授，現(xiàn)任教于東南大學(xué)。主要從事光電科學(xué)與技術(shù)方面的研究，包括顯示、照明、光伏等領(lǐng)域，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

項(xiàng)目來源:江蘇省自然科學(xué)基金項(xiàng)目(BK20130629);江蘇省科技支撐工業(yè)項(xiàng)目(BE2013087)

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

EEACC:2810