comsol太赫兹开关。 编号026
最近实验室里几个师弟在折腾太赫兹开关，整天对着COMSOL界面挠头。这玩意儿确实有点意思——用电磁波控制光路通断，关键频段还卡在0.1-10THz这个特殊位置。咱们先别急着点仿真按钮，得把模型骨架搭结实了。
先整活个波导结构。COMSOL的几何建模跟搭积木似的，但要注意单位换算。随手贴段参数设置代码：
double lambda = 300e-6; // 1THz对应波长 double wg_width = lambda*0.8; model.param().set("gap", "10e-6"); // 开关间隙 model.param().set("L_t", "lambda*5"); // 太赫兹路径长度这里lambda算的是真空波长，实际介质中还得考虑折射率。gap参数控制着开关灵敏度，就像水龙头旋钮——太小容易击穿，太大又关不严实。
重点在端口设置。太赫兹波的模式可比微波复杂，得用模式分析先探探路。看这段端口边界条件配置：
model.physics("emw").feature("port1").set("PortType", "Rectangular"); model.physics("emw").feature("port1").set("WaveExcitation", "TE"); model.physics("emw").feature("port1").set("TE_mode", new String[]{"TE10"});注意TE模式的下标从0开始计数，别手滑写成TE11。有次师弟把模式设错了，结果仿真出来的电场分布活像毕加索的画，整个实验室笑了半小时。
材料属性才是重头戏。石墨烯这货在太赫兹频段会变身，电导率随偏压变化比川剧变脸还快。用分段函数描述：
double V_bias = 3; // 偏置电压 String sigma_expr = "if(V_bias>2, 1e4*(V_bias-2), 0.1) [S/m]"; model.physics("emw").feature("mat1").set("sigma", sigma_expr);这表达式相当于给开关装了个触发器。实际做参数扫描时，建议用COMSOL的扫参功能遍历0-5V，步长别超过0.2V，太赫兹响应曲线可能有突变点。
最后说说网格。自适应网格生成虽然方便，但在关键区域还是手动加密更靠谱。比如开关间隙处的网格得细过绣花针：
model.mesh("mesh1").feature("size1").set("hauto", 3); model.mesh("mesh1").feature("size1").set("hgrad", 1.5); model.mesh("mesh1").feature("size1").set("hcurve", 0.7);这三个参数控制着网格生长率和曲率适应度。有次偷懒全用默认值，结果电场强度算出来比实际高了两个数量级，差点以为要搞出什么黑科技。
跑完仿真记得看S参数矩阵。特别关注S21在开关状态下的衰减量，-30dB算及格，-40dB能发二区。要是发现谐振峰位置飘移，八成是边界条件没设对或者材料参数飘了。搞太赫兹器件就像养猫，参数稍微不对就给你摆脸色看。
（代码示例基于COMSOL 5.6的Java API，实际操作时建议配合App开发器食用更佳）
