基于平面波展开法的三维光子晶体能带计算系统

资 源 简 介

本项目旨在利用数值计算方法分析三维光子晶体中的电磁波传播特性。核心功能是通过平面波展开法（Plane Wave Expansion Method, PWE）将麦克斯韦方程组在倒格空间中转化为线性本征值问题。系统支持多种典型晶格类型，包括简单立方（SC）、面心立方（FCC）和体心立方（BCC）等，并允许用户自定义介质构型，如球体、立方体或管状结构。 该项目详细实现了介电常数的傅里叶空间变换，通过构建哈密顿矩阵并调用高性能特征值求解器，计算在第一布里渊区沿高对称路径分布的波矢量k对应的本征频率。其实现过程涵盖

详 情 说 明

三维光子晶体能带结构计算系统

项目介绍

本项目是一个基于平面波展开法（Plane Wave Expansion Method, PWE）的高性能数值模拟工具，专门用于计算和分析三维周期性介质结构（光子晶体）中的电磁波色散关系。通过将麦克斯韦方程组转化为倒格空间中的线性本征值问题，系统能够精确求解光子能带结构，帮助研究人员识别光子禁带（Photonic Band Gap）并探索光与物质的相互作用。

主要功能特性

• 多晶格类型支持：系统内置了简单立方（SC）、面心立方（FCC）和体心立方（BCC）三种典型三维晶格的基矢定义及高对称路径规划。
• 参数化结构定义：支持自定义背景介质与填充材料的介电常数、填充球体半径以及倒格矢截断数量，便于进行收敛性分析和结构优化。
• 倒格空间自动离散化：根据截断频率自动生成倒格矢集合，并采用球形截断方案以保持倒空间的对称性。
• 矢量波本征方程求解：考虑了电磁波的横波特性，在每个平面波基底上构建两个正交的偏振分量，完整求解三维矢量波方程。
• 自动化可视化：实时生成三维单元晶格结构示意图，并自动绘制沿第一布里渊区高对称路径的能带图。
• 禁带分析：系统能够自动计算并输出能带结构中存在的频率禁带信息。

实现逻辑分析

核心算法遵循以下严谨的数值计算流程：

1. 参数初始化与坐标构建：设定晶格常数、填充率及折射率对比度。根据选定的晶格类型（如FCC），计算实空间基矢 a1, a2, a3 以及对应的倒空间基矢 b1, b2, b3。
2. 生成倒格矢集合：在三维空间内通过整数线性组合生成候选 G 矢量，并利用模值筛选（球形截断）确定参与计算的平面波数量，确保基底的完整性与计算效率。
3. 构建介电矩阵（Ho方法）：

* 计算倒格空间中介电常数的傅里叶变换 $epsilon(G-G')$。 * 对于球体结构，采用解析形状因子进行计算，以获得更高的数值精度。 * 先构建 $epsilon$ 矩阵再进行求逆运算得到 $eta$ 矩阵（$1/epsilon$），此方法在处理高对比度介质时具有更好的收敛性。

1. 路径采样与循环求解：

* 在倒空间高对称点之间进行线性插值，生成连续的波矢量 $k$ 采样序列。 * 对于每个 $k$ 点，根据横波条件（$nabla cdot D = 0$）构造偏振基矢量。 * 组建大小为 $2N_{pw} times 2N_{pw}$ 的哈密顿矩阵，其中 $N_{pw}$ 为平面波数量。

1. 特征值并行求解：调用高性能特征值求解器计算频率本征值，提取归一化频率，并剔除由于数值误差产生的零模。
2. 后处理：将计算得到的数个能带按顺序排列，进行绘图展示，并标识出高对称点位置。

核心辅助功能说明

• 晶格与路径配置：负责定义不同晶格的几何常数及其在布里渊区巡行的标准坐标点，如 FCC 的 X-U-L-Gamma-X-W 路径。
• 倒格矢生成逻辑：通过网格化搜索与距离判定，筛选出指定能量范围内的所有有效倒格矢。
• 傅里叶空间系数计算：利用解析公式计算三维形状因子的傅里叶变换，区分 $G=G'$ 和 $G neq G'$ 的不同情况。
• 偏振基矢量构造：针对每一个波矢量 $k+G$，构造两个相互正交且同时垂直于波矢方向的单位向量，以满足电磁场横波属性。
• 路径插值处理：在离散的高对称点之间自动插入计算点，确保能带结构的平滑度。

使用方法

1. 配置参数：在主函数开头部分修改晶格类型（lattice_type）、几何形状（geometry_type）、介电常数（eps）以及填充半径（r_a）。
2. 精度控制：通过调整 N_cutoff 控制平面波展开的数量，较大的数值会提高精度但显著增加计算时间。
3. 运行模拟：执行主脚本，系统将依次进行倒格空间矩阵构建、特征值求解和图形化展示。
4. 查看结果：程序左侧窗口显示晶格物理结构，右侧窗口展示能带图，命令行窗口将输出禁带计算结果。

系统要求

• 环境：MATLAB R2016b 或更高版本。
• 工具箱：基础 MATLAB 内置函数即可支持核心计算，图形显示需要绘图工具箱支持。
• 内存建议：对于 $N_{cutoff} > 5$ 的计算，建议配备 8GB 以上内存以处理大规模矩阵运算。