comsol三次谐波本征手性BIC,远场偏振图手性透射曲线二维能带图Q因子图电场图所见即所得。 大子刊nc复现在连续介质中束缚态驱动下具有最大和可调谐手征光学响应的平面手征超表面。最近在复现Nature子刊关于手性超表面的工作时发现利用连续介质中的束缚态BIC搞手性操控是真有意思。这次咱们直接甩干货聊聊怎么用COMSOL玩转三次谐波和本征手性BIC顺带看看那些酷炫的电场图和Q因子图到底怎么整出来的。先看二维能带图——这玩意儿是超表面设计的GPS导航。用COMSOL的波动光学模块跑参数化扫描记得在频域研究中把波矢k从0扫到π/aa是晶格常数。有个坑别踩边界条件里的Floquet周期端口设置必须用波矢kx,ky模式否则扫出来的能带直接跑偏。下面这段脚本能自动提取特征频率model mphload(band_structure.mph); freqs mphglobal(model, freq, dataset, dset1); plot(real(freqs), imag(freqs), o); % 实部虚部分别对应频率和损耗重点来了——本征手性BIC的实现。在结构对称性破缺的位置引入梯度扰动比如把纳米棒设计成梯度倾斜排列。这时候在电场图里能看到明显的旋涡状分布用后处理里的emw.Ex*emw.Ey交叉项能直接反映手性强度。注意看电场相位分布左右旋圆偏振光的相位差会呈现顺时针/逆时针螺旋特征。comsol三次谐波本征手性BIC,远场偏振图手性透射曲线二维能带图Q因子图电场图所见即所得。 大子刊nc复现在连续介质中束缚态驱动下具有最大和可调谐手征光学响应的平面手征超表面。远场偏振图是个直观的验证手段。用远场投影功能时记得勾选计算Stokes参数然后自己写脚本画偏振椭圆import pyvista as pv stokes model.result().export(Stokes_params).data() ellipse pv.PolarPlot(stokes[:,0], stokes[:,1]) # S1和S2参数决定椭圆形态手性透射曲线要玩双线交叉才够劲。在频扫时同时监控左右旋光的透射率用参数耦合把结构不对称度Δ设为扫描变量。当Q因子图出现尖峰时建议用对数坐标说明BIC被激活。这时候透射曲线会出现明显的手性劈裂交叉点位置随Δ值移动——这就是可调谐性的关键。最后说说三次谐波增强。在COMSOL中需要同时跑基频和三频两个物理场用非线性材料属性链接。重点观察电场增强因子的空间分布是否与基频模式重叠这时候的远场图会出现六重对称花瓣结构和线性区的圆偏振特征形成鲜明对比。搞这些图表的终极奥义是参数化扫描时把结构尺寸、入射角、材料非线性系数全设为全局参数然后直接在后处理里用LiveLink联动MATLAB生成动态图表。记住玩转超表面仿真就是要在物理直觉和数值把戏之间反复横跳毕竟——电场分布不会说谎但如何让它说出真相就得看操作了。