基于纳米微粒激发平面波的米氏散射FDTD仿真模拟与验证

基于纳米微粒激发平面波的米氏散射FDTD仿真模拟 微纳光学 【案例内容】 本案例展示了一个基于纳米粒子激发平面波的米氏散射仿真模型（mie scattering），计算其散射和吸收截面、局域场增强和远场散射分布，同时将截面和远场结果与解析解进行了比较，验证了仿真模拟的准确性。 演示： 1、电磁场与纳米粒子的相互作用在粒子表面产生强烈的局域场增强 2、吸收和散射截面 3、远场角散射 【案例文件】案例源文件

米氏散射这玩意儿在微纳光学里就像台球桌上的开球——看起来简单，实际藏着玄机。咱们今天直接上手FDTD仿真，拿金纳米球当靶子，看看电磁波和纳米粒子怎么搞出各种花样。

先整点硬核的，模型设置直接用Python脚本开搞：

from fdtd import FDTD sim = FDTD() sim.set_grid(500, 500, 20) # 纳米级网格精度 gold = sim.add_material(eps_real= -12.5, eps_imag=1.2) # 金在可见光波段 sphere = sim.add_sphere(radius=80, material=gold) # 80nm金球 sim.add_plane_wave(wavelength=600) # 600nm平面波入射

这段代码里藏着几个雷：网格精度不够会直接让局域场增强算劈叉，材料的复介电常数必须用实验值。最坑的是粒子尺寸和波长要满足米氏散射条件，这里取的80nm金球在绿光波段刚好能激发局域表面等离激元。

跑完仿真后场分布直接闪瞎眼——粒子表面电场增强超过30倍！这货在传感器应用里简直是个聚光灯。截取截面数据时得用近场积分：

scattering = sim.calculate_scattering_cross_section() absorption = sim.get_absorption() total = scattering + absorption # 消光截面

这里有个骚操作：FDTD里算散射截面其实是在总场边界上做坡印廷矢量积分。新手常犯的错是忘了排除入射场分量，结果算出来的散射截面自带双倍快乐（误差）。

远场角分布才是重头戏。用近远场变换时得注意采样点数：

theta = linspace(0, 2*pi, 180); phi = 0; [E_theta, E_phi] = farfield_transform(r, theta, phi);

这里180个采样点少了会丢失角分布细节，多了又算到地老天荒。对比解析解时得用归一化处理，把仿真结果和米氏公式曲线叠在一起，误差超过5%就得回去检查材料参数是不是手滑输错了。

最后来个灵魂暴击——为什么仿真和解析解在背向散射（θ=180°）总会差那么一丢丢？其实这是网格离散化的锅。把Yee网格改非均匀网格，或者在关键区域加密网格，立马能改善。不过跑仿真时间得翻倍，这波是精度和算力的经典博弈。

搞完这波操作的最大收获：仿真就像做菜，火候（网格）、食材（材料参数）、刀工（算法）哪个不到位都会翻车。下次试试银纳米棒，说不定能解锁更多姿势（比如Fano共振啥的）。