1. 项目背景与核心思路
最近在实验室里折腾光子晶体结构时，突然想到个有趣的玩法——用COMSOL模拟两个光子晶体薄片扭转叠加后的能带结构和透射特性。这种结构让我联想到魔角石墨烯中那个著名的1.1°扭转角现象，只不过这次的主角换成了光子晶体。当两个周期性介质以特定角度叠加时，它们的电磁模式会发生奇妙的耦合，这在光学器件设计和新型光子学材料开发中都是个值得深挖的方向。
光子晶体本质上是通过周期性介电常数变化来操控光子运动的"光学半导体"，其能带结构直接决定了哪些频率的光可以传播（导带），哪些会被禁止（禁带）。传统单层光子晶体的能带计算已经比较成熟，但双层扭转结构的模拟却有几个技术难点：首先是扭转后的超晶胞建模，其次是耦合模式的计算方法选择，最后是如何准确提取透射谱特征。
2. 模型构建关键步骤
2.1 基础参数设定
我选择的是二维三角晶格光子晶体，介质柱材料设为硅（ε=12），背景为空气。晶格常数a=500nm，圆柱半径r=0.3a。这个参数组合在近红外波段（1500-1600nm）会产生明显的带隙，便于后续观察扭转带来的影响。
在COMSOL中建立单层模型时需要注意：
- 使用"波束包络"接口而非传统频域求解，可以更好地处理周期性边界条件
- 设置Floquet周期性边界时，需要正确定义基矢方向
- 网格划分要保证每个介质柱边缘有至少5层单元
2.2 扭转结构建模技巧
当两层晶体相对旋转θ角后，系统会形成摩尔纹超晶胞。这里有个计算技巧：超晶胞尺寸由公式L=a/(2sin(θ/2))决定。例如当θ=5°时，超晶胞边长约为5.7a，即2.85μm。
实际操作中我采用了分步建模法：
- 先建立单层晶体的完美周期结构
- 通过旋转复制得到第二层晶体
- 使用"形成装配体"功能确保两层独立可控
- 用参数化扫描研究不同扭转角的影响
关键提示：COMSOL在处理大尺寸超晶胞时容易内存溢出，建议先用小角度（如θ>3°）测试模型可行性，再逐步细化。
3. 能带计算与结果分析
3.1 求解器配置要点
在能带计算阶段需要特别注意：
- 使用"频域"研究配合"布洛赫周期条件"
- 设置正确的波矢扫描路径（通常沿超晶胞的第一布里渊区边界）
- 启用"参数化扫描"同时求解多个k点
- 内存不足时可启用"集群计算"选项
我采用的典型求解器配置：
matlab复制study = mphstudy(model, "BandStructure");
study.set("pname", ["kx", "ky"]);
study.set("plistarr", {linspace(0,pi,50), linspace(0,pi,50)});
study.run();
3.2 典型能带特征
当扭转角θ=3°时，观测到三个显著现象：
- 原带隙区域出现扁平化能带（类似魔角石墨烯的Dirac点）
- 在ω=0.28(2πc/a)附近出现新的局域态
- 导带边缘出现明显的能带分裂
这些特征可以通过电场分布图进一步验证：
- 扁平能带对应着两层晶体的局域耦合模式
- 新出现的态密度峰与摩尔超晶胞的共振相关
- 能带分裂源于扭转导致的对称性破缺
4. 透射谱模拟实战
4.1 边界条件设置
透射谱计算需要建立有限尺寸模型：
- 上下设置完美匹配层(PML)吸收边界
- 左右采用周期性边界条件
- 顶部添加平面波激励端口
- 底部设置功率传输监测端口
关键参数关系：
- PML厚度 ≥ 工作波长/2
- 计算区域高度 ≥ 5倍晶格常数
- 频域扫描步长 ≤ 最小特征线宽/5
4.2 典型透射特征
在θ=5°时观测到：
- 原带隙区域出现窄带透射峰（Q因子>2000）
- 透射谱线型呈现不对称Fano共振特征
- 角度变化1°会导致透射峰位移约15nm
这种现象可以解释为：
- 窄带峰对应扭转诱导的局域共振模
- Fano线型源于连续态与离散态的干涉
- 角度敏感性来自摩尔纹周期的指数依赖
5. 常见问题排查指南
5.1 收敛性问题
症状：能带出现锯齿状波动
解决方法：
- 增加PML层数（建议8-10层）
- 细化k点采样密度（至少50点/路径）
- 检查周期性边界条件设置
5.2 内存不足报错
优化策略：
- 改用频域求解器的"直接"矩阵求解法
- 降低网格最大单元尺寸（不超过λ/10）
- 分段计算不同k点区间后合并结果
5.3 物理场耦合异常
典型表现：电场分布不连续
检查清单：
- 确认装配体接合面正确设置"连续性"条件
- 验证材料参数在各域中赋值正确
- 检查网格在介质界面处的过渡是否平滑
6. 进阶应用方向
这种扭转光子晶体结构有几个潜在的应用价值：
- 角度可调光学滤波器：通过机械旋转实现波长精准调控
- 高Q值微腔设计：利用扁平能带增强光场局域
- 新型传感器开发：摩尔纹结构对环境折射率异常敏感
我在实验中发现一个有趣现象：当在两层晶体间加入10nm厚的氧化铝间隔层时，θ=2.5°会出现特别尖锐的共振峰（线宽<0.1nm）。这可能是由于间隔层改变了界面耦合强度，具体机制还需要进一步研究。