[论文解读] 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$ 对多体波函数与场进行去对称化，降低傅里叶变换中的计算成本。
• 通过使用满足 $\boldsymbol{\Gamma}_1 \mathbf{T}_{nl} \boldsymbol{\Gamma}_1 = 0$ 或 $\boldsymbol{\Gamma}_2 \widetilde{\mathbf{T}}_{nl} \boldsymbol{\Gamma}_2 = 0$ 的零-$\mathbf{T}$ 算子 $\mathbf{T}_{nl}$ 或 $\widetilde{\mathbf{T}}_{nl}$ 对 $\mathbf{L}$ 进行平移，解决 $\mathbf{L}$ 中的退化问题，从而消除零空间并确保可逆性。
• 将该框架应用于量子系统（如N电子薛定谔方程）、波传播以及具有共振子结构的超材料。

实验结果

研究问题

• RQ1 如何将一大类时间谐振线性连续体方程统一于单一数学框架之下？
• 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$ 通过对前两个坐标的所有排列进行6项交错求和实现，从而支持高效的傅里叶空间计算。
• $\boldsymbol{\Lambda}_A$ 对矢量场的作用由一个涉及索引排列的对称6项表达式给出，保持反对称性并降低计算成本。
• 该方法通过直接应用去对称化形式，避免了在实空间与傅里叶空间之间反复迭代，仅需对 $\mathbf{k}_1$ 与 $\mathbf{k}_2$ 进行傅里叶变换，适用于 $\mathbf{L}^D$-相关系统。
• 由于零空间导致的 $\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}$ 进行平移，得以解决。