表面曲率對石英毛細(xì)管微瓶模式及傳感特性影響

張 峰，吳根柱,2*，汪成程

(1.浙江師范大學(xué) 數(shù)理與信息工程學(xué)院，金華 321004；2.浙江師范大學(xué) 行知學(xué)院，金華 321004)

引 言

近十幾年來,回音壁模式(whispering gallery mode,WGM)光學(xué)微腔憑借著高Q值、模式體積小、制作方便等特點(diǎn)引起人們青睞?；匾舯谀Ｊ轿⑶皇侵敢话愠叽鐬槲⒚琢考壍膱A對稱波導(dǎo)結(jié)構(gòu)，光在微腔內(nèi)部以全反射的形式傳播，使得光被約束在微腔內(nèi)部繞著赤道運(yùn)行，當(dāng)光程為波長的整數(shù)倍即滿足相位匹配條件時(shí)，會(huì)產(chǎn)生等間距分立的共振模，這種電磁場模式也就是回音壁模式。回音壁模式微腔在非線性光學(xué)、光力學(xué)、腔量子電動(dòng)力學(xué)、高靈敏折射率傳感等方面有廣泛應(yīng)用前景。迄今為止人們設(shè)計(jì)制作了多種結(jié)構(gòu)回音壁模式微腔，包括微球型、微環(huán)型、微盤型、微環(huán)芯型、微管型等。最近幾年人們提出一種新型微腔，即微泡型微腔[1-3]，該微腔是種混合型微腔，具有不同的介電層，倏逝光場會(huì)耦合進(jìn)核芯中的液體，而大部分能量仍留在壁中，因此,該結(jié)構(gòu)仍保持著高Q值同時(shí)又具有很高的傳感靈敏度。

微泡型微腔可由CO2激光對內(nèi)部加氣壓的石英毛細(xì)管進(jìn)行局部加熱熔融而制成[4-6]。通過調(diào)節(jié)加熱區(qū)域和加熱功率以及內(nèi)部氣壓，可制成一個(gè)特定大小和壁厚的微球或有一定曲率的微瓶形狀結(jié)構(gòu)，該結(jié)構(gòu)可具有單柄或雙柄[7]，然后利用注射泵往微泡內(nèi)注入液體，液體內(nèi)物質(zhì)分子與倏逝波發(fā)生耦合而影響內(nèi)壁內(nèi)回音壁模式特性，從而通過對微腔諧振波長變化的檢測可以間接實(shí)現(xiàn)對微腔內(nèi)待測物的檢測。微泡的幾何尺寸對其模式特征影響非常大，曾有不少學(xué)者[8-10]對球形和橢球形微泡結(jié)構(gòu)進(jìn)行了仿真模擬，共同特點(diǎn)是利用有限元方法模擬分析了微泡膨脹半徑和壁厚對其Q值、有效折射率等參量以及溫度、壓力和折射率傳感性能等影響，而較少涉及到有關(guān)微泡曲率的影響，但實(shí)際制作過程中可能無法避免形成壁內(nèi)外不同曲率的微瓶狀微泡，為此需要了解內(nèi)外表面曲率對微泡模式性能影響，基于此，本文中使用有限元方法分析研究了不同曲率對微瓶狀空芯和液芯微氣泡回音壁模式特征及其傳感特性的影響，以便達(dá)到最佳微泡結(jié)構(gòu)設(shè)計(jì)目的。

1 物理模型

典型微瓶狀微泡結(jié)構(gòu)示意圖如圖1a所示，由一段中間略鼓起的毛細(xì)管構(gòu)成。圖1b和圖1c是它的截面圖。

Fig.1 Schematic of microbottle

a—3-D geometry of microbottle b—cross section of microbottle c—longitudinal section of microbottle

過z軸的截面輪廓滿足[6,9,11]：R2(z)=R20[1-(Δkz)2/2]，可簡化成類似一條拋物線，其中R2(z)為微瓶的外徑，R20(z)為在z=0處的外徑，Δk是微瓶的曲率，Δk越大則曲率越大。為了區(qū)分壁內(nèi)外微瓶的表面曲率，定義Δk1為外壁拋物線的曲率，Δk2為內(nèi)壁拋物線的曲率，當(dāng)Δk1一定值時(shí)，通過改變?chǔ)2，可以實(shí)現(xiàn)不同曲率情況下的微瓶模型。

2 數(shù)值仿真及分析

選用λ=1.55μm通訊工作波長，首先通過有限元數(shù)值計(jì)算方法求解弱項(xiàng)亥姆霍茲方程[8,12-13]：

(1)

Fig.2 Mode field distribution of microbottle for different modes withR=25μm,t=1.4μm,Δk12=0.0035，Δk22=0.006

a—the first order mode b—the second order mode c—the second tangential mode

回音壁模式的主要損耗包括輻射損耗(倏逝波損耗)、材料吸收損耗和表面粗糙損耗。采用熔融方法能有效降低表面粗糙損耗，為了模擬實(shí)際情況，把材料吸收作為一個(gè)附加的虛部εi=-3.56×10-10引入材料介電常數(shù)。并在模擬的時(shí)候引入一個(gè)完全匹配層(perfectly matched layers，PML)，作為各向異性吸收器，能實(shí)現(xiàn)在有限計(jì)算空間內(nèi)模擬輻射到無窮遠(yuǎn)倏逝波。品質(zhì)因子Q的計(jì)算公式為:

(2)

式中,fr是本征(復(fù)數(shù))頻率。實(shí)部表示諧振頻率，虛部表示本征損失。Q值越大說明微腔儲(chǔ)存光子能量能力越強(qiáng)，反之則越弱，是衡量微腔性能優(yōu)劣的一個(gè)重要參量。先以空氣為核芯介質(zhì)研究微瓶回音壁模式特征。

由參考文獻(xiàn)[10]可知，當(dāng)微泡的壁厚為1μm時(shí)更有利于傳感，能保持較高分辨率和靈敏度，且不易破裂，易于制作。為此這里選用赤道面處壁厚為1μm，外壁曲率平方Δk12=0.0035時(shí)，內(nèi)壁曲率的平方Δk22=0.005時(shí)微瓶。在不同切向模式下,品質(zhì)因子Q隨膨脹半徑R變化趨勢如圖3所示。膨脹半徑R<19μm時(shí)，3種切向模式Q值均隨膨脹半徑R增大而增大，這是由于此時(shí)微瓶尺寸較小，輻射損耗占主導(dǎo)，材料損耗相對比較小可以忽略。當(dāng)膨脹半徑R>19μm時(shí)，3種切向模式下微瓶Q值均趨于飽和，此時(shí)膨脹半徑較大，輻射損耗較小，模式的Q值主要由材料損耗決定。切向高階模(如二階模和三階模)Q值始終高于一階模，且先達(dá)到飽和，這是由于高階模式下模場中心不全部在赤道面處，受壁厚影響，一般情況下壁厚越大Q值也越高[8]。

Fig.3 Relationship betweenQandRwhent=1μm, Δk12=0.0035, Δk22=0.005

圖4是在膨脹半徑R=12μm時(shí)，切向不同模式下Q值隨著曲率變化情況。圖4a中，赤道面處壁厚為1μm，固定外壁的曲率平方Δk12=0.0035，可以看出，相同半徑時(shí)Q值并不總是隨曲率的增大而增大。當(dāng)膨脹半徑R較小時(shí)，如R=12μm,切向二階模存在一個(gè)Q值極大的點(diǎn)，此時(shí)內(nèi)壁的曲率平方Δk22=0.008，之后Q值隨Δk22的變大而降低;切向三階模在R=12μm時(shí)同樣具有類似的規(guī)律，Δk22<0.007時(shí)Q值單調(diào)遞增，之后Q值隨Δk22變大而降低。對比3種模式可以看出，當(dāng)曲率改變時(shí)，切向高階模Q值變化大于一階模。這是由于內(nèi)壁曲率變大時(shí)，切向高階模(處于赤道面以外)所在壁厚變厚，從而模式Q值變大，而Q值受輻射損耗和材料損耗共同影響，高階模處于非赤道面處勢必導(dǎo)致模式所在半徑發(fā)生變化，從而輻射損耗和材料損耗存在一個(gè)折中產(chǎn)生模式最佳曲率，有一個(gè)極大值。而一階模處于赤道面處，所在壁厚始終保持不變，故受曲率的影響較小。圖4b中，赤道面處壁厚為1μm，固定外壁的曲率平方Δk12=0.007，顯然與圖4a有相同的趨勢，其中星號點(diǎn)代表RΔk=1即微瓶成為微氣泡[9]，一階模Q值高達(dá)5.5×106。

Fig.4 Relationship between Q and Δk22a—Δk12=0.0035 b—Δk12=0.007

液芯微泡的研究更有意義[14]，下面討論液芯下微瓶狀微泡回音壁模式特征。由參考文獻(xiàn) [8]～參考文獻(xiàn) [10]可知，曲率Δk的不同對其模場分布及諧振波長偏移有很大的影響；且液芯時(shí)，微瓶的壁厚對其模式影響很大，當(dāng)微瓶壁非常厚時(shí)，模式大部分電磁場在壁中傳播；隨著壁厚減小，更多電磁場耦合進(jìn)核芯中。這里對不同曲率不同壁厚液芯微瓶進(jìn)行數(shù)值模擬，核芯的介電常數(shù)被液體取代(例如水的介電常數(shù)εreal=1.332=1.7698)。由以下公式可計(jì)算出模式在核芯中電磁場所占的比例:

(3)

(4)

式中，w(r)是能量密度，H(r),E(r),ε(r)和μ0分別是磁場強(qiáng)度、電場強(qiáng)度、介質(zhì)介電常數(shù)和真空磁導(dǎo)率[15]。

Fig.5 Relationship between percentage in core and thickness

圖5描述了在膨脹半徑R=25μm，Δk12=0.0035，Δk22分別為0.004和0.005時(shí)，TM基模和切向二階模下微瓶核芯光場所占百分比隨赤道面處壁厚變化規(guī)律?？梢钥闯?，核芯能量百分比隨赤道面處壁厚增大而減小，壁越厚則越多能量處于壁中。內(nèi)壁表面曲率Δk2越小，核芯百分比所占能量越大。對比兩種模式，顯然在壁厚小于一個(gè)工作波長時(shí)，Δk2發(fā)生改變時(shí)，切向二階模所受影響大于一階模。

3 微瓶傳感特性

回音壁模式微腔應(yīng)用于傳感領(lǐng)域，主要由于微腔內(nèi)部或者自身與外界相互作用引起諸如頻率或者光譜改變[16-18]。特別對于液芯傳感[1,19-20]，當(dāng)外界環(huán)境(核芯物質(zhì)濃度、溫度、壓力)發(fā)生變化時(shí)，微腔內(nèi)部諧振波長會(huì)產(chǎn)生一個(gè)較小漂移，通過掃描可精確測量出該譜線漂移量。表面曲率不同對微瓶傳感特性影響很大[4,9,21]，將折射率靈敏度定義為回音壁模式的頻移比:

S=?λ(nc)/?nc

(5)

式中,nc表示核芯物質(zhì)有效折射率，λ(nc)表示微腔內(nèi)部諧振波長。通過引入對微瓶狀微泡內(nèi)核芯物質(zhì)折射率一個(gè)微小改變量(如0.001)進(jìn)行模擬，計(jì)算出模式諧振波長偏移量。選用外壁表面曲率平方Δk12=0.0035，赤道面處壁厚為1μm，不同膨脹半徑、內(nèi)壁表面曲率進(jìn)行模擬計(jì)算，得到不同切向模式諧振波長移動(dòng)變化曲線,如圖6和圖7所示。

Fig.6 Relationship between wavelength and R

Fig.7 Relationship between wavelength and Δk22

圖6中描繪了在內(nèi)表面曲率平方Δk22分別為0.004,0.005,0.006、一階模和二階切向模式下液芯微瓶核芯折射率發(fā)生0.001改變量時(shí)諧振波長移動(dòng)量?？梢钥闯?各階模式下波長變化量均隨膨脹半徑增大而增大，且內(nèi)壁表面曲率越小，各模式下諧振波長變化量越大，越有利于傳感。對比兩種模式，相同曲率時(shí)，隨膨脹半徑變化切向二階模波長變化量小于一階模，且曲率越大時(shí)效果越明顯。選用更多內(nèi)壁曲率數(shù)據(jù)進(jìn)行模擬，如圖7所示。各階模在內(nèi)壁表面曲率越大時(shí)諧振波長偏移量越小，且切向二階模波長偏移量減少得更快，表明曲率對切向高階模的影響更為明顯，即曲率越小越有利于二階切向模式折射率傳感應(yīng)用。

4 結(jié) 論

基于有限元數(shù)值模擬方法對不同表面曲率微瓶狀微泡回音壁模式光學(xué)特性進(jìn)行了研究。探討了微瓶在不同表面曲率、膨脹半徑時(shí)模式特征，并討論了不同表面曲率下不同模式核芯能量百分比，最后對其折射率傳感特性進(jìn)行了分析。通過調(diào)節(jié)表面曲率可獲得高Q值、高靈敏度的折射率傳感器。雖然只討論了小尺寸微泡和一種工作波長(1.55μm)情況，但對大尺寸和其它波長同樣存在相似結(jié)論。該傳感器除了折射率傳感外,還可以用于其它物理參量如壓強(qiáng)和溫度等傳感領(lǐng)域[22-23]。