FDTD脚本建模实战：纳米柱阵列生成与参数扫描自动化

我们平时用 FDTD 仿真（比如 Lumerical FDTD Solutions）做光学设计，绝大多数人上手都是从图形界面（GUI）拖拽结构开始的。鼠标点一点，画个矩形、圆柱，设个材料，好像也挺方便。但一旦你开始做参数扫描、批量优化、模型复用，或者需要几十上百个结构单元时，纯 GUI 操作会让人崩溃——重复拖拽、重复设置、重复导出，既慢又容易出错。这时候，脚本就是唯一靠谱的出路。这篇教程我会聚焦在“用 FDTD 脚本创建结构”这件事上，从最基础的命令语法，到一个可以直接抄作业的纳米柱阵列建模案例，再讲参数扫描和排错经验。内容面向已经能用 GUI 建简单模型、但想转向脚本自动化建模的人，也适合刚接触 FDTD 脚本、被语法和单位折腾得头疼的新手。

1. 为什么要用脚本创建结构——当 GUI 建模不够用的时候

1.1 脚本的真正应用场景

我先说说我个人的使用习惯。并不是说 GUI 建模不好，而是脚本建模解决的痛点完全不一样。GUI 适合“一次性、简单结构”的验证，比如你想快速看一个单层薄膜的透射率，拖个矩形、设个光源就能跑了。但现实中稍微实际一点的需求，GUI 都很难受。

第一种典型场景是参数扫描。比如说你想研究纳米柱直径从 100 nm 到 200 nm、步长 5 nm 的情况下，透射光谱会怎么变化。用 GUI 的话，你需要手动改 21 次直径，每次都要点开属性面板、输入数值、运行、记录结果。一次两次还能忍，规模一大或者变量变成两三个（直径、周期、高度），手动操作就完全不可行了。而脚本里只需要一个for循环就能搞定，跑完之后结果还能自动导出，去洗杯咖啡回来数据就齐了。

第二种场景是复杂周期性阵列。超表面、光子晶体、衍射光栅这类器件，往往需要在模拟区内排布成百上千个单元结构。你不可能手动拖几百个六边形、圆柱到 FDTD 里，就算你有这个耐心，后面改一个参数也要全部重新来一遍。而脚本建模时，一个嵌套循环就能轻松生成整个阵列，而且修改参数只需要改一个变量名，后边的循环和设置全部自动跟进。

第三种场景是模型复用和交接。脚本本质上就是一个可复现的记录文件。半年前建的模型，三个月前改的设置，写得规不规范，一看脚本就全都清楚了。你不需要去翻仿真文件里的哪个隐藏属性被改过，也不需要担心发给别人之后对面看不懂你的边界条件和网格设置。一份写清楚注释的脚本，就是一份最标准的建模说明书。我自己在团队里做项目时，能明确感受到脚本交付比.fsp文件交付靠谱得多——接收方哪怕没有 GUI 操作基础，也能照着脚本理解你的建模逻辑。

1.2 脚本与 GUI 建模的边界与取舍

有些人会有种误解，觉得用了脚本就等于放弃了图形界面。其实不是。Lumerical 的脚本模型和 GUI 模型是同源的，脚本里创建的每个对象都会实时出现在布局编辑器和结构树中。反过来，你在 GUI 里手动建好的模型，也能通过保存脚本或者导出命令把建模过程转成脚本。两边是联动的，关键看你在这个环节用哪个更顺手。

我的一个习惯是“GUI 看结构，脚本管流程”。建模型的过程中我通常会穿插使用：先用脚本批量生成主体结构，然后用 GUI 手工微调个别材料参数或者检查几何是否重合，再回到脚本里设置扫描和优化流程。这样做不仅比纯 GUI 快，也比纯脚本试错更直观。毕竟脚本创建完结构之后，立刻可以通过select命令选中对象，再在布局编辑器里看它的具体位置和形状，完全不耽误可视化。

另外要说清楚一个认知：FDTD 脚本不是一种需要先系统学完“编程”才能上手的语言，它本质上和 Python 有点像，但比 Python 简单得多。你只需要掌握创建对象的命令、设置属性的语法、循环和条件判断的基础写法，再加上几个常用的文件操作函数，就可以覆盖 90% 的日常建模脚本需求。真的不用怕，我最早开始学的时候也就花了半天时间，能用的核心命令二三十个封顶了。

还有一个我踩过的坑：别一上来就追求“全脚本化”。非要让一个不规则形状的 3D 模型完全用代码生成，那是跟自己过不去。Lumerical 脚本擅长的是参数化规则结构，比如矩形、圆柱、球、圆环、多边形，以及这些基本原语的组合和阵列。真正不规则的复杂模型，通常的做法是用 CAD 或者建模软件导入.step、.stl文件，这个在脚本里一条import命令就能搞定。两者结合起来用，效率和精度都能保证。

2. 从零开始：脚本创建结构的基础语法与常用命令

2.1 结构原语的创建命令

Lumerical FDTD Solution 的脚本（官方叫 Lumerical Script，简写有时候也叫 LSF）建结构使用的一组命令，核心就是“创建对象”和“设置属性”。创建原语的命令很直白：addrect加矩形、addcircle加圆、addsphere加球、addcylinder加圆柱、addpoly加多边形、addring加圆环，以及addstructuregroup加结构组。

这些命令执行之后，会在结构树里新增一个对象，紧接着就要用set来设置它的属性。例如你创建一个硅材料的矩形柱子：

addrect; set("name", "Si_pillar"); set("x", 0); set("y", 0); set("z", 0); set("x span", 0.2); set("y span", 0.2); set("z span", 0.5); set("material", "Si (Silicon) - Palik");

这段脚本执行后，一个中心位于原点、横截面 0.2 μm × 0.2 μm、高度 0.5 μm 的硅矩形柱就出现了。这里我要重点提醒一个新手最容易踩的坑——Lumerical 的默认单位是微米（μm）。如果你脑子里想的是纳米，直接填 100，那建出来的结构就是 100 μm，等于在你看不见的尺度上放了个庞然大物。所以写脚本之前一定要先明确单位，或者干脆通过set("x span", 100e-3)这么写，用科学计数法强行转成微米。

2.2 属性设置与单位陷阱

set命令的通用格式是set("属性名", 值)，属性名必须是对象属性面板里存在的标准名称。常用的几何属性包括位置（x、y、z）、尺寸（x span、y span、z span），材料属性就是material，名字和 GUI 材料库里显示的名字保持一致就行。还有index这类渲染优先级参数，之前在多层结构里需要指定哪个物体盖在哪个物体上面时，就要手动改。

另外还有一个很实用但容易被忽略的命令是set("override color", 1)和set("color", 某个RGB值)，方便你区分同一模型里的不同部件。虽然不设置也不影响仿真结果，但脚本建模的后期，结构数量一多，不同结构用不同颜色显示对排查几何重叠问题特别有帮助。

单位陷阱我再展开说。Lumerical 里位置x、y、z和尺寸x span等的单位是微米，但有些属性用的不是几何单位，比如角度用的是度（degree），材料折射率是纯数值，薄膜层层的厚度单位反而要根据介质材料属性而定。我见过有人写光栅周期时填了 0.7（其实他想表达 700 nm），结果仿真出来的衍射波长完全不符合物理预期。排查了半天，最后发现是单位少了一个数量级。所以每次运行脚本前，我都会在脚本开头写一个clear;，紧接着输出一行所有关键尺寸的微米值，确认数字没问题再继续。

2.3 用 for 循环批量创建阵列结构

脚本建模最爽的地方，就是可以用循环批量生成结构。我拿二维光栅举例。假设要在 5 μm × 5 μm 的区域里做一个周期为 0.7 μm 的正方形柱阵列，直径 0.3 μm、高度 0.6 μm。手写几十个addcylinder显然不现实，用循环写就非常简洁：

# 周期和尺寸 period = 0.7; diameter = 0.3; height = 0.6; Nx = 7; # x 方向柱子数量 Ny = 7; # y 方向柱子数量 for (i = -floor(Nx/2); i <= floor(Nx/2); i++) { for (j = -floor(Ny/2); j <= floor(Ny/2); j++) { addcylinder; set("name", "pillar_" + num2str(i) + "_" + num2str(j)); set("x", i * period); set("y", j * period); set("z", 0); set("radius", diameter / 2); set("z span", height); set("material", "Si (Silicon) - Palik"); set("override color", 1); set("color", [1, 0.5, 0]); # 橙色显示 } }

这段代码里我用了floor、嵌套循环、字符串拼接num2str，这些都是 Lumerical 脚本自带的基础功能，和 Python 的语法风格非常像。跑完之后，结构树里会出现 49 个圆柱，布局编辑器里可以看到整个阵列。这里还有个小细节：set("name", ...)时尽量把循环索引加到名字里，后面做参数扫描、选择性修改某个柱子或者单独选中检查时，会方便得多。

还有一个和循环配合很高的命令是select。循环里创建完对象之后，如果要给某些局部结构单独设置属性，可以用select("pillar_1_-1")精确选中它，再set属性。这种“按名字操作”的方式比用鼠标在结构树里点来点去快得多，尤其当结构数量多到结构树都卡的时候，命令行反而更流畅。

3. 实战案例：用脚本搭建纳米柱阵列模型

3.1 建模目标与参数设定

这一节我会完整演示一个纳米柱阵列建模案例，这个案例几乎涵盖了脚本创建结构的所有核心操作。目标：在 FDTD 模拟区内创建一个硅纳米柱阵列，中心在原点，结构基底是二氧化硅衬底，上方是周期排列的硅柱。我们需要扫描硅柱高度对透射率的影响，高度从 300 nm 到 500 nm，步长 20 nm。

先明确基本参数：

• 模拟区尺寸：x 方向 1.4 μm，y 方向 1.4 μm，z 方向 3 μm（含 PML 边界）
• 阵列周期：0.7 μm
• 柱子半径：0.15 μm（直径 0.3 μm）
• 柱子材料：Si (Silicon) - Palik
• 衬底：SiO2 (Glass) - Palik，厚度 1 μm
• 光源：平面波，从上方垂直入射，波长范围 400 nm 到 800 nm

这里要特别提一下模拟区边界条件的设置。对于周期性结构，x 和 y 方向应该用周期性边界条件（periodic），z 方向用 PML。这样模拟区只需要包含一个周期就可以等效无限阵列，极大地节省计算资源。脚本里设置边界条件的方式是set("x bc", "periodic");等命令，我们后面会用到。

3.2 完整脚本解析

这个案例完整脚本如下，每一段我都加了注释说明意图：

# 纳米柱阵列建模脚本 clear; switchtolayout; # ---------- 基本参数 ---------- period = 0.7; # 周期，单位 μm radius = 0.15; # 圆柱半径，单位 μm h_si = 0.4; # 硅柱高度，单位 μm (初始值) h_glass = 1.0; # 衬底厚度，单位 μm z_glass_top = -0.5; # 衬底上表面 z 坐标 z_si_bottom = z_glass_top; # 硅柱底部贴着衬底 # ---------- 创建衬底 ---------- addrect; set("name", "SiO2_substrate"); set("x", 0); set("y", 0); set("z", z_glass_top - h_glass / 2); set("x span", 5); # 衬底在 x、y 方向做大一些，模拟无限大衬底 set("y span", 5); set("z span", h_glass); set("material", "SiO2 (Glass) - Palik"); # ---------- 创建硅柱阵列 ---------- N = 3; # 单方向柱子数量，3×3 阵列 for (i = -floor(N/2); i <= floor(N/2); i++) { for (j = -floor(N/2); j <= floor(N/2); j++) { addcylinder; set("name", "Si_pillar_${i}_${j}"); # 注意这里的字符串模板写法 set("x", i * period); set("y", j * period); set("z", z_si_bottom + h_si / 2); set("radius", radius); set("z span", h_si); set("material", "Si (Silicon) - Palik"); set("override color", 1); set("color", [0.8, 0.4, 0.2]); } } # ---------- 设置模拟区域 ---------- select("FDTD"); set("x span", 1.4); set("y span", 1.4); set("z span", 3); set("x bc", "periodic"); set("y bc", "periodic"); set("z bc", "PML"); set("z min", -1.5); set("z max", 1.5);

这里面有几个地方值得详细解释。首先是switchtolayout;，这个命令会把 FDTD 切换到布局设计模式（layout mode），而不是仿真计算模式。建模阶段一定要先执行这个命令，否则你后续添加的结构可能不会正确显示在布局环境中。这个命令我几乎每个建模脚本都会写，算是个基础习惯。

其次是set("name", "Si_pillar_${i}_${j}")这种写法，它利用的是 Lumerical 脚本对字符串模板的支持，比"Si_pillar_" + num2str(i) + "_" + num2str(j)更简洁易读。需要注意的是，模板里的变量会自动转换类型，但有些很老版本的 Lumerical 可能不完全支持${}写法，如果不支持就退回用num2str拼接。

第三是模拟区边界的设置。set("z min", -1.5)和set("z max", 1.5)定义了 FDTD 模拟区在 z 方向的范围，配合z span的数值，这三个需要保持一致。否则你设置了z span是 3，但z min和z max的差不是 3，Lumerical 会以z span优先还是z min/z max优先，不同版本行为不太一样，容易造成边界条件设置和预期不符。我写脚本时习惯于显式设置z min和z max，避免歧义。

3.3 结构组 dry-run 与 GUI 联动

建完结构之后，千万别急着直接跑仿真。Lumerical 脚本里常用的一个技巧是“先 dry-run”——也就是只建模、不运行仿真，在布局编辑器里目视检查所有几何结构的位置是否正确。我通常是这么做的：脚本跑完之后，打开布局编辑器，从正上方（xy 平面）和侧方向（xz 平面）各看一眼，确认柱子间隔均匀、和衬底贴合、没有明显重叠。

这里就要用到overlap检查了。Lumerical 脚本里可以逐对选中物体，然后点击右键的“check overlap”功能，也可以在命令行用相关函数检查结构重叠。但我个人更倾向于用目视检查加上一个“渲染顺序”的技巧——多个物体在同一位置时，后创建的物体会遮挡先创建的。如果你想看衬底和柱子的交界面，可以把衬底的透明度调高。设置透明度的方法是在结构属性面板里找“render type”，或者通过脚本set("alpha", 0.5)调整。

另外一个重要的联动操作：当你确认模型正确后，可以执行：

save("nanopillar_array.fsp");

把当前布局保存为仿真文件。以后打开这个文件，所有用脚本创建的结构都还在，而且你可以继续用脚本修改、运行、提取数据。再配合 GDS 导入或者.stl扩展，Lumerical 可以和其他 EDA/CAD 工具无缝衔接。

还有一个小技巧：如果你需要把脚本创建的几何导出到其他软件查看，可以选中结构后，用菜单栏的 File -> Export -> Geometry 导出为.stl或者.step格式。这些文件可以直接导入到 SolidWorks、COMSOL 或者自家画图工具里做进一步处理。在高频电磁仿真和光学设计中，几何文件交换是很常见的需求。

4. 进阶技巧：用脚本驱动参数扫描与优化

4.1 参数关键字的定义与使用

结构建好了，下一步就是让脚本帮我们批量跑仿真，分析不同结构参数下的光学响应。Lumerical 脚本提供了强大的参数扫描功能（sweep），你不需要自己写循环去一次次修改结构并运行——虽然那样也是可行的——更规范的做法是定义参数关键字（parameter）和扫描任务。

参数关键字本质上是一个模型里的全局变量，脚本、结构属性和扫描任务都可以引用它。例如我们想扫描硅柱高度，可以这样定义：

addglobalsweep; setsweep("sweep_1", "sweep type", "parameter"); addglobalsweepparameter("sweep_1", "h_si_scan", 0.3, 0.5, 11);

这里addglobalsweep创建了一个参数扫描任务，addglobalsweepparameter定义了扫描的变量h_si_scan从 0.3 到 0.5、均匀取 11 个点。然后你还需要告诉 Lumerical，在扫描过程中要用这个变量值来更新结构中“z span”属性。这一步需要在set("z span", h_si_scan)里引用同名的参数关键字。

实际脚本中，更稳妥的流程是先把关键几何属性指向一个中间变量，再用这个变量去设置属性。比如：

h_si = h_si_scan; # 将扫描参数赋值给几何变量 select("Si_pillar_0_0"); set("z span", h_si);

addglobalsweepparameter中的数值会自动按等差序列生成，不需要你手动构造循环。扫描任务建好之后，在 GUI 的优化和扫描窗口里就能看到这个任务，直接运行即可。

4.2 自动批处理与结果导出

参数扫描跑完之后，我们通常需要把每个参数下的结果汇总分析。Lumerical 脚本有专门的数据提取和保存命令。比如我们只需要透射率，可以这样写：

select("T"); # 假设监视器名为 T T = getresult("T", "T"); transmission = T.T; # 保存数据到文件 write("transmission_vs_height.txt", transmission);

但这只是在一次仿真结束后的数据导出。配合扫描任务后，我们可以把数据保存放到每次仿真结束时自动触发的地方，比如在“sweep”设置里的“post processing”回调函数中写入文件。我这里建议一个更直觉的做法：在扫描之前写好一个用于保存多个结果数据的初始化文件，然后每次扫描循环里用matlabsave或者write导出，文件名带上当前参数值，例如：

write("T_data_h_" + num2str(h_si_scan*1000) + "nm.txt", transmission);

批量跑完所有仿真之后，你就能得到一个带参数的文档集合，后续用任何数据处理软件（Python、MATLAB、Origin）导入即可。相比手动记录或者跑完一次截一次图，这种方式数据可追溯、格式统一，做论文或者工程报告时写“数据可用性声明”都要方便得多。

4.3 效率优化心得

脚本自动化建模很长时间之后，我才真正意识到“效率优化”的本质不是一味缩短单次 FDTD 仿真时间，而是减少人机交互的开销。很多新手会花大把时间手动监视器调位置、等每次仿真结束再去看结果图，结果一天下来没跑几个点。而脚本化之后，人只需要在开始前把参数范围和输出路径定义好，剩下交给计算机跑，自己的时间可以拿去做物理分析或者准备下一轮的模型。

还有一个很实在的小技巧：在跑大批量扫描之前，先用一个粗网格、较小的模拟区做快速试跑，验证脚本和数据提取代码没有 bug。等确认逻辑没问题了，再恢复高精度网格和完整模拟区正式跑。这样一次就能省下好几个小时的无效等待。这个“小规模试跑”的脚本完全可以复用，只要加一个mesh_accuracy和模拟区大小的开关变量就行。

比如在我的工作流里，脚本开头常会有这种结构：

quick_test = 1; # 1 表示快速验证，0 表示正式仿真 if (quick_test) { set("x span", 0.7); set("y span", 0.7); set("mesh accuracy", 2); } else { set("x span", 1.4); set("y span", 1.4); set("mesh accuracy", 4); }

这样同样的脚本既能快速验证，又能正式跑，省去维护两套脚本的麻烦。

5. 常见问题与排查技巧实录

5.1 脚本报错排查速查表

脚本报错是家常便饭，关键是快速定位。我把这些年遇到过的高频报错整理成一个速查表，方便大家对照。

这张表里，我发现最坑的是“属性名拼写错误”，因为它不一定报错，只是设置了等于没设置。比如你想设置x span，结果输入成了xspan（少了空格），Lumerical 可能不会直接报错，但也不会生效。这时候你选中对象看属性面板，会惊讶地发现尺寸没变。所以我在脚本里凡是设置完关键属性，都会加一句查询：

select("Si_pillar_0_0"); ?get("x span");

?命令可以在脚本控制台直接输出表达式的结果，多写几行就能精确定位是哪一步的问题。

5.2 单位、精度与几何重叠问题

单位问题我在前面提过，这里再补充一个和网格相关的精度问题。FDTD 仿真是基于网格的，结构边缘会被网格阶梯化处理（staircasing）。如果你的结构尺寸和网格尺寸不匹配，仿真结果会有不小的误差。脚本建模时，可以通过设置网格精度和添加网格覆盖区域（mesh override region）来控制局部网格大小。比如对纳米柱阵列，通常会在结构所在区域加一个细网格覆盖，网格步长设置为最小特征尺寸的 1/10 到 1/20。

几何重叠的问题则是另一个经常被忽视的坑。脚本创建结构时，如果你不小心让两个不同的结构在空间上重叠，FDTD 会默认后创建的对象优先显示，但实际计算时两个物体都会参与材料平均或折射率混合（取决于网格材料算法）。这可能导致材料的等效折射率出乎意料。检查几何重叠可以这样做：

# 选中两个物体，检查重叠 select("Si_pillar_0_0"); select("SiO2_substrate"); ?overlap("Si_pillar_0_0", "SiO2_substrate");

不过 Lumerical 里的重叠检查更多依赖 GUI 中的“check overlap”按钮，脚本命令不一定在所有版本都支持。我个人常用的方式是目视检查加上显式控制 z 坐标。比如在本案例中，我设置硅柱底部z_si_bottom = z_glass_top，就让接触面严格重合（而不是交叉重叠），这样既避免了材料混合带来的意外，也让模型更符合真实加工情况。

5.3 性能和稳定性问题

最后聊一下脚本跑大批量仿真时的稳定性和效率问题。最常见的一个坑是内存泄漏或者长时间跑批导致软件卡死。我的做法是每跑完一个仿真点，主动清理不需要的数据：

clear; # 清理内存中的临时数据

但要注意clear也会把结构体定义和对象清掉，所以如果有需要保留的对象设置，可以用select先保存到变量或者直接不要频繁clear。更稳的做法的脚本里定期重启 FDTD。不过这通常涉及到外部调度，一般用户不需要做得这么极端。

还有一个大坑是循环跑仿真时不要用run;而是用run;前面的保存步骤。我见过有人每跑一次扫描点就手动点一下 run，结果一个下午就在点鼠标中过去了。正确做法是脚本里循环执行run;然后立刻读结果，读完后自动进入下一个参数点，全程不需要人工介入。我在大量实践中发现，如果中途需要人工干预，99% 是脚本没有写好——不是边界条件没设对，就是结果提取部分有 bug。

稳定性方面，一个很好的习惯是“每步写日志”。我的脚本末尾通常会加几行：

logfile = "scan_log.txt"; write(logfile, "height=" + num2str(h_si_scan) + " completed at " + time());

哪怕是简单的字符记录，也能帮你在几百个仿真点里快速定位哪一步出了问题。不要小看这个习惯，我至少两次靠这种日志文件，在两三千行的批量仿真里准确找到了异常点的参数值，省去了无数重新仿真的时间。

最后再分享一个小技巧

在我做 FDTD 建模的这些年里，脚本带给我的帮助远不止是“省时间”。它真正的价值是让建模过程变得可重复、可共享、可追溯。以前我还在 GUI 里手动添加长方体的时候，总担心某个参数忘了改导致整轮结果作废。自从切换到“参数变量 + 循环 + 扫描任务”的模式之后，这种担心的比例降低了很多，因为所有参数都在脚本顶部一目了然，改一个变量名，整个模型和所有输出数据都会跟着动，想错都难。如果你刚开始接触 FDTD 脚本，我建议你从今天这个案例入手，先跑通一个最简单的参数扫描，感受一下“一键改参、自动出图”的流程。上手之后，你会觉得以前在 GUI 里一个一个拖结构的做法，真的就像是回到了手工时代。