在光学仿真里FDTD Solutions现在挂在 Ansys 名下的 Lumerical FDTD是我用得最顺手的工具之一。但说实话绝大多数新人甚至一些老手对“入射波长模式”这件事的理解都停留在表面——以为就是在光源面板里填一个波长数字然后点运行。直到仿真结果跟实验对不上才回头怀疑光源设置。今天这篇东西我想把这个经常被忽略、却直接决定仿真成败的参数组彻底掰开揉碎讲一遍。不管你是做等离激元、光子晶体、波导器件还是超表面只要你用 FDTD 做纳米光学仿真入射波长模式的设置逻辑你都绕不开。1. 入射波长模式到底在设什么一个被严重低估的全局参数1.1 波长模式不是“填一个数”那么简单很多初学者第一次打开 FDTD Solutions会先找光源Source。界面上确实有一个很显眼的“Wavelength”标签里面可以填中心波长、波长范围、频率点数。但我见过太多人在这里栽跟头有人做一个纳米球散射随手填了个 532 nm 单波长就开始跑结果消光峰死活算不出来有人想扫一个 400-800 nm 的透射谱结果把“wavelength start”和“wavelength stop”填反了得到一条负斜率的废曲线还有人压根没注意光源有“Broadband”和“CW”两种模式仿真跑完了才发现自己算的其实只是一个频点上的场分布。这里必须先厘清一个核心概念所谓“入射波长模式”本质上是光源在时间域里以什么形式激发结构。FDTD 算法的优势就是它能用一次仿真得到宽谱响应原理是发射一个包含多个频率分量的脉冲信号通常是高斯脉冲然后在时域上追踪电磁场在结构里的演化最后通过傅里叶变换把时域信号转成频域也就是波长域结果。这就是为什么 Lumerical 里叫它“Broadband”光源——它不是一个“单一波长平面波”而是一束“包含一段波长范围的脉冲波”。反过来“CW”模式就是 Continuous Wave连续波它只包含一个频率对应单一波长。CW 光源在时域上不衰减一直振荡适合看某个固定波长下的稳态场分布。但在实际工程里我们很少直接用 CW 扫描参数因为那等于每个波长都要跑一遍完整仿真效率太低。正确的做法是用 Broadband 光源一次性把所有关心的波长都“射”进去再在监视器里按频点提取结果。所以你填入的那个波长范围本质上定义了你要激励的物理过程的范围。范围选择是否合理直接决定了你能不能捕捉到结构的共振峰、透射谷、带隙边界这些关键特征。1.2 单频与宽谱的工程取舍这里我必须说一句得罪人的话很多人用单频波长仿真是因为看教程里“金球在 532 nm 下的近场增强”就是这么做的于是照猫画虎。但单频仿真往往掩盖了结构本身的色散特性。比如一个法布里-珀罗腔在 1550 nm 处是透射峰但如果你只仿 1550 nm你根本不知道这个峰有多宽、旁边还有没有别的模式。这对设计带宽、容差分析都是致命伤。我个人的习惯是第一次摸一个新结构绝对不做单频仿真。先用宽谱粗扫把整个感兴趣区域的响应曲线拉出来看清楚了哪里有峰、哪里有谷、大概有多宽然后再决定要不要在某个特定波长做精细的单频验证。宽谱粗扫就像拿广角镜头先拍全景单频精算就是换长焦镜头对准局部细节。两个配合用既不会漏特征也不会浪费机时。当然也有一些场景必须用单频。比如你要算某个特定波长下的电位移场分布用于非线性效应估算或者你要给实验上的一个特定激光波长做严格优化这时候 CW 光源反而更干净因为它不包含额外的频率分量不会给结果带来“看不到的邻居”的干扰。2. 光源类型决定波长模式怎么用平面波、高斯束、模式源各自的门道2.1 平面波与高斯束宽谱扫描的主力与陷阱在 FDTD Solutions 里最常见的宽谱入射是平面波Plane Wave和高斯束Gaussian Beam。很多教程一上来就默认你用平面波但平面波有个使用前提它覆盖整个模拟区域边界且入射方向垂直于光源平面。如果你研究的是周期性结构比如光栅、超表面平面波加周期性边界是标配。可一旦结构是有限尺寸的孤立粒子平面波宽谱扫出来的结果就会包含衍射效应带来的“伪共振”。高斯基模光束在实际实验中更容易对应激光光源它的束腰半径、发散角都不是随便填的。很多人把高斯束的束腰设得远远大于结构尺寸结果算出来的响应和平面的区别不大还白白增加计算量。我一般建议如果结构尺寸远小于光束束腰直接上平面波节省资源如果你的目的是模拟真实激光聚焦照射微纳结构那用高斯束并把束腰设在结构中心附近这样波长模式的展宽效应尤其是短波长端束腰变窄才不会被忽略。这里要特别提醒一个坑高斯束的命名里包含波长相关性。在 FDTD 里拉一个 Gaussian Beam 光源它默认的束腰半径是一个固定几何量但实际高斯光束的聚焦能力跟波长直接相关——同样一个束腰对 400 nm 的光已经是衍射极限附近的紧聚焦对 1600 nm 的光可能根本聚焦不到那么小。如果你用同一个物理参数去扫宽谱短波端和长波端的入射场空间分布其实发生了明显变化这时候你观察到的“波长依赖性”里混入了“光源聚焦状态变化”的贡献需要在后处理时小心甄别。2.2 模式源和偶极子什么时候才适合宽谱扫描做完波导和微环类结构的朋友一定没少用 Mode Source。模式源是把某个波导模式作为激励源注入到结构里它跟平面波最大的区别是模式源的横向分布由波导本征模决定入射“波长模式”的含义就变成了“你要注入哪个模式参数随波长如何变化”。用模式源扫宽谱有一个典型误区模式源在射入结构之前会根据你给定的模式位置和折射率分布计算一个基模或指定阶数模式。不同波长下这个模式的模场分布、有效折射率都不一样。FDTD 在宽谱模式下计算模式源时通常会自动在各个波长点重新求解本征模式所以它给出的注入模式是“随波长变化”的这本身没问题。问题是如果你在结构里设置了一段弯曲波导或锥形过渡模式之间的耦合会强依赖入射模式的有效折射率匹配情况宽谱扫描结果里出现的振荡很可能不是结构本身的光谱响应而是模式失配造成的干涉伪影。偶极子源我这边用得较少主要是在发光体近场增强、荧光增强这类涉及“点光源”的场景里。偶极子背后的物理是自发辐射偶极矩它天然是宽谱激励。做偶极子源时光源位置、取向要仔细考虑因为偶极子距离结构表面的距离会指数级影响近场耦合强度。比如你要仿真一个荧光分子靠近金纳米颗粒通常会在距离颗粒表面 5-10 nm 的地方放一个偶极子源。这个“入射位置”虽然不是波长但它的影响比波长更敏感——偶极子发出的倏逝波分量在短波长端更强频谱不均匀性更明显。3. 网格精度、材料色散与波长范围三个容易连环翻车的隐形变量3.1 每波长网格数固定网格精度为什么不够FDTD 方法是把麦克斯韦方程组在时间和空间上离散化空间步长网格尺寸直接决定了数值色散和稳定性。很多人对“mesh accuracy”的理解就是“越大越准”但实际上精度上去了、计算量也上去了一个三维结构很容易从几分钟算几个小时的台词变成“计算一晚上还在跑”。更关键的是网格精度跟波长是挂钩的。Lumerical 的 auto non-uniform mesh 默认会保证每个材料区域里大约有 10 个网格覆盖最短波长实际上受 mesh accuracy 设置控制。如果你用宽谱光源起始波长 400 nm、终止波长 1000 nm那么 FDTD 会根据 400 nm 这个最短波长来确定全局网格基准。很多人不修改 mesh accuracy认为默认精度够了然后在 900 nm 附近看到一个莫名其妙的振荡怎么检查逻辑都对不上。其实原因往往就是在较长波长区域网格相对过于细密导致数值色散和不稳定性溢出或者反过来如果 mesh accuracy 设太低长波段虽然覆盖了足够多网格但短波段已经离散化到“面目全非”。我建议做宽谱扫描时先用一个较低的 mesh accuracy比如 2配合宽谱粗扫目的是找特征位置。找到关键特征之后再在那个峰/谷周围切出一个小波长范围提升 mesh accuracy 到 4 甚至 5同时加一两个 mesh override 把结构表面或尖角处的网格局部加密。这样既节省时间又保证关键波长下的数值收敛性。记住FDTD 的每一段网格都对应一个波长范围的数值分辨能力网格精度不是“越大越好”而是“跟波长和目标特征匹配才好”。3.2 材料色散在宽谱模拟中经常被忽略我一直觉得材料模型比光源设置本身更坑人。做宽谱扫描时材料折射率的色散折射率随波长变化是必须考虑的。比如金、银这类等离激元材料在可见光到近红外区间折射率变化极其剧烈。如果你只用某个单一波长下的折射率常数去跑整个宽谱结果基本是废的——共振峰位置、谱线宽度都会漂移得离谱。Lumerical 的材料库里提供了多种色散模型选项比如多系数模型Multi-coefficient model、Drude 模型、Lorentz 模型、Sellmeier 模型等。对于金属默认的 AuGold- Johnson and Christy 数据是常见的经典选择对于介质比如 Si、SiO2Sellmeier 模型也比较准。问题是你在光源波长范围设定好后有没有检查过所用材料在该范围内的折射率行为是不是合理的有次我做一个 Si 波导结构模拟 1260-1360 nm 的传输发现结果里出现了异常的吸收峰排查了很久才发现材料库里我误选了“Si (Palik)”的损失数据这组数据在近红外区有一个杂质吸收峰跟我的结构完全无关。解决办法其实很简单在跑正式仿真之前先画一个材料的折射率-波长曲线肉眼确认它在你的波长区间内没有异常跳变。这步成本极低但真的能救你于水火。3.3 PML 边界与入射角度的配合问题还有一个容易被甩锅给“入射波长模式”的变量是模拟区域边界的设置。Lumerical 里常用的边界是 PMLPerfectly Matched Layer中文叫完美匹配层用来吸收向外传播的电磁波。PML 的吸收性能跟入射角度和波长都有关系。在宽谱扫描中如果结构散射特别强大量的光会以掠射角几乎平行于边界的角度打到 PML 上这时候默认 PML 配置的吸收效果会大打折扣表现为边界反射率升高结果里出现明显的驻波振荡。做粒子散射、等离激元增强这类模拟时我一般会把模拟区域边界留得离结构表面足够远至少 1-2 个中心波长并且把边界类型从默认的“PML”改成“Steep Angle PML”专门应对大角度入射的场景。这同样适用于斜入射平面波当你给平面波设置了入射角度边界处的场分布已经不是垂直入射那么简单了PML 必须跟着调整否则你的“angle response”曲线里会混入大量来自边界反射的伪结构。4. 监视器与频谱提取怎么把波长模式变成可用的数据4.1 2D 监视器与透射/反射计算光有宽谱光源还不够你要在结构后面放监视器Monitor才能把“入射波长模式”转换成可读的光谱。Lumerical 里最常用的是 2D 监视器它可以记录通过某个平面的电场/磁场分量。计算透射率或反射率时关键是定义归一化基准——一般要跑一次“没有结构”的空白模拟作为参考然后把有结构的透过功率除以空白参考的透过功率才能得到真正的透射率。这里有个细节经常被忽略光源是宽谱脉冲时监视器记录的是时域场后处理时做傅里叶变换的频率范围必须跟光源的波长范围对齐。如果监视器面板里设置的 frequency points 比你实际想要的分辨率低你会看到谱线有锯齿如果 frequency profile 设成了 single wavelength那你扫宽谱光源等于白扫。我见过一个同事光源列表里明明设了 400-900 nm 宽谱扫描但监视器 Profiles 里只勾选了一个 532 nm 的频率点出来的结果自然只有一个波长他还在那纳闷为什么“宽谱扫描”输出不是连续曲线。4.2 频率点数与波长谱的锯齿问题“Frequency points”这个参数是很多人调参时的盲区。Lumerical 监视器里的频率点数决定了对时域信号做傅里叶变换时的输出频率采样密度。频率点是均匀分布在线性频率轴上的而波长是频率的反比所以如果你在 400-800 nm 范围均匀取了 100 个频点换算到波长域后短波端的波长间隔比长波端更密。这在提取光谱时就意味着波长长端的分辨率其实低于波长短端。遇到这种情况我会选择把频率点数提高到 200-500或者直接在后处理脚本里用均匀波长数组重新插值。否则当你绘制“透射率 vs 波长”曲线时长波端可能因为采样太稀疏而漏掉窄谷。尤其是等离激元共振峰常常只有几十纳米宽如果长波端采样不够峰值位置和半高宽都会测不准。频谱锯齿问题也同样跟频率点数太少有直接关系当时域信号还没衰减完全傅里叶变换会在频域产生振荡频率点数多可以缓解“看到”锯齿但不能根治根治的方法是延长仿真时间增加 auto shutoff 的宽容度或者加窗函数。还有一点透明材料里的共振峰往往很窄时域场衰减得很慢默认的 auto shutoff level通常是 1e-5 或者 -80 dB可能不够。你要留意一下模拟运行结束的信息如果提示“auto shutoff reached during simulation”确保是在时域场确实衰减完才停的如果是“simulation time reached limit”强制结束那说明还有一些尾场没衰减干净频谱会有轻微波纹。遇到这种情况把 simulation time 加长或者把 auto shutoff level 设得更严格比如 1e-8是两张最好用的牌。5. 实战推演从金纳米球案例看波长扫描的完整链路5.1 第一步先粗扫别急着精算为了把上面这些抽象理论落到具体操作上我拿一个经典例子来讲直径 60 nm 的金纳米球放在折射率 1.33 的水环境里仿它的消光光谱。物理直觉告诉我们金球的局域表面等离激元共振LSPR应该在 520-560 nm 附近有一个吸收峰散射峰稍微红移一点。这种结构的共振峰宽度不算太窄适合演示。我的第一步永远是在 Simulation 区域里设一个大一点的盒子把背景折射率设为 1.33光源用平面波 Broadband波长范围直接拉到 400-800 nm网格先给一个松散的精度mesh accuracy 2监视器在球后方放一个 2D 透射监视器同时在球的周围放一个 transmission box 用于计算吸收和散射。这样跑完一次之后我不用去细究某个点的场分布只是看消光曲线的大形貌峰大概在哪个位置、宽度多宽、有没有多峰结构。这个步骤通常在几分钟内能跑完三维结构如果网格太粗可能也就 10-20 分钟。它的意义在于给你整个参数空间画了一张“地图”你知道后面该往哪儿使劲。5.2 第二步锁定共振区加密频点粗扫结果出来假设峰值出现在 540 nm 附近但我发现波长 500-600 nm 之间的采样点太少曲线有些毛糙。第二步我会做两件事第一把光源的 wavelength range 缩窄到 480-620 nm把 frequency points 设置到 200 甚至 300第二把 mesh accuracy 提升到 4并在金球表面加一个 mesh override设定表面网格尺寸为 2 nm 左右。为什么光源范围要跟着缩窄因为 FDTD 的自动网格是根据整个光源范围里的最短波长来定标准网格的范围越宽最短波长越短全局网格越细、耗时越长。缩窄到共振区附近后最短波长也变长了网格可以相应放宽计算效率大幅提升而你在共振区的频率分辨率反而因为 frequency points 堆叠而更高了。这一步做完出来的消光曲线应该是光滑的共振峰位置和半高宽都稳定了。这一步还经常暴露一个问题粗扫时的共振峰位置在加密后发生了移动。这很正常尤其是金属等离激元结构对网格敏感。我的建议是如果加密前后峰位移动超过 5-10 nm那你还得继续把 mesh accuracy 往上调或者调整 override 网格直到峰位对网格不再敏感为止。这就是所谓的“网格收敛性测试”是光学仿真里过不了审就会被打回的一步。5.3 第三步单频验证与场分布检查宽谱扫描确认了谱线形状之后我会回到单频。选定共振峰对应的波长比如 538 nm把光源模式切成 CW 单波长模式重新跑一次短时仿真然后提取电场分布图、电荷分布图。这一步是为了确认共振模式的近场增强位置在哪里、热点是不是你预期的那个界面。为什么宽谱模式不能直接替代这个单频验证因为宽谱模式的场分布是所有波长贡献的叠加。如果你在宽谱模式下查看某一时刻的电场图看到的实际上是多个波长的干涉混合场没法干净地展示某一特定共振模式的形态。切到单频之后场分布图像干净得多物理图像也更容易汇报给同事或写进论文。单频验证时要注意CW 光源不衰减仿真不能靠 auto shutoff 自动结束时间步数到了就停。如果你只看场分布一般跑 100-200 个光学周期就足够了。我在这一步最常见的问题是光源距离结构太近导致高次倏逝波还没衰减就撞上结构让近场图出现一些假的高频振荡条纹。处理方式是把平面波光源放在模拟区域边界附近用一块平整的空白区过渡到结构区域。6. 我把这些坑记成了自己的检查清单6.1 每次新建模拟我都会先过一遍的清单做 FDTD 仿真这几年我总结了一套自己的排查清单。每次项目出了问题我都会从头到尾按顺序过一遍80% 的情况能在前五个条目里找到根因检查项具体内容对应常见问题光源波长范围是否覆盖目标响应区间start stop漏掉共振峰、曲线方向颠倒光源类型Broadband 还是 CW跟目标匹配单频结果误当宽谱、场分布叠加混叠光源位置距离结构是否足够远是否在 PML 安全距离外倏逝波扰动、边界反射干扰材料色散折射率-波长曲线是否平滑模型是否适用伪吸收峰、峰位偏移网格精度是否做过收敛性测试override 是否正确峰位漂移、谐振展宽或变窄监视器频率点频率点数够不够波长长端是否采样不足频谱锯齿、漏窄峰PML 边界是否需要 steep angle距离结构是否够远驻波伪影、反射干扰这张表不是死的但每次排查时我都按这个顺序走。有一次结果里出现了一个完全不在物理预期内的小峰我把材料、网格、监视器翻了个遍都没发现问题最后发现是光源的“amplitude”被我误设成了一个带频率依赖的表达式。所以后来我在清单里又加了一条检查光源参数表达式里有没有隐性的波长或频率变量。6.2 一些省时间的习惯最后分享几个我自己觉得很省时间的小习惯。第一个习惯是“先粗后精两次成型”。不管结构多复杂永远不要一上来就把 mesh quality 拉到 7、频点设到 1000。先用粗网格宽谱扫全貌再用细网格锁细节最后用单频验证场分布。这个流程不仅省时间还能让你在中间步骤中发现一些之前预期之外的物理现象。第二个习惯是善用脚本批处理。Lumerical 有比较成熟的脚本接口。我会写一个简单的脚本循环自动改变入射角或环境折射率批量跑一组宽谱扫描然后把结果统一导出。这样比自己手动改参数、一个一个点运行要省太多时间更不容易出错。第三个习惯是归档光源设置截图。我是个很依赖经验记录的人每次项目到一个里程碑我会把光源面板、监视器设置、网格设置的截图存到项目文件夹里。这样过了两个月后回看数据我还能清楚地回忆起当时的仿真条件而不是看着一堆数字发呆。FDTD Solutions 的入射波长模式说到底是仿真者与物理模型之间的一次约定你在时间域里发一个特定形状的脉冲它回你一份宽谱响应。你把约定理顺了它就是最趁手的工具理顺不好它就是吞掉你大量时间的坑。希望这篇文章能让你少踩几个我踩过的雷。