1. 先說清楚“off-gamma BIC”到底算什么物理1.1 BIC不是玄學連續(xù)譜里的“隱身態(tài)”我最早聽到Bound State in the Continuum第一反應是“這不就是一個永遠不輻射的共振嗎跟理想模式有什么區(qū)別”。后來做超表面高Q諧振設計才發(fā)現(xiàn)事情沒那么簡單。BIC指的是這樣一個態(tài)它的頻率落在輻射連續(xù)譜內理論上可以和所有平面波通道耦合但因為對稱性或參數(shù)配平的原因它對外輻射的分量剛好為零能量就老老實實待在結構里。換個說法它是一間漏風的房子里的某個房間外面再怎么刮風下雨這個房間里的能量就是出不去。在COMSOL里算BIC核心工作其實不是“找不輻射的模式”而是把輻射通道在數(shù)值上完整地呈現(xiàn)出來然后找到輻射阻尼消失的那個參數(shù)點。BIC在倒空間中通常是一個孤立點只有一個布洛赫波矢kx對應無限高Q。實際樣品總有尺寸、入射角展寬你不可能精確落在那個k點上所以近幾年大家都在做“Merging BIC”也就是把多個孤立BIC通過調節(jié)幾何參數(shù)合并到一起讓高Q區(qū)域從“一個尖峰”變成“一片平臺”。1.2 Γ點BIC和off-Γ BIC是兩種完全不同的存在Γ點kx0, ky0出現(xiàn)的BIC一般是對稱性保護的。結構里有C2或C4旋轉對稱性某些模式在這些對稱操作下的變換性質和所有輻射通道都正交對應模式就無法輻射。這種BIC只要結構對稱性不被破壞它就一直存在不太需要專門調參數(shù)。off-Γ的BIC就麻煩得多。它在kx不為零的某個位置出現(xiàn)沒有嚴格的對稱性保護純粹靠結構的幾何參數(shù)在某組特定數(shù)值下讓偶然輻射通道的系數(shù)碰巧為零。這類BIC也被稱為偶然BIC或參數(shù)型BIC它的位置、存在性和Q值對幾何參數(shù)極其敏感。橢圓柱的橢圓度、柱高、周期只要微調一點點off-Γ BIC的位置就會在k空間里移動甚至消失?！癕erging off-gamma BIC”就是利用這種敏感性把某個幾何參數(shù)掃過去讓兩個原本分開的off-Γ BIC在k空間里逐漸靠近最后合并。合并成功之后Q值在很大的k區(qū)間內都保持高數(shù)值。這個概念對做傳感、窄帶濾波、低閾值激光的人特別有用因為實際光斑尺寸有限角譜展寬不可避免合并后的BIC對入射角不再苛刻。1.3 Merging解決的真實痛點Q值對環(huán)境誤差的敏感我見過不少剛接觸這個方向的人跑COMSOL找到某個kx處的Q峰興奮得不行結果稍微改一下周期或者厚度那個峰就消失了。這不是算法的問題而是單一off-Γ BIC本質上就是“參數(shù)懸崖”上的平衡存在條件非常脆弱。Merging BIC的意義在于通過把BIC合并讓極端高Q區(qū)域從零維變成一個有一定寬度的區(qū)域。這樣樣品即便存在幾十納米的尺寸誤差或者激發(fā)光有一個小的角度偏差Q值仍然能保持較高水平。所以我建議所有打算在COMSOL里復現(xiàn)Merging BIC計算的人動手之前先把這三個概念分清對稱保護BIC在Γ點、偶然BIC在off-Γ、通過參數(shù)掃描讓多個偶然BIC匯合。章節(jié)后面講的所有建模和掃描策略都是圍繞這三個對象的辨別與追蹤展開的。2. COMSOL建模之前的選型決定了后面所有計算的成敗2.1 用RF模塊還是波動光學模塊很多人在這一步猶豫。我的習慣是直接用RF模塊的“電磁波頻域”接口做特征頻率研究。波動光學模塊也能做特征頻率但它的邊界條件和端口設置更偏向波導器件對周期性超表面來說RF模塊的Floquet周期性邊界條件用起來更順手材料折射率設置也直白。尤其在近紅外或可見光波段用RF模塊求解本征模問題時介質柱的折射率用實數(shù)先不考慮吸收損耗這樣算出來的Q值反映的是純輻射Q方便你判斷BIC是否真的“不輻射”了。如果你最終關心的是Q_total也就是實驗上直接測到的品質因子那要把材料吸收加回來的話就采用1/Q_total 1/Q_rad 1/Q_abs這種思路分兩步算。這個拆分極其重要后文踩坑部分會詳細說為什么不能一開始就開吸收損耗。2.2 幾何參數(shù)化為什么橢圓柱是主流選擇如果只是為了得到某個高Q模式隨便一個矩形柱或者圓柱都行。但你要跟蹤的是off-Γ BIC隨參數(shù)移動的過程幾何參數(shù)就必須選得均勻、單調、不引起模式身份突變。我建議用橢圓柱結構短半軸設為r1長半軸設為r2橢圓度定義為γ r2 / r1。γ從1變到0.6結構對稱性從圓對稱逐漸過渡到橢圓對稱。圓對稱時有些模式在Γ點簡并打破圓對稱后簡并消除模式沿k空間的色散會重新排列而off-Γ BIC在這種“由對稱操作導致模式重組”的過程中出現(xiàn)或消失。用橢圓度當掃描參數(shù)比直接改柱高或周期更平滑。拿TiO2介質超表面舉例常見初始量級是周期P700 nm柱高H300 nm長半軸r2200 nm短半軸r1150 nm襯底用熔石英折射率約1.45TiO2在近紅外折射率約2.4。這些數(shù)值只是起跑點真正的合并條件要靠掃描來確定但以這個量級開始掃描計算量比較友好物理畫面也比較清晰。2.3 Floquet邊界和倒空間掃描路徑的設置思路要在COMSOL里掃k向量常見做法是給周期性超表面單元設置“周期性邊界條件”并啟用Floquet子類型然后指定布洛赫波矢分量kx和ky。很多教程會讓你直接在“Floquet周期”里填具體的k_x和k_y數(shù)值。但在實際掃描時應該把kx定義成全局參數(shù)kx s * (2*pi / P)其中s是輔助參數(shù)從0掃到0.5。這樣做的好處是掃描的每一步都對應倒空間里的規(guī)則路徑不會因為手動填數(shù)出現(xiàn)k點分布不均的問題。掃描路徑我一般沿著Γ-X方向也就是固定ky0掃kx如果要找更精細的合并行為再沿Γ-M方向補掃。這方面的設置還有一個坑不要把“布洛赫波矢”和“入射角”混在一起。特征頻率研究里求解的是本征頻率ω同一個本征頻率對應自由空間波矢k0 ω/c而布洛赫波矢kx是倒空間里的獨立坐標。只要你在邊界條件里指定kx求解器就能算出對應的模式頻率不需要再通過k0 kx/sinθ去反推入射角否則會給自己制造不必要的非線性。3. 特征頻率研究里如何把Q因子“逼”出來3.1 從復本征頻率到Q值COMSOL特征頻率研究輸出的本征頻率是復數(shù)。物理上實部代表模式的諧振頻率虛部代表能量衰減速率。Q值按下面這個公式從本征頻率提取Q Re(f) / (2 * |Im(f)|)如果是無損耗介質Im(f)的來源就是輻射。BIC點對應的理想情況是Im(f) → 0但數(shù)值上永遠不可能得到無窮大的Q因為離散化網(wǎng)格、求解器誤差、邊界條件近似都會引入虛假泄漏。所以判斷BIC不能等Q變成無窮大而是看Q在某個k點是不是出現(xiàn)尖銳的極大值同時這個極大值對網(wǎng)格加密是否穩(wěn)定。我在實際操作時習慣用log10(Q)來畫圖而不是線性Q。原因很直接普通共振的Q可能在幾千到幾萬BIC附近的Q可能到10^5甚至10^7線性坐標會把低Q部分壓成一條直線根本看不出off-Γ峰的位置。取對數(shù)后Q峰的形狀和位置一目了然。3.2 用輔助參數(shù)掃描k空間的兩種做法在COMSOL里讓kx變化有兩種常用方式。第一種直接在研究設置的參數(shù)列表里把kx設為參數(shù)選擇“參數(shù)掃描”把所有想算的k點列出來。這種方式直觀、結果清晰但每個k點都要重新做一次特征頻率求解比較費時間。第二種用輔助掃描。在研究設置的“掃描”里加一個輔助參數(shù)s從0以固定步長掃到0.5同時設定kx s*(2π/P)。這兩種方式本質沒區(qū)別我推薦先粗掃。第一步用20個k點、較粗網(wǎng)格把Q曲線趨勢摸清楚確認哪些區(qū)域有可疑的Q峰再在峰附近加密取點。一次特征頻率掃描不要貪多寧可多跑幾輪。3.3 模式追蹤怎么確保你盯的就是同一個模式特征頻率問題最大的麻煩在于隨著kx變化本征值的排序會跳。同一個物理模式在不同k點可能排到第幾個完全不一樣尤其是頻率接近的模式密集區(qū)COMSOL列出的第5個特征值可能在下一個k點變成了第3個。我有兩個辦法應對。一是把頻率范圍縮窄在“特征頻率”設置里選擇“搜索特征頻率附近”給定一個預期頻率把無關的低頻模式和高頻偽模式擋在外面。二是每次算完一個k點就輸出電場分布確認模式剖面和你要跟蹤的模式一致。我一般用兩個判據(jù)諧振頻率的連續(xù)性以及模式在z方向的電場剖面或模式偏振特征是否保持一致。如果看到頻率出現(xiàn)不連續(xù)跳變或者同一支色散上突然冒出一個Q完全異常的數(shù)值大概率是模式跟蹤跳線了需要回到那個k點檢查場分布?？梢宰鲆粋€小建議在參數(shù)掃描時開啟COMSOL的模式跟蹤功能它會嘗試從上一個參數(shù)點的解來找當前點的對應模式。并不是萬能但在小步長掃描下很好用一旦步長太大模式跟蹤照樣會跟丟。4. Merging off-γ BIC的完整跟蹤流程4.1 第一步固定幾何參數(shù)摸清Q-kx曲線拿到模型之后不要急著掃橢圓度γ先把一個固定幾何下的Q曲線做出來。選擇γ1.0或接近1的值沿Γ-X方向掃描kx從0到π/P。此時輸出模式的Q值曲線重點關注兩個位置Γ點是否有對稱保護的尖峰以及kx≠0的非Γ區(qū)域是否有第二個尖峰。這一步能幫你確認幾件事Floquet邊界設置是否正確、目標模式在k空間里的連續(xù)色散長什么樣、哪些區(qū)域有off-Γ BIC的雛形。如果Γ點沒有出現(xiàn)預想中的高Q尖峰先別急著調物理參數(shù)檢查邊界條件和網(wǎng)格Γ點的對稱保護BIC在理想條件下應該非常明顯。4.2 第二步把橢圓度γ加進來觀察兩個峰之間的相對運動確定一個基礎幾何能出off-Γ峰后把γ變成掃描參數(shù)。一般來說γ偏離1越遠結構的面內對稱性調整越劇烈off-Γ BIC在倒空間的位置會移動Q峰值也會變化。掃描時可以先做γ從1.0到0.6的粗掃步長0.05每個γ下用20個k點跑趨勢。你會看到Q峰的位置隨γ連續(xù)移動有些γ下峰朝Γ點方向走有些γ下兩個峰之間靠攏。這個掃描階段的關鍵是同時記錄三個量kx、頻率、log10(Q)。不要只記Q峰值因為Merging BIC的判據(jù)是整個k空間里高Q區(qū)域的寬度不是單個孤峰。建議在二維掃描結束時直接畫出以γ為縱軸、kx為橫軸的Q熱力圖能很直觀地看到兩個高Q軌跡如何相交和匯合。4.3 第三步合并判定的客觀標準很多人以為看到兩個Q峰扎到一塊就算Merging BIC成立。我的經(jīng)驗是至少再加兩個判據(jù)。第一個判據(jù)是Q^{-1}的展寬。合并后在kx區(qū)域中間的一段范圍內Q^{-1}應該比合并前低一個量級以上。實際操作中你可以在某個特定γ參數(shù)下統(tǒng)計Q^{-1}低于某個閾值比如10^-5的kx區(qū)間寬度再對比單BIC情況下的寬度。如果合并后寬度明顯擴展說明高Q區(qū)域確實變大了。第二個判據(jù)是遠場或PB分量檢查。Merging BIC在物理上對應兩個輻射通道同時關閉因此模式下某個遠場方向的輻射應顯著下降。如果有條件在候選參數(shù)下加一層PML并計算模式在結構周圍的漏泄場看看漏泄波是不是真的被抑制了。COMSOL里加PML之后Im(f)會變成一個有限的較小值此時Q值才是可信的“開放結構下的數(shù)值Q”。4.4 計算資源不夠時的掃描策略完整二維掃描γ×kx非常耗時。我通常分兩輪第一輪用二維模型、粗網(wǎng)格、20個k點×8個γ值找出靶心區(qū)域第二輪只在靶心區(qū)域加密比如γ掃描步長收窄到0.01kx點加密到100個。如果機器內存吃緊可以把頻率求解范圍進一步收緊并且把COMSOL的并行線程設置成和核心數(shù)一致。計算能跑通是一回事跑通之后會不會判斷結果是另一回事省下的時間要用在驗證上。5. 那些會讓你誤判BIC結果的坑我都替你踩過5.1 網(wǎng)格疏密直接決定Q值“假峰”的位置BIC計算對網(wǎng)格極其敏感這點必須放在最前面說。我遇到過不止一次粗網(wǎng)格下在某個kx處掃出了一個Q很高的峰表面上看是off-Γ BIC加密網(wǎng)格之后峰位置居然移動了甚至消失。原因很簡單網(wǎng)格離散引入的人工泄漏在k空間分布并不均勻它可能在局部與真實模式的能量分布“耦合”出虛假的高Q尖峰。應對方法很笨但有效對每個候選峰位置至少做三套網(wǎng)格密度對比最大單元尺寸分別取λ/10、λ/20、λ/30。如果Q峰值隨著網(wǎng)格加密持續(xù)上升且峰位幾乎不動說明是物理的BIC如果峰位漂移或Q值先升后降那基本是數(shù)值假象。這一條是所有BIC計算的第一守則。5.2 絕對不要一開始就開材料吸收材料吸收和輻射損耗都會讓Im(f)非零但它們的物理來源完全不同。BIC的判斷依據(jù)是輻射損耗趨近于零。如果你在模型里用了復折射率那算出來的Q是總Q根本無法判斷該模式的輻射是否關閉。更麻煩的是在參數(shù)掃描里吸收項會掩蓋BIC帶來的Q尖峰讓你誤以為沒有BIC。推薦流程是無損耗折射率→找到輻射Q尖峰→加吸收項→對比Q_total。如果你想在無損耗條件下判斷合并是否發(fā)生就用無損耗折射率。這也是為什么我前面一直強調用實折射率起步的原因。5.3 周期邊界方向填反導致Γ點對稱保護BIC消失Floquet周期邊界條件里的布洛赫波矢分量順序千萬不能填錯。對于二維晶格第一個分量對應x方向的相移第二個分量對應y方向的相移。如果你調換了方向在kx0時應該受對稱保護的BIC會莫名其妙地出現(xiàn)有限輻射而且Γ點的Q值會低得不正常。驗證方法很簡單在kx0、ky0處算一個理想對稱結構理論上這個模式Q應該趨近于一個很大的數(shù)值。如果算出來只有幾千強烈提示邊界條件或網(wǎng)格有問題先解決這個問題再繼續(xù)掃描。另外周期單元的相對邊界之間相位要一致否則結構的平移會引入額外相位誤差看似細微實際會讓off-Γ BIC位置偏移好幾個k點。5.4 偽模式、局域模式和真實模式的辨別特征頻率一次能吐出一堆本征值里面混著很多沒有物理意義的偽模式常見的有兩種網(wǎng)格模式場分布集中在幾個網(wǎng)格節(jié)點附近波動劇烈和邊角模式能量集中在結構邊緣或周期邊界角落。這類模式通常Q值極端要么極低要么異常高頻率也不按色散規(guī)律變化。我的識別手段是掃一幅k空間場剖面圖看電場在單元內部是否呈現(xiàn)規(guī)整的駐波形式。如果是繼續(xù)跟蹤如果場分布里出現(xiàn)強烈的高頻波紋或邊界尖點直接放棄不要在這種模式上浪費時間。判斷模式“干凈”與否的經(jīng)驗是所有BIC計算里最容易被人忽略但又最關鍵的一環(huán)。6. 用PML和遠場結果驗證Merging BIC是否成立6.1 加PML對開放結構的必要性周期性結構在水平方向用Floquet邊界垂直方向是開放的能量可以從上下表面輻射進自由空間。如果上下邊界直接用完美磁導體或完美電導體去截斷輻射會被反射回來造成Q值虛高。這種虛高在BIC點附近尤其陰險因為BIC本來輻射就小反射回來會掩蓋真實的有限輻射。建議在結構上下加PML層厚度一般取一個波長量級材料用PML類型并與結構之間留一段空氣或襯底材料的間隔間隔不小于半個波長。加PML之后Im(f)會包含真實的輻射通道Q值才代表開放結構里的實際Q。做Merging BIC判斷時至少在最終候選參數(shù)下用PML模型驗證一次不要只在無PML模型里看峰。6.2 遠場方向圖和Q熱力圖怎么用來“下結論”Merging BIC驗證的最后一公里是看遠場輻射有沒有在某個輻射方向被壓掉。在COMSOL中如果研究里包含遠場域可以在后處理里加一個遠場輻射方向圖。合并前后的對比非常明顯合并前遠場圖在某幾個角度上有明顯的漏泄瓣合并后這些瓣變窄消失說明輻射通道徹底關閉。此外把所有掃描結果整理成Q熱力圖時強烈建議對Q值取對數(shù)并在色標上封頂。比如把超過10^6的Q一律顯示為同一顏色否則數(shù)值噪聲會造成彩色糊斑你根本看不清合并區(qū)域在哪。這一手處理也是寫論文出圖時經(jīng)常用到的對讀者、對合作者都更友好。6.3 從仿真走向實驗對照時要注意的現(xiàn)實差異如果你后面打算做實驗驗證COMSOL模型最好從一開始就包含襯底和上方覆蓋介質而不是懸空薄膜結構。加入襯底后模式的等效折射率分布會變化off-Γ BIC的合并參數(shù)也會跟著移動。實驗上最常用的驗證手段是角分辨反射光譜固定波長掃描入射角觀察窄帶共振在角譜上的展寬行為。Merging BIC在實驗上的標志是在一個較寬的角度范圍內都能看到高Q窄帶反射谷而不是只在某個特定角度出現(xiàn)。就算只做仿真不發(fā)文章我也建議你把襯底加進去。因為襯底的引入會讓某些在懸空結構里看起來完美的BIC直接消失提前發(fā)現(xiàn)這個問題遠比你投稿或被審稿人問到再補算要舒服得多。