comsol三次谐波，本征手性BIC,远场偏振图，手性透射曲线，二维能带图，Q因子图，电场图

comsol三次谐波本征手性BIC,远场偏振图手性透射曲线二维能带图Q因子图电场图所见即所得 大子刊nc复现在连续介质中束缚态驱动下具有最大和可调谐手征光学响应的平面手征超表面当平面超表面玩起手性魔术BIC驱动的光学极限探索comsol三次谐波本征手性BIC,远场偏振图手性透射曲线二维能带图Q因子图电场图所见即所得 大子刊nc复现在连续介质中束缚态驱动下具有最大和可调谐手征光学响应的平面手征超表面最近在玩COMSOL的时候偶然复现了大子刊NC上那个关于本征手性连续介质束缚态BIC的超表面设计。这玩意儿简单来说就是让二维材料在特定波长下产生“光学拧麻花”效应——左手和右旋光的透射率差异直接拉满甚至还能通过结构微调实现动态控制。下面直接上干货边跑数据边唠嗑。1. 电场图与BIC的“幽灵模式”先看电场分布图1用COMSOL的ewfd.Ez切个剖面。在BIC频率点附近电场会在纳米柱间隙形成漩涡状分布像被无形的手拧了一把。这种模式本质上是结构对称性破缺和辐射损耗抵消的结果。代码里有个关键操作model.study(std1).feature(freq).set(plist, 0.98*BIC_freq); // 故意偏移0.98倍频率 model.sol(sol1).runAll; // 暴力求解为什么要故意偏移频率因为纯BIC其实是数学上的奇点实际仿真得在邻近频点观察泄漏态。这时候电场图会暴露隐藏的对称性破缺——原本该完美束缚的场开始“漏气”但漏得很有规律为后续手性响应埋下伏笔。2. 二维能带图寻找那个完美的Q因子能带计算用了Floquet边界条件脚本里这个循环特别实用for kx in np.linspace(-0.5,0.5,50): # 扫第一布里渊区 model.param.set(kx, f{kx*pi/period}) model.study(std1).run() Q abs(real(freq))/2/imag(freq) # 手动算Q因子扫完k-path后能带图会在Gamma点kx0附近出现一个陡然上升的Q因子图2。这个尖峰对应的就是BIC模式——Q值理论上无限大实际仿真中我最高跑到1e6级别。不过要注意网格划分用extra-fine模式反而可能漏掉高Q模式适当粗化网格physics-controlled mesh有时更稳定。3. 手性透射曲线左右旋光的掰手腕大赛重点来了——手性透射差图3。这里需要同时算左旋和右旋圆偏光的透射% 圆偏光入射设置 s_param [1, 1i]; // 右旋 p_param [1, -1i]; // 左旋 excit model.physics(ewfd).feature(port1); excit.set(Sparam, {s_param, p_param});当频率扫过BIC点时左右旋透射率会出现断崖式差异。有意思的是在三次谐波频段约3倍BIC频率这个差异还会二次放大——因为非线性效应和模式耦合产生了协同作用。不过需要打开Nonlinear Optics模块并且要设置density of states修正才能准确捕捉。4. 远场偏振图光学莫尔条纹秀远场偏振分布图4用farfield数据集导出配合mphpolarplot画图。当结构处于临界对称性破缺状态时远场会呈现花瓣状分布且左右旋分量在特定角度上完全正交。有个技巧用mphinterp函数提取特定方向的Stokes参数theta 45; phi 0; S1 mphinterp(model, S1, coord, [theta, phi], dataset, dset2)这组数据可以直接用来反推结构的手性强度。当调整纳米柱的旋转角度比如从30度到45度远场图案会像万花筒一样旋转但主瓣方向始终锁定在特定偏振态——这种“所见即所得”的特性让实验党狂喜。结语调参侠的自我修养这个设计的核心在于利用BIC的极端场局域能力把原本微弱的手性效应放大到离谱的程度。但实操中最抓狂的是参数敏感度——纳米柱的旋转角度差1度Q因子可能掉一个量级。后来发现用L-BFGS优化算法配合Global Sensitivity分析模块能救命不过那就是另一个故事了...所有配图请脑补电场漩涡图/Q因子悬崖/透射率剪刀差/偏振花瓣图——毕竟咱的COMSOL截图还在渲染队列里排着队呢