周振婷 楊理 姚潔 葉燃 徐歡歡 葉永紅

(南京師范大學(xué)物理科學(xué)與技術(shù)學(xué)院，南京 210046)

(2012年11月23日收到;2013年6月3日收到修改稿)

1 引言

隨著微納技術(shù)的發(fā)展，金屬納米材料作為一種人工制備的具有多種電磁學(xué)性質(zhì)的光學(xué)材料，引起了科學(xué)界的廣泛興趣.金屬顆粒在入射光的照射下，其表面的傳導(dǎo)電子會隨入射光場做受迫振動，當金屬顆粒中傳導(dǎo)電子和入射光發(fā)生共振時，金屬表面的局域等離子波被激發(fā)[1]，對入射光發(fā)生很強的吸收或散射，產(chǎn)生一個明顯的共振峰，這就是局域表面等離子體共振(localized surface plasmon resonance，LSPR)現(xiàn)象.顆粒的幾何形狀和周圍的介質(zhì)對于LSPR的峰位影響很明顯[2-4].這一高靈敏性，使得LSPR被廣泛應(yīng)用在納米光學(xué)[5]、生物化學(xué)傳感[6]、超快的光速開關(guān)[7-9]、非線性光學(xué)[10]等領(lǐng)域.

近年來，科學(xué)家研究了不同形貌的Ag，Cu金屬顆粒的光學(xué)特性，如金屬納米球[11]、納米環(huán)[12]、納米棱鏡[13]等.有一部分研究針對的是金屬/介質(zhì)交替層疊的多層納米顆粒，當堆迭層大小一致時，理論和實驗研究均表明，納米顆粒的共振峰位置隨著金屬/介質(zhì)層數(shù)的增加而紅移[14].對于金屬/介質(zhì)層上下不一致的多層納米顆粒，目前已有的研究工作主要是針對其磁響應(yīng)和負折射效應(yīng)等[15-17]，還沒有表面等離激元共振峰位隨層數(shù)增加而發(fā)生變化的研究.

本文結(jié)合納米球刻蝕法(nanosphere lithography，NSL)和熱蒸發(fā)鍍膜的技術(shù)，在玻璃襯底上制備了銀和氧化硅交替層疊的周期排列的錐狀三角納米顆粒陣列.實驗表明，透過率曲線在近紅外波段有一個明顯的透射谷，該透射谷的位置隨著金屬/介質(zhì)層數(shù)的增加發(fā)生藍移.我們采用Ansoft HFSS(high frequency structure simulator)商業(yè)軟件對該結(jié)構(gòu)進行了模擬仿真，利用局域電場場強分布的情況對該實驗現(xiàn)象進行了解釋和討論.

2 實驗

樣品的制備分為三步.第一步，二維膠體晶體模板的制備.將玻璃片切割成2.5 cm×4 cm的基片，對玻璃片做親水性處理后將其浸入重鉻酸鉀溶液中保持12 h，再用去離子水反復(fù)沖洗后烘干備用.將兩塊玻璃片相互接觸組成夾角為2°楔形的三面開口的微腔，選用直徑為960 nm，體積分數(shù)為10%的聚苯乙烯微球懸浮液，用注射器將配置好的溶液注入微腔[18].將微腔的斜面向下放在恒溫恒濕箱中，直至樣品干燥.拆開微腔，得到在玻璃基片上的大面積六角密堆排列的單層膠體晶體.第二步，用熱蒸發(fā)的方法在二維有序的膠體晶體模板表面交替沉積20 nm的銀和20 nm的氧化硅.最后一步，通過超聲的方法去除二維膠體晶體以及膠體晶體表面的金屬/介質(zhì)半球殼陣列，在襯底上形成一種相互分離的三角形金屬/介質(zhì)納米“島”組成的有序陣列.我們可以得到1到3層的非對稱金屬/介質(zhì)多層納米陣列，即樣品A是Ag/SiO層疊的點陣，樣品B是(Ag/SiO)2層疊的點陣，樣品C是(Ag/SiO)3層疊的點陣.

3 結(jié)果與討論

圖1(a)是膠體晶體掩模板的掃描電鏡圖，圖1(b)，(c)和(d)是樣品A、B和C的電鏡圖，三角形的納米顆粒清晰可見.圖1(d)中圈出來的部分可以比較清楚地看出顆粒是具有一定高度的上小下大的三維結(jié)構(gòu).由于金屬/介質(zhì)是沿著聚苯乙烯小球邊緣的空隙沉積在基板上，故而剩余在玻璃基板上的三角形點陣的邊緣呈現(xiàn)沿球表面的弧形向上的形態(tài).

圖1 (a)NSL制備的掩模板的電鏡圖;(b)，(c)，(d)分別是樣品A，B和C的電鏡圖

圖2 樣品A，B和C的透過率 (a)樣品的實驗值;(b)樣品的模擬值

我們利用SHIMADZU UV-3600紫外可見分光光度計測量樣品的透過率.圖2(a)是樣品A，B和C的透過率光譜曲線.從圖2(a)中可以看到，樣品A，B和C的LSPR共振峰的峰位分別位于1320，1276和1004 nm.我們發(fā)現(xiàn)，隨著金屬/介質(zhì)層數(shù)的增多，透射谷的位置發(fā)生明顯的藍移.透過率曲線在850 nm左右有微小震蕩，是在測量過程中光度計換探測器導(dǎo)致的.已有文獻中研究的金屬/介質(zhì)層疊的結(jié)構(gòu)均為上下大小一致的對稱結(jié)構(gòu)，共振峰位隨著金屬/介質(zhì)層的增多而紅移[14].我們制備了上下層大小不同的錐狀結(jié)構(gòu)，其透射谷的位置隨著金屬/介質(zhì)層的增多而發(fā)生藍移.為了研究藍移的原因，我們采用商業(yè)軟件Ansoft HFSS軟件對樣品進行了數(shù)值模擬，并通過局域的電場分布分析了該現(xiàn)象產(chǎn)生的原因.

4 HFSS模擬計算及討論

在采用HFSS對三個樣品的透過率曲線進行數(shù)值模擬的時候，模型的尺寸依據(jù)樣品電鏡圖測量尺寸和微球的大小.如電鏡圖所示，三角形納米顆粒圍繞黑色圓形呈“領(lǐng)結(jié)”[19]狀排列的周期性陣列，圖中黑色圓形部分的直徑大小和納米球的直徑一致，單個三角顆粒的大小由顆粒的中垂線長度a和高度b決定.球徑大小為960 nm，計算得三角陣列橫截面的中垂線長度約為317 nm.從電鏡圖還可以看出顆粒沿微球周邊排列，且隨著金屬/介質(zhì)層的增多顆粒頂端面積不斷減小，所以金屬/介質(zhì)層的幾何參數(shù)采用圖3所示的錐狀顆粒結(jié)構(gòu).模型計算點陣的一個周期，采用垂直于y軸的一對電場邊界和垂直于x軸的一對磁場邊界來擬合連續(xù)周期邊界條件，入射光沿z軸負方向垂直樣品入射.模型中銀的材料參數(shù)采用Drude模型進行擬合，εAg=1-ωp/[ω(ω+iωx)]，其中ωp=1.37×1016Hz是Ag的體等離子振蕩頻率，ωx=8.5×1013Hz是共振頻率.氧化硅的折射率采用1.46，點陣周圍均設(shè)為空氣.

樣品A，B和C透過率的模擬值如圖2(b)所示，每個樣品都有一個明顯的透射谷，樣品A，B和C的共振峰位分別為1322，1185和1050 nm，且隨著金屬/介質(zhì)層的增加，共振峰位發(fā)生藍移.比較圖2(a)和(b)中不同層數(shù)的金屬/介質(zhì)點陣的透過率曲線的實驗值和仿真值，兩者的峰位及其移動方向基本符合.仿真結(jié)果和實驗結(jié)果的差別主要來源于實際顆粒表面不平整和樣品周期性沒有達到仿真模型中的那么完美.從電鏡照片可以看出，在有些地方，納米顆粒存在粘連的現(xiàn)象，大小也不完全一致，這會導(dǎo)致整體的透射率下降而且透射谷變寬.此外，為了驗證模型的正確性，我們還用該軟件模擬了上下層大小一致的金屬/介質(zhì)納米顆粒陣列，得到透射谷隨層數(shù)增多而紅移的結(jié)論，這與已報道的結(jié)果是符合的.

圖3 模擬計算的結(jié)構(gòu)圖

為了分析透射谷隨層數(shù)增加而藍移的原因，對模擬結(jié)果中的電場強度的分布情況進行了研究.圖4是樣品A，B和C在各自透射谷位置處的電場圖，圖4(a)，(b)和(c)分別是樣品A，B和C在x-y平面內(nèi)三角納米顆粒陣列頂端銀層的電場分布圖，圖4(d)，(e)和(f)分別是樣品A，B和C的三角側(cè)邊電場分布圖.當金屬納米顆粒發(fā)生LSPR共振時，金屬表面局域電場會大大增強[20].圖4(a)，(b)，(c)顯示，發(fā)生共振時，三角形的三個尖端電場強度比其他邊界的電場強度強很多，這種現(xiàn)象和以往文獻中三角納米顆粒的結(jié)果是一致的[21].當考慮一層到三層金屬/介質(zhì)納米顆粒時，圖 4(d)，(e)，(f)表明，電場增強被局域在納米顆粒最上面的銀層，在三角錐狀陣列的其他金屬/介質(zhì)交界面的電場都非常弱，甚至沒有電場分布.仿真結(jié)果表明金屬/介質(zhì)層增多時，電場增強的主要是頂層金屬層，隨著層數(shù)的增多，起主要作用的頂層金屬層的面積減小，故而共振峰位藍移.

為了確定下層金屬層對頂層金屬層的影響，通過數(shù)值計算來觀察介質(zhì)層厚度變化對共振峰位的影響.對樣品B的結(jié)構(gòu)(二層金屬/介質(zhì)納米顆粒)，在前面的仿真基礎(chǔ)上，保持金屬層厚度20nm不變，將介質(zhì)層的厚度依次取10，20，30和40 nm，得到介質(zhì)層不同厚度的三角納米陣列的透過率曲線如圖5所示.

從圖5可以看出，當金屬層厚度保持20 nm不變，介質(zhì)層的厚度設(shè)為10，20，30和40 nm時，點陣共振峰位分別為1192，1189，1161和1135 nm.可以得到這樣的規(guī)律：隨著中間介質(zhì)層厚度的增加，三角納米點陣的共振峰位發(fā)生藍移.介質(zhì)層厚度增加后，下層對頂層的作用更弱且頂層的面積減小，LSPR峰位發(fā)生藍移.對比對稱結(jié)構(gòu)，介質(zhì)層變厚時金屬層之間偶合作用減弱，共振峰位發(fā)生紅移.

圖4 (a)，(b)，(c)是樣品A，B和C共振峰位波長處的x-y平面的三角納米顆粒陣列頂端銀層的電場分布圖;(d)，(e)，(f)是樣品A，B和C共振峰位波長處三角側(cè)邊電場分布圖

圖5 介質(zhì)層厚度為10，20，30和40 nm時的透過率曲線

5 結(jié)論

本文利用納米球刻蝕和熱蒸發(fā)的方法制備了不同層數(shù)的銀和氧化硅交替層疊的三角納米顆粒陣列.采用掃描電鏡觀測其表面形貌，可以看出納米顆粒呈以球徑為960 nm周期性有序排列的錐狀三角陣列.樣品的透過率曲線測量表明，該非對稱多層納米陣列的透過率有一個非常明顯的透射谷.隨著金屬/介質(zhì)層的增加，透射谷的位置藍移.利用仿真軟件HFSS理論計算了不同層數(shù)的樣品和介質(zhì)層厚度變化的(Ag/SiO)2納米陣列的透過率曲線，并結(jié)合電場強度的分布，解釋了藍移現(xiàn)象產(chǎn)生的原因.

[1]Willets K A，van Duyne R P 2007Annu.Rev.Phys.Chem.58 267

[2]Jensen T R，Malinsky M D，Haynes C L，van Duyue R P，Haglund Jr R F，Yang L，Magruder III R H，Wittig J E 2000J.Phys.Chem.B 104 10549

[3]Wang K，Yang G，Long H，Li Y H，Dai N L，Lu P X 2008Acta Phys.Sin.57 386(in Chinese)[王凱，楊光，龍華，李玉華，戴能利，陸培祥2008物理學(xué)報57 386]

[4]Fukumi K，Chayahara A，Kadono K，Sakaguchi T，Horino Y，Miya M，F(xiàn)ujii K，Hayakawa J，Satou M 1994J.Appl.Phys.75 3075

[5]Huang W，Qian W，E-l Sayed M A 2005J.Phys.Chem.B 109 18881

[6]Haes A J，Hall W P，Chang L，Klein W L，van Duyne R P 2004Nano Lett.4 1029

[7]Feldstein M J，Keating C D，Liau Y H，Natan M J，Scherer N F 1997J.Am.Chem.Soc.119 6638

[8]Jensen Fukumi K，Chayahara A，Kadono K，Sakaguchi T，Horino Y，Miya M，F(xiàn)ujii K，Hayakawa J，Satou M 1994J.Appl.Phys.75 3075

[9]Haglund Jr R F，Yang L，Magruder III R H，Wittig J E，Becker K，Zuhr R A 1993Opt.Lett.18 373

[10]Wang W，Wang Y，Sun Y 2007Appl.Sur.Sci.253 4673

[11]Chien W Y，Szkopek T 2008Opt.Express16 1820

[12]Aizpurua J，Hanarp P，Sutherland D S，K¨all M，Garnett W B，Garc′?a de Abajo F J 2003Phys.Rev.Lett.90 057401

[13]Jin R C，Cao Y，Mirkin C A，Kelly K L，Schatz G C，Zheng J G 2001Science294 1901

[14]Su K H，Wei Q H，Zhang X 2006Appl.Phys.Lett.88 063118

[15]Jeyaram Y，Jha S K，Agio M，L¨of fler J F，Ekinci Y 2010Opt.Lett.35 10

[16]Ma W Y，Yao J，Yang H，Liu J Y 2010Nanotechnology and PrecisionEngineering8 240

[17]Yang X Y，Yao J，Junsuk Rho，Yin X B，Zhang X 2012Nature Photon.6 450

[18]Sun J，Tang C J，Zhan P，Han Z L，Cao Z S，Wang Z L 2010Langmuir26 7859

[19]Iwami K，Ono T，Esashi M 2006J.Microelectromech.Syst.15 1201-8

[20]Huang Q，Zhang X D，Wang S，Cao L R，Sun J D，Geng W D，Xiong S Z，Zhao Y 2009Acta Phys.Sin.58 2731(in Chinese)[黃茜，張曉丹，王爍，曹麗冉，孫建東，耿衛(wèi)東，熊紹珍，趙穎2009物理學(xué)報58 2731]

[21]Ma W Y，Yang H，Hilton J P，Lin Q，Liu J Y，Huang L X，Yao J 2010Opt.Express18 843