企业官网建设流程全解析

1. 仿真目标与整体思路：一套模型串起线性手性和非线性增强

1.1 这套仿真到底在算什么

先给还没完全进入状态的同学交个底：标题里这些东西，其实是在讨论同一个物理问题——一个介质超表面，在什么几何条件下能同时具备两个性质：一个是线性光学层面的“本征手性BIC”，另一个是非线性光学层面的“三次谐波增强”。

如果你最近在调研手性超表面的文献，大概率会看到“本征手性连续谱束缚态”这个说法。这里的“本征”两个字，是指结构本身在二维平面上就是手性的，不需要斜入射，不需要外加磁场，更不需要把结构做成三维螺旋，只要圆偏振光垂直打上去，透射谱里就会自然出现左旋圆偏光和右旋圆偏光的差异，也就是CD信号。而BIC呢，通俗讲是一种被“锁”在结构里的模式：理论上它和外部辐射通道完全解耦，所以没有辐射损耗，Q因子可以趋近无穷，但由于它不往外漏光，你直接测远场是看不见它的。这时候“手性”的作用就体现出来了——它通过破坏结构的面内镜像对称性，把部分BIC从“完全暗态”变成“准BIC”，让能量既能被高度局域增强，又能在远场被激发和读取。

这一个逻辑链条，直接决定了整套仿真的思路。我见过不少刚入门的同学一上来就把全部目标堆到一起，想一次性得到所有图，结果参数乱成一锅粥。我做这套模型的经验是：先把它拆成线性本征问题、线性透射问题、非线性谐波问题三个层次，每个层次对应特定的计算方法和后处理手段。本文后面所有章节，基本就是沿着这条主线展开的。

1.2 为什么把三次谐波和手性BIC放在同一个体系里

三次谐波THG是非线性光学里最经典的“探针”效应之一。它的物理来源是三阶非线性极化，辐射强度理论上正比于局域基频电场强度的六次方——注意，是六次方。这意味着如果你能把结构里的局域场增强10倍，THG信号理论上能增强10^6倍量级，这是一个非常恐怖的非线性增强杠杆。

普通超表面的局域场增强，靠的是等离激元或者米氏共振，增强倍数往往只有几倍到几十倍。但BIC不一样，特别是近BIC模式下，理论Q因子可达10^5甚至更高，对应的局域场增强可以达到两个数量级以上。所以一个自然的想法是：把“高Q的准BIC模式”当作非线性增强的平台，而“本征手性”则给这个平台增加了偏振选择性——也就是说，这种增强只对某一旋向的圆偏光有效。

从应用角度看，这个组合能解决一个实际问题：手性非线性器件通常要求结构同时具备“强手性响应”和“强非线性转换效率”。单独做手性超表面，共振峰处的CD信号能做到很高，但普通共振Q因子低，非线性效率上不去；单独做高Q非线性超表面，结构往往是对称的，对手性没有响应。本征手性BIC走的是另一条路——在手性单元结构里构造BIC，等于是在同一套几何里同时调用对称性破坏和高Q束缚态两套物理机制。

对仿真来说，这意味着你必须同时获得五类数据：二维能带图用来确认BIC的存在与位置，Q因子图用来量化模式束缚强度，手性透射曲线用来验证线性的偏振选择性，电场的近场分布图用来定位增强热点和分析模式对称性，远场偏振图则用来确认辐射光的偏振状态并指导后续实验测量。

1.3 一个容易理解的结构设计出发点：既有“手性旋钮”，又有“泄漏阀门”

在我参考过的多篇介质超表面文献里，构造本征手性BIC最常用的手段，是取一个原本具有镜像对称性的高Q超表面单元，然后在几何上引入一个小扰动，比如把方柱的角削掉一块、在圆盘上挖一个斜向椭圆孔、把矩形柱改为Z形或S形等。这个“削角”或者“开孔”的动作干了两件事：它破坏了面内镜像对称，使模式对圆偏振产生手性选择性；同时让原本和外部连续谱完全解耦的BIC产生了一个微弱的辐射通道，把Q因子从无穷拉回到有限但依然很高的值。

这个逻辑可以用一个很直白的类比来说明：把BIC想象成一栋只有一扇门的密室，BIC状态就是门反锁、里面的人死活出不来。你如果只想看到这个人，永远看不到，因为他不透气。现在你在墙上开了一条很小的缝，这门缝就是“手性扰动”——人终于能往外走一点，Q因子仍然很高，因为他大部分时间还是被关着的，但你已经有办法感知到他的存在和活动状态了。门缝开得越大，信号越强，但Q下降得也越快，共振增强打折扣。所以几乎所有相关研究都会做一个参数扫描，比如从0度到30度去扫这个结构几何参数，观察Q因子、CD信号和THG增强之间的权衡。

因此，仿真时我给读者的第一个建议就是：不要急着复现某个具体论文的全部结果，先把你的基础结构建好，确保能在一或两个参数变化下系统性地从BIC过渡到准BIC，再做后续优化。

2. 几何建模、材料选型与COMSOL参数化设置

2.1 晶胞参数和材料参数怎么起步

如果你手里没有现成的器件尺寸，从一个相对通用的近红外介质超表面配置起步是比较稳妥的。我用过比较多的是非晶硅或者晶体硅纳米柱，波长范围放在1300到1700纳米。这个区间最方便的地方在于，硅的折射率约3.4到3.5，损耗比较低，单元的尺寸大概在亚波长量级，网格量和计算时间都比较可控。

举个例子，假设周期P=900纳米，柱高H=600纳米，基础的方形柱或者矩形柱边长约450纳米，当你需要它进入本征手性区时，最直接的做法是把方形柱的一组对角切掉，或者在一个圆盘上引入一个“S形”的镂空。这些几何操作都可以在COMSOL的几何序列里用布尔运算一次性完成。

材料设置上，记得硅的折射率要按波长变化来设定。如果直接填固定折射率，在宽谱透射扫描中会丢掉很多色散细节，尤其是BIC共振本身带有很强的模式色散，折射率不对会导致共振位置和Q因子都出现偏移。折射率数据可以从常见的光学常数库中提取，再在COMSOL材料节点里用插值函数引入。

2.2 在COMSOL里建立手性单元并定义扫参

进入COMSOL波动光学模块之前，先想清楚你要用的是二维还是三维模型。说句实话，二维能带图、手性透射曲线和Q因子图，很多情况下可以在合理的二维化模型里快速出结果；但如果你的结构是介质柱阵列，电场在高度方向上有显著分布，那我建议直接上三维，因为二维模型会丢掉很多垂直方向的物理细节。我实测下来，三维模型配合周期性边界条件，在网格数量控制得当的情况下，计算负担并不像想象中那么可怕。

三维建模的关键步骤是几何操作和参数化：

1. 建立一个长方体单元，长宽均为P，高度为H。
2. 在中心添加一个需要做成手性形状的介质柱体。最简单的做法是建立一个长方体，然后在两个对角处添加两个小圆柱或三角块，用差集把它们“切掉”。
3. 把用于切角的块的尺寸或者旋转角设为全局参数，比如d_corner。这个参数最好独立于柱子的整体尺寸，方便后面扫参时只改变“不对称度”。
4. 在单元上下面添加空气层，厚度至少为工作波长的1.5倍以上，给倏逝场留有足够衰减空间；同时设置周期性端口。
5. 模型的四周，Z方向用周期性边界条件，顶面和底面用端口边界条件。

第一次跑的时候可以把周期P和柱高H设为扫描参数，但我经验是优先扫“不对称度”这一个参数。因为它才是控制扰动强度、决定Q因子和手性大小的核心旋钮。等你能观察到随扰动参数增大、Q因子下降且CD增强的趋势后，再去优化周期和柱高放大效果，效率会高很多。

2.3 网格和求解器设置上的几个关键点

BIC模式之所以难算，是因为它的辐射损耗非常小，透射谱上的共振峰极窄，需要很高的频率分辨率才能解析。COMSOL的频域扫描如果只是均匀步长扫描，经常会漏掉共振点；更常见的问题是网格太粗时，有限元数值泄露会人为增加辐射损耗，让BIC模式根本看不出高Q特征。

我的建议是分两次求解。第一次用稍微均匀的网格和中等频率步长，比如2纳米间隔，扫出透射谱上的大概共振位置；确定共振位置后，再在共振波长附近做局部细化扫描，步长缩到0.1到0.01纳米量级，同时把网格在介质柱内部和上下表面的边界层加密。注意在BIC模式里，电场能量主要集中在介质内部和表面附近，空气区网格密度过高意义不大，反而会拖慢速度。

这里有一个特别容易被忽略的点——边界条件的取向。手性结构对左右旋圆偏振光的响应不同，因此你的端口中必须能独立区分LCP和RCP。在COMSOL中使用周期性端口时，可以设置两个正交线偏振的本征模式，然后通过设置相位关系来合成圆偏振入射。如果你不希望手工操作端口模式，也可以直接用背景场替代——在散射场公式下，把入射平面波设为Ex分量相位0度、Ey分量相位正负90度，就能得到两种不同旋向。第一次算的时候最好把两种入射的几何都跑一遍，确保CD信号的符号和手性结构的方向自洽。

3. 二维能带、Q因子提取与BIC模式判定

3.1 能带图不要只扫Γ点

很多人在验证BIC的时候，习惯性地只在布里渊区中心点即Γ点附近做频率扫描，看到某个模式Q因子很高就下结论说这是BIC。这其实不够严谨，尤其当结构引入手性扰动之后，原本在Γ点的对称保护BIC可能会移动到布里渊区内部成为驻波状态，而不是仍然严格在Γ点。

更稳妥的做法是做一个“二维能带图”，也就是沿布里渊区边界走一圈路径，通常对正方晶格取Γ->X->M->Γ，甚至完整遍历第一布里渊区内的波矢。在这个路径上求解特征频率，就能看出模式是如何随波矢变化的。一个真正的BIC，在能带图上的表现是：在特定的波矢点附近出现一条和其他能带没有交叉的平坦带或者高Q模式带，且该点在频率上不与任何辐射连续谱重合。

在COMSOL里做能带图，方法是周期性边界条件里定义波矢分量kx和ky，用参数化扫描改变它们，对每个波矢做一次特征频率求解，然后把得到的本征频率和波矢位置画出来。生成能带图的脚本思路其实很简单，难的是确保每条带上的模式都是物理真实的模式。我的排查方式是在每个波矢点保存一个电场分布快照，检查模式是否为周期性边界下的布洛赫模式，而非局部边界伪影。

3.2 Q因子怎么从仿真里提取

BIC本身Q是无穷大的，但你实际能计算到的只是准BIC。Q因子的提取有两种常用途径。

第一种，透射谱拟合。在频域扫描得到透射曲线后，把共振峰用Fano线型或者Lorentz线型做拟合，从线宽估Q = λ₀ / Δλ。这个方法直观，但对网格要求高，窄峰下如果频率扫描步长不够小，很容易低估Q。

第二种，本征值法，也是我更推荐的做法。在COMSOL的特征频率求解器中设置包含复数本征值，也就是考虑阻尼情况，每得到一个模式的复数频率ω = ω_r + iω_i，就可以得到Q = ω_r / (2|ω_i|)。这个数值在有限的网格精度下，会有一个人为上限，一般是10^5到10^7量级。如果你把扰动参数调小，Q计算结果不再上升，甚至出现振荡，那就说明网格已经到极限了，需要加密或者换用更高阶的形函数。

我一开始做这类仿真时，习惯直接看透射谱的峰宽，后来写论文才发现审稿人更认本征值提取的Q曲线。两种方法同步使用很有必要，因为它们能从不同角度校验模式是否真实。

3.3 怎么判断你找的模式真的是手性BIC

判断方法可以从两个维度去做。第一个维度是场分布：真正的BIC模式通常具有独特的近场对称性，其远场辐射分量为零，你在近场电场图中会看到一个高局域的涡旋状或节点状结构。当加入手性扰动后，准BIC的电场图中会出现不对称化趋势，某些热点向切角方向偏移。第二个维度是Q因子随扰动参数的幂律变化：理论上一维扰动导致的准BIC，其Q和扰动参数的平方成反比，即Q ∝ 1/δ²。你可以把不同切角尺寸对应的Q画在双对数坐标里，看斜率是否有趋近2的线性区间。

有了这两个判断，你才能比较有把握地说：“我看到的这个高Q模式，确实来源于BIC泄漏后的准BIC，而不是普通的米氏共振。”

4. 手性透射曲线与远场偏振图

4.1 圆偏振入射设置的前后一致性问题

手性透射曲线本身不难算，但初次接触的人最大的困惑是圆偏振的定义。我在仿真代码里习惯把迎着光传播方向看的电场旋转方向作为判据。假如光沿Z方向传播，LCP的定义是电场的x分量超前y分量90度，数学上写作E = E0(x̂ + iŷ)。RCP则对应E = E0(x̂ - iŷ)。

但是要特别提到的是，不同文献的坐标定义可能反过来，因此你仿真得到的CD正负号并不是物理的第一要素，关键是你在透射曲线、远场偏振和后续THG曲线里要保持同一套定义。否则你可能会看到一套矛盾的结果——比如透射谱显示右手圆偏透过率高，而远场偏振图却显示共振辐射是左旋圆偏，这种冲突九成不是物理问题，而是你的坐标定义在某一环节翻转了。

实际计算时，一次性在模型中加两种旋向的入射背景场可能会占用更多内存。我通常分两次扫描，一次LCP、一次RCP，用脚本把S参数导出后，在绘图阶段再合并成手性透射图。

4.2 透射T_LCP和T_RCP能读出什么

在手性透射曲线图里，完整的物理信息不仅包括同偏振透射T_LL和T_RR，还包括交叉偏振透射T_LR和T_RL。对纯本征手性超表面且不考虑非对称反射的情况下，交叉偏振项和同偏振项的差异可以用于计算椭偏参数和圆二色性。仿真中，如果只设置了单一端口模式，一些软件会自动忽略交叉偏振分量，这是不够的。我都是用两个正交模式和它们的相位叠加来定义端口矢量，这样S参数矩阵可以包含全部四个透射分量。

一个典型的本征手性BIC透射谱长什么样？在共振波长附近，你会看到某一旋向的透射谱中出现一个非常窄的类Fano凹陷或者尖峰，而另一旋向的透射谱几乎不变。CD曲线在共振处出现一个符号明确、峰值显著的信号。这里有个容易被误导的地方：很多新手会把CD峰位直接等同于BIC共振峰位，其实由于Fano干涉背景的存在，CD峰的极值位置通常会和透射凹陷位置有微小偏移。不要慌，这是离散态和连续背景干涉的正常结果，不是计算错误。

4.3 远场偏振图到底怎么画

从仿真的角度看，远场偏振图通常有两种画法。一种是直接在COMSOL的远场节点里，计算某个频率下超表面整个周期阵列在远场球面上的电场角分布。另一种是手动推导：取近场平面的数据，做平面波展开，然后对每个衍射级次或者每个k方向的辐射分量提取偏振态。

如果你使用COMSOL远场节点，最需要注意的问题是：远场计算需要一个封闭的积分表面，而且电流源J和磁流源M的提取需要保证在结构外部没有其他源。对周期结构来说，我习惯在一个完整周期单元周围建立一个立方体积分表面，把实际的切向场都提取出来，再在该表面上计算远场贡献。

提取偏振图的具体操作，可以把远场球面上的电场写成球坐标分量E_theta和E_phi，然后投影到传播方向的垂直平面，分解成左旋和右旋分量的振幅。可以用归一化的Stokes参数S1/S0和S2/S0来定量描述辐射场的偏振态。在准BIC共振频率处，你会观察到辐射场以某一旋向为主，偏振椭圆度接近1，而偏离共振时背景辐射主要是线偏振成分，椭圆度降下来。这个现象非常直观，也算是本征手性BIC在远场留下的“指纹”。

5. 三次谐波增强的原理与建模仿真流程

5.1 电场图为什么是三步走里的第一步

计算THG之前，我强烈建议你先把基频共振处的电场图认认真真看几遍。THG信号正比于局域电场的六次方，这个关系决定了并不是结构所有区域都对远场THG有贡献，只有那些场强极高的“热点”才是真正贡献非线性信号的区域。

画电场图时有几个实用建议：对三维模型，至少画两个截面，一个是z= H/2处的水平截面，观察面内模式分布和手性扰动导致的电场不对称；另一个是沿着x方向的垂直截面，观察模式在垂直方向上的局域化程度。颜色标尺用对数坐标，因为BIC模式下的场增强可能达到几十到几百倍，线性坐标下低场区域会全部变成蓝色，看不出热点位置。

BIC的电场图有明显特点：模式能量高度局域，电场强度在柱体内部远大于外部空气区。当你加入手性扰动后，场分布中会出现明显的旋转特征。我常见的现象是，在某些相角下，电场的旋向会切角附近形成“涡旋结构”。这些涡旋往往就是非线性增强的源点。如果这种区域恰好和材料的高非线性系数区域重叠，THG效率就会非常好。

5.2 三次谐波仿真的“两趟计算”

严格全波模拟THG，可以理解为一次三阶非线性极化过程：基频电场E(ω)在介质中产生非线性极化P_NL(3ω) = ε₀χ(3)E(ω)³，这个极化在三倍频处充当一个等效电流源，向外辐射三倍频光。

工程上我推荐把计算分成两步走。第一步是纯线性的基频场求解：在COMSOL的频域求解器中，将一个周期性超表面结构放在端口激励下，频点在基频波长，获得整个结构的电场分布Eω(r)。第二步是把非线性极化当作一个已知的“等效源”，在三倍频频率处再做一次线性电磁场求解——此时介质的折射率要用三倍频波长处的实际折射率，并且把源项设置为从第一步电场推导出的P_NL。

直接在COMSOL中一次性同时求解基频和三倍频是一套相当复杂的多物理场耦合，多数时候收敛困难和内存开销都会让你怀疑人生。分两步做，每一步物理目标单一，排查问题也容易。如果你只是想要一个相对THG效率估算，而且是非磁性的非手性材料体系，也可以走捷径：把每个像素的局部增强因子|E_ω(r)/E₀|取六次方再在体积内积分，用这个量来表征THG的远场变化趋势。很多文献里的THG增强倍数图，其实都源于这个简化模型，在缺少高精度实验对比时完全够用。

5.3 手性THG曲线怎么从两次仿真中提取

当你做完了基频的线性仿真后，会很自然地拿到两个量：LCP入射下的电场分布E_L(r)和RCP入射下的电场分布E_R(r)。用这两组电场去计算THG源的强度对比，可以得到一个估算：手性非线性响应强度差。

如果需要得到“手性THG转换效率曲线”这种更完整的结果，就需要在上述两步法的基础上做一个波长扫描：在基频波长附近选一组合适的波长点，比如从共振中心向两边各取3到5个点，每个波长都做一次线性求解和一次三倍频源辐射求解，最后把三倍频频率处的透射功率画成一条曲线。这组曲线做两个旋向各跑一遍，就能拿到漂亮的LCP-THG和RCP-THG对比图。

这里提醒一个取点问题：因为Q因子很高时共振极窄，你选的波长扫描点在远离共振时THG效率会快速低到数值噪声水平，这会给结果带来很大的误差。此时需要注意增加频率点密度，或者在计算THG源时使用相对场增强，而不是直接用绝对功率值。另外一个坑是三次谐波频率处材料的折射率。很多材料在近红外基频处损耗很低，但在三倍频处可能落入吸收带，此时必须使用带虚部的复折射率，否则你会得到偏高的THG功率估算。

实际操作中，我还见过有人直接在基频求解时不加三阶项，然后在后处理里用“E^3”的公式去画非线性源分布，这本身没问题；但如果你的结构里电场增强幅度太大，局部功率密度可能已经超过了非微扰阈值，那理论预估本身就值得怀疑了。

6. 常见问题与排查技巧

6.1 这些问题我基本都踩过

下面我挑几个反复出现的坑位总结一下，按自检优先级整理成表，方便你遇到问题时按图索骥。

6.2 三个提升仿真效率的小技巧

第一，手性BIC结构本身并不需要非常复杂的网格。共振模式高度局域在介质内部，这时你要做的是在柱体内部使用比空气层密得多的“局部细化”。可以创建一个只覆盖介质柱的域，然后在网格序列中单独给它设置最大单元尺寸，把它设为波长的十分之一以内。空气区的网格则可以设为波长的四分之一甚至更粗。这样在不牺牲精度的前提下，计算时间能压缩到原来的三分之一以下。

第二，二维能带图和Q因子扫描的计算量通常不大，但每个起点都手动设置扫描会很痛苦。写一个简单的循环脚本，把k路径、扰动参数和特征频率搜索范围都参数化。脚本输出后，你就可以一次跑出所有需要的数据，睡一觉起来直接收结果。

第三，如果电场的后处理结果异常、热点位置在每次网格调整后变化明显，这往往意味着你的结构处于“临界状态”，不是BIC发生了偏移，就是某个参数刚好落在模式简并点上。这时不要单纯追求加密网格，而是去扫描一个控制参数，比如增大一点切角尺寸，看模式是否能稳定。模式不稳定通常意味着你需要调整模型本身，而不是网格。

6.3 怎样保证结果的物理可解释性

仿真做得再漂亮，最后也需要能和实验或者文献对比。在这类超表面计算里，实际器件总有制造误差、衬底折射率的影响、入射光的发散角效应，这些都会让测量到的Q因子远低于理论值。所以有条件的话，在仿真后处理阶段加一个参数化分析，稍微改变一下周期大小或者柱高尺寸，记录Q因子的变化范围，这样能帮助理解实验上的容差范围。

如果你准备把这套方法复用在其他波段或者其他材料体系，比如把近红外硅换成可见光区的氮化硅，或者中红外的砷化镓，结构尺寸需要按比例缩放。核心逻辑是保持单元尺寸和波长比基本不变，同时注意新材料的折射率实部和虚部即损耗特性是否满足高Q条件。材料折射率虚部只要达到10的负三次方量级，就已经可能把Q因子压在几千以内，这会直接改变你的“非线性增强杠杆”效果。