梅 菲，楊治杰，許 磊，徐進霞，陳志強

(1.湖北工業(yè)大學，湖北 武漢 430212;2.華北水利水電學院，河南 鄭州 450011;3.武漢大學，湖北 武漢 430072)

基于FDTD的三角形銀納米顆粒陣列等離子共振吸收特性的模擬

梅 菲1，楊治杰1，許 磊2，徐進霞3，陳志強1

(1.湖北工業(yè)大學，湖北 武漢 430212;2.華北水利水電學院，河南 鄭州 450011;3.武漢大學，湖北 武漢 430072)

利用時域有限差分法，數(shù)值模擬了銀納米三角陣列的等離子共振吸收特性，仿真計算了銀納米顆粒尺寸大小對陣列等離子體共振的影響.數(shù)值計算結果表明，三角形銀納米顆粒陣列的等離子共振吸收峰隨顆粒尺寸的增大而紅移，可通過改變顆粒大小調節(jié)共振吸收波長數(shù)值以用于不同的應用研究.基于FDTD理論模擬的計算結果與實驗規(guī)律基本相符，并較好地闡釋了銀納米三角有序陣列表面等離子共振吸收峰的頻移與其納米粒子長徑比之間的一些依賴關系，此結論將有利于納米光傳感器應用的進一步探索研究.

銀納米顆粒;時域有限差分法;等離子共振吸收

表面等離子體學的發(fā)展與納米材料制備技術的結合使光和電的協(xié)同作用進入了一個嶄新的階段.表面等離子體共振(Surface Plasmon Resonance，SPR)是金屬納米結構非常獨特的光學特性，基于表面等離子體共振的納米結構體系的研究已成為國際上迅猛發(fā)展的熱點研究領域之一.金屬納米結構中被激發(fā)的表面等離激元由于其聚光特性可以在亞波長的近場范圍內顯著提高激子源與表面等離激元結構之間的耦合與能量傳遞效率，可以使許多光學過程的效率得到顯著的提高［1］.

特別是基于有序納米顆粒陣列的表面等離子體共振性能，人們可通過調節(jié)顆粒大小、間距乃至成分來大范圍地控制其特性.隨著表面等離子體理論研究的深入以及各種結構器件的成功制作，其在光學各領域應用具有巨大的潛力，尤其在一些經(jīng)典光學長期不能解決的問題上均有所突破，有序顆粒陣列在等離激元增強［2］、數(shù)據(jù)存儲［3］、生物化學傳感［4］、增強光催化性能［5］、單電子晶體管［6］和量子激光器［7］等多個領域有著重要的應用前景.2004年Okamoto等［8］首次在Natute Materials報道通過表面等離子體激元與InGaN/GaN量子阱之間的能量遷移，以增強半導體內的態(tài)密度并提高自發(fā)輻射效率，其熒光發(fā)射量子效率可提高約17倍，從而發(fā)光二極管的內量子效率得以增強;Atwater等［9］利用時域有限差分法(FDTD)的理論分析表明，在僅考慮散射效應的情況下，利用金屬納米顆粒的表面等離子體效應可使光在晶硅太陽電池中的傳播路程增大114倍.

由于表面等離子體共振是一種物理光學現(xiàn)象，需要滿足入射光波的頻率與金屬中自由電子的固有振動的頻率相等.因此文章主要研究銀納米三角有序陣列結構的表面等離子共振吸收特性，利用時域有限差分方法(Finite Difference Time Domain，F(xiàn)DTD)模擬計算不同尺寸的銀納米顆粒表面的電場強度和分布，分析入射激發(fā)光的波長、等離子體的共振頻率和納米顆粒大小等參數(shù)對等離子共振吸收波長和強度的影響.

1 數(shù)值模擬計算方法

目前，研究表面等離子體特性的數(shù)值模擬方法主要有:時域有限差分算法［10］，嚴格耦合波分析方法［11］，轉移矩陣方法［12］.其中時域有限差分算法是目前最成功的電磁數(shù)值方法之一.麥克斯韋方程組是分析研究金屬表面等離子體共振性能的理論基礎，應用于FDTD方法計算的Maxwell旋度方程和本構關系為:

1.1 實驗原型以及數(shù)值模擬的結構模型

該實驗中的主要材料是由PS球模板法制備的金屬Ag納米顆粒陣列，其陣列結構如圖1所示.

圖1 三角形銀納米顆粒陣列的掃描電鏡圖

隨著現(xiàn)代納米技術的發(fā)展，不同形貌的貴金屬納米顆粒陣列的制備及對它們相應光學特性的研究引起了人們廣泛的興趣.由于SPR光譜峰位對納米顆粒的形狀、大小、分布、外部介質環(huán)境的變化非常敏感，可通過顆粒大小、間距乃至成分來大范圍地控制其特性.實驗采用磁控濺射技術作為銀薄膜沉積方法，以自組裝制備出的按六角密排結構排列的PS球為掩膜板，在二氧化硅襯底上成功制備出三角形銀納米顆粒陣列結構.襯底為高純二氧化硅基底.基于FDTD方法模擬三角銀納米顆粒陳列的光學特性，計算了三角形銀鈉米顆粒陣列對光場的增強效應以及詳細討論了三角形銀納米顆粒的大小等參數(shù)的變化對表面等離子體共振場強分布的影響.

1.2 銀三角納米陣列的結構模型及FDTD參數(shù)

表面等離子體波是一種電磁表面波，與塊狀金屬對比，金屬納米粒子具有大得多的比表面積(單位質量物質所具有的總面積)，導體中傳導電子的運動往往受到粒子邊界的影響比較強烈，所以這種納米結構的光學性質可以認為由表面等離子體激元的性質來決定.由于金屬表面等離子體激元的響應，金屬納米粒子的表面能夠產(chǎn)生極大的電磁場增強，并且能夠對一定頻率的入射電磁波產(chǎn)生透射增強.

文章中設定三角形顆粒模型為三角柱形狀，按照圖1掃描電鏡的實驗結果設定模型具體參數(shù)，模擬環(huán)境中的三角形銀納米顆粒陣列如圖2所示.實驗分別制備了3種尺寸大小的三角形銀納米顆粒陣列，特征尺寸見表1.其中D為PS球的大小，L為三角形的邊長，H為三角形的高度.可以通過改變PS球的大小靈活地調整三角形的間距，通過改變沉積時間來控制三角性顆粒的高度，改變入射角可以控制形狀.

圖2 仿真模型橫截面三角形銀納米顆粒陣列的結構圖

表1 仿真模型橫截面三角形銀納米顆粒陣列的幾何尺寸大小 nm

通過在Y和X方向(垂直與紙面方向)加上周期性的邊界條件，就可以將這樣的一個金屬鏈沿著X，Y方向無限擴展，在Z軸方向采取PML邊界條件，從而模擬一個有限厚度的金屬三角形顆粒陣列的光學特性.運用FDTD solution軟件模擬了2個三角形型銀納米粒子靠近時的電耦合效應，計算中采用的模型如圖3所示，其中，納米陣列M在模擬中心.由于采取周期邊界，因此模型只需要覆蓋到FDTD的計算邊緣即可.R，T代表監(jiān)視器，分別用來記錄反射率、透射率.實驗采用的是平面波光源(TFSF光源)，平面波的范圍正好覆蓋Ag，波長范圍為300"800 nm.光源S的傳播方向(k vector)向下，而電場的偏振方向(E field vector)用AB雙箭頭表示.首先要分析的是傳輸特性和波長之間的關系.要分析它們之間的關系必須選擇合適的觀測點，并對相關的電磁場分量進行分析.

圖3 FDTD模擬界面及對三角形銀納米粒子模型

2 典型模擬結果及分析

實驗結合仿真的主要目的是探討三角形銀納米粒子的大小變化對表面等離子體共振場強分布的影響并利用仿真的結論為今后實驗提供理論依據(jù)和實驗指導.為了研究鄰近極化粒子對場增強作用的影響，選取了序號1中距離為260 nm的實驗作為例子，在Z方向距離樣品襯底h/2處放置監(jiān)視器，所觀察到的場增強效果如圖4所示.

圖4 銀納米陣列模型電場增強模擬圖

260 nm直徑的實驗中吸收峰位于波長478 nm附近，圖4是入射光波長在478 nm時三角形銀納米陣列的表面電場強度圖.從圖中可以看到沿Y軸方向2個三角形的對角尖銳處都出現(xiàn)了較強的電場增強效應，最大場增強放大到了1倍，這種增強只出現(xiàn)在與光源的偏振方向垂直的三角形對角處.文獻中將這些能夠在局域處獲得場增強的位置成為“熱點”［13］.模擬結果表明，在外電場的激發(fā)下具有尖銳角形狀的金屬納米粒子的表面等離子體局域場增強效應較為顯著，這是由于曲率半徑小的位置電荷密度相對較高，形成了較強的局域電磁場分布，而且由于其形狀的不對稱性，產(chǎn)生了多處近場增強的熱點.這個結論與文獻中運用有限元法模擬平面波激勵下三角形粒子所得到的場增強模擬結果相符合［14－15］.

圖5給出了不同顆粒尺寸的三角形銀納米顆粒的吸收譜線圖.數(shù)值計算結果表明，銀納米三角有序陣列的等離子共振吸收峰隨顆粒尺寸的增大而紅移，因此可通過改變顆粒大小調節(jié)共振吸收波長數(shù)值以用于不同的應用研究.

圖5 不同顆粒尺寸的三角形銀納米粒子吸收譜

3 結語

通過基于FDTD方法數(shù)值模擬了2個三角形型銀納米粒子靠近時的等離子共振吸收特性.從仿真結果可以看出，這種三角形型銀納米結構具有明顯的等離子共振增強效應，同時共振吸收波長受銀納米顆粒尺寸的影響非常大，其等離子共振吸收峰隨顆粒尺寸的增大而出現(xiàn)紅移現(xiàn)象.此結論不僅為今后實驗提供理論依據(jù)和實驗指導，還可通過模擬調節(jié)獲取不同共振吸收波長以用于不同的應用研究.

［1］吳民耀，劉威志.表面電漿子理論與模擬［J］.物理雙月刊，2006，28(2):486 －496.

［2］Xu W L，Zheng M J，Ding G Q，et al.Fabrication and optical properties of highly ordered ZnO nanodot arrays［J］.Chem Phys Lett，2005，411:37 －42.

［3］Chou S Y，Krauss P R，Zhang W，et al.Sub － 10 nm imprint lithography and applications［J］.J Vac Sci Technol，1997，B 15(6)，2897 －2904.

［4］Elghanian R，Storhoff J J，Mucic R C，et al.Selective colorimetric detection of polynucleotides based on the distancedependent optical properties of gold nanoparticles［J］.Science，1997，277:1078 －1081.

［5］Mathur A，Roy S S，Hazra Kiran S，et al.Growth of carbon nanotube arrays using nanosphere lithography and their application in field emission devices［J］.Diamond Relat Mater，2010，19:914 －917.

［6］Kurihara K，Namatsu H，Nagase M，et al.Room temperature operated single electron transistor fabricated by electron beam nanolithography［J］.Microelectronic Engineering，1997，35(1－4):261－264.

［7］Kim J，Benson O，Kan H，et al.A single-photon turnstile device［J］.Nature，1999，397(6719):500 －503.

［8］Okamoto K，Niki I，Shvartser A，et al.Surface-Plasmon-Enhanced Light Emitters Based on InGaN Quantum Wells［J］.Nat Mater，2004，3:601 －605.

［9］Atwater H A，Polman A.Plasmonics for improved photovoltaic devices［J］.Nature Mater，2010，9:205 －213.

［10］高本慶.時域有限差分法FDTD Method［M］.北京:國防工業(yè)出版社，1995.

［11］Moharam M G，Gaylord T K.Rigorous coupled-wave analysis of metallic surface-relief gratings［J］.Opt Soc Am，1986，A3:1780 －1787.

［12］曹莊琪.導波光學中的轉移矩陣方法［M］.上海:上海交通大學出版社，2000.

［13］Fang N，Zhang X.Imaging properties of a metamaterial superlens［J］.Appl Phys Lett，2003，82(2):161 － 163.

［14］Liu Z，F(xiàn)ang N，Yen T J，et al.Rapid growth of evanescent wave by a silver superlens［J］.Appl Phys Lett.，2003，83(25):5184－5186.

［15］Fischer U C，Zingsheim H P.Sub-microscopic pattern replication with visible-light［J］.J Vac Sci Technol，1984，19(4):881－885.

Simulation on Plasma Resonance Absorption Properties of Triangular Ag Nanoparticle Array Based on FDTD Method

MEI Fei1，YANG Zhi-jie1，XU Lei2，XU Jin-xia3，CHEN Zhi-qiang1

(1.Hubei University of Technology，Wuhan 430068，China;2.North China Institute of Water Conservancy and Hydroelectric Power，Zhengzhou 450011，China;3.Wuhan University，Wuhan 430072，China)

Plasma resonance absorption properties of the triangular Ag nanoparticle array are simulated by Finite Difference Time Domain(FDTD)method and the influence of the size of the triangular Ag nanoparticle on the plasma resonance are discussed.The results show that the absorption peak of the plasma resonance shifts to the longer wavelength with increasing the size of the triangular Ag nanoparticles，and vice versa.The theoretical analysis shows that the calculated results agree well with the experimental，and this well illustrates the relationship between the nanoparticle size and surface plasmon resonance absorption peaks.These approaches will surely facilitate further exploration of metal nanostructures for applications in nanoscale optical sensors.

Ag nanoparticle;FDTD;plasma resonance absorption

1002－5634(2012)03－0121－04

2012－03－18

梅 菲 (1980—)，女，湖北武漢人，講師，博士，主要從事半導體材料生長方面的研究.

(責任編輯:蔡洪濤)