On a semi-discrete model of Maxwell's equations in three and two dimensions

技術要約：マクスウェル方程式の幾何学的構造保存型半離散モデル

問題提起
数学的物理学の問題における固有の幾何学的および位相的構造を保存する離散モデルの構築は、信頼性が高く、物理的に一貫した数値シミュレーションを実現するために不可欠である。外延的計算（exterior calculus）を用いた電磁理論の離散化に関する研究は数多く存在するが、その多くは格子スキームやホイットニー形式（Whitney forms）に依拠している。本論文は、空間変数は離散化される一方で時間は連続的なままとなる、半離散的な定式化の必要性に対処するものである。目的は、ホイットニー形式や余チェイン（cochains）と微分形式間のド・ラーム写像を用いずに、離散外延的計算（DEC）の代数的および位相的構造を厳密に遵守する、3次元および2次元におけるマクスウェル方程式のフレームワークを開発することである。

手法
著者は、2次元モデルに関する先行研究を3次元へと拡張した、離散外延的計算に基づく幾何学的離散化フレームワークを採用している。手法の手順は以下の通りである：

1. 組合せ論的モデルの構築： $\mathbb{R}^3$ の組合せ論的モデルを、0次元および1次元の基底要素のテンソル積によって生成される3次元鎖複体 $C^{(3)}$ として定義する。離散形式（0形式、1形式、2形式、3形式）を表現するために、実数値の係数を持つ双対余チェイン複体 $K^{(3)}$ を導入する。
2. 離散演算子： 外延的計算の基本操作をこれらの複体上で定義する：
   • 余境界演算子（coboundary operator） $d_c$ は、境界演算子 $\partial$ との双対性を通じて定義され、外微分（exterior derivative）の離散的な類似体として機能する。
   • $\cup$-積（$\cup$-product） および 離散ホッジ・スター演算子（$*$） は、$K^{(3)}$ の基底要素上で定義される。特筆すべき点として、離散ホッジ・スターは二重適用時に指数のシフトが生じる（$** \neq \text{Id}$）。これは連続的な場合とは異なる重要な相違点である。
   • 余微分（codifferential） $\delta_c$ は、特定の内的積に関する $d_c$ の随伴として定義され、離散ラプラシアン $\Delta_c$ は $d_c\delta_c + \delta_cd_c$ として構成される。
3. 半離散定式化： 空間微分を離散演算子 $d_c$ で置き換える一方で、連続的な時間微分 $\frac{d}{dt}$ は保持することで、マクスウェル方程式を定式化する。構成関係（constitutive relations）は、定義されたホッジ・スターを用いて離散設定に適応させ、差分微分方程式系を導出する。
4. 2次元への還元とトーラス解析： 3次元モデルを、組合せ論的なトーラス上の2次元設定へと還元する。連続的な偏微分方程式は、解析解を得るために、定数係数を持つ一次線形常微分方程式（ODE）系へと変換される。

主要な貢献および結果

• 幾何学的保存： 本論文は、外延的計算の本質的な構造を保存する半離散モデルの構築に成功している。離散的なファラデーの法則、アンペールの法則、ガウスの法則、および構成関係の離散的な類似体を導出した。
• ポインティングの定理： ポインティングの定理の離散的な対応物が導出され（命題 3.1）、離散代数的枠組み内での電磁エネルギーの保存が示された。恒等式 $d_c(E \cup H) = -\frac{1}{2}\frac{d}{dt}(E \cup D + B \cup H) - E \cup J$ が確立されている。
• 波動方程式： 電場、磁場、および電磁ポテンシャルのための半離散的な波動方程式を導出した。これらは、離散ラプラシアン $\Delta_c$ と二階の時間微分を含む差分微分方程式の形式をとり、離散ホッジ・スターに固有のシフト演算子を組み込んでいる。
• ゲージ不変性： 半離散的なゲージ変換の類似体が定義されており、スカラーポテンシャルおよびベクトルポテンシャルの特定の変換の下で、電場1形式が不変であることが示されている。
• トーラス上での明示的解： 具体的な応用として、モデルを組合せ論的なトーラス（周期境界条件を持つ $2 \times 2$ グリッド）に適用する。源（$Q=0, J=0$）が存在しない場合、系は線形同次常微分方程式系へと簡約される。本論文は、系の行列の固有値および固有ベクトルを計算することにより、一般解の明示的な表現を提供する。固有値は $0$（重複度3）、$\pm 2$（重複度2）、および$\pm 2\sqrt{2}i$（重複度1）であり、定数項、指数項、および振動項からなる解が得られる。

意義
本論文は、連続理論の幾何学的および位相的な完全性を維持する、マクスウェル方程式の一貫した空間離散化を提供しているという点で意義がある。ホイットニー形式を回避し、離散外延的計算に基づく純粋に組合せ論的なアプローチを利用することで、本モデルは標準的な有限要素法や有限差分法に対する明確な代替案を提示している。半離散方程式を、組合せ論的なトーラス上での解けるODE系へと変換できる能力は、この手法の解析的な扱いやすさを実証している。本研究は、数学的物理学の基本方程式の離散的な類似物を開発するという著者の一連の研究の継続であり、構造保存型の数値シミュレーションのための厳密なフレームワークを提供している。