[论文解读] The Perturbative Resolvent Method: spectral densities of random matrix ensembles via perturbation theory

Wenping Cui, Jason W. Rocks+1 more|arXiv (Cornell University)|Dec 1, 2020

Random Matrices and Applications参考文献 5被引用 6

一句话总结

Perturbative Resolvent Method (PRM) 是一种新颖的微扰方法，通过分析随机线性方程组解对微小扰动的响应，在热力学极限下计算随机矩阵系综的谱密度。该方法利用零温腔体方法，通过响应矩阵计算格林函数，从而实现对 Wigner 半圆律、Marchenko-Pastur 律以及实矩阵和非厄米矩阵的圆律/椭圆律等关键定律的简洁推导。

ABSTRACT

We present a simple, perturbative approach for calculating spectral densities for random matrix ensembles in the thermodynamic limit we call the Perturbative Resolvent Method (PRM). The PRM is based on constructing a linear system of equations and calculating how the solutions to these equation change in response to a small perturbation using the zero-temperature cavity method. We illustrate the power of the method by providing simple analytic derivations of the Wigner Semi-circle Law for symmetric matrices, the Marchenko-Pastur Law for Wishart matrices, the spectral density for a product Wishart matrix composed of two square matrices, and the Circle and elliptic laws for real random matrices.

研究动机与目标

• 开发一种简单、微扰的框架，用于在热力学极限下计算随机矩阵系综的谱密度。
• 通过响应矩阵和格林函数，将零温腔体方法与谱密度计算联系起来。
• 使用统一且易于理解的方法，推广现有结果（如 Wigner 半圆律和 Marchenko-Pastur 律）。
• 将该方法扩展至非厄米矩阵和实随机矩阵，包括圆律和椭圆律的推导。
• 提供一种灵活且可解析处理的传统 RMT 技术（如复制法或图解法）的替代方案。

提出的方法

• 该方法构建一组随机线性方程，并分析其解在小的恒定扰动下的变化。
• 将响应矩阵定义为解变量对扰动的导数，证明其等于预解式（格林函数）。
• 应用零温腔体方法计算这些响应，将问题简化为求解平均响应矩阵和响应场的自洽方程。
• 对于许多系综，腔体方程可简化为求解单个关于响应的多项式方程，避免求解完整的平均场方程。
• 通过厄米化约化，将该方法推广至非厄米矩阵，有效将系统尺寸加倍，将非厄米问题映射为等价的厄米问题。
• 通过标准解析延拓技术，在虚频率趋于零的极限下，从预解式的虚部提取谱密度。

实验结果

研究问题

• RQ1 基于响应理论的微扰方法能否在热力学极限下精确导出随机矩阵系综的谱密度？
• RQ2 如何系统地应用零温腔体方法，推导出如 Wigner 半圆律和 Marchenko-Pastur 分布等已知谱律？
• RQ3 该方法能否扩展至非厄米矩阵和实非对称随机矩阵，包括圆律和椭圆律的推导？
• RQ4 响应矩阵在通过格林函数将微扰响应与谱密度联系起来的过程中起什么作用？
• RQ5 在何种条件下，腔体方法可简化为求解关于响应的单个多项式方程？

主要发现

• PRM 通过在腔体框架下分析微扰响应，成功推导出对称随机矩阵的 Wigner 半圆律。
• 该方法通过系统性地处理腔体方程，重现了 Wishart 矩阵的 Marchenko-Pastur 律。
• 对于由两个方阵构成的乘积 Wishart 矩阵，PRM 在热力学极限下给出了谱密度的闭式表达式。
• PRM 通过证明实 i.i.d. 随机矩阵的特征值在半径为 σ 的圆盘内均匀分布，提供了 Girko 圆律的简洁推导。
• 该方法推广至实非对称矩阵，导出了椭圆律，其特征值分布为 ρ(u,v) = 1/(πσ²(1−ζ²))，分布在椭圆 u²/(1+ζ)² + v²/(1−ζ)² ≤ σ² 内。
• PRM 建立了响应矩阵与格林函数之间的直接联系，使得在许多情况下无需求解完整的自洽方程即可高效计算谱密度。