FDTD Solutions新手必看：从零开始搭建你的第一个纳米光学仿真模型（附完整脚本）

FDTD Solutions新手必看：从零开始搭建你的第一个纳米光学仿真模型（附完整脚本）

第一次打开FDTD Solutions时，面对密密麻麻的工具栏和参数面板，大多数新手都会感到无从下手。作为计算电磁学领域的黄金标准工具，它确实需要一定的学习曲线——但别担心，我们今天要做的不是系统学习所有功能，而是用20分钟完成一个能立即看到光学现象的"最小可行案例"。想象一下，当你第一次看到自己设计的纳米结构周围出现绚丽的电磁场分布时，那种成就感会瞬间点燃继续探索的热情。

我们将以金纳米球的光散射仿真为例，这个案例包含了材料定义、几何建模、光源设置和结果分析的全流程。选择它有三个原因：首先，金属纳米颗粒是光子学研究的经典模型；其次，其对称结构能简化网格划分；最重要的是，你能直观看到著名的局域表面等离激元共振现象。以下是完成本案例前需要准备的：

• 已安装FDTD Solutions 2023R2或更新版本
• 8GB以上内存的计算机（集成显卡即可）
• 约500MB的可用磁盘空间

1. 软件初始化与项目设置

启动软件后，首先在File菜单中选择New Project，建议命名为Nanosphere_Scattering.fsp。此时你会看到三个主要窗口：左侧的Object Tree（对象树）、中央的Visualizer（可视化区）和右侧的Property Editor（属性编辑器）。对于我们的纳米球仿真，需要先做两项基础配置：

# 设置仿真区域大小 set("x span", 1e-6); # X方向1微米 set("y span", 1e-6); # Y方向1微米 set("z span", 1e-6); # Z方向1微米

在Simulation标签页下找到Mesh Settings，将网格类型设为"auto non-uniform"。这种智能网格会在场强变化剧烈的区域自动加密，对于纳米结构特别重要。接着配置边界条件：

提示：PML(完美匹配层)能有效吸收 outgoing 波，防止虚假反射。对于300-800nm的光学仿真，默认的8层PML完全够用。

2. 构建金纳米球模型

在Object Tree中右键点击Structures，选择Add Sphere。这时你会看到一个红色球体出现在可视化窗口中央。我们需要调整两个关键参数：

1. 半径设为50nm（在Property Editor中修改radius值）
2. 材料选择"Au (Gold) - Johnson & Christy"

# 通过脚本添加纳米球 addsphere( "name", "gold_nanosphere", "radius", 50e-9, "material", "Au (Gold) - Johnson & Christy", "x", 0, "y", 0, "z", 0 );

金属材料的色散特性对仿真结果影响巨大。FDTD Solutions内置了多种金材料的实验数据模型，其中"Johnson & Christy"在可见光波段最为准确。你可以在Materials标签页查看其折射率随波长的变化曲线——注意在约520nm处会出现明显的虚部峰值，这正是等离激元共振的特征。

3. 光源与监视器配置

点击Sources标签添加平面波光源，关键参数设置如下：

• 光源类型：Plane Wave
• 注入方向：Z轴负方向
• 波长范围：400-700nm（覆盖可见光波段）
• 偏振方向：X轴（电场振动方向）

addplane( "name", "source", "injection axis", "z", "direction", "backward", "wavelength start", 400e-9, "wavelength stop", 700e-9, "polarization angle", 0 );

接下来添加两个关键监视器：

1. 频率域场分布监视器（在XY平面截取）
   • 位置：Z=0（纳米球赤道平面）
   • 记录分量：|E|²（电场强度平方）
2. 散射截面监视器（包围整个纳米球）
   • 类型：Power
   • 记录模式：散射功率

注意：场监视器要放在共振可能发生的区域附近，而功率监视器需要完全包围散射体才能准确计算截面。

4. 运行仿真与结果可视化

点击Run按钮开始计算。在1μm³的仿真区域内，普通笔记本大约需要3-5分钟完成运算。当状态栏显示"Simulation completed"后，右键点击场监视器选择Visualize，你将看到：

最令人兴奋的部分来了——在520nm附近扫描波长时，你会观察到：

• 散射截面出现明显峰值
• 纳米球两侧形成强烈的电场增强
• 特征性的偶极子辐射模式

这些现象正是金纳米颗粒等离激元共振的典型表现。为了量化分析，我们可以提取散射截面数据：

# 导出散射截面数据 f = getdata("scattering","f"); lambda = 3e8./f; scattering = transmission("scattering"); plot(lambda*1e9, scattering); xlabel("Wavelength (nm)"); ylabel("Scattering cross section");

5. 参数优化与扩展思路

现在你已经完成了基础仿真，接下来可以尝试以下进阶操作：

• 尺寸效应：修改纳米球半径（20-100nm），观察共振峰位移
• 材料对比：将金换成银或铝，比较共振品质因数的差异
• 环境介质：添加背景折射率（如n=1.33模拟水中情况）
• 阵列结构：复制多个纳米球研究耦合效应

下表展示了不同尺寸金纳米球的共振峰位置变化：

当需要保存结果时，推荐使用脚本批量导出数据而非手动截图：

# 批量保存场分布图 for wl = [400e-9, 520e-9, 700e-9] select("field_monitor"); set("wavelength", wl); image("E", "component", "x"); saveimage(sprintf("Field_at_%dnm.png", wl*1e9)); end

完成第一个仿真后，建议尝试修改纳米球形状为立方体或棒状，观察模式变化。你会发现，仅仅改变几何对称性就会产生完全不同的电磁响应——这正是纳米光子学最迷人的地方。