comsol二次谐波SHG计算，单介质柱，展示图片即为仿真文件截图，所见即所得

今天在实验室折腾COMSOL的二次谐波（SHG）仿真时，发现单介质柱的非线性效应比想象中更有意思。特别是当基频波穿过圆柱结构时，表面极化产生的二次谐波场分布简直像科幻片里的能量场。直接上模型截图（图1），咱们边看边聊怎么从零开始搭这个模型。

先搞几何结构：半径500nm的二氧化硅圆柱，基底用硅。建模时手滑设置错了三次参数，最后发现用COMSOL的圆柱体工具直接画最省事。注意介质柱高度设为1μm，周围包裹空气域（图2红框）。这里有个坑——空气域边界至少离结构3倍波长，否则仿真结果会抽风。

!图1 模型结构示意图

（此处插入模型结构截图）

comsol二次谐波SHG计算，单介质柱，展示图片即为仿真文件截图，所见即所得

材料库里的二氧化硅默认没有非线性特性，得手动加戏。在材料属性里找到"Nonlinear Optics"标签，把二次非线性极化率χ²设置为1pm/V（实际材料参数看具体文献）。代码层面其实是在后台生成：

model.param.set('chi2', '1e-12[m/V]'); 

这个参数直接关联到极化强度计算项P=ε0χ²E²，直接影响谐波强度。但别急着改大数值——收敛性问题会让你想砸键盘。

物理场设置最烧脑：需要同时存在基波（1550nm）和二次谐波（775nm）。用两个"Electromagnetic Waves, Frequency Domain"接口嵌套，在频域设置里勾选"Second harmonic generation"。关键代码藏在研究步骤里：

study.step('freq').set('plist', {'1.94e14','3.88e14'});

这相当于告诉求解器：先算基频，再用基频场作为激励源计算二次谐波。记得在频域扫描里设置好两个频率点的计算顺序，顺序反了绝对报错。

网格剖分建议用物理场控制网格（图3箭头处）。手动调整过三次，发现边界层数设3层时，电场强度计算误差能控制在5%以内。用代码查看网格参数：

mphmesh(model,'mesh1').feature('size').set('hauto',3);

这个'hauto'参数控制自动细化等级，3级刚好平衡精度和计算量。算力不够的同志建议降到2级，虽然边缘场会有点锯齿感。

最后看结果（图4），基波在介质柱边缘形成环形热点，二次谐波则在轴线方向形成明显增强。后处理时用这个积分表达式抓取转换效率：

int2_solid = postint(fem,'esq2.normE','unit','W/m^2','dl','solid');

结果约0.3%的转换效率虽然不高，但作为微纳结构已经达标。重点观察场分布是否符合相位匹配条件——看截图中谐波场是否呈现预期干涉条纹。

调试时发现三个常见翻车点：

1. 极化方向设置错误导致无谐波产生（检查材料坐标系）
2. 光源偏振方向与晶体对称轴不匹配（改线偏振角度）
3. 扫频范围覆盖不全（用参数扫描找共振峰）

建议先做波长参数扫描锁定共振点（图5），能提升两个数量级的转换效率。这个单柱模型虽然简单，但改造成阵列或异质结构后，说不定能搞出些有趣的非线性光学现象。