精读笔记

Problem Setting

论文标题：Lehner's operator norm formulas, semidefinite programming, and spiked matrix models（arXiv preprint / 2026）。

这篇论文处理的不是一个新的随机矩阵模型本身，而是 Lehner 公式作为工具的“可计算性与可解释性”问题。Lehner 已经把 A0⊗1+∑Ai⊗mi 的谱边写成有限维正定矩阵优化；但这个优化仍然是非线性的，半圆情形虽可看作凸结构，原形式却不显式给出 dual geometry，Rademacher 情形甚至原目标非凸。

真正困难点在于：自由概率给出了谱边的 variational formula，但随机矩阵应用通常需要更强的东西——稳定计算、上下界证书、灵敏度、以及对 eigenvector 的二阶信息。原始 Lehner 公式主要给 value，不自然给 certificate。已有 BBP 应用中，上界来自构造 Z，下界则通过排除更小 Z 的反证完成，证明逻辑不对称，也难以迁移到 fluctuation 问题。

关键矛盾是：自由算子谱边是无限维对象，但可计算公式是有限维非线性对象；如果不能把后者变成凸优化，就很难系统地利用 duality。本文的核心就是把这个矛盾改写为 SDP primal-dual geometry。

Motivation

已有路线不够的地方有两层。

第一，Lehner 公式本身看起来已经“finite-dimensional”，但还不是一个好的 computational primitive。原文的迭代数值方法沿着 K(Z)=cI 的隐式曲线移动，接近最优时 derivative 可能奇异；半圆情形还能忍，Rademacher 情形的非凸性则更根本地阻碍了稳定求解。

第二，BCSvH24 等工作把 Lehner 公式用于 Gaussian matrix series 和 spiked models，但主要利用的是谱边值。对于 spiked model，最有统计意义的是 leading eigenvector 的行为；而 eigenvector second moment / covariance 在原始 Lehner 公式中并不显式。作者看到的缺口是：如果能获得 Lehner optimization 的 dual optimizer，它按凸优化 sensitivity 原理应当对应谱边对 A0 的导数，而这个导数在有限维随机矩阵中正是 leading eigenprojection 的期望。

所以动机不是“换一种写法更漂亮”，而是：把 Lehner 公式变成 SDP 后，value、certificate、sensitivity、numerics 和 eigenvector fluctuation 都进入同一个框架。

Core Idea

核心思想是把 Lehner 公式里的矩阵非线性项视为 non-commutative perspective，然后用 Schur complement / auxiliary PSD variables 进行精确 SDP 化。半圆情形的 f(t)=t^2 是 operator convex，AiZ^{-1}Ai 可直接被 block PSD 约束线性化；Rademacher 情形的 frad(t)=(-1+√(1+4t^2))/2 不是 operator convex，原目标确实非凸，但仍可通过引入 Yi 和 Riccati-like 结构得到精确 SDP。

本质区别在于，prior 把 Lehner 公式当成一个 variational evaluation formula；本文把它改造成一个 primal-dual certificate machine。这个转变带来的 inductive bias 是“谱边由 PSD cone 的支撑函数/对偶变量刻画”，而不是由某个隐式 fixed point 单独刻画。对随机矩阵应用而言，这重新组织了信息流：primal Z 控制 eigenvalue upper bound，dual S 控制 lower bound，同时 S 还携带 eigenvector second moment 的信息。

Method

半圆 small SDP：从 Λsc=inf_Z λmax(A0+Z^{-1}+Φ(Z)) 出发，用 Schur complement 得到 [[cI-A0-Φ(Z), I],[I,Z]]⪰0。它解决的是非线性 inverse 项不可直接对偶的问题。核心变化是把 Z^{-1} 从目标中移出，变成 PSD cone 的线性约束，从而暴露 dual：max_S,T ⟨A0,S⟩+2ReTr(T)，约束 [[S,-T],[-T*,Φ(S)]]⪰0, TrS=1。

半圆 large SDP：从原 Lehner 形式 A0+Z+∑AiZ^{-1}Ai 出发，对 ∑AiZ^{-1}Ai 整体做 block Schur complement。它解决 small SDP 可能不取到 primal optimum 的边界问题，同时给出 Λsc 作为固定 feasible set 的 support function 的解释。代价是维度扩大到 d(n+1)。

互补松弛：核心关系是 S=ZΦ(S)Z，以及在满秩情况下 A0+Φ(Z)+Z^{-1}=cI。这不是形式推导，而是解释 Lehner 原公式中“可限制到 scalar multiple of identity”现象的一部分。更重要的是，作者把它量化为近最优性 lemma：给定 feasible Z，只要找到 ΨZ(S)=ρS 的 PSD eigenmatrix，就能自动得到 lower bound。

Rademacher SDP：先确认 Lehner Rademacher 目标非凸，避免误以为只是 operator convex perspective 的直接应用。然后用约束 Z⪰0, [[Yi,Ai],[Ai,Yi+Z]]⪰0, Z+A0+∑Yi⪯cI 给出精确 primal SDP。证明精确性有两种机制：一个用自由 Rademacher 在 free product Hilbert space 上的 projection decomposition；另一个用 P=A(P+W)^{-1}A 的 Riccati-like characterization 和单调迭代。

Spiked model：对 A0=θvv^T, Φ(I)=I 的 isotropic Gaussian noise，构造 Z=αvv^T+β(I-vv^T) 给上界；下界用 ΨZ 的 leading PSD eigenmatrix 生成 Sest。这里 Sest 不只是证明工具，而被解释为 dual optimizer 的近似，进一步预测 Cov(\hat v)。

Key Insight / Why It Works

最重要的 insight 是：Lehner 公式隐藏的是 SDP geometry，而不是单纯的 nonlinear fixed point。半圆情形中，这几乎由 Schur complement 决定；真正有意思的是作者系统地开发了 dual side，并把 dual variable 解释成 sensitivity。这个 sensitivity 连接到 eigenvector second moment，是本文最值得迁移的想法。

方法有效的原因有三点。

第一，Φ(S)=∑AiSAi 是 completely positive map，这保证 PSD cone 与 eigenmatrix 构造相容。Lemma 4.7 之所以能工作，是因为互补松弛 S=ZΦ(S)Z 可以被看作 CP operator ΨZ 的 Perron-Frobenius 型 eigenproblem。也就是说，下界证书不是硬猜出来的，而是由 primal Z 自动诱导出来的。

第二，SDP dual 把“谱边对 A0 的导数”显式化。对于有限维 X=A0+W，∂λmax/∂A0=\hat v\hat v^T；如果 Eλmax≈Λsc 且可交换微分，则 dual optimizer S⋆≈E[\hat v\hat v^T]。这个逻辑不是定理，但非常有解释力，也解释了为什么 eigenvector covariance 会被 Lehner dual 捕捉。

第三，Rademacher 部分的实质贡献是证明非凸原问题仍有精确 SDP lift。这里不是 better engineering，而是隐藏 convexification。frad 不是 operator convex，按常规 non-commutative perspective 路线会失败；作者绕过失败点，利用特定函数满足的 Riccati 方程和自由 Rademacher 的投影结构。这是本文最不显然的技术点。

哪些可能只是辅助？数值稳定性讨论主要是 SDP 相比 Lehner 原迭代的自然收益，属于 convex optimization engineering，但合理。BBP eigenvalue limit 的新证明不是结果层面突破，因为已有 BCSvH24；它的价值在于证书结构，而非更强 theorem。eigenvector fluctuation 的数值实验是 hypothesis generation，不应被当作已验证理论。

Relation To Prior Work

最接近的谱系有三条。

第一是 Lehner 1999 及 operator-valued free probability 中关于 free operator norm 的计算。本文不是改进自由概率极限定理，而是对 Lehner variational formula 做 convex-analytic reorganization。实质新增信息是 SDP primal/dual、互补松弛、Rademacher exact lift。

第二是 BB21、BBvH23、BCSvH24 关于 Gaussian matrix series 与自由概率比较的路线。那些工作把 Λsc 当成 deterministic equivalent，用于 matrix concentration 和 spiked model。本文与它们的差异在于：不只估计 value，而是把 Λsc 的 dual optimizer 作为统计对象。BBP eigenvalue transition 本身不是新结论，新的是 dual certificate proof 和 eigenvector fluctuation conjecture。

第三是 SDP / matrix convexity / non-commutative perspective。半圆 small SDP 从技术上看是 Schur complement 的标准操作，不应过度神化；但把它与 Lehner 公式、free probability 和 spiked eigenvector sensitivity 串起来是新的。Rademacher SDP 更实质，因为它越过了 non-convex objective 的障碍。

看似新的部分中，半圆 SDP 本身可能已有相近想法，文中也提到 van Handel 私下得到过 dual form；真正新增的是完整 primal-dual 对、large SDP、近最优性机制、Rademacher 精确 SDP，以及把 dual optimizer 用于 eigenvector fluctuation 的解释。

Dataset / Evaluation

这篇不是 ML 论文，没有 dataset 意义上的评测。evaluation 分两类。

理论 evaluation：主要结论是精确等价和新证明。SDP reformulation 由定理支撑，强度很高；BBP eigenvalue limit 的新证明验证了 certificate framework 至少能重现已有重要应用。

数值 evaluation：作者比较 SDP solver 与 Lehner 原迭代，结论是 SDP 更稳定、更准确；这符合预期，但文中没有系统复杂度/规模实验，不能据此断言大规模上更优。Rademacher SDP 变量随 n 增加，实际可扩展性未充分说明。

Eigenvector fluctuation 部分的实验覆盖了 GOE、独立非均匀方差、平滑方差 profile、相关 Gaussian noise，能说明 conjecture 不只是 GOE symmetry 的产物。实验真正支持的是：anisotropic covariance 可由 Sest 预测，有限维投影近似 Gaussian。但它没有验证高维极限的严格条件，也没有覆盖非 Gaussian、强局部化噪声或 σ∗不消失场景。核心 claim 中 SDP 等价已由理论支持；eigenvector fluctuation claim 仍只是数值证据。

Limitation

最核心的限制是 eigenvector fluctuation 理论尚未闭合。S⋆≈E[\hat v\hat v^T] 的推理依赖 value approximation 与 differentiation 的交换，文中未充分说明；这一步在随机矩阵中可能需要 eigenvalue simplicity、uniform differentiability、误差对 perturbation 的稳定控制等条件。当前只是合理 heuristic。

第二，Sest 与真正 dual optimizer S⋆的关系没有一般误差界。Definition 7.8 构造的是 CP operator leading eigenmatrix，能用于 BBP 下界，但把它当 covariance predictor 仍缺少理论保证。若 dual optimizer 不唯一，sensitivity 解释也会变弱；文中给了一些 uniqueness 条件，但不是普适。

第三，方法的可扩展性存在上限。small semicircular SDP 变量是 d×d，尚可；large SDP 是 d(n+1) 维；Rademacher SDP 有 n 个 Yi，n 大时可能非常重。论文说 SDP 比 Lehner 迭代稳定，但这主要是小规模数值和 convex solver 的优势，不能自动推广到大规模 operator series 或 graph lift 问题。

第四，BBP 应用依赖 isotropy Φ(I)=I 和 σ∗→0。这个条件本质上排除了沿 spike 方向有强 bias 的噪声；作者也指出 Gaussianity 可能在 σ∗=Ω(1) 时破裂。因此所谓非 universal anisotropic Gaussian fluctuation 的适用边界还不清楚。

第五，v 方向涨落没有被解释。作者明确说明需要 E⟨\hat v,v⟩ 的二阶项，而当前 SDP sensitivity 只自然给 E[\hat v\hat v^T] 的一阶投影信息。这说明 dual optimizer 并没有编码全部 eigenvector fluctuation 结构。

Takeaway

• 1. Lehner 公式应被看作 SDP value，而不是单纯的自由概率 nonlinear variational formula；一旦这样看，dual certificate 和 sensitivity 会自然出现。
• 2. 对 Gaussian matrix series / spiked models，dual optimizer 是比谱边本身更有统计信息的对象：它很可能编码 leading eigenvector 的 second moment 与 anisotropic fluctuations。
• 3. Rademacher 情形的精确 SDP lift 说明“原目标非凸”不等于“问题非凸”；某些自由算子范数公式可能存在隐藏的 convex representation，值得系统寻找。
• 4. 后续真正值得做的是把 conjectural sensitivity-to-fluctuation 链条定理化，并开发利用 Φ/CP map 结构的 scalable algorithms，而不是依赖通用 SDP solver。

一句话总结

这篇论文把 Lehner 自由算子范数公式从非线性 variational tool 重构为 SDP primal-dual certificate framework，并由此打开了用 dual optimizer 预测 spiked matrix eigenvector 二阶涨落的路线。