|
|
计算平面波激发的纳米粒子的散射和吸收截面、局部场增强和远场散射分布（Mie 散射）。将截面和远场结果与解析解进行比较，以验证仿真的准确性。（联系我们获取文章附件） Z.rR) H". [&VP5Z 概述 CE;J`; &2\^S+4  Zxv Bo4>tH rIJv(&l 纳米粒子的散射特性通常用场增强、横截面和远场分布来描述。本例展示了如何从单个 FDTD 仿真中获得这些结果。 yHIZpU|( j |fhYft 运行和结果 W34_@,GD I5mtr 1.打开仿真文件，然后单击“运行”按钮。 h}SP` 2.可以通过右键单击监视器或分析组并选择感兴趣的参量来手动浏览结果。 DL]\dD 3.关联的脚本文件可用于绘制如下所示的代表性结果。 0P7sMCY u 本地字段增强 Tw+V$:$$ 电磁场与纳米粒子的相互作用可以在粒子表面产生强烈的场增强。频域场监测仪直接测量局部场增强。下图显示|E|的平方在XY，XZ和YZ平面中，在最接近脚本中指定的“目标波长”的波长点穿过粒子中心。 -}$mv V~! lY\ 可以注意到，TFSF 源的边缘在图中可见，因为图像颜色的突然变化。源内的字段是“总计”字段（即事件字段 + 分散字段），而只有“分散”字段在源外部可见。 pf[m"t6G~ (N
0kTi]b  $ l+DkR+ ! Pf6UNN' 吸收和散射截面 -gK*&n~ 吸收截面（总吸收功率除以入射光束每单位面积的功率）由位于 TFSF 源内的分析组计算。分析组测量流入颗粒的净功率，并通过将其归一化为源强度，返回吸收截面。同样，散射截面由位于 TFSF 源外部的分析组计算。 Iq["(!7E5 @r=v*hu 根据定义，横截面以m的平方用于 3D 模拟和m用于 2D 模拟。 < 2,NWn . |Ta-D++]' 横截面测量通常被标准化为散射物体的大小，如下图所示。Mie 效率定义为横截面与几何面积的比， πr2对于球体（3D）和对2r于圆（2D），并且通常相对于尺寸参数 （2πn1/λ)，其中n1是 FDTD 区域的背景指数，对于空气为1。 TcKt 将 FDTD 结果与从 mie3d 脚本获得的分析溶液进行比较。两个结果之间的差异很明显，希望对模拟设置进行一些改进。这将是下一节关于收敛检验的主题。 vn(j i= fs,]%g^  c[<>e#s+;  xixdv{M<FF
远场角散射 )p!.V( , 在大多数散射实验中，散射场（辐射图）的测量相对于所考虑的波长尺度远离散射体。“scat_ff” 监视器返回远场中的散射场分布。以下极坐标图显示了 X-Y、X-Z 和 Y-Z 平面中远场中的散射场。每个图都包含两种颜色的线条：蓝色表示 FDTD 仿真结果，绿色表示 mie3ds12 脚本命令的分析结果。第一个图显示了如何在每个平面中定义极角。 T;`2t; Kd;Iu\4hv  poU1Q#+4p*
重要模型设置 jL]Y;T8 .4O~a 模型设置脚本 r~q 3nIe/, 模型对象中的设置脚本用于设置网格大小、仿真跨度和粒子位置。该脚本是确保模拟区域、网格覆盖区域、源、scat 和 abs 监视器位置正确的便捷方法。例如，TFSF 源必须位于 scat 和 abs 监视器之间，对象之间至少有两个网格单元。这些对象的位置必须通过安装脚本进行设置。其他属性（如模拟时间）可以直接在对象中修改。 2PTAIm Rq TnXx;v  vp |'Yy(9z >O$ JS, TFSF 来源 3]wV 1<K TFSF 光源是专门为这种情况而设计的，其中非周期性物体被平面波照亮。通过将散射场与入射场分开，使纳米粒子的散射分析变得简单明了。为了使散射分析正常工作，确保散射体完全在 TFSF 源内至关重要。 >_'0 s 使用 TFSF 源进行电源归一化 rJKac"{ TFSF 源的电源规范化可能会令人困惑。与其将结果归一化为源功率（对于理想平面波来说，这是无限的，因为它具有无限的范围），不如按源强度进行归一化。这导致功率测量值以横截面型单位返回。 "OWW -m “abs”和“scat”分析组 %yPjPUHy 由六个2D监视器组成，形成一个封闭的盒子，测量流入/流出盒子的净功率。这些分析组的位置非常重要。测量吸收功率的 “abs” 分析组必须完全在 TFSF 源内，但在粒子之外。“scat” 监视器必须完全位于 TFSF 源之外。 oo-O>M#5 H zc5BC 网格覆盖区域 R8bKE(*rxj 对于金属仿真，网格覆盖区域通常用于更准确地解析金属界面的位置，尤其是曲面。在此仿真中，网格覆盖区域设置得足够大，不仅包括金球，还包含整个TFSF区域。这是有意为之，因为 TFSF 源在均匀网格化区域中效果最佳。 J5 "d|i 另请注意，网格大小会影响总监视器和 scat 监视器与源的距离。最好在源和监视器之间保持至少两个网格单元间距，以避免放置在灰色阴影源注入区域中的监视器返回非物理结果。请注意，这些条件由“模型”设置脚本强制执行。 ;m#_Rj6 wmB_)`QNP 对称 ?R":"*eu 此模拟在 X 和 Z 维度上都具有对称平面。为了将仿真时间和内存减少 4 倍，将 X min 边界条件设置为对称，将 Z min 边界条件设置为反对称。请注意，只有当粒子和源都具有必要的对称性时，才能使用对称性。 ,Z$!:U 使用参数更新模型 p>9|JMk ;}KT 3Q<^ 对仿真文件进行参数化，以便更轻松地设置仿真。该模板目前使用球形粒子，但它可以与任意形状的粒子或多个粒子一起使用。在“模型”中指定参数后，其余仿真对象的大小将自动调整。 A~M .v0 · 设置源波长范围和偏振。 ?DgeKA"A · 设置纳米颗粒的材料或索引。 C/#?S=w`4 · 在“模型”中设置纳米颗粒的跨度和位置、网格覆盖的网格大小以及仿真跨度。源和 “abs”/“scat” 分析组将自动由最多两个网格单元分开，纳米颗粒被 “abs” 分析组完全包围。 X+[h]A · 在模拟非球形粒子或多个粒子时，可能需要更新边界条件以匹配新结构的对称性。还需要修改关联的脚本文件，以校正散射体的几何面积和大小参数。 7xh91EU:4 w:nLm,  KD^n7+w% 进一步推广模型 t>} (` 0 )K!!Zq3;| 基板上的颗粒 L}= t"y 此示例使用被均匀材料包围的粒子。如果颗粒在基板上，则必须修改分析的远场部分。此示例中使用的技术（从封闭的监视器盒投影）仅在所有监视器都位于向外延伸到无穷大的单一均匀材料中时才有效。当存在基板时，计算远场散射模式的最佳方法是使用一个位于粒子上方或下方的监视器（取决于散射的主要方向）。然后，您可以使用标准的 farfield3d 函数。使用单个监视器时，必须使仿真跨度足够大，以使大多数散射光在到达 PML 吸收边界之前可以通过监视器。此问题仅适用于远场分析。无需更改横截面和近场测量的分析。 zGaqYbQD ^I:f4RWo  Plj >+XRO 非偏振照明 X 4CiVV 对于具有非相干非偏振照明的系统，运行第二次仿真，将源偏振旋转 90 度，然后对结果求平均值。这可以通过对源偏振角进行 2 点参数扫描轻松实现。 noZ!j>f{@l 收敛 ~c`%k> $
Av>xgfX 使用当前设置（模拟范围为 1x1x1 um3，网格精度3，5nm网格附近粒子）仿真需要大约150 MB的内存，运行时间约为1分钟。这些设置提供了合理的精度水平，同时最大限度地减少了仿真时间。以下更改将提供更高的准确性。 r H$M6S 网格细化 * -X`^R 将网格细化设置为“共形变体 1”，以实现金颗粒边界的子单元分辨率。如果网格很粗糙，并且在目标频率下金属和周围介质之间的介电常数差异很大，则选择此设置时必须小心。最好执行一些收敛测试。 *?rO@sQy] 网孔尺寸 tSVU,m 将网格覆盖网格尺寸设置为 0.8nm ^GnR1.ux 模拟跨度 ?
J /NYV 在所有方向上将模拟跨度设置为 2um。当模拟区域太小时，共振表面等离子体模式的倏逝尾部将与 PML 边界条件相互作用。 ib%'{?Q. PML 反射 9c{T|+ ] 从 PML 边界条件反射的任何光都可能影响结果。更多的 PML 层将减少反射。但是，如果您使用默认 8 个图层的“拉伸坐标 pml”，则无需更改它，除非您需要更高的精度。 7G=Q9^J.H DGTD 求解器 3FG'A[x3O 考虑使用米氏散射 （DGTD）获得金属纳米颗粒的高精度结果。DGTD 求解器中有限元网格的性质可以实现更好的收敛，并且不易出现阶梯和热点问题。 6%Pvh- ~_ 下图显示了更高精度 FDTD 仿真的横截面。FDTD 与理论结果之间的一致性显然要好得多。此外，较小的网格会产生更高分辨率的场轮廓，从而更好地解析金属界面附近的场。 e/b
|
s l Zx Ak 
|
