COMSOL复现金螺旋超材料：圆偏振控制建模实战

1. 金螺旋超材料建模实战：从COMSOL复现到原理深挖

十年前发表在Science上的金螺旋超材料论文，至今仍是圆偏振控制的经典案例。作为光学超材料领域的标杆设计，这种结构展现出的宽带圆二色性特性令人着迷——金属结构居然能像光学滤波器一样选择性透过特定旋向的光。今天我们就用COMSOL完整复现这个模型，过程中会揭示许多教科书上不会讲的实操细节。

1.1 模型物理背景与核心挑战

金螺旋结构的精妙之处在于其手性几何特征与电磁波的相互作用。当圆偏振光通过螺旋阵列时，由于金属-介质界面的等离激元共振和结构本身的法布里-珀罗效应，会产生强烈的旋光选择性。这种效应可以用等效介质理论来描述：

\[
\epsilon_{eff} = \epsilon_0 \begin{bmatrix}
\epsilon_{xx} & i\epsilon_{xy} & 0 \
-i\epsilon_{xy} & \epsilon_{xx} & 0 \
0 & 0 & \epsilon_{zz}
\end{bmatrix}
\]

其中非对角项ε_xy正是产生圆二色性的关键。在建模过程中，我们需要特别注意三个核心问题：

1. 如何准确描述金属（金）在光学频段的色散特性
2. 怎样设置边界条件才能正确模拟无限周期阵列
3. 端口激励的圆偏振态定义与螺旋手性的匹配关系

1.2 几何建模技巧与参数化实现

论文中的螺旋参数非常精确：螺距2.7μm，半径0.9μm，线宽0.2μm，总圈数4圈。在COMSOL中，用参数化曲线建模比传统CAD软件更高效：

java复制double pitch = 2.7e-6;  // 单位：米
double radius = 0.9e-6;
double wireWidth = 0.2e-6; 
int turns = 4;

// 螺旋中心线参数方程
curve.parametric(
    radius*cos(t*2*pi*turns),
    radius*sin(t*2*pi*turns),
    pitch*t,
    t, 0, 1
);

// 螺旋横截面处理
wire = curve.sweep(circle(wireWidth/2));

这里有几个关键细节：

1. 使用归一化参数t∈[0,