[論文レビュー] A unifying perspective on linear continuum equations prevalent in science. Part II: Canonical forms for time-harmonic equations

Graeme W. Milton|arXiv (Cornell University)|Jun 3, 2020

Thermoelastic and Magnetoelastic Phenomena参考文献 113被引用数 5

ひとこと要約

この論文は、拡張された複合材料理論に再定式化することで、物理学および工学における時間調和線形連続体方程式を統一的な数学的枠組みで提示する。高階の勾配を含まない方程式に対して標準形を導出し、フーリエ空間への射影作用素とリゾルベント法を用いた体系的解析と解法を可能にする。応用分野には波動伝播、メタマテリアル、量子系が含まれる。

ABSTRACT

Following some past advances, we reformulate a large class of linear continuum science equations in the format of the extended abstract theory of composites so that we can apply this theory to better understand and efficiently solve those equations. Here in part II we elucidate the form for many time-harmonic equations that do not involve higher order gradients.

研究の動機と目的

• 物理学および工学分野の多様な線形時間調和方程式を、共通の数学的枠組みに統一すること。
• 高階の勾配を含まない方程式を、拡張複合材料理論の標準形に再定式化すること。
• フーリエ空間への射影とリゾルベント技術を用いて、波動方程式、メタマテリアル、量子系の体系的解析と解法を可能にすること。
• 材料応答テンソルの退化を解消するために、作用素のシフトを用いてゼロラグランジアン寄与を特定・除去すること。
• 非局所的および隠れ変数系、特にメタマテリアルや量子多体系における潜在的微細構造を含む系の取り扱いを一般化すること。

提案手法

• フーリエ空間における自己随伴射影作用素 $\boldsymbol{\Gamma}_1$ を用いて、時間調和方程式を $\mathbf{J}(\mathbf{x}) = \mathbf{L}(\mathbf{x})\mathbf{E}(\mathbf{x}) - \mathbf{s}(\mathbf{x})$ の形に再定式化し、制約条件 $\boldsymbol{\Gamma}_1\mathbf{E} = \mathbf{E}$ および $\boldsymbol{\Gamma}_1\mathbf{J} = 0$ を課す。
• 支配的作用素を $\mathbf{D}(\nabla)^\dagger \mathbf{L} \mathbf{D}(\nabla) = \underline{z}\mathbf{I} - \mathbf{D}(\nabla)^\dagger \underline{\mathbf{B}} \mathbf{D}(\nabla)$ と表現し、$\underline{z}$ を $\omega^2$ またはエネルギー $E$ とすることで、リゾルベントに基づく解法を可能にする。
• フーリエ空間における微分的制約を満たすために、射影 $\boldsymbol{\Gamma}_1(\mathbf{k}) = \mathbf{D}(i\mathbf{k})[\mathbf{D}(i\mathbf{k})^\dagger \mathbf{D}(i\mathbf{k})]^{-1} \mathbf{D}(i\mathbf{k})^\dagger$ を適用する。
• 対称および反対称射影作用素 $\Lambda_a$ および $\boldsymbol{\Lambda}_A$ を用いて、多体波動関数および場を非対称化し、フーリエ変換における計算コストを低減する。
• $\mathbf{L}$ の退化を解消するために、$\mathbf{T}_{nl}$ または $\widetilde{\mathbf{T}}_{nl}$ と呼ばれるゼロ-$\mathbf{T}$ 作用素によるシフトを施し、零空間を除去し、可逆性を保証する。
• この枠組みを、N電子系（例：N電子シュレーディンガー方程式）、波動伝播、共振子構造を有するメタマテリアルに適用する。

実験結果

リサーチクエスチョン

• RQ1 広範な時間調和線形連続体方程式のクラスを、1つの数学的枠組みに統一することは可能か？
• RQ2 拡張複合材料理論に表現した場合、高階の空間勾配を含まない方程式がどのような標準形をとるか？
• RQ3 フーリエ空間への射影作用素 $\boldsymbol{\Gamma}_1(\mathbf{k})$ は、波動および連続体方程式においてどのように微分的制約を満たすか？
• RQ4 $\mathbf{L}$ におけるゼロラグランジアン寄与は可解性にどのように影響し、作用素のシフトによってどのように除去できるか？
• RQ5 多体波動関数の $\Lambda_a$ および $\boldsymbol{\Lambda}_A$ を用いた非対称化は、シュレーディンガー型方程式の解法における計算複雑性をどのように低減するか？

主な発見

• 関心の対象となる時間調和方程式は、フーリエ空間における自己随伴射影を用いて、体系的に $\mathbf{J} = \mathbf{L}\mathbf{E} - \mathbf{s}$ の標準形に書き換えることができ、$\boldsymbol{\Gamma}_1\mathbf{E} = \mathbf{E}$ および $\boldsymbol{\Gamma}_1\mathbf{J} = 0$ を満たす。
• 支配的作用素は $\underline{z}\mathbf{I} - \mathbf{D}(\nabla)^\dagger \underline{\mathbf{B}} \mathbf{D}(\nabla)$ の形を取り、リゾルベント $[\underline{z}\mathbf{I} - \mathbf{D}(\nabla)^\dagger \underline{\mathbf{B}} \mathbf{D}(\nabla)]^{-1} \mathbf{f}$ を用いた解法が可能になる。
• N電子系において、反対称化作用素 $\Lambda_a$ は最初の2つの座標に関する置換の6項の交代和として作用し、フーリエ空間での計算を効率化する。
• $\boldsymbol{\Lambda}_A$ がベクトル場に作用する際の式は、インデックスの置換を含む対称的な6項の表現で与えられ、反対称性を保持するとともに計算コストを低減する。
• 反復的な実空間／フーリエ空間の往復を回避するため、直接的に非対称化された形を適用し、$\mathbf{L}^D$-依存系では $\mathbf{k}_1$ および $\mathbf{k}_2$ のみにフーリエ変換を施すことで十分である。
• $\mathbf{L}$ の零空間に起因する退化は、$\boldsymbol{\Gamma}_1 \mathbf{T}_{nl} \boldsymbol{\Gamma}_1 = 0$ または $\boldsymbol{\Gamma}_2 \widetilde{\mathbf{T}}_{nl} \boldsymbol{\Gamma}_2 = 0$ を満たすゼロ-$\mathbf{T}$ 作用素による $\mathbf{L}$ または $\mathbf{L}^{-1}$ のシフトによって解消される。