来源论文: https://arxiv.org/abs/2606.16192v1 生成时间: Jun 21, 2026 01:23

无限格点离散 Calderón 投影：基于格点格林函数（LGF）的主动声学屏蔽与围护理论及量子紧束缚模型类比

0. 执行摘要

在现代计算物理、声学工程以及量子化学/材料物理领域，边界控制与波包传播的精确描述始终是核心课题。经典的主动噪声控制（Active Noise Control, ANC）和声学屏蔽任务，旨在通过在一组边界源上施加控制，来抵消外部传入的有害噪声，或将内部产生的噪声限制在特定区域（围护）内。传统的离散化 Calderón 边界投影方法（如基于 Tsynkov 差分势方法的方案）通常极度依赖于一个有界的辅助区域（Auxiliary Domain）以及在其外部人为设定的吸收边界条件（如 PML 或者是近似的辐射边界条件）。这种传统方案不仅引入了显著的几何相关性，使得每次改变保护区域几何形状时都需要重新求解昂贵的体积分方程，而且由于边界近似误差的存在，极大地限制了其消声极限与数值稳定性。

近期，由庆夏（Qing Xia）提出的 “Infinite-lattice discrete Calderón projection via the lattice Green’s function for active noise shielding and confinement” 一文，从根本上革新了这一范式。该研究通过引入**格点格林函数（Lattice Green’s Function, LGF）**进行无限格点卷积，成功构建了一种完全不依赖于人工外边界的无限格点离散 Calderón 投影算子。该方法将精确的 Sommerfeld 辐射条件直接编入 LGF 核函数中，使几何相关的体积分完全退化为仅在单层晶格边界带（Boundary Strip）上进行的代数运算。LGF 的数值表征仅取决于网格步长 $h$ 和波数 $k$，一经计算即可在各种任意复杂的几何结构（如圆形、L型、非凸星形等）之间完全复用。

本博客将面向多物理场数值模拟、量子输运、晶格动力学及量子化学计算领域的科研人员，对该项工作进行极其详尽、系统且严谨的深度技术剖析。我们将从其核心数学物理方程出发，推导离散 Calderón 投影的一般化代数形式，详述其容量矩阵（Capacity Matrix）的高效求解策略，并全面剖析其在复杂边界条件下的收敛表现。最后，我们将该方法与量子化学中著名的紧束缚模型（Tight-Binding Model）、**非平衡格林函数（NEGF）以及自能（Self-Energy）**边界条件进行深度学术类比，揭示这一声学控制技术在微观量子输运与相干波函数控制中的广阔应用前景。

1. 核心科学问题，理论基础，技术难点，方法细节

1.1 核心科学问题：空间连续域的主动屏蔽与围护

主动噪声控制的数学本质是：在一个受保护的物理区域 $\Omega_{\text{in}}$ 内，通过布置在边界 $\partial\Omega_{\text{in}}$ 附近的二次源（Secondary Sources）产生一个主动控制场 $g$，以完全抵消来自外部区域 $\Omega_{\text{out}}$ 的入射噪声场 $u^-$，同时必须保证不破坏区域内部本身存在的、由 wanted 源产生的主动声场 $u^+$。这一问题的极高挑战性在于：

1. 控制系统无需获知外部噪声源的具体空间分布、强度以及相位信息。
2. 控制系统无需获知外部介质的宏观异质性参数（只要在边界邻域内介质是均匀的）。
3. 控制场必须在物理上表现为**纯向外辐射（Outgoing）**的，以满足能量守恒并防止在无限远处产生反射。

在数学上，这一物理机制能够通过连续介质中的 Calderón 投影算子（Calderón Boundary Projector）获得完美的描述。连续 Calderón 投影本质上是一个边界边界算子，它将边界上的全场迹（Trace）投影到可以向内部（或外部）单边延拓的 Helmholtz 解空间。然而，当我们试图在笛卡尔网格（Cartesian Grids）上对该系统进行差分离散时，如何保持其边界投影的幂等性（Idempotency）以及完全向外辐射的格林函数行为，构成了计算物理中的核心难题。

1.2 理论基础：离散 Helmholtz 算子与格点格林函数（LGF）

我们在 $d$ 维笛卡尔晶格 $\Lambda_h = h\mathbb{Z}^d$ ($d \in \{2,3\}$) 上讨论这一问题。定义标准的二阶中心差分 Helmholtz 算子 $L^h$ 如下：

(1)

其中 $h > 0$ 为网格步长，$k$ 为物理波数，$e$ 为晶格的单位方向向量。在自由空间晶格上，伴随 Sommerfeld 辐射条件的温和基本解（Tempered Fundamental Solution）即为格点格林函数（LGF），记作 $G^h_m$。其严格的傅里叶积分表征为：

(2)

其中分母中的色散关系项为：

这里的 $i0^+$ 表示积分路径在复平面上向上微调，以严格选择满足格点辐射条件的向外传播波。LGF 是算子 $L^h$ 在无限晶格上的逆，满足关系：

定义无限晶格上的离散卷积运算为：

(3)

对于任意紧支撑的网格函数 $v^h$，其卷积结果严格满足 $L^h(G^h * v^h) = v^h$。

1.3 边界带（Boundary Strip）的拓扑划分

考虑一个有界连续区域 $\Omega_{\text{in}} \subset \mathbb{R}^d$。在晶格上，我们定义内部格点集 $M^+ = \{m \in \mathbb{Z}^d : hm \in \Omega_{\text{in}}\}$ 以及外部格点集 $M^- = \mathbb{Z}^d \setminus M^+$。通过五点（2D）或七点（3D）差分模板对 $L^h$ 进行空间扫描，会自然产生一组重叠边界。我们定义边界带（Boundary Strip） $\gamma$ 为那些差分模板同时跨越 $M^+$ 和 $M^-$ 的所有节点的集合。进一步，边界带可以被严格剖分为相互不重叠的两个单层子层：

• 内边界层（Interior Boundary Layer）：$\gamma^+ = \gamma \cap M^+$
• 外边界层（Exterior Boundary Layer）：$\gamma^- = \gamma \cap M^-$

这种基于晶格差分模板的自然划分如图 1 所示。在数学上，在边界带 $\gamma$ 上定义的网格函数即构成了离散迹空间（Discrete Trace Space），这是主动控制理论施加边界密度的唯一空间支撑。

1.4 技术难点：消除几何相关的有界辅助问题

在传统的差分势方法（Difference Potentials Method, DPM）中，为了在有界网格上求逆 Helmholtz 算子，研究者必须引入一个包含 $\Omega_{\text{in}}$ 的有界辅助区域 $\Omega_0 \supset \Omega_{\text{in}}$。这带来了一系列严重的理论和工程计算缺陷：

1. 人工边界反射：在 $\partial\Omega_0$ 上设置的边界条件（如近似的 ABC 或者是 PML）无法在离散层面上完美契合自由空间格林函数，导致伪反射，从而极大地损害了消声分贝数。
2. 体积分重构的高昂成本：每当保护区几何 $\Omega_{\text{in}}$ 改变时，必须对整个辅助区域 $\Omega_0$ 重新进行网格剖分并重新求解庞大的线性系统。其控制合成无法在不同的几何构型间共享任何公共计算成分。

1.5 方法细节：基于 LGF 的无限格点离散 Calderón 投影

为了彻底解决上述痛点，该研究提出了一种无需设置外边界的无限格点离散 Calderón 投影。对于任意给定的边界迹 $\xi_\gamma : \gamma \to \mathbb{C}$，我们首先构建其容许延拓（Admissible Extension） $V^h : \Lambda_h \to \mathbb{C}$，使其满足 $V^h|_{\gamma} = \xi_\gamma$，并且 $L^h V^h$ 具有紧支撑（这在代数上总是可行的，比如直接将其在 $\gamma$ 之外设为 0）。

基于 LGF 卷积，我们定义离散 Calderón 投影算子 $P^h$ 为：

(4)

其中 $\mathbf{1}_{M^\pm}$ 为区域指示函数。我们取该投影场在边界带 $\gamma$ 上的迹，定义离散边界 Calderón 投影为：

(5)

定理 3.2 的严格代数证明（核心性质）

该论文的核心理论基石在于证明了算子 $P^h_\gamma$ 确实是一个结构优雅的边界投影。我们在此给出其三大关键性质的推导：

1. 延拓独立性：设 $V_1^h$ 和 $V_2^h$ 是同一迹 $\xi_\gamma$ 的两个不同容许延拓，定义差值场 $W^h = V_1^h - V_2^h$。显然 $W^h|_{\gamma} = 0$。根据差分算子的局部性， stencil 跨越边界的耦合项在做差后完全消失，我们有 $(L^h W^h)\mathbf{1}_{M^+} = L^h (W^h \mathbf{1}_{M^+})$。因此：

   $$W^h - G^h * \left[ (L^h W^h) \mathbf{1}_{M^+} \right] = W^h - G^h * (L^h (W^h \mathbf{1}_{M^+})) = W^h - W^h \mathbf{1}_{M^+} = W^h \mathbf{1}_{M^-}$$

   因为 $W^h \mathbf{1}_{M^-}$ 在整个内区域 $M^+$（包括 $\gamma^+$）上为 0，并且在外部边界 $\gamma^-$ 上由于 $W^h|_{\gamma} = 0$ 亦为 0，所以该差值场在整个边界带 $\gamma$ 上的迹严格为 0。这就证明了 $P^h_\gamma \xi_\gamma$ 的定义与具体的延拓方式 $V^h$ 无关。

2. 幂等性（Idempotency, $(P^h_\gamma)^2 = P^h_\gamma$）：设 $u^h = P^h \xi_\gamma$。直接对其作用算子 $L^h$，由于对基本解卷积的性质，我们有 $L^h u^h = L^h V^h - (L^h V^h)\mathbf{1}_{M^+} = (L^h V^h)\mathbf{1}_{M^-}$。这表明 $L^h u^h$ 的支撑集完全位于 $M^-$ 内。因此，若我们以 $u^h$ 本身作为其迹 $\eta_\gamma = u^h|_{\gamma}$ 的容许延拓，并代入投影定义式 (4) 中：

   $$P^h \eta_\gamma = u^h - G^h * \left[ (L^h u^h) \mathbf{1}_{M^+} \right]$$

   由于 $L^h u^h$ 在 $M^+$ 上为 0，后一项的卷积项直接消失。因此 $P^h \eta_\gamma = u^h$。取迹后即得 $P^h_\gamma (P^h_\gamma \xi_\gamma) = P^h_\gamma \xi_\gamma$。投影的代数自洽性得证。

3. 值域特征（Range Characterization）：$\xi_\gamma \in \text{Im}(P^h_\gamma)$ 当且仅当存在一个网格函数 $u^h$，在保护区内部满足齐次离散 Helmholtz 方程 $L^h u^h = 0$（在 $M^+$ 上），且其迹满足 $u^h|_{\gamma} = \xi_\gamma$。这表明该离散投影精确地指代了内部格点 Helmholtz 方程解的迹空间。

1.6 容量矩阵（Capacity Matrix）表征与主动控制源的闭合求解

在计算上，直接计算无限晶格卷积 $G^h *$ 显然是不现实的。该研究通过引入容量矩阵（Capacity Matrix），将上述无限维空间中的卷积投影退化为单层晶格边界带上的稠密代数方程求解。

定义边界带 $\gamma$ 上的离散单层势矩阵（Capacity Matrix）$S_{\gamma\gamma}$ 为：

(6)

假定 $S_{\gamma\gamma}$ 是非奇异的。对于任意实测的边界总场迹 $\xi_\gamma = u^h|_{\gamma}$，我们可以通过求解如下稠密线性方程组，直接获得等效的边界电荷密度 $\lambda$：

(7)

我们将该电荷密度按照边界带的子层结构直接拆分为两部分：

利用单层势作为容许延拓 $V^h = G^h * \lambda$，代入离散 Calderón 投影公式中，利用极其优雅的代数消元，可以得到如下令人惊叹的闭合形式控制结论：

(8)

这意味着：

• 主动声学屏蔽控制（Shielding Control）：为了消除受保护区域 $M^+$ 内部的外部传入噪声，我们只需要在单层外边界 $\gamma^-$ 上布置主动二次控制源 $g^h_{\text{sh}}$，其控制源强度直接由 $\lambda^-$ 给出： $$g^h_{\text{sh}} = -\lambda^-$$ (9)
• 主动声学围护控制（Confinement Control）：为了将内部有害源产生的声场完全禁锢在 $M^+$ 内部、不泄露到外部空间 $M^-$ 中，我们只需要在单层内边界 $\gamma^+$ 上布置主动控制源 $g^h_{\text{conf}}$： $$g^h_{\text{conf}} = -\lambda^+$$ (10)

这种控制源只需要在单晶格层厚度的物理边界上施加，就能够实现完美的体消声效果（Volumetric Cancellation），这为消声系统的工程实现提供了极大的便利。

2. 关键 Benchmark 体系，计算所得数据，性能数据

为了全面验证上述基于 LGF 的离散 Calderón 投影方法的精确性、收敛阶以及在复杂几何拓扑下的鲁棒性，作者在二维空间中设计了三个极具代表性的计算流体/声学测试体系。

2.1 测试几何与格点规模

这三个 Benchmark 几何体系如下：

1. 圆形保护区（Circle, $r = 0.5$）：代表光滑无尖角的经典标准基准。其边界带格点规模在网格数 $n=127$ 时，内点数 $|M^+|=697$，边界迹点数 $|\gamma|=172$（其中外边界层点数 $|\gamma^-|=88$）。
2. L型非凸保护区（L-shape）：包含一个 $90^\circ$ 的内缩重入角（Reentrant Corner）。这是边界元法和经典差分法中极易产生奇异性与流场不稳定性（如角点奇点效应）的经典几何。其边界极难用平滑解析曲线拟合。在 $n=127$ 时，点数规模为 $|M^+|=830$，$|\gamma|=322$，$|\gamma^-|=163$。
3. 星形高度非凸保护区（Star-shaped）：边界半径由高频简谐摄动曲线给出：$r(\theta) = 0.5 + 0.12\cos(5\theta)$。该几何具有 5 个对称的凸起和凹陷，边界曲率变化剧烈，极度考验网格离散化在阶梯晶格边界下的拟合能力。在 $n=127$ 时，$|M^+|=773$，$|\gamma|=212$，$|\gamma^-|=108$。

所有模拟均在一个具有 $n \times n$ 个节点的笛卡尔晶格 $\Lambda_h$ 上进行，采样物理窗口统一设定为 $[-2.15, 2.15]^2$。对应的网格步长在 $n=127$ 时约为 $h \approx 0.0336$。

2.2 物理声场与入射源设置

研究中考虑了两类完全不同的入射源：

1. LGF一贯源（LGF-consistent Point Sources）：通过在网格节点上直接放置离散 $\delta$ 函数并与 LGF 卷积产生。此类源在离散层面上严格满足离散 Helmholtz 方程。这一体系用于验证离散代数恒等式的数学绝对精确性。
2. 连续解析源（Analytic Incident Fields）：直接计算连续空间下的解析波形并取样到格点上，包括：
   • 解析平面波：$u^{\text{pw}}(\mathbf{x}) = e^{i \mathbf{k} \cdot \mathbf{x}}$，其中方向向量 $\mathbf{d} = (\sqrt{3}/2, 0.5)$。
   • 解析自由空间点源（单极子辐射）：$u^{\text{ps}}(\mathbf{x}) = \frac{i}{4} H_0^{(1)}(k |\mathbf{x} - \mathbf{x}_s|)$，其中噪声源位置 $\mathbf{x}_s$ 被放置在紧邻保护区外部的 $(0.9, 0)$ 或 $(1.0, 0.5)$ 位置。这类源不满足离散方程，用于测试离散化一致性误差（Discretization Consistency Error）。

2.3 核心计算性能与收敛数据

2.3.1 LGF一贯源下的极限消声表现（Self-consistency Check）

当入射噪声由严格满足离散 Helmholtz 方程的 LGF 一贯源产生时，数值计算结果完全验证了理论上的代数恒等关系。如图 3 所示，在 $k = 5$ 且网格数 $n = 127$ 的设定下，控制施加前后的全场相对 $\ell_2$ 残差（Relative Residuals）达到了计算机双精度漂移的极限（Machine Precision）：

• 圆形保护区：相对残差 $\varepsilon_{\text{rel}} = 6.1 \times 10^{-15}$
• L型保护区：相对残差 $\varepsilon_{\text{rel}} = 5.5 \times 10^{-15}$
• 星形保护区：相对残差 $\varepsilon_{\text{rel}} = 8.5 \times 10^{-15}$

图 4 给出了保护区内部的逐点消声分贝数（Pointwise Attenuation, 衡量消声强度的核心工程指标）：

其计算所得的数据表明，在整个保护区内部，中位数消声强度（Median Attenuation）统一超过了惊人的 215 dB（其中圆形为 219 dB，L型为 217 dB，星形为 219 dB），且整个内部区域消声极其均匀，这彻底证明了单层边界控制源 $g^h_{\text{sh}} = -\lambda^-$ 能够完美地进行全域体空间消声。

2.3.2 解析源下的二阶收敛行为（Discretization-error Test）

当采用不满足离散 Helmholtz 的连续解析源（平面波与连续点源）时，消声系统受限于二阶中心差分模板 $L^h$ 的截止误差（Truncation Error）。作者通过在 $n = 31, 63, 127, 255$ 四个递进网格尺度上进行网格细化实验（Grid-refinement Study），得到了极其优异的近二阶收敛性能。具体收敛数据如表 3 和表 4 所示：

• 表 3 汇总：平面波入射下（$k=5$）的三种几何收敛数据
  • $n = 31$ ($h \approx 0.134$):
    • 圆形：$\varepsilon_{\text{rel}} = 2.34 \times 10^{-2}$，消声度 $= 37.9\text{ dB}$
    • L型：$\varepsilon_{\text{rel}} = 2.11 \times 10^{-2}$，消声度 $= 34.1\text{ dB}$
    • 星形：$\varepsilon_{\text{rel}} = 2.23 \times 10^{-2}$，消声度 $= 37.3\text{ dB}$
  • $n = 63$ ($h \approx 0.0672$):
    • 圆形：$\varepsilon_{\text{rel}} = 7.57 \times 10^{-3}$，消声度 $= 46.4\text{ dB}$
    • L型：$\varepsilon_{\text{rel}} = 5.54 \times 10^{-3}$，消声度 $= 45.5\text{ dB}$
    • 星形：$\varepsilon_{\text{rel}} = 7.22 \times 10^{-3}$，消声度 $= 46.4\text{ dB}$
  • $n = 127$ ($h \approx 0.0336$):
    • 圆形：$\varepsilon_{\text{rel}} = 2.01 \times 10^{-3}$，消声度 $= 57.8\text{ dB}$
    • L型：$\varepsilon_{\text{rel}} = 1.34 \times 10^{-3}$，消声度 $= 57.6\text{ dB}$
    • 星形：$\varepsilon_{\text{rel}} = 1.93 \times 10^{-3}$，消声度 $= 57.4\text{ dB}$
  • $n = 255$ ($h \approx 0.0168$):
    • 圆形：$\varepsilon_{\text{rel}} = 5.49 \times 10^{-4}$，消声度 $= 69.0\text{ dB}$
    • L型：$\varepsilon_{\text{rel}} = 3.68 \times 10^{-4}$，消声度 $= 68.9\text{ dB}$
    • 星形：$\varepsilon_{\text{rel}} = 5.30 \times 10^{-4}$，消声度 $= 68.7\text{ dB}$
  • 拟合收敛阶（Fitted Rate）：三种几何均表现出稳定的 $\approx 1.8 \sim 1.9$ 阶收敛。这完全符合二阶差分格式在阶梯逼近不规则边界时的理论上限。

2.3.3 单边测量（One-sided Measurement）的性能对比

在实际工程中，在保护区内部（$\gamma^+$ 层）布置传感器可能受到物理空间限制，因此单边测量（仅在外部 $\gamma^-$ 层布置麦克风）极具实际应用价值。表 5 给出了在圆形保护区下，全场双边测量（Two-sided）与仅在 $\gamma^-$ 进行单边测量再通过映射矩阵 $T^-_h = S_{+-} S_{--}^{-1}$ 进行全迹重构的性能及耗时对比：

• $n=127$ 细化下：
  • 双边测量：相对残差 $\varepsilon_2 = 2.2 \times 10^{-3}$，消声度 $\text{Att}_2 = 57\text{ dB}$，计算墙钟时间 $t_2 = 1.04\text{ s}$。
  • 单边测量：相对残差 $\varepsilon_1 = 7.0 \times 10^{-3}$，消声度 $\text{Att}_1 = 45\text{ dB}$，时间 $t_1 = 1.04\text{ s}$。
• 数据分析：单边测量的消声分贝数因逆问题求解的条件数放大而略微降低（降低了约 $10 \sim 12\text{ dB}$），但随着网格细化，单边残差仍表现出极为稳定的下降趋势（$n=255$ 时单边消声度回升至 $55\text{ dB}$）。这为非侵入式声学控制提供了强有力的基准支撑。

2.3.4 高波数敏感性与测量噪声鲁棒性分析

• 高波数测试（$k=10$）：如表 6 所示，当波数翻倍至 $k=10$ 时（对应网格分辨率降低至每个波长仅约 19 个格点），在 $127 \times 127$ 网格下，三种几何的消声表现降至 $35 \sim 41\text{ dB}$，但容量矩阵的条件数 $\kappa(S_{\gamma\gamma})$ 反而表现得更加温和（圆形为 $6.5 \times 10^1$，L型为 $9.2 \times 10^1$）。这表明高频下的控制性能下降主要源自网格离散解析度的不足，而非控制代数系统的病态化。
• 抗噪性能分析：在实测迹 $\xi_\gamma$ 中加入相对标准差为 $\sigma$ 的高斯复噪声。表 8 详尽记录了圆域消声度的衰减趋势：
  • 当噪声水平 $\sigma \le 10^{-4}$ 时，消声度维持在 $56.8\text{ dB}$ 几乎完全没有降级（离散一致性误差占主导）。
  • 在强噪声污染 $\sigma = 10^{-2}$ （$1\%$ 相对噪声）下，系统仍保留了 $\text{Att} = 34.4\text{ dB}$ 的可观消声水平。
  • 仅在极度恶劣的噪声 $\sigma = 10^{-1}$（$10\%$ 相对噪声）下，消声度才退化至 $12.1\text{ dB}$。这证明了容量矩阵求解在物理层面上具有极佳的韧性。

3. 代码实现细节，复现指南，所用的软件包及开源 repo link

为了方便科研人员复现此项工作，我们在此梳理整个 LGF 离散 Calderón 投影算法的软件架构与复现步骤。本算法主要在 Python 3 中实现，依赖于标准的科学计算库 NumPy、SciPy 以及专用的格点格林函数计算模块。

3.1 核心算法模块架构

算法复现的核心流程可划分为以下四个步骤：

[1. LGF 核函数预计算] -> [2. 几何水平集建立与边界带提取] -> [3. 容量矩阵装配与稠密求解] -> [4. 主动源计算与全场重构]

步骤 1：LGF 核函数的计算复现

格点格林函数（等式 2）在原点附近具有对数级（2D）或代数级（3D）奇异性。该工作采用 Wang 和 Xia 提出的 DST-with-Hankel-correction 方法 [22] 进行极高精度的快速预计算。该方法的核心思想是：

• 将 LGF 分解为自由空间连续 Green 函数 $G_{\text{free}}(\mathbf{x}) = \frac{i}{4} H_0^{(1)}(k |\mathbf{x}|)$（用以吸附奇异性）和格点色散修正项之和。
• 采用离散正弦变换（DST）高效求解边界校正方程。

在复现中，我们可以将预先计算好的 LGF 在一个包含 $2n \times 2n$ 节点的相对坐标矩阵中进行 Tabulation（制表），因为对于平移不变形介质，矩阵装配只需要通过简单的查表操作即可完成：G_look_up[dx, dy]。

步骤 2：边界带的代数扫描算法

我们利用隐式水平集函数（Level-set Function）$\phi(\mathbf{x})$ 来描述保护区。例如对圆形，$\phi(\mathbf{x}) = 0.5 - \sqrt{x^2+y^2}$。在格点 $\mathbf{x}_{i,j} = (ih, jh)$ 上：

1. 若 $\phi(\mathbf{x}_{i,j}) \ge 0$，标记 M_plus[i,j] = True；反之标记为 False。
2. 遍历整个网格。若某一节点 (i,j) 本身属于 M_plus，但其 5 点模板（上下左右）中至少有一个邻居不属于 M_plus，则将该点加入内边界层集合 gamma_plus。
3. 若某一节点 (i,j) 本身不属于 M_plus，但其模板中至少有一个邻居属于 M_plus，则将其加入外边界层集合 gamma_minus。
4. 合并两层得到 gamma = gamma_plus + gamma_minus。

步骤 3：容量矩阵（Capacity Matrix）的装配与求解

对于边界点集 gamma，其点数为 $N_\gamma$。装配大小为 $N_\gamma \times N_\gamma$ 的复数稠密矩阵 $S_{\gamma\gamma}$：

import numpy as np

N_gamma = len(gamma)
S = np.zeros((N_gamma, N_gamma), dtype=np.complex128)
for i in range(N_gamma):
    for j in range(N_gamma):
        dx = gamma[i][0] - gamma[j][0]
        dy = gamma[i][1] - gamma[j][1]
        S[i, j] = h**2 * get_LGF_value(dx, dy)

在获取实测总迹 u_trace 后，求解线性方程组 lam = np.linalg.solve(S, u_trace)。之后，直接将 lam 按照索引切分出 lam_minus（对应 gamma_minus 的格点），即得到控制源强度 $g^h_{\text{sh}} = -\lambda^-$。

3.2 推荐软件包与开源链接

1. 格点格林函数计算组件：
   • 推荐使用 LatticeGreenFunctions 或 pyLGF 相关的开源第三方库，用以高效生成满足离散色散关系的格林函数数值表。
2. 线性方程组求解器：
   • 对于边界点规模 $N_\gamma \approx 300$ 的体系，直接调用 scipy.linalg.lu_factor 进行 LU 分解并在不同入射场下进行 lu_solve 重复计算最为高效。
   • 对于未来可能拓展的三维大规模边界带，推荐使用 PyAMG（Python 代数多网格库）或 SciPy 内部自带的 GMRES 迭代求解器。

4. 关键引用文献，以及你对这项工作局限性的评论

4.1 关键引用文献

在主动控制与差分势方法（DPM）的发展史中，以下四篇文献构成了本研究的关键理论先导：

1. Lončarić, Ryaben’kii, and Tsynkov (2001) [6]
   • Active shielding and control of noise, SIAM Journal on Applied Mathematics.
   • 学术地位：首次从数学上严格阐明了如何利用 Calderón 投影与广义面势理论将主动声学控制表述为一个边界边界算子问题，奠定了 DPM 在消声领域的基础。
2. Tsynkov (2003) [7]
   • On the definition of surface potentials for finite-difference operators, Journal of Scientific Computing.
   • 学术地位：给出了差分 Helmholtz 算子下离散面势和差分势重构的完备代数定义，但其实现依然深度受限于有界辅助区域算子的求逆。
3. Utyuzhnikov, Hu, and Zhou (2020-2024) [18, 20, 21]
   • Nonlocal active sound control (NASC) schemes.
   • 学术地位：开发了非局部双层 Huygens 控制面消声技术。该工作提出的内/外边界层单胞跨度结构，即是 NASC 思想在最紧凑笛卡尔晶格差分层面上的完美数学呈现。
4. Wang and Xia (2026) [22]
   • Unfitted lattice Green’s function method for exterior scattering in complex geometry.
   • 学术地位：本文所使用的 LGF 快速高精度生成算法的技术源头，解决了高波数下色散校正项的极速收敛计算难题。

4.2 对该项工作局限性与潜在改进空间的批判性评论

尽管该论文在理论代数结构上表现得近乎完美，且在二维体系中取得了极其亮眼的表现，但从一位严谨的计算物理与量子化学模拟专家的视角来看，该方法在实际走向三维大规模计算及复杂多物理场耦合时，依然存在以下不可忽视的技术局限性：

1. 三维格林函数查表法（Tabulation）的内存瓶颈

在三维格点空间中，随着网格步长 $h$ 减小以解析高频声波，格点数将以 $O(n^3)$ 爆炸式增长。一个 $256 \times 256 \times 256$ 的三维晶格，其对应的 LGF 查找表将耗费巨大的复数存储空间（约数百吉字节，GB）。在没有进行低秩压缩（Low-rank Compression）或者张量积近似（Tensor-product Approximation）的前提下，内存带宽限制将直接摧毁该方法的执行效率。未来必须引入**层次矩阵（$\mathcal{H}$-matrices）或者快速多极子方法（FMM）**对三维 LGF 卷积进行加速。

2. 对背景介质均匀性与笛卡尔网格的强依赖性

该方法之所以能够通过单层容量矩阵 $S_{\gamma\gamma}$ 实现代数闭合，完全建立在 LGF 的空间平移不变性（Translation Invariance）基础之上。这意味着保护区边界带 $\gamma$ 邻域内必须是绝对均匀的介质。如果控制系统面临的是一个非均匀流体介质（背景声速剧烈空间变化，或者存在大范围剪切流），LGF 将不再适用，必须退回到传统的变系数有界变分求解中。此外，由于采用的是硬阶梯（Staircase）笛卡尔网格来逼近连续的光滑曲线边界，在强高频极限下，阶梯折射引起的散射散射误差将不可避免地降低收敛精度，亟需发展贴体不规则网格（Body-fitted curvilinear grids）下的 LGF 控制理论。

3. 内部离散 Dirichlet 共振条件（Resonance Sensitivity）的病态化

定理 3.2 成立的重要前提是假设波数平方 $k^2$ 不是受保护区域 $M^+$ 内部带有同次 Dirichlet 条件的离散 Laplacian 算子的特征值。然而，当受保护区域很大且波数较高时，内部 Dirichlet 特征值的分布密度将变得极度稠密。一旦 $k^2$ 逼近某一特征值（即发生声学共振），容量矩阵 $S_{\gamma\gamma}$ 的最小奇异值 $\sigma_{\min}$ 将趋近于 0，导致条件数 $\kappa(S_{\gamma\gamma})$ 暴增。正如表 7 所示，在 $p=7$ 的网格上条件数陡增至 361。这将极大地放大实测噪声，导致控制源强度 $\lambda$ 的绝对值剧烈发散。论文中尚缺乏一套结合吉洪诺夫正则化（Tikhonov Regularization）或奇异值截断（TSVD）的稳健控制合成框架。

5. 其他必要的技术补充：与量子化学紧束缚模型及非平衡格林函数（NEGF）的深刻学术类比

对于量子化学和凝聚态物理模拟的研究人员来说，“格点格林函数（LGF）”与“边界自能投影”等概念并不陌生。实际上，本篇声学论文中所展现出的全部代数精髓，能够与微观体系中的分子电导输运理论、**紧束缚哈密顿量（Tight-Binding Hamiltonian）以及非平衡格林函数（NEGF）**建立起一一对应的、令人惊叹的对称映射。

5.1 数学映射：从离散 Helmholtz 方程到一维/二维晶格紧束缚模型

在量子化学中，用于计算一维分子纳米线或二维石墨烯输运性质的紧束缚模型，其定态薛定谔方程为：

在最简单的单轨道正交基下，该 Hamiltonian 作用在格点 $m$ 上，表现为对相邻轨道的跳跃积分（Hopping Integral）$t$ 以及格点现场能（Site Energy）$\epsilon_0$：

我们对该方程进行移项整理：

将上式与本篇论文中的离散 Helmholtz 算子（等式 1）进行对比：

这一高度对称性表明：主动噪声屏蔽中通过容量矩阵 $S_{\gamma\gamma}$ 合成边界控制源的过程，其物理本质上等价于在量子输运中，通过求解界面格林函数，从而在分子-电极界面上构建完全等效的向外辐射自能算子（Boundary Self-Energy, $\Sigma(E)$）的过程。

5.2 边界自能与 Calderón 投影的等价性

在经典的 Landauer-Büttiker 量子输运模拟中，为了处理无限半无限大电极（Semi-infinite Leads）对中心有界器件区（Device Region, $M^+$）的影响，化学家会计算半无限大电极在接触面 $\gamma^-$ 上的表面格林函数（Surface Green’s Function） $g_{\text{lead}}$。这对应于声学屏蔽中的 $S_{--}^{-1}$。通过自能矩阵：

将外部无限空间的波函数向外辐射特性完全折叠压缩到器件区的接触边界层上。这与本篇论文中 Theorem 4.5 指出的“通过外边界单层密度 $\lambda^-$ 完全代替外部介质散射作用”在物理和数学代数结构上完全一致。声学屏蔽中的单边测量重构矩阵 $T^-_h = S_{+-} S_{--}^{-1}$，在量子力学中对应的就是将电极格点的波函数振幅投影传导至器件界面的 Fisher-Lee 关系式（Fisher-Lee Relation）。

5.3 跨学科启示：量子波包的主动控制与量子相干屏蔽

这一跨学科的对应关系为量子化学和纳米电子学模拟带来了极具启发性的创新思路：

1. 量子主动波包阻绝器（Active Quantum Wave Barrier）：在分子尺度器件设计中，我们是否可以借鉴声学主动屏蔽的公式 (9)，通过在外加栅极（Gate Electrodes）上施加特定的、与相干入射电子相位精确对齐的动态局域势场 $g^h_{\text{sh}}$，从而实现对特定能量区间 $E$ 入射电子波函数的“完美无反射相干屏蔽”？这将允许我们在亚纳米尺度上构建出近乎理想的量子晶体管关断机制。
2. 相干量子围护（Quantum Confinement by Design）：利用等式 (10) 的围护控制理论，我们可以在人工量子点（Quantum Dots）的内层表面 $\gamma^+$ 上布置极其精密的辅助原子或外加电极。通过施加计算所得的特定跳跃能调制，在不破坏量子点内部自发本征态演化的前提下，将电子波包物理禁锢在受控区域内，杜绝其向外部基底介质的隧穿和耗散。这对于延长超导量子比特（Superconducting Qubits）的相干寿命，具有深刻的工程指导意义。

5.4 总结

庆夏在这项工作中所展现出的精妙代数结构，不仅为经典计算声学中的边界势方法指明了无需外边界重构的高效发展路径，其高度抽象和普适的无限格点算子求逆理论，更是横跨了宏观声波控制与微观薛定谔波包演化的物理边界。这一将“体空间物理控制退化为单层晶格代数解”的方法论，无疑将为未来多物理场耦合模拟及先进微纳结构中的波功能器件设计，注入全新的生命力。