基于FDTD的二维表面等离子体波导仿真系统

资 源 简 介

本项目利用二维有限差分时域法（FDTD）对表面等离子体波导结构的电磁特性进行深度模拟与分析。系统核心实现了基于Yee网格的麦克斯韦方程组离散化迭代过程，能够精确刻画光波在亚波长金属结构表面激发的表面等离子体激元（SPP）的传播特性。功能涵盖了针对频率色散材料（如金、银等贵金属）的Drude模型处理，采用辅助差分方程（ADE）技术将材料的频率相关性实时耦合进时域计算中，解决了金属材料在光频段的复杂介电常数仿真难题。该项目可以模拟金属-介质-金属（MDM）或介质-金属-介质（DMD）等多种波导拓扑结构，支持自

详 情 说 明

基于FDTD算法的二维表面等离子体波导电磁场仿真系统

项目介绍

本项目是一个基于有限差分时域法（FDTD）的数值仿真系统，专门用于研究二维空间内表面等离子体激元（SPP）的传播特性。系统聚焦于金属-介质-金属（MDM）型波导结构，通过模拟亚波长尺度下电磁波与金属表面自由电子的相互作用，揭示光场在纳米结构中的局域与传输规律。该仿真系统能够处理光频段下具有显著色散特性的金属材料，并提供完整的时域更新逻辑与频域分析功能。

功能特性

• 色散材料建模：内置Drude模型以模拟真实金属（如金）在光频段的物理响应，通过处理等离子体频率与碰撞频率，精确还原复介电常数。
• 先进边界条件：实现了非分裂卷积完美匹配层（CPML）吸收边界，有效抑制计算域边缘的数值反射，保证开放区域仿真的准确性。
• 波导结构定制：支持MDM拓扑结构设计，允许用户自定义介质通道宽度、背景介质折射率等几何与物理参数。
• 实时场量动态可视化：仿真循环中集成实时绘图功能，可同步观察磁场（Hz）与电场（Ex）的波动演化过程以及透射端的时域波形。
• 频谱分析能力：利用快速傅里叶变换（FFT）将时域监测数据转化为频域传输谱，支持对波导能量透过率进行定量分析。

使用方法

1. 在环境配置完成后，直接运行主函数。
2. 仿真开始后，程序会弹出图形窗口，左侧展示Hz和Ex的空间场分布，下方展示透射监测点随时间变化的振幅。
3. 系统将按照设定的1500个时间步长进行迭代更新。
4. 迭代完成后，程序会自动执行频谱转换，在窗口右下方绘制出归一化的能量透射频谱图。
5. 命令行终端将同步输出仿真摘要，包括波导结构信息、通道宽度以及捕捉到的最大场强。

系统要求

• 环境支撑：MATLAB及其基础工具箱（推荐R2018a及以上版本）。
• 硬件要求：建议起步内存8GB，主频2.5GHz以上的处理器以保证计算效率，因为FDTD涉及大量矩阵循环运算。

主程序功能实现逻辑说明

系统的核心流程严格遵循FDTD的时间步迭代算法，具体逻辑步骤如下：

1. 物理与网格参数初始化：定义真空光速、介电常数等基础常数。建立空间长度为4微米、宽度为1微米的二维仿真网格，步长精确到10纳米，并根据Courant稳定条件计算采样时间步长。
2. 材料与结构定义：通过矩阵操作构建MDM波导。中间100纳米宽度被设定为介电常数为1的空气层，上下两侧被定义为满足Drude特性的金属区域。程序利用属性标记矩阵来区分介质和金属。
3. CPML参数预处理：在迭代开始前，针对计算域的四周边际计算CPML系数。通过指数建模和电导率梯度设置，生成用于更新psi卷积项的辅助系数。
4. 辅助差分方程（ADE）实现：针对金属部分的色散特性，引入电流密度J作为中间变量。利用双线性变换将色散关系从频域转换至时域，计算出金属内部电场与感应电流的更新系数。
5. 时域迭代更新循环：

* 磁场Hz更新：利用电场的旋度更新磁场值，同时叠加对应于x和y方向的CPML卷积项以实现边界吸收。 * 激励源注入：在预设的源位置引入中心波长1.55微米的高斯脉冲软源，脉冲在空间上呈现高斯分布以匹配波导基模。 * 电流密度J更新：在金属网格点处，根据当前的电场值和前一时刻的电流状态，更新此时的电流密度Jx和Jy。 * 电场Ex和Ey更新：利用磁场的旋度更新电场，对于金属格点，需额外扣除计算得到的电流密度项，以模拟金属支路的电流损耗与色散效应。

1. 监测与后处理：在仿真的特定位置设立监测点，记录每个步长的场强数据。最后对时域序列进行FFT运算，得到特定频率范围内的能量分布。

算法与关键逻辑细节分析

• TM波型离散化：代码实现了TM（横磁）模式的计算，即场分量包括磁场Hz以及电场Ex、Ey。这种模式是激发表面等离子体激元的必要前提。
• CPML卷积项更新：为了避免传统的PML分裂场方法，代码采用了卷积形式。通过引入psi_hzx、psi_hzy等辅助场，将复坐标拉伸效应实时耦合进麦克斯韦更新方程，这显著提高了吸收边界在掠射角下的表现。
• Drude模型时域耦合：使用了ADE-FDTD方案。 c1和c2系数是基于中心差分法推导出的，它们将复杂的频率相关性转化为了当前时刻场与上一时刻电流密度的加权组合，使得仿真能够实时处理金属的自由电子震荡响应。
• 软源（Soft Source）技术：电磁波的引入采用场量叠加而非直接赋值的方式，这允许从激励源后方反射回来的波无损地穿过源位置，避免了对仿真结果的干扰。
• 频谱归一化处理：在最终计算透射率时，通过对FFT后的分量求模平方，并寻找特定频率范围内的最大值进行归一化，使得结果能够直观反映能量随频率的变化趋势。