Comsol光学仿真模型：包括纳米球/柱 Mie散射多级分解

---

折腾纳米颗粒的光学仿真？先搞懂Mie散射多级分解

最近在COMSOL里搞纳米颗粒的光学仿真，发现Mie散射这玩意儿真是让人又爱又恨。尤其是当你要拆解散射场的多级分量（比如电偶极、磁四极这些）时，参数设置一不小心就翻车。今天就来唠唠，怎么用COMSOL把纳米球/纳米柱的散射场拆得明明白白。

---

第一步：先搭个纳米球的几何模型

COMSOL建模的核心就两个字——别急。

1. 几何建模：直接拉个球体，半径设成50nm（典型纳米颗粒尺寸）。
2. 材料赋值：选金（Au）或者银（Ag），从内置材料库调取介电常数，记得勾选“分散材料”选项，否则波长变化时参数不会自动更新。
3. 背景场设置：用“散射边界条件”包裹整个结构，再套一层“完美匹配层”（PML），防止反射干扰。

// COMSOL的模型树操作（类似脚本）
model.geom.create("geom1", 3);
model.geom("geom1").create("sph1", "Sphere");
model.geom("geom1").feature("sph1").set("radius", "50e-9");
model.material.create("mat1");
model.material("mat1").propertyGroup("def").func("eps").set("au");

为啥要用PML？

简单说，PML就是个“吸波黑洞”，让出去的电磁波别反射回来捣乱。如果不加，仿真结果可能会像在镜子里照鬼影——全是数值误差。

---

关键一步：电磁场频域求解

这里重点在波动光学模块的“电磁波，频域”接口。

1. 入射波设置：选平面波，极化方向沿x轴，波长范围建议400-800nm（可见光区）。
2. 扫频技巧：别傻傻地一个波长一个波长跑，用参数化扫描，直接输入range(400e-9, 50e-9, 800e-9)，自动跳波长。

model.study.create("std1");
model.study("std1").create("param", "Parametric");
model.study("std1").feature("param").set("pname", "lambda");
model.study("std1").feature("param").set("plist", "range(400e-9, 50e-9, 800e-9)");

注意坑点：

Comsol光学仿真模型：包括纳米球/柱 Mie散射多级分解

如果纳米颗粒尺寸接近波长（比如100nm以上），别偷懒只用“散射场”变量，得用“总场-散射场”公式，否则远场计算会崩。

---

重头戏：Mie散射的多级分解

Mie散射的理论公式能把散射场拆成无限级数，但仿真时我们只取前几级（比如l=1到4）。

COMSOL怎么操作？

1. 在“结果”里新建“远场投影”数据集，选球坐标系。
2. 用“球谐函数”展开，输入公式：
   
sum((2l + 1)/(2k^2) * (|al|^2 + |bl|^2), l=1..4)

   这里al和bl对应电多极和磁多极系数。

代码分析：

COMSOL其实内置了Mie系数计算，但需要手动提取。比如用emw.al和emw.bl调用。

// 提取电偶极(l=1)的散射效率
a1 = emw.a_l(1);
scattering_efficiency = (abs(a1)^2) * (2*1 + 1)/(2*emw.k^2);

翻译成人话：al和bl对应不同级次，l=1是偶极，l=2是四极…数值越大，对应更复杂的场分布。

---

结果验证：纳米球 vs 纳米柱

1. 纳米球：散射谱通常有一个主峰（偶极共振），高阶峰可能藏在短波长区域。
2. 纳米柱（直径50nm，高100nm）：各向异性导致偏振依赖，TM波可能出现磁偶极模式（b_l系数飙升）。

实战技巧：

在“派生值”里画emw.al(1)和emw.bl(1)随波长的变化曲线，如果两条线交叉，说明该波长点存在模式竞争。

---

最后说点人话

Mie散射多级分解的本质，就是把复杂的场分布拆成“乐高积木块”。COMSOL的优势在于能直接调用球谐函数系数，不用自己手撕公式。但记住：

• 别迷信默认设置，PML厚度和散射边界的位置都得试错；
• 高阶模式（l>3）对网格精度敏感，加密到λ/10以下再谈收敛；
• 纳米柱仿真比球麻烦十倍，先搞定球再去挑战柱子。

仿真跑崩了？正常。回头检查材料参数和边界条件，多半是这两个在搞鬼。