[论文解读] Decay properties for functions of matrices over $C^*$-algebras

Michele Benzi, Paola Boito|arXiv (Cornell University)|Jul 26, 2013

Quantum many-body systems参考文献 16被引用 8

一句话总结

本文通过全纯函数演算，建立了属于 $C^*$-代数的带状矩阵与稀疏矩阵的解析函数的指数非对角线衰减界。关键贡献是将以往仅限于实数或复数矩阵的经典衰减结果，推广至非交换及无限维 $C^*$-代数，包括希尔伯特空间上有界算子、紧空间上的连续函数以及四元数，且给出了与矩阵规模 $n$ 无关的可计算衰减界。

ABSTRACT

We extend previous results on the exponential off-diagonal decay of the entries of analytic functions of banded and sparse matrices to the case where the matrix entries are elements of a $C^*$-algebra.

研究动机与目标

• 将已知的矩阵函数指数衰减结果从实数或复数矩阵推广至一般 $C^*$-代数取值的矩阵。
• 回应 P.-L. Giscard 提出的关于参数化与非交换设定下衰减性质的问题。
• 推导出 $f(A)$ 的非对角线条目范数的可计算、与 $n$ 无关的界，其中 $A$ 是 $C^*$-代数上带状或稀疏的厄米特矩阵，$f$ 为解析函数。
• 通过适当的代数嵌入，将理论扩展至非交换代数，如希尔伯特空间上有界算子与四元数。
• 为量子化学、控制理论与数值分析中涉及结构化矩阵函数的应用提供统一框架。

提出的方法

• 利用全纯函数演算定义 $f(A)$，其中 $f$ 为在包含 $A$ 的谱的区域上解析的函数，$A$ 为 $n\times n$ 厄米特矩阵，其条目属于 $C^*$-代数 $\mathcal{A}_0$。
• 运用谱论与 $C^*$-代数性质，特别是 $C^*$-恒等式 $\|a^*a\| = \|a\|^2$，推导 $f(A)$ 的矩阵条目范数界。
• 推导出形如 $\|f(A)_{ij}\| \leq C \cdot e^{-\gamma \cdot d(i,j)}$ 的衰减估计，其中 $d(i,j)$ 为稀疏结构中的图距离，常数 $C$ 与 $\gamma$ 与 $n$ 无关。
• 将四元数矩阵表示为 $2n \times 2n$ 复数块矩阵，以应用复 $C^*$-代数技术，但对允许的函数类施加限制。
• 通过切比雪夫展开的逼近理论验证衰减行为的数值表现，使用理论框架导出的界。
• 利用 $\mathcal{B}(\mathcal{H})$、$C(\mathcal{X})$ 与 $\mathbb{H}$ 上的矩阵数值示例，验证指数衰减并对比理论界。

实验结果

研究问题

• RQ1 能否将矩阵函数的指数非对角线衰减特性从复数矩阵推广至一般 $C^*$-代数取值的矩阵？
• RQ2 当 $A$ 的条目属于非交换或无限维 $C^*$-代数时，$f(A)$ 的衰减界是否仍保持一致且与矩阵规模 $n$ 无关？
• RQ3 如何在 $C^*$-代数框架下处理四元数矩阵函数？这会对允许的函数类施加何种限制？
• RQ4 全纯函数演算在推导 $C^*$-代数上结构化矩阵的可计算、与 $n$ 无关的衰减界中起何作用？
• RQ5 理论衰减界能否通过数值方法验证？其与实际矩阵条目中观察到的衰减相比如何？

主要发现

• 对于任意 $C^*$-代数（包括非交换与无限维）中带状或稀疏的厄米特矩阵的解析函数 $f(A)$，其非对角线条目均建立了指数衰减。
• 衰减界的形式为 $\|f(A)_{ij}\| \leq C \cdot e^{-\gamma \cdot d(i,j)}$，其中常数 $C$ 与 $\gamma$ 与 $n$ 无关，确保了在不同矩阵规模下的统一性。
• 对于 $50 \times 50$ 厄米特三对角四元数矩阵，$\exp(A)$ 与 $\log(A)$ 的指数衰减被数值确认，非对角线块的范数迅速衰减。
• 在 $A \in \mathbb{H}^{50 \times 50}$ 且条目为随机的情况下，$\exp(A)$ 的衰减符合理论界，第 $k$ 个非对角线块的范数随 $k$ 呈指数衰减。
• 对于 $20 \times 20$ 块矩阵，其块大小为 $5 \times 5$ 复矩阵条目，$\exp(A)$ 中块的谱范数随距对角线距离的增加而指数衰减，与理论预测一致。
• 使用 15 项切比雪夫展开近似 $f(A) = \ln(A + 5I)$ 时，逼近误差被控制在 $1.7509 \times 10^{-14}$ 以内，验证了函数演算方法的精度。