[论文解读] Solving Three Dimensional Maxwell Eigenvalue Problem with Fourteen Bravais Lattices

Tsung‐Ming Huang, Tiexiang Li+4 more|arXiv (Cornell University)|Jun 28, 2018

Photonic Crystals and Applications参考文献 32被引用 7

一句话总结

本论文通过在正交坐标系中重新表述布洛赫边界条件并利用离散旋度算子的可交换性，提出了一种统一的、无零空间的三维麦克斯韦特征值问题求解方法，适用于全部14种布拉维晶格。关键贡献在于通过酉合同变换对双旋度算子进行系统性特征分解，得到标准的3×3复反对称形式，从而实现对三维各向同性光子晶体能带结构的高效且精确计算。

ABSTRACT

Calculation of band structure of three dimensional photonic crystals amounts to solving large-scale Maxwell eigenvalue problems, which are notoriously challenging due to high multiplicity of zero eigenvalue. In this paper, we try to address this problem in such a broad context that band structure of three dimensional isotropic photonic crystals with all 14 Bravais lattices can be efficiently computed in a unified framework. We uncover the delicate machinery behind several key results of our work and on the basis of this new understanding we drastically simplify the derivations, proofs and arguments in our framework. In this work particular effort is made on reformulating the Bloch boundary condition for all 14 Bravais lattices in the redefined orthogonal coordinate system, and establishing eigen-decomposition of discrete partial derivative operators by systematic use of commutativity among them, which has been overlooked previously, and reducing eigen-decomposition of double-curl operator to the canonical form of a 3x3 complex skew-symmetric matrix under unitary congruence. With the validity of the novel nullspace free method in the broad context, we perform some calculations on one benchmark system to demonstrate the accuracy and efficiency of our algorithm.

研究动机与目标

• 开发一种统一的计算框架，用于求解全部14种布拉维晶格的三维麦克斯韦特征值问题。
• 克服由无散度约束引起的广义特征值问题中大规模零空间的挑战。
• 在正交坐标系中对所有晶格重新表述布洛赫边界条件，以确保与数值离散化的兼容性与一致性。
• 利用离散偏导数与旋度算子之间的相互可交换性，建立其系统性特征分解。
• 通过酉合同变换将双旋度算子约化为标准形式，从而实现能带结构问题的高效且精确求解。

提出的方法

• 在由原始斜交晶格矢量导出的新定义的正交坐标系中，对全部14种布拉维晶格重新表述布洛赫边界条件。
• 采用类似Yee的有限差分格式将旋度算子∇×E和∇×H离散化为矩阵-向量乘积，确保与重新表述后的边界条件兼容。
• 证明离散偏导数算子C1、C2、C3彼此可交换，从而实现同时对角化。
• 证明双旋度算子A = C*C可通过酉合同变换变换为3×3复反对称矩阵的标准形式。
• 利用C1、C2、C3的特征分解，系统性地识别旋度算子C的零空间与值空间的正交基。
• 通过投影到值空间构造无零空间的广义特征值问题，确保仅保留物理模态。

实验结果

研究问题

• RQ1 能否开发一个单一的统一框架，用于求解所有14种布拉维晶格的三维麦克斯韦特征值问题？
• RQ2 如何在正交坐标系中对任意布拉维晶格（尤其是三斜晶系）一致地表述布洛赫边界条件？
• RQ3 离散偏导数算子之间的可交换性在简化双旋度算子特征分解中起到何种作用？
• RQ4 双旋度算子能否在酉合同变换下约化为标准形式？这一约化如何实现广义特征值问题的高效求解？
• RQ5 如何在所有晶格中系统性地分析并消除旋度算子的零空间？

主要发现

• 作者成功地在正交坐标系中对全部14种布拉维晶格重新表述了布洛赫边界条件，实现了离散化的一致性。
• 证明了离散偏导数算子C1、C2、C3彼此可交换，从而支持同时对角化与解析特征分解。
• 通过酉合同变换，双旋度算子A = C*C被约化为标准的3×3复反对称矩阵形式，显著简化了特征值问题。
• 旋度算子的零空间被系统性地识别并消除，实现了无零空间的公式化，避免了虚假模态的出现。
• 在基准系统上的验证表明，该方法在计算三维光子晶体能带结构时具有高度的准确性和效率。
• 该框架具有普适性，适用于全部14种布拉维晶格，其中三斜晶系作为通用情形。