简介《FDTD Solutions软件教程》面向微纳光学、超表面与光电器件仿真方向的学习者和科研人员围绕Lumerical时域有限差分法FDTD的建模、求解与分析展开适合希望系统掌握软件操作并理解仿真物理设置的中级读者。资源为1个PDF文件压缩包约11.61MB内容按新功能、求解物理问题的方法、单位和规范化、计算机辅助设计模拟编辑、模拟实体对象、材料数据库、模拟计算和分析、脚本语言等模块组织便于按主题检索。教程重点覆盖8.6版用户可定义材料模型、非对角各向异性介质、模式扩展监视器、可旋转模式光源、场分析工具及TFSF光源任意角度导入等特性并对脚本语言中的系统控制、变量操作、函数、循环条件、绘图、近远场投影和光栅投影作了分节说明。已有10743人学习可作为微纳光学器件与超表面仿真入门至进阶的参考手册也便于查漏补缺。1. 从一次超表面仿真跑崩说起这套 8.6 教程解决的是什么问题去年帮人看一个超表面的案子结构画完、光源摆好、点下 run二十多分钟后场图一片糊透射率算出来大于 1。问题不在模型在没人提醒他监视器要做归一化、金属附近网格要加密、PML 要完整包住结构。这类翻车在微纳光学里太常见了——时域有限差分法本身不玄学玄学的是参数。FDTD Solutions 干的事很直白把麦克斯韦旋度方程离散到 Yee 网格上做时间步进光学器件、超表面、波导耦合器、光栅都能建。这套 8.6 版教程的价值不在教你点按钮而在把求解方法选型、网格预算、材料拟合、脚本控制四件事串成一条能复现的链路。刚接手仿真任务的研究生能照着走已经会用但总在收敛性上打转的熟手也能从里面翻到网格与归一化的边界条件。2. 求解器怎么选FDTD、本征模与 2.5D 传播算法的分工动手之前先想清楚用哪个引擎这一步选错后面网格调得再细也白搭。这套文档第二部分把五种求解路径摊开讲了时域有限差分法、本征模计算、2.5D 传播、特征值扩展传播以及光子线路级的时域/频域求解和半导体器件级的漂移-扩散求解。前三种是微纳光学日常用得最多的理解它们的假设边界比记住菜单在哪重要得多。2.1 四种方法的适用边界与代价对照FDTD 是直接对时间和空间同时离散宽带响应一次跑完靠傅里叶变换拿到频域结果复坡印亭矢量、反射透射、远场投影都能派生出来。代价是三维仿真对内存和时间都很贪网格点数按波长三次方涨。本征模计算解决的是另一个问题它不跑时间直接把波导横截面上的麦克斯韦方程组装成矩阵特征值问题用稀疏矩阵技术解出模场分布和有效折射率。器件沿传播方向不变、你只关心支持哪些模式时它比 FDTD 快一到两个数量级。2.5D 传播算法是把三维几何压成一组二维有效指标再跑 2D FDTD适合平板波导体系的长距离传播几百微米量级的环形谐振器、光子晶体腔用它比全 3D 划算。方法核心假设典型场景主要代价3D FDTD无近似全波超表面单元、散射体、任意三维结构内存与时间随网格三次方增长2D FDTD结构在第三维无限延伸光栅、平板波导截面、快速验证忽略纵向有限尺寸带来的损耗本征模求解横截面不变求模式波导模式、有效折射率、模式重叠不解瞬态与非线性2.5D 传播支撑板模之间弱耦合长距离平板波导器件强耦合结构误差放大选型上我一般这么走先问器件在传播方向上是不是基本不变是就优先本征模或传播算法再问是不是平板体系、长度超过几十微米是就压成 2.5D只有当结构三维特征明显、或者涉及非线性与超短脉冲时才老老实实上 3D FDTD。反过来拿 3D FDTD 去算一根几十微米长的锥形波导纯属给自己找罪受。2.2 把三维压成二维2.5D 传播算法的四步操作文档里 2.3.1 节把流程写得很清楚四步标定核心波导的垂直平板模、压扁结构算等效二维指标、跑 2D FDTD、需要时把场重新扩展到三维。中间那步的等效指标计算有两条路——变分法和基于互惠理论的算法输出都是带色散的有效材料色散同时来自原始材料属性和波导几何。这里最容易翻车的地方是生成的有效折射率可能带上原始材料没有的物理特性比如人为的负虚部折射率所以脚本里要给指标范围加钳位。# 设定波长范围标定垂直方向平板模 lambda_min 1.5e-6; lambda_max 1.6e-6; setanalysis(wavelength, linspace(lambda_min, lambda_max, 11)); # 指定传播算法使用的等效指标计算方法 # 1 变分法2 互惠理论 set(effective index method, 1); # 对生成的有效指标做上下限钳位防止非物理的负虚部 set(min index, 1.0); set(max index, 3.6); # 压扁三维结构生成二维等效材料 switchtolayout; set(dimension, 2); run;这段脚本的关键在三个参数波长采样点数决定了等效材料色散拟合的精度采样太稀会让拟合出来的多系数模型在带外发散等效指标方法的选择影响强限制波导的精度变分法在折射率对比大的体系里通常更稳上下限钳位是保命的尤其当原始材料本身带增益或强吸收时压扁过程很容易算出虚部符号反掉的指标跑出来的场会莫名其妙地增长。压扁完成之后我会先看一遍有效折射率表确认基模的实部随波长单调、虚部数量级合理再进 2D 计算。跳过这一步后面所有结果都不值得信。2.3 网格预算共形网格能救回多少精度FDTD 的网格是矩形笛卡尔的材料属性、几何信息、电磁场都在网格顶点上算。网格越细越准成本和内存跟着涨。软件默认会自动生成网格并且对高折射率材料和高吸收材料自动加密——前者是为了保证单位波长内的网格点数固定后者是为了解析穿透深度。问题在于自动网格判断不了你的结构里哪块区域场变化最剧烈比如金属尖角、纳米间隙、表面等离激元的界面。共形网格是 7.0 版本就引入的能力作用是在不整体加密的前提下让曲面和斜面的台阶误差降下来。对薄层、曲面、高折射率波导里的硅这类结构共形网格几乎是必开项否则你只能靠全局加密去换精度成本翻几倍。我的习惯是全局网格精度先给个中等档然后在关键区域叠一个网格覆盖区只在那块加密。# 全局网格精度1 最粗 5 最细日常从 2 或 3 起步 setglobalmonitor(mesh accuracy, 2); # 打开共形网格处理曲面与斜壁 set(conformal variant 1, 1); # 在金属间隙处叠加局部网格覆盖区 addmesh; set(x, 0); set(x span, 200e-9); set(y, 0); set(y span, 200e-9); set(z, 0); set(z span, 200e-9); set(override x mesh, 1); set(dx, 2e-9); set(override y mesh, 1); set(dy, 2e-9); set(override z mesh, 1); set(dz, 2e-9);mesh accuracy是个打包参数改动它会连带改最小网格步长和网格点数上限所以别在扫描循环里动态调它否则不同扫描点的离散误差不一致结果曲线会出现假的抖动。conformal variant有两个档变体 1 更保守变体 2 在强弯曲结构上精度更高但更吃内存超表面 pillar 这类结构我通常先用 1收敛性对不上再换 2 做交叉验证。局部覆盖区里的dx/dy/dz是硬步长写死之后不受全局精度影响这才是控制关键区域误差的正确姿势。边界条件这边PML 从 7.0 起可以直接指定目标反射率而不用去猜层数——层数多了拖慢仿真少了反射压不住。常见的做法是设一个 1e-5 量级的反射目标让软件自己决定层数。另外文档里反复强调的一点PML 区域要完整覆盖结构周边默认行为会自动把材料属性延伸到 PML 里这对吸收是好事但如果你手动改过结构范围得回头确认一遍延伸是否还成立。周期和 Bloch 边界则用于单元结构超表面用周期边界时源也必须是对应的平面波或 TFSF混用会导致归一化错位。3. 脚本语言实战从变量、循环到 run 与远场投影图形界面能干的活脚本都能干反过来不成立。这套文档第八部分把脚本拆成十六节从系统变量、运算符、函数、循环条件一直到添加对象、运行仿真、量度归一化、近远场投影、材料数据库、自定义 GUI 和自定义命令。真正决定效率的是前三节之外的东西怎么组织参数、怎么把扫描跑起来、怎么从结果里把数捞出来。3.1 脚本骨架先建结构再挂监视器最后才 run一段能复用的脚本顺序永远是先清空、再定义参数、再建对象、再挂监视器、最后执行。参数放最上面是为了后面改成扫描时只动一处。结构用addrect、addcircle这类基元拼复杂组群用addgroup打包组内对象用groupscope切换作用域——这一条新手最容易忽略导致set改错对象。# 清空当前工程保证脚本可重复执行 newproject; # 参数区所有可变量集中在这里 L 400e-9; # 单元周期 H 600e-9; # 柱高 R 90e-9; # 柱半径 n_mat 3.48; # 柱材料折射率 # 仿真区 addfdtd; set(dimension, 3); set(x, 0); set(x span, L); set(y, 0); set(y span, L); set(z, 0); set(z span, 2e-6); # 周期边界单元结构仿真用 set(x min bc, periodic); set(y min bc, periodic); # 衬底与柱子 addrect; set(name, substrate); set(z max, 0); set(z min, -500e-9); set(material, SiO2 (Glass) - Palik); addcircle; set(name, pillar); set(x, 0); set(y, 0); set(z min, 0); set(z max, H); set(radius, R); set(material, Si (Silicon) - Palik); # 监视器反射与透射功率 addpower; set(name, R); set(monitor type, 2D Z-normal); set(z, -600e-9); addpower; set(name, T); set(monitor type, 2D Z-normal); set(z, 900e-9); run;参数区用变量而不是直接写死数字是为了后面接扫描。addfdtd建的是求解区域dimension设 3 就是三维周期边界只作用于 x、y 两个方向z 方向仍然是 PML这是超表面单元仿真的标准配置。两个功率监视器一个放衬底下方、一个放柱子上方注意都要留在 PML 之外否则监视器会落在吸收层里读数直接失真。run是同步执行脚本会阻塞到计算结束如果要做批量任务后面换成异步的作业接口更合适。3.2 参数扫描与优化把 for 循环换成 runsweep早期的参数扫描得自己写 for 循环每次改参数、run、存结果。7.0 之后图形界面能直接做扫描和基于粒子群的优化7.5 又补了并行能力8.6 加了良率分析。脚本层面runsweep、addjob、runjobs、clearjobs这几个命令是把任务丢进队列的关键getsweepresult、getsweepdata负责把扫描结果捞回来。# 定义扫描变量柱半径从 70nm 到 110nm取 9 个点 sweep_name radius_sweep; r_vals linspace(70e-9, 110e-9, 9); # 建立扫描任务指定要采集的结果 runsweep(sweep_name); # 逐个扫描点修改参数并入队 for (i 1:length(r_vals)) { switchtolayout; select(pillar); set(radius, r_vals(i)); addjob(sweep_name); } # 提交执行多核或许可足够时会自动并行 runjobs; # 回收透射率随半径的变化 T_data getsweepresult(sweep_name, T); T_curve T_data.T;runsweep建的是一个任务容器本身不计算addjob每调一次就往里塞一个待算配置runjobs才真正分发。这里有个容易被忽略的点switchtolayout必须写在每次改参数之前否则你改的是上一次仿真完成后的状态参数根本没进到排队的任务里最后九个点算出来一模一样排查半天以为是物理问题。并行分发依赖引擎许可数量文档里提到 7.5 起可以在局域网内多机并行三台机器就是三倍吞吐但每个并行实例都要占许可许可不够时会退化成串行队列看着卡住其实是在排队。getsweepresult返回的是数据集对象取字段用点号别用索引硬取不同监视器类型的字段名不一样。优化这边图形界面提供的是粒子群算法也可以自己写优化流程。我的经验是参数少于三个、且你大致知道最优区间时网格扫描比优化更可控参数多、目标函数多峰时再上优化并且要给参数加边界否则粒子群很容易跑到非物理区域比如半径大于周期结构直接重叠。3.3 近场到远场投影与光栅投影监视器类型决定投影命令超表面和光栅最常要的输出不是近场是远场方向图、衍射级次效率、或者到某个模式的耦合系数。文档 8.12 和 8.13 分别讲近远场投影和光栅投影路径不一样前者用普通功率或场监视器加远场分析后者需要周期结构专用的光栅投影直接给出各级衍射效率。# 在结构上方加一个二维场监视器用于远场投影 addprofile; set(name, near_field); set(monitor type, 2D Z-normal); set(z, 800e-9); set(x span, L); set(y span, L); # 远场投影指定投影面与频率点 farfield_data farfield(near_field); theta farfield_data.theta; E_theta farfield_data.Etheta; # 周期结构的光栅投影直接取衍射级次 grating_data gratingff(near_field); n_orders grating_data.n; eff grating_data.T;farfield做的是近场到远场的变换输入监视器必须完整覆盖一个投影面且要在结构近场之外、PML 之内位置太靠近结构会截断倏逝波太靠近 PML 会把吸收层的衰减算进去。gratingff是针对周期结构的它依赖周期边界和单元尺寸所以你改了L就必须重算不能拿旧结果套新周期。衍射级次效率的归一化基准是入射功率如果源用了 TFSF归一化会自动扣掉散射场这点和普通平面波源不一样混用会出现效率总和大于 1 的情况。3.4 自定义 GUI 与自定义命令把重复劳动封成按钮文档 8.15 和 8.16 讲的是用户自定义图形接口和自己的脚本命令这两节很多人跳过其实对团队协作价值最大。把一套固定的建模仿真流程封成带输入框的小面板新同事不用理解网格参数也能跑出一致的结果把常用的后处理写成自定义命令比如取某个监视器的峰值透射率并返回波长可以省掉每次手写索引。常见做法是把自定义命令写在一个单独的脚本文件里通过脚本编辑器加载然后在主脚本里像内置命令一样调用。需要注意的是作用域自定义命令内部如果用了set而没有明确select它会作用到当前选中的对象上行为变得不可预测。稳妥的写法是在命令内部先groupscope或显式select把上下文锁死。4. 材料与光源设置自定义模型、各向异性张量与斜入射 TFSF材料这一块是 8.6 版本改动最大的地方。用户可定义的材料模型允许直接改更新方程非线性材料、负折射率材料、任意形式的场相关介电常数都能挂进去材料数据库里补了顺磁性材料、拉曼-可尔模型和四能级二电子激光模型硅的拉曼效应、孤子传播、激光动力学这类问题有了现成的起点。非对角各向异性介质支持完整的 9 元素介电常数张量磁光效应和液晶分子轴向可以建模了。4.1 材料拟合先建 FDTD 模型再用材料浏览器检查材料数据库里给的是采样数据FDTD 时域求解必须把它变成可迭代的多系数模型这个转换叫拟合。7.0 的材料分析模型允许输入介电常数或反射率实虚部的公式但文档特别强调了一句仿真过程并不直接使用这个公式必须先生成 FDTD 材料模型。这句话翻成操作就是——你写完公式不算完得点生成然后用材料浏览器把拟合曲线和原始数据叠在一起看。# 新建材料并导入采样数据 mat addmaterial(MyDispersion); setmaterial(mat, sampled data, 1); importdataset(mat, nk_data.txt); # 生成多系数时域拟合模型 setmaterial(mat, fit, 1); setmaterial(mat, max coefficients, 6); setmaterial(mat, tolerance, 1e-3); # 强制在材料浏览器中检查拟合误差 materialanalysis(mat);max coefficients是拟合用的系数上限数值越大拟合越准但时域更新方程越重内存和时间都跟着涨tolerance是拟合残差目标低损耗色散材料在 8.6 里拟合效果有改进但如果你把容差压到 1e-5 以下很可能拟合不收敛反而生成一个在带外疯狂振荡的模型。materialanalysis是必须走的一步看两个东西拟合区间是否覆盖了你实际仿真的波长范围以及虚部在带内是否为正——虚部符号错了材料会变成增益介质场会指数增长。我见过最隐蔽的一次翻车就是拟合区间只覆盖了 1.5 到 1.6 微米但扫描范围设到了 1.7 微米边界外的结果全是噪声。4.2 九元素介电常数张量先给对角再靠 Grid Attribute 旋转各向异性介质的本构关系里介电常数是个 3×3 张量。文档给的操作顺序是先指定张量的对角元素再通过格子特性实体对象对张量做任意旋转。这个顺序不能反因为旋转是作用在已定义的张量上的你先把角度设了再填对角元素旋转结果会被覆盖。# 定义各向异性材料先写对角元素 mat_aniso addmaterial(AnisoMedium); setmaterial(mat_aniso, anisotropy, 1); setmaterial(mat_aniso, epsilon xx, 2.1); setmaterial(mat_aniso, epsilon yy, 2.1); setmaterial(mat_aniso, epsilon zz, 2.5); # 添加格子特性对象指定旋转角 addgridattribute; set(name, rotate_attr); set(attribute type, grid attribute); set(rotation 1, 30); # 绕第一个轴的旋转角单位度 set(rotation 2, 0); set(rotation 3, 0); # 把旋转属性绑定到目标结构 select(waveguide_core); set(grid attribute name, rotate_attr);对角元素决定主轴上的三个介电常数旋转则由格子特性对象在网格层面完成所以它不像几何旋转那样需要重建结构。三个rotation参数的顺序是按轴依次施加的不是欧拉角的任意约定文档里给了明确说明照抄顺序就行。这里最常见的错误是把旋转属性挂到了错误的组上——格子特性对象是按空间位置生效的如果你在一个组里放了多个结构属性会作用到所有落在该网格区域的对象上结果就是不该转的地方也转了。液晶和磁光材料的仿真里我一般给每个需要独立取向的区域单独建一个格子特性对象宁可多建几个也别指望一个对象管全部。4.3 任意角度 TFSF两个前提没满足结果一定错8.5 起支持任意角度入射的 TFSF 光源这对斜入射超表面和光栅的仿真意义很大。但文档里给的前提条件很硬要保证波导或光纤沿 PML 边界扩展出去同时要关闭高级选项里的extend structure through pml。这两条其实是矛盾的——默认行为是自动把材料属性延伸到 PML 里对普通结构是好事对 TFSF 斜入射就会让散射场和总场边界打架。操作上先在光源的 General 标签里设入射角再进 FDTD 求解区域的高级设置取消结构延伸选项然后手动把衬底或波导延伸到 PML 外侧。判断有没有做对看两个地方TFSF 边界上的场是否连续以及反射监视器在无结构时是否读到一个接近零的本底。本底压不下去说明边界附近的材料分布破坏了总场/散射场的分离这时候调角度、调网格都没用得回去改结构延伸。角度越大需要的网格也越细。斜入射时场在网格上的相位变化更快单位波长网格点数不够会出现数值色散表现为谐振峰偏移。经验值是入射角超过 30 度后把关键区域的网格步长再压三分之一然后做一次角度收敛性对照。5. 避坑排查五条踩过的记录每条都有现象、原因和解法仿真跑不出来八成不是算法问题是设置问题。下面五条是我和周围人反复踩到的按现象、原因、解决三段写遇到同类症状可以直接对号入座。5.1 透射率算出大于 1PML 没完整包住结构现象很直观无结构时反射监视器本底不接近零加結構后透射率加反射率超过 100%场图在边界附近有明显条纹。原因是 PML 吸收层没有完整覆盖结构的横向范围或者结构被手动截断在 PML 内侧导致部分散射光被边界反射回来和入射场叠加。文档里提过默认行为会自动延伸材料属性到 PML但只要你手动改过对象范围这个延伸就可能失效。解决方法是先跑一次空跑确认反射本底在 1e-5 量级以下再把结构范围放大到超出 PML 内边界最后才看正式结果。这一步我现在的习惯是写进脚本里空跑不过就不允许继续。5.2 谐振峰随网格跳变网格粗加上材料拟合没过关现象是同一条透射谱换个网格精度峰位偏移几十纳米或者扫描到某个波长点结果突然发散。原因通常有两个叠在一起结构最小特征尺寸处的网格步长不够以及色散材料的拟合区间没覆盖扫描范围带外拟合发散。排查顺序是先看材料浏览器里的拟合残差和虚部符号再在关键区域加网格覆盖区做收敛性对照把步长减半看峰位是否稳定。别一上来就全局加密那样只是把问题拖慢不会消失。5.3 扫描结果全一样switchtolayout 漏写或参数没进任务现象是九次扫描的输出曲线完全重合或者最后一组参数覆盖了所有结果。原因是addjob之前没有switchtolayout改参数改在了仿真完成后的状态上或者扫描变量名和getsweepresult里的名字不一致。解决是把切回布局、改参数、入队这三步固定成一个顺序中间不插别的操作并且在取结果前先打印一次扫描点的参数列表确认每个点确实不同。这个错犯一次会浪费一整轮计算时间非常亏。5.4 并行任务排队不动许可数量和任务粒度没算清现象是runjobs之后进度条长期停在第一个任务或者多机环境下只有本机在算。原因是并行实例数量受引擎许可限制许可不够时任务会退化成串行排队看着像卡死另一种情况是单个任务太大调度器分不出空闲槽位。解决是先确认可用许可数把扫描拆成不超过许可数的批次提交单批任务规模控制在单机能吃下的内存以内。另外 7.5 起支持并行时的动态监视器但并行环境下监视器刷新会拖慢吞吐长任务建议关掉实时刷新。5.5 脚本改错对象groupscope 和单位制这两处最容易出问题现象是明明改了某根波导的宽度结果另一根变了或者几何尺寸差了一千倍。原因是set作用在当前选中对象上组内对象没切换作用域时命令会落到组或上一个选中对象单位方面除特殊说明外全部按国际单位写 400 却想表达 400 纳米结构直接变成 400 米。解决是在每次select之后紧跟一次get打印关键参数做回读校验几何量统一在脚本开头用带单位的变量定义杜绝裸数字。回读校验多花两行代码能省掉几小时的排查。6. 结果验证把良率分析和模式扩展监视器真正用起来收敛性对照只能说明数值稳了说明不了设计稳。8.6 在优化和扫描窗口里加的良率分析工具本质是扩展的蒙特卡罗分析给加工误差范围内的多个参数各设一个分布跑多参数扫描统计整条透射谱或某个性能指标的变化。做法上先把关键尺寸参数化给每个参数设正态分布和容差比如柱半径 ±5 纳米、高度 ±10 纳米再指定蒙特卡罗样本数。样本数不是越多越好二十到五十次通常能看出统计趋势上百次就是拿时间换小数点后第三位。# 良率分析多参数蒙特卡罗 yield_setup yield_run; addsweep(yield_setup, montecarlo); # 参数分布均值 ± 标准差 setparameter(yield_setup, radius, R, 5e-9, normal); setparameter(yield_setup, height, H, 10e-9, normal); # 样本数与目标指标 set(number of trials, 30); set(figure of merit, T_peak); runjobs; yield_result getsweepresult(yield_setup, yield);setparameter里的标准差要和实际工艺能力对齐拍脑袋给一个很小的值良率当然接近 100%但没有参考意义。figure of merit要选一个能代表设计意图的标量比如峰值透射率或谐振波长偏移量选错指标统计出来的是噪声。模式扩展监视器是另一个值得单独拎出来的验证手段。它把监视器上的任意场分布投影到指定的一组模式上返回的是重叠积分和 S 参数这些 S 参数可以直接喂给线路级仿真做系统联调。验证波导耦合器时我用它同时看基模和二阶模的功率占比如果高阶模占比随长度剧烈振荡说明结构里有寄生反射得回去查网格或边界。它的一个前提是投影用的模式必须和监视器所在截面的材料分布一致截面跨了不同材料层还硬套单模重叠系数会失真。验证项手段通过标准数值收敛网格步长减半对照峰位偏移小于线宽 1/10边界吸收空跑反射本底本底低于 1e-5材料模型材料浏览器拟合残差带内虚部符号正确、残差平滑工艺容差蒙特卡罗良率分析目标指标分布集中在设计窗口内模式纯净度模式扩展监视器目标模占比高且随长度稳定从那以后我每次交付仿真结果前都强制走一遍空跑、网格对照、材料回读这三步一样不过就不往下走。这套教程真正值钱的地方也在这——它把 8.6 的材料模型、张量、TFSF、脚本接口和 7.x 时代的并行、共形网格、参数扫描放在同一条流程里讲你照着搭一遍坑基本都能提前撞上。希望帮到你。本文还有配套的精品资源点击获取
