亞波長微納光學的前沿研究 （一）＊

莊松林 王 琦 朱亦鳴 耿 滔 張大偉

①中國工程院院士，②博士，③⑤教授，④副教授，教育部光學儀器與系統(tǒng)工程研究中心；上海市現(xiàn)代光學系統(tǒng)重點實驗室；上海理工大學光電信息與計算機工程學院，上海200093＊上海市重點學科項目第三期項目（S30502）；上海市超精密光學加工與檢測服務平臺建設（11DZ2290301）

亞波長微納光學的前沿研究 （一）＊

莊松林①王 琦②朱亦鳴③耿 滔④張大偉⑤

①中國工程院院士，②博士，③⑤教授，④副教授，教育部光學儀器與系統(tǒng)工程研究中心；上海市現(xiàn)代光學系統(tǒng)重點實驗室；上海理工大學光電信息與計算機工程學院，上海200093
＊上海市重點學科項目第三期項目（S30502）；上海市超精密光學加工與檢測服務平臺建設（11DZ2290301）

負折射 反常多普勒效應 黑硅

材料與結(jié)構(gòu)在微納米尺度展現(xiàn)了許多不同于宏觀尺度的新特征，微納加工技術(shù)已經(jīng)成為當前科學研究與工業(yè)開發(fā)的熱門領(lǐng)域之一。筆者簡要介紹了負折射材料和黑硅這兩種微納光學材料的制備及其特性，展示了微納光學材料在新技術(shù)中的廣闊前景和科技創(chuàng)新中的重要作用。

1 概述

微納光學是指微納米尺度下的光學現(xiàn)象。隨著光學器件和光學系統(tǒng)不斷的微型化，其光學特性也會發(fā)生變化，當特征尺寸達到微納米量級時，會出現(xiàn)很多宏觀條件下沒有的新現(xiàn)象。負折射材料的出現(xiàn)，為人們提供控制光和電磁波行為的全新手段；而利用飛秒激光脈沖與硅材料相互作用獲得的黑硅材料，在太陽能電池等應用研究領(lǐng)域有著廣泛的應用前景。筆者簡要介紹其所在研究小組在特殊的光學材料——負折射材料和黑硅材料制備和特性研究方面的相關(guān)進展。

2 負折射材料及研究現(xiàn)狀

在自然界中，物質(zhì)的電磁性質(zhì)可以使用介電常數(shù)ε和磁導率μ來描述。根據(jù)麥克斯韋方程可知，平面電磁波的波矢k，電場強度E和磁場強度H三者兩兩垂直，遵守右手定則k＝E×H。電磁波的能量傳播方向由坡印廷矢量S決定，在一般電磁介質(zhì)中，坡印廷矢量S的方向和波矢k的方向相同。那么，是否存在一種特殊的介質(zhì)，當平面電磁波在其中傳播時，波矢k，電場強度E和磁場強度H三者不再遵守右手定則k＝E×H，而是遵守左手定則k＝－E×H，使得坡印廷矢量S的方向和波矢k的方向相反？蘇聯(lián)物理學家Veselago［1］最早對這一問題進行了研究，他假想了一種具有負的介電常數(shù)ε和負的磁導率μ的雙負材料，來實現(xiàn)這一構(gòu)想。由于這種材料遵守左手定則k＝－E×H，因此這種材料被稱為“左手材料”（left-handed material，LHM）。Veselago發(fā)現(xiàn)當平面電磁波從一般的右手材料入射到左手材料中時，在左手材料內(nèi)部會發(fā)生入射光線和折射光線位于法線同側(cè)的現(xiàn)象。此時，如果將Snell定律做一下推廣，認為左手材料的折射率為負值，這樣使得這種情形也滿足Snell定律的形式；因此，左手材料又稱為“負折射率材料”（negative index material，NIM）。

盡管負折射率材料體現(xiàn)出諸多的奇特性質(zhì)，但由于人們沒有在自然界中發(fā)現(xiàn)介電常數(shù)ε和磁導率μ同時為負值的材料，所以，Veselago的研究成果并沒有引起當時學術(shù)界的重視。直到20世紀末，英國物理學家Pendry重新開啟了這方面的工作。他在理論上仔細研究了導線陣列［2］和有缺口的環(huán)形共振器［3］（split-ring resonator，SRR）陣列的電磁特性，極大地推動了這一領(lǐng)域的發(fā)展。其后，Pendry又提出了“完美透鏡”（perfect lens）的概念［4］，即利用介電常數(shù)和磁導率同時為－1的負折射率平板，不僅可以聚焦光場的傳播場部分，更重要的是可以恢復光場的倏逝場部分，突破瑞利衍射極限成像，這一無可比擬的特性引起學術(shù)界極大的興趣和重視［5－10］。

在實驗上，美國物理學家Smith等［11］走出了關(guān)鍵的一步，他們把導線陣列和有缺口的環(huán)形共振器這兩種結(jié)構(gòu)做在一起，在微波段實現(xiàn)了在同一材料中介電常數(shù)ε和磁導率μ同時為負值。接下來，他們改進陣列結(jié)構(gòu)，并首次在實驗中觀察到負折射現(xiàn)象［12］。2003年，研究者對Smith等人的實驗進行了改進［13－14］，證明在遠場也能觀察到負折射現(xiàn)象，且不同入射角下測量到的負折射率是一致的，完全符合Snell定律。這一系列的實驗工作使負折射率是否真的存在這一爭論暫時告一段落，并被Science雜志評為當年十大科學進展之一。當然，由于實際的材料必然存在色散和吸收，理想的完美透鏡并不存在；但利用負折射率平板成像，仍然能夠突破傳統(tǒng)的衍射極限。因此，光波段的負折射率材料的應用會對高分辨率成像、光通訊、高密度光存儲、超大規(guī)模集成電路的光刻技術(shù)以及集成光學等領(lǐng)域產(chǎn)生深遠的影響。另外，負折射率的LC網(wǎng)絡可應用于耦合器、諧振器、天線和導波方面，能極大地改善這些器件的性能?？傊?，負折射率材料的出現(xiàn)，為人們提供了控制光和電磁波行為的全新手段。

在一般材料中，當波源和觀察者發(fā)生相對移動時，會發(fā)生多普勒效應：兩者相背而行，觀察者接收到的頻率低于波源頻率，如圖1（a）所示。但在負折率材料中，情況正好相反。如前所述，群速的方向（即能流傳播的方向）和相速的方向（即相位傳播的方向）相反，所以如果波源和觀察者相背而行，觀察者接收到的頻率會升高，發(fā)生反常多普勒效應，如圖1（b）所示。由于在實驗上很難將探測器放在人造的負折射率材料內(nèi)部測量頻移，所以大部分研究者都采用理論研究和數(shù)值模擬、仿真相結(jié)合的方法［15－18］。到目前為止，只有在音頻波段［19］和GHz波段［20］的間接測量報道，但一些科學家對實驗工作提出了不同看法［21－23］。因為多普勒效應現(xiàn)象早已被廣泛應用于工農(nóng)業(yè)生產(chǎn)、醫(yī)療衛(wèi)生、科技國防等等社會生活的許多方面，可以預見，反常多普勒效應的存在將會給這些應用領(lǐng)域帶來根本性的影響。另一方面，目前雖然在很多人工合成材料中觀察到了負折射現(xiàn)象，但這并不能作為評判材料折射率為負的依據(jù)。因為發(fā)生負折射現(xiàn)象不代表材料的折射率為負，例如全息光柵的衍射、光子晶體等頻面形狀的內(nèi)凹［24－25］等原因都會發(fā)生負折射現(xiàn)象，但該材料折射率并不為負。因此，反過來，反常多普勒效應的實驗測量將有力地證明負折射率材料的存在性，驗證相關(guān)理論的正確性，為后續(xù)工作鋪平道路。為此，我們搭建了自主設計的實驗裝置，利用自制的實驗平臺，在國際上首次觀測到了光頻段的反常多普勒現(xiàn)象［26］。

我們實際制備了負折射光子晶體棱鏡，Si介質(zhì)柱的截面半徑r＝0.2a，其中晶格常數(shù)a＝5μm。為了接近理想的二維光子晶體，在當前工藝許可的條件下，盡量增加介質(zhì)柱的高度，使得高度h＝50μm。為了測試光路的需要，光子晶體棱鏡樣品制作成菱形，頂角為60°，棱邊長5 mm，Si介質(zhì)柱分布在高度為550μm的Si襯底上。圖2（a）和2（b）分別給出實際樣品的俯視圖和側(cè)視圖。接下來對樣品的負折射效應做了實驗驗證，實驗入射光源選擇E偏振的CO2激光器，工作波長10.6 μm。實驗結(jié)果表明，所有具有相同晶格結(jié)構(gòu)的不同樣品（超過10塊樣品）都體現(xiàn)出明顯的負折射特性。在這里需要注意的是，根據(jù)實驗測量結(jié)果，樣品的等效折射率在－0.524～－0.376之間變化，絕對值小于－0.7的理論絕對值，我們認為這是由于Si柱的有限高度造成的。另外，我們制備的樣品下覆層材料為Si，沒有上覆層，Si柱暴露在空氣中，因此光子晶體平板上下覆層材料的不對稱，也會對光子晶體的能帶結(jié)構(gòu)產(chǎn)生影響［27－28］。

圖2 負折射光子晶體棱鏡的掃描電鏡圖（a）俯視圖；（b）側(cè)視圖

反常多普勒效應的測量有兩大困難：①探測器不能置于光子晶體介質(zhì)內(nèi)部，且相對于光子晶體介質(zhì)運動；②需要找到一個明確的判據(jù)來判斷所測頻移的正、負，以區(qū)分正常和反常多普勒效應。我們設計如圖3的實驗光路，CO2激光光束垂直于光子晶體棱鏡的第一個棱邊入射，即沿ΓM方向。光子晶體棱鏡固定在一個平移臺上，平移臺可沿x方向勻速運動。通過計算分析可知探測器所測拍頻為

其中f0為光源頻率；f′2為探測器接受到的參考光頻率；f1為測量光束經(jīng)過光子晶體棱鏡，到達棱鏡第二條棱邊時，棱鏡在界面上接受到的光波頻率；k為略小于1的正常量。式（1）中第一項f′2－f0k的值在同一次實驗過程中是一個常量，且與光子晶體棱鏡的等效折射率無關(guān)，給定平移速度后可以很容易算出。當實驗參數(shù)θ＝－26°，β＝45°，平移速率為0.012 3 mm／s，0.024 5 mm／s，0.048 8 mm／s和0.073 2 mm／s時，f′2－f0k均為正值，其值分別等于1.90 Hz，3.77 Hz，7.51 Hz和11.27 Hz。注意到參量k恒為正值，那么當測量到的拍頻Δf＜f′2－f0k時，f1－f0能且只能為正，即f1＞f0；也就是說，當光程在光子晶體棱鏡中增加時，在光子晶體棱鏡第二個界面上接受到的光波頻率大于CO2激光器的光源頻率，即在光子晶體棱鏡中發(fā)生的多普勒效應為反常。

圖3 實驗光路示意圖

圖4 實驗測量的頻差Δf與f′2－f0k值在不同平移速度下的對比

圖4中的插圖給出了在上述四個平移速度下，探測器測得的信號。通過快速傅立葉變換得到拍頻Δf分別為0.89 Hz，1.83 Hz，3.65 Hz和5.14 Hz，實驗測量值與理論值符合得很好，實驗值與理論值的誤差小于5%。將實驗測量值Δf與f′2－f0k值相比較，如圖4所示，實際測量的頻差Δf均小于f′2－f0k，這一結(jié)果清楚無誤地證明實驗觀測到了反常多普勒效應。使用具有相同晶格結(jié)構(gòu)的不同光子晶體棱鏡樣品，我們多次重復了上述實驗，實驗結(jié)果類似。將光子晶體棱鏡替換為普通的ZnSe棱鏡，做了同樣的實驗，實驗結(jié)果表明，ZnSe棱鏡中的多普勒效應是正常的，實驗測量值與理論值吻合得很好，誤差小于3%。這一驗證性實驗證明我們的實驗系統(tǒng)可以很好地區(qū)分正常和反常多普勒效應，且實驗結(jié)果可信，與理論值符合得很好。

3 黑硅材料的研究進展

近年來，利用飛秒激光在一定氣壓的氣體環(huán)境里照射硅材料表面，可以生成一種新型的“黑硅”材料。這些微結(jié)構(gòu)可以顯著改變硅材料的屬性［29－30］，特別是在很寬的波長范圍內(nèi)（0.2～2.5μm），光吸收效率可被顯著提高到90%以上［31－32］。因此這一特性可被廣泛應用于太陽能電池［33］、太赫茲輻射［34］、傳感器和光電探測器［35－38］的研究。為獲得更大的光電吸收系數(shù)，科學家們不斷致力于不同表面微結(jié)構(gòu)形貌的硅材料的研究，包括激光通量、脈 沖 寬 度［39－40］、偏 振 方 向［41］、脈 沖 數(shù)［42］、氣 體 種類［43］、氣體壓力［44］等。然而，由于很多參數(shù)的研究范圍不夠大，有些特殊的現(xiàn)象未被發(fā)現(xiàn)。我們注意到，不同波長的激光脈沖在硅表面上可以形成不同的波紋圖案，且對相同材料的穿透深度也不同，所以最后形成的表面微觀結(jié)構(gòu)也應有所不同。然而，不同波長飛秒激光脈沖蝕刻的微結(jié)構(gòu)的具體區(qū)別至今沒有完善的研究結(jié)果。我們研究小組詳細研究并討論了在兩種不同波長的飛秒激光脈沖的驅(qū)動下，通過改變實驗條件，包括激光功率、脈沖數(shù)，以及在相同的激光通量下它們之間的比例關(guān)系，硅表面微觀結(jié)構(gòu)變化的差異。

3.1 黑硅材料的實驗裝置

圖5 實驗裝置圖

制備表面微結(jié)構(gòu)硅材料的實驗裝置如圖5所示。由鈦寶石再生放大系統(tǒng)輸出飛秒光源，其輸出脈沖的中心波長為800 nm，重復頻率1 kHz，脈沖寬度為45 fs，激光光斑的空間分布為近高斯形。激光脈沖經(jīng)過焦距為f＝1 m的聚焦透鏡會聚到真空腔中，真空腔的入射窗口為厚度0.4 mm的窗片，腔內(nèi)充入氣壓可調(diào)的SF6氣體，真空腔被固定在一個三軸移動平臺上，實現(xiàn)三維精密調(diào)節(jié)。一個可翻轉(zhuǎn)的反射鏡則用于兩種不同光源之間的轉(zhuǎn)換，使其分別與樣品相互作用：首先，當該反射鏡向上翻轉(zhuǎn)，800 nm波長的激光束將直接入射到樣品上；將該反射鏡向下翻轉(zhuǎn)，800 nm波長的激光束先入射到0.2 mm厚的BBO倍頻晶體上，可以得到波長為400nm的飛秒激光脈沖。晶向為100的硅片垂直于光源的入射方向放置于真空腔內(nèi)，且位于光源焦點前方，以避免激光過強損傷窗片。硅片表面的激光光斑大小通過光學成像CCD來監(jiān)測確定。通過控制硅片和激光焦點之間的距離，始終保持樣品表面處光斑的直徑大小約為300 μm。為實現(xiàn)激光功率和脈沖數(shù)的可調(diào)性，在腔前還放置了可調(diào)衰減片和脈沖計數(shù)器。

3.2 實驗結(jié)果及討論

在背景氣體SF6氣壓為500 Torr（注：1 Torr＝1 mmHg＝133.322 Pa）時，我們研究了隨800 nm和400 nm激光功率的變化，激光強度、脈沖數(shù)以及它們的比例關(guān)系與產(chǎn)生的不同表面形貌硅的依賴關(guān)系。對應的硅表面上形成的平均尖峰高度的變化如圖6所示。

3.2.1 激光功率與激光波長的關(guān)系

圖6 波長（a）800 nm和（b）400 nm的激光光源蝕刻硅材料表面：（a1和b1）分別給出當脈沖數(shù)固定為1 000時，硅的平均峰高與激光功率之間的關(guān)系；（a2和b2）分別給出當800 nm與400 nm激光光源的功率分別為400 mW和200 mW時，平均峰高與脈沖數(shù)之間的關(guān)系；（a3和b3）給出在相同的激光通量2.8 kJ／m2下，激光功率與脈沖數(shù)的不同比例關(guān)系與微結(jié)構(gòu)硅的平均尖峰高度之間的關(guān)系

在實驗中，脈沖數(shù)均固定為1 000，改變激光功率觀測微結(jié)構(gòu)平均峰值高度的變化。從圖6（a1）和6（b1）我們可以看到，800 nm和400 nm激光脈沖所蝕刻的硅微觀結(jié)構(gòu)在峰高上有類似的變化趨勢，即尖峰高度隨脈沖能量的升高而先升高后降低。已經(jīng)知道，微觀結(jié)構(gòu)形成的整個過程可以分為三個階段［45］：首先，入射光和初始類似鏡面的硅表面的反射光之間的干涉，導致能量在硅表面的不均勻沉積。當入射能量超過硅材料燒蝕和熔化的閾值，將在非均勻的深度產(chǎn)生燒蝕和熔化，由此產(chǎn)生以激光波長為周期的張力波；其次，隨著燒蝕和熔化時間的增加，這些張力波逐漸成為波紋圖樣，然后變?yōu)闇手芷陉嚵械乃螤?；最后，水滴邊緣的物質(zhì)繼續(xù)燒蝕，形成錐形的微觀結(jié)構(gòu)。因此，通過適當增加激光的能量，可促進水滴邊緣材料的燒蝕，繼而形成更高的錐形結(jié)構(gòu)。

此外，從另一個角度可以看到，激光照射在硅表面，如果每個脈沖的能量均被樣品有效吸收，那么硅表面的能量密度將為最大，然后逐漸向深層降低。該密度分布導致頂層燒蝕剝離以及下面的材料熔化。在這個過程中，脈沖數(shù)代表激光與硅之間的反應時間，這就決定了能量能夠?qū)氲讲牧蟽?nèi)部的深度，而激光功率決定了硅表面材料的燒蝕和揮發(fā)程度。因此，對于固定的脈沖數(shù)，單位時間內(nèi)轉(zhuǎn)移到材料內(nèi)部的能量是固定的。如果激光的能量非常大，激光能量無法全部導入到硅材料的深層，有一部分能量可能會累積在硅表面。當積累的能量達到閾值時，表面材料會燒蝕過度變成液態(tài)硅。從而導致在初始的幾百個脈沖階段，不能有效形成水滴狀結(jié)構(gòu)。但是隨著時間的增加，能量逐漸導入材料內(nèi)部，那么后期的激光脈沖可以再次產(chǎn)生錐形結(jié)構(gòu)?？v觀整個過程，前一部分脈沖對尖峰結(jié)構(gòu)的形成并無貢獻，也就意味著形成尖峰的有效脈沖個數(shù)減少了。因此，當脈沖數(shù)固定時，隨著激光功率的增加，可以在尖峰高度隨激光功率的變化曲線上觀察到一個轉(zhuǎn)折點。此外，當激光功率超過脈沖數(shù)的能量傳輸限制時，激光功率越高，能量來不及傳導而導致的熔化效果會越明顯，有效的脈沖數(shù)越少，平均峰高也就越低。圖7為相應的掃描電子顯微鏡（SEM）照片。

圖7 脈沖數(shù)為1 000，波長為（a）800 nm 和（b）400 nm 的激光脈沖在不同激光功率（a1）400，（a2）1 000，（a3）1 400，（a4）1 800，（a5）2 000，（b1）50，（b2）125，（b3）175，（b4）275和（b5）325 mW 下，所形成的微結(jié)構(gòu)硅材料的掃描電子顯微鏡圖像。觀察角度均為傾斜45°

通過比較由800 nm和400 nm的激光脈沖形成的兩種微結(jié)構(gòu)硅材料，如圖6（a1）和6（b1），我們發(fā)現(xiàn)，在入射激光波長、激光功率和形成的峰高之間存在一個特殊關(guān)系：當800 nm的激光脈沖的能量是400 nm激光脈沖的8倍時，其形成的微結(jié)構(gòu)硅材料的平均峰高為400 nm激光脈沖形成的3倍。我們認為這個有趣的現(xiàn)象也可以通過微結(jié)構(gòu)形成的過程來分析。能量沉積不均勻形成以激光波長為周期的張力波，接著形成波紋圖樣、水滴的準周期陣列，以及隨后形成錐形結(jié)構(gòu)。在這些過程中，可以發(fā)現(xiàn)，圓錐尖峰形是由于波紋點間的微小空間的燒蝕和熔化產(chǎn)生的。根據(jù)激光輻照能量密度分布［46］，每個小空間可以近似為一個圓錐體，則這些圓錐體的單位體積的能量可以寫成：

另一方面，對于相同的硅材料形成的微結(jié)構(gòu)，其燒蝕和熔化閾值應該是相同的。因此，對于兩個波長不同的激光脈沖（JVλ＝JVλ），當其單位體積的能量相同，并

12且超過燒蝕熔化閾值時，硅表面可以同時形成微結(jié)構(gòu)。則式（1）可以寫成：

這個方程式表明，如果不考慮其他因素影響，對于一定波長和吸收效率的激光脈沖，在輻射激光能量與形成峰高之間有一個比例關(guān)系。利用這個確定的比例關(guān)系，可以預測不同波長下的激光脈沖形成的峰高。實驗中，我們用800 nm和400 nm激光脈沖來形成微結(jié)構(gòu)的硅材料。因此，通過代入?yún)?shù)，A800≈60%，A400≈40%［39］，R800＝2R400，則式（2）可以寫成：

其中，J800和J400分別表示800 nm 和400 nm 激光脈沖的入射能量，d800和d400分別表示800 nm和400 nm激光脈沖的穿透深度。由式（3）可證明，當800 nm激光脈沖的能量是400 nm激光脈沖能量的8倍時，所形成的平均峰高將為400 nm的3倍。此數(shù)值分析與我們的實驗結(jié)果非常吻合［49］。

3.2.2 脈沖數(shù)與激光波長的關(guān)系

用800 nm和400 nm激光脈沖在硅材料表面蝕刻，所形成的峰高變化的平均值與脈沖數(shù)也將構(gòu)成一個函數(shù)關(guān)系，如圖6（a2）和圖6（b2）所示。我們可以看到，800 nm激光脈沖所形成的平均峰高隨著脈沖數(shù)的增加而單調(diào)增加，這是由硅表面微結(jié)構(gòu)的形成機制決定的。我們知道，隨著脈沖數(shù)的增加，激光和硅材料的相互作用時間也增加，并且更多的能量導入到材料內(nèi)部，這會導致表面以下的材料熔化，促進形成錐形尖峰［46］。因此，峰高隨著脈沖數(shù)的增加而單調(diào)增加。對于400 nm激光脈沖，不同的是，其形成峰高與脈沖數(shù)之間構(gòu)成一個非線性關(guān)系。通過對比800 nm的激光脈沖，400 nm激光脈沖蝕刻尖峰的形成過程可以分為兩階段：當脈沖數(shù)小于1 000時，峰高的增長速度更快些；當脈沖數(shù)大于1 000時，其增長速度將更慢些。

至于使用400 nm激光脈沖形成尖峰的第一階段，我們推測其原因是在硅材料內(nèi)入射激光脈沖單位體積的能量與激光波長成反比關(guān)系［式（2）］。當激光脈沖入射到硅材料表面，能量被硅材料吸收并聚集在一個錐形區(qū)域。對于波長較短的激光脈沖，其錐形區(qū)域更小，因此更容易達到燒蝕和熔化閾值并且形成尖峰微觀結(jié)構(gòu)的速度越快。

至于第二階段，我們通常認為，當光照射到錐形微結(jié)構(gòu)上時，光被尖峰斜邊和頂上的球反射進尖峰之間的峰谷，提高了入射光通量和燒蝕速率［46］；然后，尖峰可變得更加銳利，并且尖峰高度也相應地增加。正如在第一部分所討論的，分布能量的體積與激光波長的三次方成反比［式（2）］，因此，400 nm激光脈沖蝕刻得到的微結(jié)構(gòu)硅材料的谷間燒蝕速度會比800 nm激光脈沖的慢。當尖峰有效形成之后，400 nm激光脈沖的峰高增長的速率將比800 nm的激光脈沖的慢，與我們實驗中所用的脈沖數(shù)超過1 000時的情況相吻合。

3.3.3 相同激光通量下激光功率和脈沖數(shù)的比例關(guān)系與激光波長的關(guān)系

激光通量Φ的表達式可以寫為Φ＝（P·t）／S＝（P·m）／（f·S），其中，P為激光平均功率，t為激光脈沖的作用時間，S為照射在硅表面的激光光斑面積，m為入射到硅表面的脈沖數(shù)目，f為激光的重復頻率，即激光通量的數(shù)值正比于功率和脈沖數(shù)兩者的乘積。此式表明激光通量與激光功率和脈沖數(shù)成正比。此外，根據(jù)上面的討論，我們知道，無論激光功率還是脈沖數(shù)，都可以有效地影響硅表面微結(jié)構(gòu)的形成。然而，決定硅材料的表面形貌和吸收屬性的主要因素是激光功率還是脈沖數(shù)，仍然不能確定。因此，我們在相同的激光通量2.8 kJ／m2下，按照不同激光功率和脈沖數(shù)的組合蝕刻硅表面，測量和分析結(jié)果如圖6（a3）和圖6（b3），相應的硅表面微結(jié)構(gòu)SEM圖見圖8。顯然可以看出，當激光功率在不超過600 mW的情況下，隨著激光功率增加，脈沖數(shù)減少，尖峰高度和尖峰間隔逐漸增加。當激光功率超過800 mW時，在激光脈沖刻蝕區(qū)域的中央部分出現(xiàn)了一個平坦的區(qū)域［圖8（a2）］。當增大激光功率并且進一步降低脈沖數(shù)時，平坦的區(qū)域范圍逐漸增大，然后在激光功率為1 200 mW，脈沖數(shù)為167時平坦區(qū)域達到最大［圖8（a3）］。有趣的是，在相同的激光通量下，當激光功率大于1 200 mW，并且脈沖數(shù)小于167時，平坦區(qū)域范圍逐步縮小并最終在激光功率為1 800 mW和脈沖數(shù)為111時消失［圖8（a4）］。

圖8 波長（a）800 nm 和（b）400 nm 的激光光源蝕刻硅表面，在相同激光通量下（2.8 kJ／m2）；SEM 所觀察到的硅表面微結(jié)構(gòu)圖。觀察角度均為傾斜45°。（a1）功率400 mW，脈沖數(shù)500；（a2）功率800 mW，脈沖數(shù)250；（a3）功率1 200 mW，脈沖數(shù)167；（a4）功率1 800 mW，脈沖數(shù)111；（b1）功率50 mW，脈沖數(shù)4 000；（b2）功率100 mW，脈沖數(shù)2 000；（b3）功率200 mW，脈沖數(shù)1 000；（b4）功率300 mW，脈沖數(shù)667

從硅表面微結(jié)構(gòu)的形成過程中，我們知道，脈沖數(shù)表征激光和硅材料相互反應的時間，這就決定了能量能夠轉(zhuǎn)移到材料內(nèi)部的深度；而激光功率決定了硅表面材料的燒蝕度和揮發(fā)度。只有合適的燒蝕和足夠的相互作用時間才可以形成高尖峰的微結(jié)構(gòu)。因此，在相同的激光通量下，我們在實驗中觀察到的現(xiàn)象可以從兩個方面理解：第一，脈沖的數(shù)量足夠大時，激光能量可以及時轉(zhuǎn)移到硅材料內(nèi)部；因此，在對表面材料的激光燒蝕和揮發(fā)作用下，可以有效形成尖峰，并且尖峰高度可以隨著激光功率的增加而增加。第二，如果脈沖數(shù)不足，則難以將所有的能量導入硅材料內(nèi)部，部分能量將積累在表面上；由于積累的能量仍然可以熔化表面附近的材料，則在硅表面上也能觀察到一些平坦的區(qū)域。當脈沖數(shù)減少到一定值，熔化的量很少，則平坦的區(qū)域?qū)⒃俅蜗?。同時，脈沖數(shù)的下降也引起尖峰高度的下降。

在如圖8的實驗中，根據(jù)上述的討論，只有在激光功率不變的情況下只增加脈沖數(shù)，積聚在硅表面的能量才能被很好地導入到材料內(nèi)部，從而平坦區(qū)域消失并且尖峰高度比相同的激光功率下脈沖數(shù)較少時蝕刻的尖峰高度高。如果激光功率很高并且脈沖數(shù)過多，則材料可能因過度燒蝕而揮發(fā)，這可能會在蝕刻區(qū)域形成一個大洞。為了證明這一點，我們做了相應的實驗，結(jié)果如圖9所示。在800 nm激光脈沖的情況下，以激光功率為1 000 mW 為例：當脈沖數(shù)從167增加到1 000的過程中，平坦的區(qū)域完全消失［圖9（a1）和9（a2）］。這直接證明了脈沖數(shù)影響能量轉(zhuǎn)移的深度，進而決定了硅表面上形成微結(jié)構(gòu)的形態(tài)。另一方面，當脈沖數(shù)增加到4 000，超過了可以轉(zhuǎn)移到材料內(nèi)部的最大能量，高的激光功率導致材料表面的不斷燒蝕和揮發(fā)，最終在蝕刻區(qū)域形成一個洞［圖9（a3）］。作為對比，對于400 nm的激光脈沖，激光功率選200 mW，當脈沖數(shù)從1 000增加到2 000的過程中，平均尖峰高度增加。當脈沖數(shù)增加至6 000時，由于過度燒蝕和揮發(fā)，類似的孔再次出現(xiàn)?？紤]所有這些實驗的結(jié)果，可以得出這樣的結(jié)論，在激光功率和脈沖數(shù)之間存在著一種最理想的比例關(guān)系，以形成微結(jié)構(gòu)表面的硅材料［50］。

圖9 波長（a）800 nm 和（b）400 nm 的激光光源，對應的激光功率分別為1 000 mW 和200 mW，在不同脈沖數(shù)（a1）167，（a2）1 000，（3）4 000，（b1）1 000，（b2）2 000，（b3）6 000下蝕刻硅材料，用SEM所觀察到的硅表面微結(jié)構(gòu)圖。觀察角度均為傾斜45°

為證明脈沖激光功率和脈沖數(shù)的比例關(guān)系對于硅微結(jié)構(gòu)的光吸收能力的影響，我們測量了在各種激光功率和脈沖數(shù)的不同組合下蝕刻的樣品的吸收率（Lambda 750S，PerkinElmer）［51］，如圖10所示?？梢钥吹?，由激光功率與脈沖數(shù)的不同組合所制備的微結(jié)構(gòu)硅材料，其光吸收能力也有所不同，其中平均峰值最高的材料對應的吸收效率也最高。這一結(jié)果進一步證明在相同激光通量下，恰當?shù)剡x擇激光功率和脈沖數(shù)能有效形成表面具有較高峰值高度的微結(jié)構(gòu)硅，從而有利于高吸收效率的微結(jié)構(gòu)硅材料的獲得［52］。

總的來說，本文通過實驗，研究了兩個波長800 nm和400 nm的激光脈沖在一系列的形成條件的變化范圍內(nèi)，包括激光功率、脈沖數(shù)和在相同激光通量下它們的比例關(guān)系，分別蝕刻出不同微結(jié)構(gòu)的硅材料。實驗證明：脈沖數(shù)表征激光與硅之間的作用時間長短，即能量傳遞到材料內(nèi)部的深淺度；而激光功率表征硅表面材料的蝕刻燒蝕程度。此外，不同波長的激光脈沖有不同的能量沉積結(jié)構(gòu)并且能量穿透的深度也不相同?；谶@些，利用波長為800 nm和400 nm的激光脈沖蝕刻硅材料，在保持相同的激光通量下，分別通過改變激光功率、脈沖數(shù)以及它們的比例關(guān)系，所形成的微結(jié)構(gòu)硅材料有類似的變化趨勢，但是其微結(jié)構(gòu)密度和平均尖峰高度都不相同。

圖10 波長為（a）800 nm 和（b）400 nm 的入射飛秒激光脈沖，保持相同的激光通量，在激光功率和脈沖數(shù)的比例不同的組合下，形成的微結(jié)構(gòu)硅材料的吸收效率曲線圖

4 小結(jié)

光波段的負折射材料的應用會對高分辨成像、光通訊、高密度光存儲、超大規(guī)模集成電路的光刻技術(shù)以及集成光學等領(lǐng)域產(chǎn)生深遠的影響。負折射材料的出現(xiàn)，為人們提供控制光和電磁波行為的全新手段。利用飛秒激光脈沖分別與硅材料相互作用可以獲得在寬光譜范圍內(nèi)均具有高吸收效率的黑硅材料，這一結(jié)果對表面微結(jié)構(gòu)材料的有效制備、表面形態(tài)控制，及其相應的光電屬性等都有著重要的意義，在太陽能電池等應用研究領(lǐng)域有著廣泛的應用前景。

（2012年3月26日收到）

［1］VESELAGO V G.The electrodynamics of substances with simultaneously negative values of permittivity and permeability［J］.Sov Phys Usp，1968，10：509-514.

［2］PENDRY J B，HOLDEN A J，STEWART W J，et al.Extremely low frequency plasmons in metallic mesostructures［J］.Phys Rev Lett，1996，76：4773.

［3］PENDRY J B，HOLDEN A J，ROBBINS D J，et al.Magnetism from conductors and enhanced nonlinear phenomena［J］.IEEE Transactions on Microwave Theory and Techniques，1999，47：2075-2084.

［4］PENDRY J B.Negative refraction makes a perfect lens［J］.Phys Rev Lett，2000，85：3966-3969.

［5］WILLIAMS J M.Some problems with negative refraction［J］.Phys Rev Lett，2001，87：249703（4）.

［6］GARCIA N，NIETO-VESPERINAS M.Left-handed materials do not make a perfect lens［J］.Phys Rev Lett，2002，88：207403（4）.

［7］SHEN J T，PLATZMAN P M.Near field imaging with negative dielectric constant lenses［J］.Appl Phys Lett，2002，80：3286-3288.

［8］YE Z.Optical transmission and reflection of perfect lenses by left handed materials［J］.Phys Rev B，2003，67：193106.

［9］TAUBNERT，KOROBKIN D，URZHUMOV Y，et al.Nearfield microscopy through a SiC superlens［J］.Science，2006，313：1595.

［10］LARKIN I A，STOCKMAN M I.Imperfect perfect lens［J］.Nano Lett，2005，5：339-343.

［11］SMITH D R，PADILLA W J，VIER D C，et al.Composite medium with simultaneously negative permeability and permittivity［J］.Phys Rev Lett，2000，84：4184-4187.

［12］SHELBY R A，SMITH D R，SCHULTZ S.Experimental verification of a negative index of refraction ［J］.Science，2001，292：77-79.

［13］PARAZZOLI C G，GREEGOR R B，LI K，et al.Experimental verification and simulation of negative index of refraction using Snell’s law ［J］.Phys Rev Lett，2003，90：107401（4）.

［14］HOUCK A A，BROCK J B，CHUANG I L.Experimental observations of a left-handed material that obeys Snell’s law［J］.Phys Rev Lett，2003，90（13）：137401（4）.

［15］LEONG K M K H，LAI A，ITOH T.Demonstration of reverse Doppler effect using a left-handed transmission line［J］.Micro Optic Tech Lett，2006，48（3）：545-547.

［16］KOZYREV A B，VAN DER WEIDE D W.Explanation of the inverse Doppler effect observed in nonlinear transmission lines［J］.Phys Rev Lett，2005，94：203902.

［17］LUO C，IBANESCU M，REED E J，et al.Doppler radiation emitted by an oscillating dipole moving inside a photonic band-Gap crystal［J］.Phys Rev Lett，2006，96：043903

［18］KATS A V，SAVEL’EV S，YAMPOL’SKII V A，et al.Lefthanded interfaces for electromagnetic surface waves ［J］.Phys Rev Lett，2007，98：073901.

［19］SEDDON N，BEARPARK T.Observation of the inverse Doppler effect［J］.Science，2003，302：1537-1540.

［20］HU X，HANG Z，LI J，et al.Anomalous Doppler effects in photonic band gaps［J］.Phys Rev E，2006，73：015602.

［21］REED E J，SOLJACIC M，IBANESCU M，et al.Comment on“Observation of the inverse Doppler effect”［J］.Science，2004，305：778b.

［22］KOZYREV A B，VAN DER WEIDE D W.Explanation of the inverse Doppler effect observed in nonlinear transmission lines［J］.Phys Rev Lett，2005，94：203902.

［23］REED E J，SOLJACIC M，JOANNOPOULOS J D.Comment on“Explanation of the inverse Doppler effect observed in nonlinear transmission lines”［J］.Phys Rev Lett，2006，96：069402.

［24］NOTOMI M.Theory of light propagation in strongly modulated photonic crystals：refractionlike behavior in the vicinity of the photonic band gap［J］.Phys Rev B，2000，62：10696-10705.

［25］LUO C，JOHNSON S G，JOANNOPOULOS J D，et al.Allangle negative refraction without negative effective index［J］.Phys Rev B，2002，65：201104（4）.

［26］CHEN J B，WANG Y，JIA B H，et al.Observation of the inverse Doppler effect in negative-index materials at optical frequencies［J］.Nature Photonics，2011，5：239-245.

［27］JOHNSON S G，F(xiàn)AN S H，VILLENEUVE P R，et al.Guided modes in photonic crystal slabs ［J］.Phys Rev B，1999，60：5751-5758.

［28］SФNDERGAARD T，BJARKLEV A，KRISTENSEN M，et al.Designing finite-height two-dimensional photonic crystal waveguides［J］.App Phys Lett，2000，77（6）：785-787.

［29］HER T，F(xiàn)INLAY R J，WU C，et al.Microstructuring of silicon with femtosecond laser pulses［J］.Appl Phys Lett，1998，73（12）：1673-1675.

［30］STUBENRAUCH M，F(xiàn)ISCHER M，KREMIN C，et al.Black silicon—new functionalities in microsystems［J］.J Micromech Microeng，2006，16，S82-S87.

［31］SERPENGüZEL A，KURT A，INAN?I，et al.Luminescence of black silicon［J］.Journal of Nanophotonics，2008，2：021770.

［32］YOUNKIN R，CAREY J E，MAZUR E，et al.Infrared absorption by conical silicon microstructures made in a variety of background gases using femtosecond-laser pulses ［J］.J Appl Phys，2003，93（5）：2626-2629.

［33］YUAN H，YOST V E，PAGE M R，et al.Efficient black silicon solar cell with a density－graded nanoporous surface：optical properties，performance limitations，and design rules［J］.Appl Phys Lett，2009，95（12）：123501.

［34］HOYER P，THEUER M，BEIGANG R，et al.Terahertz emission from black silicon［J］.Appl Phys Lett，2008，93（9）：091106.

［35］HUANG Z，CAREY J E，LIU M，et al.Microstructured silicon photodetector［J］.Appl Phys Lett，2006，89（3）：033506.

［36］BRANZ H M，YOST V E，WARD S，et al.Nanostructured black silicon and the optical reflectance of graded-density surfaces［J］.Appl Phys Lett，2009，94（23）：231121.

［37］MALONEY P G，SMITH P，KING V，et al.Emissivity of microstructured silicon ［J］.Appl Opt，2010，49（7）：1065-1068.

［38］MYERS R A，F(xiàn)ARRELL R，KARGER A M，et al.Enhancing near-infrared avalanche photodiode performance by femtosecond laser microstructuring ［J］.Appl Opt，2006，45（35）：8825-8831.

［39］CROUCH C H，CAREY J E，SHEN M，et al.Infrared absorption by sulfur-doped silicon formed by femtosecond laser irradiation［J］.Appl Phys A，2004，79（7）：1635-1641.

［40］PENG Y，WEN Y，ZHANG D S，et al.Optimal proportional relation between laser power and pulse number for the fabrication of surface-microstructured silicon ［J］.Appl Opt，2011，50（24）：4765-4768.

［41］CROUCH C H，CAREY J E，WARRENDER J M，et al.Comparison of structure and properties of femtosecond and nanosecond laser-structured silicon ［J］.Appl Phys Lett，2004，84（11）：1850-1852.

［42］ZORBA V，BOUKOS N，ZERGIOTI I，et al.Ultraviolet femtosecond，picosecond and nanosecond laser microstructuring of silicon：structural and optical properties ［J］.Appl Opt，2008，47（11）：1846-1850.

［43］ZHU J，SHEN Y，LI W，et al.Effect of polarization on femtosecond laser pulses structuring silicon surface［J］.Appl Surf Sci，2006，252：2752-2756.

［44］FOWLKES J D，PEDRAZA A J，LOWNDES D H.Microstructural evolution of laser-exposed silicon targets in SF6 atmospheres［J］.Appl Phys Lett，2000，77（11）：1629-1631.

［45］SHEN M Y，CROUCH C H，CAREY J E，et al.Femtosecond laser-induced formation of submicrometer spikes on silicon in water［J］.Appl Phys Lett，2004，85（23）：5694-5696.

［46］SHEEHY M A，WINSTON L，CAREY J E，et al.Role of the background gas in the morphology and optical properties of laser-microstructured silicon ［J］.Chem Mater，2005，17（14）：3582-3586.

［47］TULL B R，CAREY J E，MAZUR E，et al.Silicon surface morphologies after femtosecond laser irradiation ［J］.Mrs Bulletin，2006，31：626-633.

［48］DASH W C，NEWMAN R.Surface barrier analysis for metals by means of Schottky deviations［J］.Phys Rev，1955，99（4）：1151-1160.

［49］SARITAS M，MCKELL H D.Absorption coefficient of Si in the wavelength region between 0.80-1.16μm ［J］.J Appl Phys，2001，61（10）：4923-4925.

［50］PENG Y，ZHANG D S，CHEN H Y，et al.Optimal proportional relation between laser power and pulse number for the fabrication of surface-microstructured silicon［J］.Appl Opt，2012，51（5）：635-639.

［51］WU C，CROUCH C H，ZHAO L，et al.Near-unity belowband-gap absorption by microstructured silicon ［J］.Appl Phys Lett，2001，78（13）：1850-1852.

［52］彭滟，溫雅，張冬生，等.飛秒激光功率與脈沖數(shù)的比例關(guān)系對制備硅表面微結(jié)構(gòu)的影響 ［J］.中國激光，2011，38（12）：1203005.

Research on Subwavelength Micro-nano Optics（Ⅰ）

ZHUANG Song-lin①， WANG Qi②，ZHU Yiming③，GENG Tao④，ZHANG Da-wei⑤
①CAEMember，②Ph.D.，③⑤Professor，④AssociateProfessor，EngineeringResearchCenterofOpticalInstrumentandSystem，Ministry ofEducation；ShanghaiKeyLabofModernOpticalSystem；Schoolof Optical-ElectricalandComputerEngineering，UniversityofShanghaifor ScienceandTechnology，Shanghai200093，China

Materials and structures in micro and nano scale present many new characteristics different from the macroscopic scale.Micro-nano processing technology has been one of the most popular fields in the current scientific research and industrial development.The preparation and special characteristics of negative refractive material and black silicon are introduced.Micro-nano optical devices play important role in new technology and broad prospects of science and technology innovation.

negative refractive，inverse Doppler effect，black silicon

10.3969／j.issn.0253-9608.2012.04.001

（編輯：沈美芳）