在COMSOL里跑彎月型結構的BIC仿真幾何建模永遠不是最耗時的環(huán)節(jié)真正讓人頭疼的是斜入射條件下的線偏振設置。入射角度一變E0和k向量之間的垂直關系被破壞SBC散射邊界條件對平面波的吸收性能也跟著波動結果看起來花花綠綠其實很可能沒有一個數(shù)可靠。這篇文章就是我最近調(diào)試這個場景的完整記錄把散射邊界條件配合斜入射線偏振的關鍵設置思路交個底給正在做納米光子諧振器、超表面或者BIC方向仿真的朋友當個參考。1. 先想明白彎月結構里的BIC斜入射為什么牽一發(fā)動全身1.1 用馬路環(huán)島來理解連續(xù)譜束縛態(tài)BIC全稱是Bound State in the Continuum連續(xù)譜束縛態(tài)。普通光模式只有兩種歸宿要么被結構束縛住要么耦合到外部輻射通道里跑向遠方。BIC的特殊之處在于它的能量處在可輻射連續(xù)譜之中但模式卻始終不向外部輻射理論壽命無限長、Q值趨向無窮大。我習慣用一個生活場景類比一條開放高速路相當于連續(xù)輻射通道普通模式就是路上的車沿著大路駛向遠方而BIC是路中央一個完全封閉的環(huán)島里面的車只能在環(huán)島內(nèi)繞圈完全不會匯入主路對觀察者來說它既不放行車輛也不產(chǎn)生擾動。這類模式在傳感、激光和超表面領域都非常有價值因為稍微改變一點點環(huán)境條件Q值就會從無限大掉落到幾千幾萬對周圍參數(shù)異常敏感。BIC大體分兩類來源。一類是對稱保護型模式場分布的對稱性和外部所有輻射平面波的對稱性都不匹配輻射通道被幾何對稱性直接鎖死另一類是參數(shù)調(diào)諧型兩個泄漏通道的輻射分量在遠場恰好幅度相等、相位相反互相抵消之后一個有泄漏的模式被洗成了束縛態(tài)。理解這兩類來源是后面做參數(shù)掃描的基礎。1.2 彎月幾何到底帶來了什么額外自由度為什么不用圓盤、方柱這類規(guī)規(guī)矩矩的結構圓盤和方柱也能承載對稱保護BIC但這類BIC往往只在正入射條件下嚴格成立。一旦改成斜入射原有的鏡像對稱等對稱元素被破壞BIC會很快退化成一個低Q的Fano共振可供調(diào)節(jié)的參數(shù)維度很有限。彎月結構不一樣兩個表面的曲率半徑、中心厚度、底面半徑都可以分別控制結構天然帶有上下和內(nèi)外的不對稱性。模式分布不再依賴嚴格的全局對稱性而是被曲率比、高寬比等局部幾何參數(shù)牢牢控制。我在模擬項目X里做過一組參數(shù)掃描固定材料折射率只改變底部曲率半徑發(fā)現(xiàn)BIC出現(xiàn)的入射角從大約46度移動到了37度附近譜線形態(tài)也從對稱的Lorentz型過渡到了明顯的Fano型。這意味著彎月結構的BIC很適合把入射角當做一個外置旋鈕斜入射仿真不是可做可不做的推廣而是這個課題本身的核心研究手段。還有一點容易被忽略彎月表面是強彎曲界面外行波在曲面上的傳播方向并不單一結構散射出去的場可能會同時朝很多角度擴散。如果邊界條件只處理某個特定方向的平面波其他方向的反射就會悄悄溜進數(shù)值結果最后表現(xiàn)成一條虛假的泄漏通道把本應很高的Q值硬生生壓低。2. 邊界條件取舍斜面入射場景我為什么棄PML改SBC2.1 SBC到底吸收了什么以及它能吸收到什么程度SBC的數(shù)學本質是在邊界面上設置阻抗匹配關系讓外行電磁波穿過邊界時盡量不生成反射波。COMSOL里的SBC分為一階和二階一階形式默認邊界處的波近似垂直邊界傳播處理斜入射時反射率明顯升高二階形式在阻抗關系里加入了切向導數(shù)項對斜入射的吸收穩(wěn)定性大幅提升。BIC仿真對邊界條件的要求比普通諧振器高得多。高Q模式意味著光子在結構內(nèi)停留成千上萬個光學周期哪怕SBC每個周期只反射0.1%的能量累積起來也足以變成一個可觀的數(shù)值損耗通道Q值上限被人為壓低。所以彎月BIC這個場景里我直接上二階SBC并關閉后向反射選項——那個功能主要用于保留駐波效應對純吸收場景不僅沒有幫助反而會引入額外的數(shù)值反射路徑。2.2 為什么我不首選PML提到吸收邊界不少人的第一反應就是PML。PML正入射時效果確實很好但它本質上是針對特定傳播方向設計的高損耗層入射角一大吸收性能會顯著滑落而且需要反復調(diào)節(jié)層數(shù)與縮放系數(shù)。彎月結構的散射場是發(fā)散且多方向的單一取向的PML很難兼顧各個角度的外行波更不用說PML還要額外占用空氣域和網(wǎng)格資源。我做過實際對比入射角在15度以下時認真配置的SBC和精心調(diào)優(yōu)的PML共振峰位置的差別不到0.01 nm但SBC模型的體量和調(diào)參周期都更短。當然不要因此把SBC當成萬能方案入射角超過40度或者散射場里有明顯的掠射分量時SBC仍然會力不從心。我個人的經(jīng)驗判斷大致如下表場景二階SBCPML入射角0-10度推薦設置最簡可用但需要額外吸收層幾何入射角10-30度推薦能穩(wěn)定工作需要細致調(diào)節(jié)PML參數(shù)入射角30-45度謹慎使用必須空盒驗證更穩(wěn)妥可靠大角度/多方向散射場不推薦通常更合適3. 斜入射線偏振的三個地基E0、k向量、相位參考點3.1 先讓入射波滿足橫波條件在SBC設置面板里入射電場參數(shù)分兩部分E0矢量決定線偏振方向和幅度k向量決定傳播方向。斜入射時最常見的錯誤是隨手填一個E0比如E0(1,1,0)覺得只要電場同時有x和y分量就行。正入射時可能僥幸蒙對斜入射一出現(xiàn)就露餡平面波要求電場嚴格垂直于傳播方向只要E0和k的內(nèi)積不為零入射場就多了一個非物理的縱向分量結構附近會產(chǎn)生虛假激勵。標準做法是先把k確定下來再用叉乘構造E0。如果入射面取xz平面入射角為theta則有k k0 * (sin(theta), 0, cos(theta)) s偏振TE電場垂直入射面 E0s (0, 1, 0) p偏振TM電場在入射面內(nèi) E0p (cos(theta), 0, -sin(theta))E0p更嚴格的寫法是直接取 y_hat 與 k 的叉乘再做歸一化這樣無論入射角的正負如何變化偏振方向都不會出錯。3.2 把入射參數(shù)做成參數(shù)表而不是寫死成數(shù)字斜入射BIC必然要做入射角掃描每次手動改設置面板里的數(shù)值既低效又容易出錯。我習慣在全局定義-參數(shù)里一次性建好一批命名參數(shù)theta_inc入射角phi_inc入射方位角kx k0*sin(theta_inc)kz k0*cos(theta_inc)E_sx、E_sy、E_szs偏振的電場分量E_px、E_py、E_pzp偏振的電場分量SBC面板里直接引用這些參數(shù)名。后續(xù)只需要把入射角定義成參數(shù)化掃描COMSOL就會自動在每個角度上重新計算SBC的入射場不用任何額外的腳本干預。這條工作流節(jié)省的時間比想象中多得多。3.3 相位參考點為什么關鍵SBC的入射場表達式通常會帶一個相位因子 exp(i*k·(r-r0))r0就是相位參考點。多數(shù)教程對此一筆帶過但斜入射時它非常關鍵。如果參考點設在結構中心入射角變化之后每個網(wǎng)格點相對參考點的相位會整體平移雖然在理論上不改變最終物理結果但在強彎曲邊界上容易引發(fā)相位跳躍影響SBC在曲面處的吸收方向。我的建議是參考點固定在全局坐標原點每次更換入射面時同步檢查k向量的符號和E0方向是否仍構成右手關系。最可靠的自檢方法是空盒子測試把彎月結構臨時隱藏只保留均勻介質和空氣盒子施加同樣的SBC平面波跑一次頻域然后后處理看體積內(nèi)總電場減去解析平面波的差值。如果差值趨近于零說明入射場和SBC的相位配置完全正確才可以放心繼續(xù)。4. 從幾何到掃描彎月BIC仿真的完整實操鏈路4.1 彎月幾何、空氣盒子和網(wǎng)格的思路彎月結構可以用球冠差集來生成。一個大曲率半徑的球冠向下延伸底部與一個小曲率半徑的球冠或圓柱面求差就能得到彎月輪廓;想獲得更精確的形狀控制可以直接用超橢圓表達式做旋轉體。建模參數(shù)至少四個頂部曲率半徑、底部曲率半徑、中心厚度、底面半徑。我建議把參數(shù)先定義好再建幾何這樣后續(xù)做幾何掃描時只需改參數(shù)、不需要重新建模?？諝夂凶臃矫嫱鈬呴L至少取工作波長的1.5到2倍。盒子太小即使SBC本身是好的結構輻射的倏逝場在邊界處還沒有衰減到足夠低截斷誤差直接污染諧振頻率。彎月表面曲率變化劇烈曲面附近網(wǎng)格要加密到波長/30左右遠處空氣域按波長/10劃分使用自由四面體過渡即可。4.2 物理場和SBC配置的標準步驟物理場選擇電磁波、頻域研究類型選頻域。材料上調(diào)試階段建議讓結構介質不帶損耗因為哪怕虛部折射率只有一點也會給BIC引入真實衰減Q值上限瞬間被鎖死后邊很難區(qū)分物理效應和數(shù)值效應。邊界配置按這個順序執(zhí)行選中所有外部邊界添加散射邊界條件階數(shù)選2類型選平面波激活平面波入射場的設置E0和k替換為全局參數(shù)確認相位參考點是坐標原點關閉后向反射如果你的COMSOL版本有曲面曲率修正一類的選項建議一并打開二階SBC在強曲面邊界上會得到進一步改善。4.3 掃描策略先波長后角度不要同時在一張表里掃一堆參數(shù)那個組合爆炸很浪費算力。先固定入射角可以是0度或預期的BIC角度附近掃描波長用先粗后細的步長把共振位置鎖定然后固定共振波長掃描入射角重點觀察Q值隨角度怎么演化。兩階段都用參數(shù)化掃描完成每段掃描讓共振線寬內(nèi)至少分布5個采樣點高Q峰如果窄得離譜就改用自適應步長。求解器我直接用直接求解器比如PARDISO或MUMPS。BIC仿真別碰迭代求解器高Q模式讓矩陣條件數(shù)差到離譜迭代器經(jīng)常長時間徘徊不收斂省下的內(nèi)存遠遠彌補不了浪費的時間。網(wǎng)格規(guī)模上說一個2λ量級的空氣盒子配合局部加密3D全模型自由度通常在幾十萬到一兩百萬普通工作站完全跑得動。5. 怎么確認你看到的是真BIC而不是數(shù)值假象5.1 三個相互獨立的證據(jù)鏈看到窄峰不等于看到了BIC我通常要求三個證據(jù)同時成立第一譜線寬度要足夠窄并且隨著逼近BIC條件線寬持續(xù)收窄Q λ/Δλ朝一個穩(wěn)定的大值收斂。第二遠場輻射功率在BIC條件附近急劇下降普通共振的輻射功率比BIC鄰域高好幾個數(shù)量級。第三近場分布呈現(xiàn)束縛態(tài)特征場極大值局域在彎月核心區(qū)域離開結構后迅速衰減而不是一個正在向外拖拽的輻射波前。5.2 網(wǎng)格與邊界對Q值的污染測試高Q仿真里數(shù)值誤差最喜歡偽裝成物理信號。網(wǎng)格太粗模式場被削掉一部分Q值可能虛高也可能虛低SBC反射會把Q值上限壓低。要分辨這兩類誤差用雙測試策略網(wǎng)格測試把彎月曲面附近的網(wǎng)格從波長/30加密到波長/50重新跑同一點看共振頻率移動量。真BIC對網(wǎng)格應當非常穩(wěn)定頻率移動小于千分之一如果明顯移動繼續(xù)加密。盒子測試把空氣盒子邊長從1.5λ擴大到3λ對比Q值和共振頻率。Q值隨著盒子變大而明顯增長說明結果被邊界卡住基本不變說明SBC吸收足夠干凈。兩個測試一起做排查偽峰的效率非常高。5.3 求解器收斂度的診察最后還要看求解器本身。直接求解器跑完后我會查看殘差和收斂報告。高Q頻點附近矩陣接近奇異容差設得不夠小的話峰形會被削平或者出現(xiàn)抖動。把相對容差壓到1e-8以下再檢查頻率掃描曲線在高Q峰附近是否光滑。如果出現(xiàn)非物理的鋸齒優(yōu)先查求解器設置不要急著去改物理參數(shù)。6. 實戰(zhàn)踩坑記錄這些細節(jié)最容易讓斜入射仿真報廢6.1 E0沒有歸一化角度掃描整體失真最早一版模型里我的E0直接填了1完全沒有考慮偏振方向歸一化。結果入射角變大時投影到傳播方向的電場有效幅度在悄悄變化透射率曲線像整體疊加了一個余弦因子和物理規(guī)律明顯不符。把所有E0向量歸一化到單位長度并確認入射總功率密度在任意角度下都一致之后問題徹底消失。6.2 s偏振和p偏振的定義一定要擺到臺面上斜入射時不同的人說水平偏振可能指的是完全不同的方向。s偏振的定義是電場垂直于入射面p偏振是電場平行于入射面。彎月結構在xy平面具有軸對稱性時s和p會激勵到不同的模式分支如果把兩者搞混譜線特征會張冠李戴。我現(xiàn)在會在模型里明確寫下incident_plane_normal、s_polarization、p_polarization三個參數(shù)后處理再畫電場矢量確認方向確保每一步都看得見。6.3 從PML模型改SBC時留下的影子吸收層有次我在一個舊PML模型基礎上直接改成SBC波形異常、Q值低了將近一個數(shù)量級。排查了很久才發(fā)現(xiàn)PML子域還留在幾何里只是邊界條件換了那層高損耗材料仍在默默吸收場。如果你也打算從PML模型改造務必把原來的PML域重新改回普通空氣域或者干脆在新模型里重新搭幾何別讓這種低級殘留浪費半天時間。6.4 空盒子自檢幾十秒的高性價比照妖鏡所有經(jīng)驗里最值得保留的一條每次換了入射角、換了偏振類型、換了邊界條件之后先別急著跑含結構的完整模型先用只有均勻介質的空盒子跑基準測試?？蘸凶幽Ｐ妥杂啥鹊蛶资刖湍芘芡?。如果空盒子里SBC都吸不干凈、或者平面波定義就有問題含結構的結果必然全錯空盒子過了再去跑正事能減少一大半無效仿真。最后再分享一個工作流習慣。每次調(diào)整完幾何或者SBC參數(shù)我會把空盒子基準、網(wǎng)格測試和盒子尺寸測試的結果存進一張對照表隨模型一起歸檔。后續(xù)再改任何參數(shù)只要對照表里對應的基準數(shù)值沒有漂移我就敢繼續(xù)往下跑。BIC仿真的坑大多埋在那些看不見的邊界里而基準測試就是專門照出這些坑的手電筒。