采用低波数提取宽带格林函数（BBGFL）对周期性结构能带图进行宽带计算

1：实验设计与方法选择：

本研究采用BBGFL方法为含不可穿透物体的二维周期结构建立表面积分方程。该方法通过低波数提取加速周期格林函数的收敛，并结合矩量法推导出适用于宽带求解的线性特征值方程。

2：样本选择与数据来源：

数值算例聚焦于半径b=0.2和b=0.4的圆形不可穿透圆柱构成的二维六方晶格，晶格常数a=1，背景介电常数ε=8.9ε0。波矢选自第一布里渊区。

3：2和b=4的圆形不可穿透圆柱构成的二维六方晶格，晶格常数a=1，背景介电常数ε=9ε0。波矢选自第一布里渊区。 实验设备与材料清单：

3. 实验设备与材料清单：计算工具包括用于数值模拟的MATLAB软件，在配备Intel Core i7-4790 CPU和32GB内存的HP ProDesk 600 G1台式机上运行。因属理论数值研究未使用实体设备。

4：实验流程与操作步骤：

流程包含：(a)选取低波数kL并计算该波数下的周期格林函数；(b)应用BBGFL方法表达具有加速收敛特性的格林函数；(c)采用脉冲基函数和点匹配法通过矩量法离散表面积分方程；(d)构建线性特征值问题同步求解多频带特征值；(e)通过与KKR方法对比验证结果。

5：数据分析方法：

计算不同Floquet模式截断数（如M=49、121、1681）对应的特征值。通过将归一化频率fN与KKR结果对比评估精度，利用CPU时间测量评估计算效率。