2025-11-11T08:10:09.207739

Weakly universal dynamical correlations between eigenvalues of large random matrices

Mallick, TÃ©llez, van Wijland

It was shown roughly thirty years ago that the density correlations of eigenvalues of large random matrices display a universal form, independent of most of the details of the distribution of the random matrix itself. We show that when the matrix elements evolve according to a Dyson Brownian motion, dynamical correlations retain a large degree of the universality found at equal times when expressed in terms of the characteristics of some partial differential equation in the complex plane.

academic

Weakly universal dynamical correlations between eigenvalues of large random matrices

基本信息

论文ID: 2511.05727
标题: Weakly universal dynamical correlations between eigenvalues of large random matrices
作者: Kirone Mallick, Gabriel Téllez, Frédéric van Wijland
分类: cond-mat.stat-mech cond-mat.dis-nn math-ph math.MP
发表时间: 2025年11月7日（arXiv预印本）
论文链接: https://arxiv.org/abs/2511.05727

摘要

大约三十年前的研究表明，大随机矩阵特征值的密度相关性呈现出一种普适形式，基本不依赖于随机矩阵分布的具体细节。本文证明了当矩阵元素按照Dyson布朗运动演化时，动力学相关性在用复平面中某些偏微分方程的特征线表示时，仍保持了等时相关性中发现的大部分普适性。

研究背景与动机

问题背景

静态普适性的发现：Brézin-Zee和Beenakker在三十年前发现，大随机矩阵特征值的静态密度相关性具有普适形式，几乎不依赖于势函数V的具体形式，这是随机矩阵理论的重要发现。
动力学扩展的需求：虽然静态情况下的普适性已被充分理解，但当矩阵元素具有动力学演化（如Dyson布朗运动）时，这种普适性是否仍然保持是一个重要的未解决问题。
理论挑战：动力学情况下的分析比静态情况复杂得多，因为需要处理时间相关的随机偏微分方程，而现有的正交多项式方法在动力学情况下难以直接应用。

研究动机

探索动力学演化下随机矩阵特征值相关性的普适性质
建立连接静态和动力学行为的理论框架
为理解强相互作用系统的动力学提供可解析的模型

核心贡献

建立了动力学相关性的普适表达式：证明了动力学resolvent相关性可以表示为： $\Gamma(z, z'; t) = -\frac{1}{\beta N^2} \frac{Q(\zeta)}{Q(z)} \frac{a^2 - \zeta z' + \sqrt{\zeta^2 - a^2}\sqrt{z'^2 - a^2}}{\sqrt{z^2 - a^2}\sqrt{z'^2 - a^2}(\zeta - z')^2}$
发现了特征线方程的关键作用：动力学演化由特征方程 $\frac{d\zeta}{dt} = Q(\zeta)\sqrt{\zeta^2 - a^2}$ 控制，其中 $\zeta(0) = z$ 。
揭示了"弱普适性"现象：虽然势函数V通过多项式Q和参数a显式出现，但相关性的整体结构仍保持普适形式。
开发了基于resolvent的宏观涨落理论：将McKean-Vlasov方程转化为含噪Burgers方程，为处理动力学问题提供了新的理论工具。

方法详解

任务定义

研究在Dyson布朗运动驱动下，大随机矩阵特征值的动力学密度相关性： $C(\lambda, \lambda'; t-t') = \langle\hat{\mu}(\lambda, t)\hat{\mu}(\lambda', t')\rangle - \mu_W(\lambda)\mu_W(\lambda')$

其中 $\hat{\mu}(\lambda, t) = \frac{1}{N}\sum_i \delta(\lambda - \lambda_i(t))$ 是经验密度， $\mu_W(\lambda)$ 是平衡密度。

理论框架

1. Dyson布朗运动

特征值演化方程： $\frac{d\lambda_i}{dt} = -V'(\lambda_i) + \frac{1}{N}\sum_{j \neq i} \frac{1}{\lambda_i - \lambda_j} + \sqrt{\frac{1}{\beta N}}\eta_i$

2. Dean-Kawasaki方程

经验密度的随机偏微分方程： $\partial_t \hat{\mu} = -\partial_\lambda j$ 其中电流包含确定性和随机部分： $j(\lambda, t) = j_{\text{det}}(\lambda, t) + \sqrt{\frac{2\hat{\mu}}{\beta N^2}}\xi(\lambda, t)$

3. Resolvent变换

引入Stieltjes变换： $\hat{G}(z, t) = \int d\lambda \frac{\hat{\mu}(\lambda, t)}{z - \lambda}$

转化为含噪Burgers方程： $\partial_t \hat{G} = -\partial_z J$ 其中： $J = -V'(z)\hat{G} + \frac{1}{2}\hat{G}^2 + \hat{P} + \Xi$

技术创新点

1. 非局域噪声的处理

噪声项的相关性： $\langle\Xi(z, t)\Xi(z', t')\rangle = -\frac{2}{\beta N^2}\delta(t-t')\partial_z\partial_{z'}\frac{\hat{G}(z,t) - \hat{G}(z',t)}{z-z'}$

2. 特征线方法

通过特征方程： $\frac{d\zeta}{dt} = Q(\zeta)\sqrt{\zeta^2 - a^2}$ 解决响应函数和相关函数。

3. 宏观涨落理论的应用

使用Janssen-De Dominicis路径积分形式主义，action为： $S[\bar{G}, G] = \int dtdz\bar{G}\left(\partial_t G + \partial_z[-V'(z)G + P + \frac{1}{2}G^2]\right) - \text{噪声项}$

实验设置

理论验证案例

1. 二次势 $V(\lambda) = \frac{\lambda^2}{2}$

$Q(\lambda) = 1$ , $a^2 = 2$
特征方程简化为 $\frac{d\xi}{dt} = 1$
可获得解析解

2. 四次势 $V(\lambda) = \frac{\lambda^4}{4}$

$Q(\lambda) = \lambda^2 + a^2/2$ , $a^2 = 2\sqrt{2}/3$
展现有限时间爆炸现象
多个松弛时间尺度

3. 高次势 $V(\lambda) = \frac{\lambda^{2n}}{2n}$

复杂的特征线流动模式
螺旋运动和振荡衰减

分析方法

复变量技术（ $\zeta = a\cosh\xi$ 变换）
留数定理计算积分
渐近分析（长时间行为）

实验结果

主要结果

1. 二次势的精确解

长时间渐近行为： $C(\lambda, \lambda'; t) \simeq \frac{1}{\pi^2\beta N^2} \frac{\lambda\lambda'}{\sqrt{2-\lambda^2}\sqrt{2-\lambda'^2}}e^{-t}$

完整表达式： $C(\lambda, \lambda'; t) = \frac{-4 - 4\lambda^2\lambda'^2 + \lambda\lambda'(7 + 2(\lambda^2 + \lambda'^2))\cosh(t) - 4(\lambda^2 + \lambda'^2 - 1)\cosh(2t) + \lambda\lambda'\cosh(3t)}{2\pi^2\beta N^2\sqrt{(2-\lambda^2)(2-\lambda'^2)}[\lambda^2 + \lambda'^2 - 1 - 2\lambda\lambda'\cosh(t) + \cosh(2t)]^2}$