[论文解读] Calculating Determinants of Block Matrices

Philip D. Powell|arXiv (Cornell University)|Dec 16, 2011

Scientific Research and Discoveries参考文献 7被引用 53

一句话总结

本文提出了一种新方法，通过递归地将 $N\times N$ 块矩阵的行列式简化为 $N$ 个较小的、经变换的块的行列式的乘积，利用广义舒尔补（generalized Schur complement）来计算其行列式。该方法可对大型、物理上相关的矩阵（如夸克物质中的 $48\times 48$ 矩阵）进行解析求解，通过系统性地消除索引，得到与已知结果一致的精确本征能级。

ABSTRACT

This paper presents a method for expressing the determinant of an N { imes} N complex block matrix in terms of its constituent blocks. The result allows one to reduce the determinant of a matrix with N^2 blocks to the product of the determinants of N distinct combinations of single blocks. This procedure proves useful in the analytic description of physical systems with multiple discrete variables, as it provides a systematic method for evaluating determinants which might otherwise be analytically intractable.

研究动机与目标

• 开发一种通用方法，用于计算 $N\times N$ 块矩阵的行列式，该矩阵包含 $N^2$ 个块，每个块大小为 $n\times n$，以简化原本在解析上难以处理的系统。
• 克服标准 $2\times 2$ 块矩阵行列式公式的局限性，后者对索引处理不对称，不适用于对称的物理分区结构。
• 提供一种系统性的递归过程，将大块矩阵的行列式简化为较小、经变换的块的行列式的乘积。
• 在复杂物理系统（如高密度夸克物质）中实现行列式的精确解析计算，其中多个离散索引（如颜色、味、狄拉克指数）非平凡地耦合。

提出的方法

• 该方法通过使用下三角辅助矩阵进行递归变换，将对角线以下的块归零，同时保持行列式的值不变。
• 在每一步 $k$，通过广义舒尔补更新块：$\bm{\alpha}^{(k)}_{ij} = \mathbf{S}_{ij} - \bm{\sigma}^{T}_{i,N-k+1}\tilde{\mathbf{S}}^{-1}_{k}\mathbf{s}_{N-k+1,j}$，其中 $\mathbf{s}_{ij}$ 和 $\bm{\sigma}^{T}_{ij}$ 是原始矩阵中块的向量。
• 该过程迭代地消除块矩阵的最后一直行和一列，使矩阵大小每步减少一个块，直到仅剩一个块。
• 最终的行列式是 $N$ 个结果块 $\bm{\alpha}^{(N-k)}_{kk}$ 的行列式的乘积，最后一个块即为原始 $\mathbf{S}_{NN}$。
• 该方法假设中间块可逆，这对依赖连续变量的矩阵通常成立，奇点为孤立点且测度为零。
• 该方法在来自夸克物质的 $48\times 48$ 矩阵上得到验证，成功将其简化为更少索引上的行列式计算，并得到精确的本征能级。

实验结果

研究问题

• RQ1 能否将一般 $N\times N$ 块矩阵的行列式表示为较小、经变换的块的行列式的乘积？
• RQ2 标准 $2\times 2$ 块矩阵行列式公式如何推广，以处理对称的、多索引的物理系统？
• RQ3 该方法能否在高密度夸克物质中出现的 $48\times 48$ 矩阵上实现行列式的解析简化，该矩阵包含多个离散索引？
• RQ4 实现此类简化的递归舒尔补变换的结构是什么？

主要发现

• 一个 $N\times N$ 块矩阵的行列式由 $\det(\mathbf{S}) = \prod_{k=1}^{N} \det(\bm{\alpha}^{(N-k)}_{kk})$ 给出，其中 $\bm{\alpha}^{(k)}$ 块通过广义舒尔补递归定义。
• 对于一个 $48\times 48$ 矩阵，若 $N=6$ 个块（每个 $8\times 8$），该方法将问题简化为仅在狄拉克和味索引上计算行列式，消除了颜色和其他量子数。
• 最终的行列式表达式为 $\det(\mathbf{S}) = \left[E+\sqrt{(E_k+\mu)^2+|\Delta|^2}\right]^8 \left[E+\sqrt{(E_k-\mu)^2+|\Delta|^2}\right]^8 \times \left[E-\sqrt{(E_k+\mu)^2+|\Delta|^2}\right]^8 \left[E-\sqrt{(E_k-\mu)^2+|\Delta|^2}\right]^8 \times (E+E_k+\mu)^4(E-E_k-\mu)^4(E+E_k-\mu)^4(E-E_k+\mu)^4$。
• 由 $\det(\mathbf{S}) = 0$ 推导出的本征能级为 $E_1 = |E_k + \mu|$（简并度 8），$E_2 = |E_k - \mu|$（简并度 8），$E_3 = \sqrt{(E_k+\mu)^2 + |\Delta|^2}$（简并度 16），$E_4 = \sqrt{(E_k-\mu)^2 + |\Delta|^2}$（简并度 16），与已知结果一致。
• 该方法成功消除了原始 48 个索引中的 6 个，将问题简化为仅在狄拉克和味自由度上的可处理行列式。