2025-11-22T05:04:16.626039

On the injective norm of random fermionic states and skew-symmetric tensors

Dartois, Radpay

We study the injective norm of random skew-symmetric tensors and the associated fermionic quantum states, a natural measure of multipartite entanglement for systems of indistinguishable particles. Extending recent advances on random quantum states, we analyze both real and complex skew-symmetric Gaussian ensembles in two asymptotic regimes: fixed particle number with increasing one-particle Hilbert space dimension, and joint scaling with fixed filling fraction. Using the Kac--Rice formula on the Grassmann manifold, we derive high-probability upper bounds on the injective norm and establish sharp asymptotics in both regimes. Interestingly, a duality relation under particle--hole transformation is uncovered, revealing a symmetry of the injective norm under the action of the Hodge star operator. We complement our analytical results with numerical simulations for low fermion numbers, which match the predicted bounds.

academic

On the injective norm of random fermionic states and skew-symmetric tensors

基本信息

论文ID: 2510.25474
标题: On the injective norm of random fermionic states and skew-symmetric tensors
作者: Stéphane Dartois (École Polytechnique), Parham Radpay (Université Paris-Saclay, CEA)
分类: math.PR (概率论), math-ph (数学物理), math.MP (数学物理), quant-ph (量子物理)
发表时间: 2025年10月29日 (arXiv预印本)
论文链接: https://arxiv.org/abs/2510.25474v1

摘要

本文研究随机反对称张量的内射范数及其相关的费米子量子态，这是不可区分粒子系统多体纠缠的自然度量。扩展了关于随机量子态的最新进展，作者分析了实数和复数反对称高斯系综在两个渐近区域的行为：固定粒子数且单粒子希尔伯特空间维度趋于无穷，以及保持固定填充分数的双尺度极限。利用Grassmann流形上的Kac-Rice公式，推导出内射范数的高概率上界，并在两种区域建立了尖锐的渐近行为。有趣的是，发现了粒子-空穴变换下的对偶关系，揭示了Hodge星算子作用下内射范数的对称性。数值模拟验证了低费米子数情况下的理论预测。

研究背景与动机

要解决的问题

本文研究随机费米子态的几何纠缠度量问题，具体关注反对称张量的内射范数及其渐近行为。这是量子信息理论中多体纠缠度量的核心问题，特别是针对不可区分粒子系统。

问题的重要性

量子信息理论基础：几何纠缠度量（Geometric Measure of Entanglement, GME）是量子信息中刻画多体纠缠的重要工具，对理解量子计算和量子通信至关重要。
费米子系统的特殊性：对于费米子（遵循泡利不相容原理的粒子），由于粒子不可区分性，传统的纠缠度量需要修正。反对称张量自然刻画了费米子态。
随机矩阵理论的推广：内射范数是随机矩阵算子范数在张量情况的推广，研究其在不同对称类下的行为是随机张量理论的重要方向。
与无序系统的联系：这一问题与自旋玻璃等无序系统的基态研究密切相关，具有统计物理意义。

现有方法的局限性

对称性类别不完整：现有随机张量研究主要集中在对称或无对称性张量，对反对称张量（对应Sp的符号表示）研究较少。
技术挑战：反对称张量的结构零元素使得直接应用现有技术（如Subag的TAP方法）困难。
费米子纠缠的特殊性：严格可分离的反对称态不存在，需要重新定义可分离态集合（使用楔积）。

研究动机

本文旨在：

将DM24关于可区分粒子随机态的分析推广到费米子情况
系统研究反对称高斯张量系综的内射范数
建立两种渐近区域（固定p、d→∞和双尺度p,d→∞）的精确上界
揭示粒子-空穴对偶性等深层对称性

核心贡献

建立了反对称张量内射范数的尖锐渐近上界：在实数和复数情况下，针对两种渐近区域（固定费米子数p与d→∞；双尺度p=⌊αd⌋, d→∞）得到高概率上界（定理1.2）。
引入了新的随机矩阵系综：定义了Block Hollow Gaussian Antisymmetric Ensemble (BHGAE)及其复数版本(cBHGAE)，证明了其算子范数界（命题2.8, 2.10）。
发现了粒子-空穴对偶性：证明了费米子态的几何纠缠在Hodge对偶（粒子-空穴变换）下保持不变（命题2.5），这在双尺度极限下表现为α ↔ 1-α的对称性。
在Grassmann流形上应用Kac-Rice公式：这是该技术的新应用场景，克服了反对称张量结构零带来的技术困难。
数值验证：对p=2,3,4的情况进行数值模拟，验证了理论预测的上界。

方法详解

任务定义

输入：随机反对称张量 $T \in \bigwedge^p \mathbb{K}^d$ （ $\mathbb{K}=\mathbb{R}$ 或 $\mathbb{C}$ ），其独立元素服从标准高斯分布。

输出：内射范数的渐近行为 $\|T\|_{inj} := \max_{x_1,\ldots,x_p \in \mathbb{K}^d, \|x_1\wedge\cdots\wedge x_p\|_2=1} |\langle T, x_1\wedge\cdots\wedge x_p\rangle|$

关键约束：

反对称性： $T_{i_{\sigma(1)},\ldots,i_{\sigma(p)}} = \text{sg}(\sigma)T_{i_1,\ldots,i_p}$
归一化费米子态： $|\psi_f\rangle = T/\|T\|_2$
优化域为Grassmann流形 $\text{Gr}(p,d)$

模型架构

1. 内射范数的等价刻画

利用Grassmann流形的几何结构，将优化问题转化为： $\|T\|_{inj} = \max_{x^{(1)}\wedge\cdots\wedge x^{(p)} \in \text{Gr}(p,d)} \frac{|\langle T, x^{(1)}\wedge\cdots\wedge x^{(p)}\rangle|}{\sqrt{|\det(X^tX)|}}$ 其中 $X$ 是以 $x^{(i)}$ 为列的 $d\times p$ 矩阵。

2. Kac-Rice公式应用

定义随机场： $f_T(x^{(1)}\wedge\cdots\wedge x^{(p)}) = \frac{1}{\sqrt{p(d-p)}} \frac{\langle T, x^{(1)}\wedge\cdots\wedge x^{(p)}\rangle}{\|x^{(1)}\wedge\cdots\wedge x^{(p)}\|_2}$

应用Kac-Rice公式估计临界点数： $\mathbb{E}[\text{Crt}_{f_T}(B)] \leq K(p,d) \int_B \sqrt{p(d-p)} e^{-\frac{p(d-p)}{2}u^2} \mathbb{E}_{\text{BHGAE}}[|\det(u-W_N)|] du$

其中 $W_N$ 是BHGAE随机矩阵， $N=p(d-p)$ 。

3. Block Hollow Gaussian Antisymmetric Ensemble (BHGAE)

定义 $N\times N$ 随机矩阵 $W_N$ ，分块结构为 $p\times p$ 块，每块大小 $d\times d$ ：

对角块为零
非对角块 $(i,j)$ 为反对称高斯随机矩阵 $G_{ij}$ （方差 $\sigma^2=\frac{1}{p(d-p)}$ ）
满足块间反对称性： $(W_N)_{(i,\alpha),(j,\beta)} = -(W_N)_{(j,\beta),(i,\alpha)}$

关键性质：

算子范数界（命题2.8）：对固定p、d→∞， $\mathbb{P}\left(\|W_N\| > 2\sqrt{\frac{p-1}{p}} + \varepsilon\right) \leq e^{-N^{1-\delta}}$
经验谱测度收敛到半圆律： $\rho_p(x) = \frac{p}{2\pi(p-1)}\sqrt{4\frac{p-1}{p} - x^2}$

技术创新点

1. Grassmann流形上的正规坐标

在参考点 $P_0 = \begin{pmatrix} I_p \\ 0 \end{pmatrix}$ 附近，使用切空间 $\mathbb{R}^{(d-p)\times p}$ 的坐标 $B$ ： $\phi(B) = \tilde{B}(\tilde{B}^t\tilde{B})^{-1}\tilde{B}^t, \quad \tilde{B} = \begin{pmatrix} I_p \\ B \end{pmatrix}$

这使得可以在欧氏空间中计算协方差函数（附录B）。

2. 矩阵Dyson方程

利用Kronecker结构： $W_N = \frac{1}{\sqrt{2p(d-p)}} \sum_{1\leq i<j\leq p} \beta_{ij}^{(p)} \otimes G_{ij} + (\beta_{ij}^{(p)})^t \otimes G_{ij}^t$ 其中 $\beta_{ij}^{(p)} = E_{ij}^{(p)} - E_{ji}^{(p)}$ 。

应用AEKN19的局部律，得到Stieltjes变换： $m_p(z) = \frac{-z + \sqrt{z^2 - 4(p-1)/p}}{2(p-1)/p}$

3. Laplace方法与复杂度函数

定义复杂度函数（固定p情况）： $\Sigma_R^{(p)}(u) = \frac{1+\log p}{2} - \frac{u^2}{2} + \Omega\left(u\sqrt{\frac{p-1}{p}}\right) + \frac{1}{2}\log\frac{p-1}{p}$

其中 $\Omega$ 是半圆律的对数势。解 $\Sigma_R^{(p)}(E_0(p))=0$ 给出渐近上界。

双尺度情况（p=⌊αd⌋）： $\gamma_\alpha(d) = \sqrt{\log d} + \frac{\log\log d}{2\sqrt{\log d}} + \frac{\beta(\alpha)}{\sqrt{\log d}} + o\left(\frac{1}{\sqrt{\log d}}\right)$ 其中 $\beta(\alpha) = \frac{3}{4} + \frac{1}{4}\frac{\alpha}{1-\alpha}\log\alpha + \frac{1}{4}\frac{1-\alpha}{\alpha}\log(1-\alpha) + \frac{1}{2}\log(\alpha(1-\alpha))$

4. 粒子-空穴对偶性的体现

Hodge对偶 $*: \bigwedge^p \mathbb{R}^d \to \bigwedge^{d-p} \mathbb{R}^d$ 保持内射范数（命题2.5）。在双尺度极限中， $\beta(\alpha)$ 关于 $\alpha=1/2$ 对称（图1），体现了 $\alpha \leftrightarrow 1-\alpha$ 的对称性。

实验设置

数值方法

使用FLN22的方法进行数值验证：

p=2情况：内射范数等于对应反对称矩阵的最大奇异值，使用SVD精确计算。

p≥3情况：使用梯度下降法近似优化： $\max_{x_1,\ldots,x_p} \frac{|\langle T, x_1\wedge\cdots\wedge x_p\rangle|}{\|x_1\wedge\cdots\wedge x_p\|_2}$

参数设置

张量阶数：p = 2, 3, 4
维度范围：
- p=2: d ∈ 10, 500
- p=3: d ∈ 10, 350（实数），d ∈ 10, 300（复数）
- p=4: d ∈ 10, 120
样本数量：
- p=2: 200个样本（所有d）
- p=3: d≤250时100个样本，d>250时20个样本
- p=4: d≤70时100个样本，d>70时10个样本

观测量

归一化内射范数： $\frac{\|T\|_{inj}}{\sqrt{p(d-p)}}$

理论预测上界：

p=2: $\alpha(2) = 2$
p=3: $\alpha(3) \approx 2.870$
p=4: $\alpha(4) \approx 3.588$

实验结果

主要结果

1. 固定p、d→∞的收敛性

图2(a-c)展示了实数情况下p=2,3,4的结果：

p=2：数值结果快速收敛到理论值2，在d=500时几乎完全吻合
p=3：收敛到2.87，d=350时仍有小幅波动但趋势明确
p=4：收敛到3.588，由于计算限制仅到d=120，但趋势清晰

关键观察：

所有情况下归一化内射范数随d增长保持有界
收敛速度随p增大而减慢（需要更大的d达到渐近区域）
样本间方差随d增大显著减小（自平均效应）

2. 复数情况的验证

图2(d)展示p=3复数情况：

与实数情况收敛到相同的上界 2.87
证实了定理4.3：实数和复数情况在归一化后具有相同的渐近行为
计算复杂度更高，仅到d=300

3. 与其他对称类的比较

图3比较了p=3时不同对称性：

反对称：收敛到2.87
对称：显著更大的值（约3.2附近）
无对称性：介于两者之间

图3(b)显示比率随d的演化：

反对称/无对称性比率约0.85-0.90
反对称/对称比率约0.75-0.80
体现了对称性对内射范数的显著影响

实验发现

理论上界的紧性：数值结果表明理论上界是渐近紧的，没有系统性的差距。
收敛速度：
- p=2最快（d~100即达到渐近区域）
- p=3需要d~300
- p=4需要更大的d（>120）
自平均现象：高维情况下，不同样本间的方差明显减小，表明内射范数具有强自平均性质。
计算瓶颈：p≥3时梯度下降法的计算成本随d和p快速增长，限制了大规模验证。

结论与讨论

主要结论

尖锐上界的建立：
- 固定p、d→∞： $\|T\|_{inj}/\sqrt{p(d-p)} \lesssim E_0(p)$ （方程2.11的解）
- 双尺度p=⌊αd⌋、d→∞： $\|T\|_{inj}/\sqrt{d^2\alpha(1-\alpha)} \lesssim \gamma_\alpha(d) \sim \sqrt{\log d}$
对偶性的发现：几何纠缠在粒子-空穴变换下不变， $\text{GME}(|\psi_f\rangle) = \text{GME}(|\psi_h\rangle)$
实数与复数的统一：归一化后两者具有相同的渐近行为
数值验证：p=2,3,4的模拟确认了理论预测

局限性

仅有上界：与DM24不同，本文尚未建立匹配的下界。BS25的技术不直接适用于反对称情况。
结构零的挑战：BHGAE矩阵的对角块为零，使得谱间隙性质（性质3）仅在 $|u|$ 足够大时成立，需要条件事件 $\mathcal{E}_K$ 。
双尺度区域的技术限制：
- 仅证明了 $\gamma_\alpha(d) = o(\exp(d^{1-\delta}))$ 情况
- 对更一般序列的Laplace方法适用性未完全建立
数值计算的限制：
- p≥3时依赖梯度下降，可能陷入局部极值
- 高阶张量的内存和计算成本限制了大规模验证
对称类的不完整性：仅研究了完全反对称情况， $S_p$ 的其他不可约表示对应的张量对称类尚未涉及

未来方向

下界的建立：
- 使用第二矩方法或其他技术证明匹配的下界
- 确立 $E_0(p)$ 和 $\gamma_\alpha(d)$ 的最优性
其他对称类：
- 研究 $S_p$ 的一般Young表对应的张量
- 建立对称类与内射范数的系统关系
物理应用：
- 将结果应用于费米子系统的实际纠缠度量
- 连接到凝聚态物理中的拓扑序等问题
计算方法改进：
- 开发更高效的内射范数计算算法
- 扩展数值验证到更高阶和更大维度
谱测度的精细结构：
- 受DLSST25启发，研究反对称张量"本征值"分布的全局性质
- 建立与自由概率论的更深联系
量子纠错码的应用：
- 如DZ25所示，内射范数与CSS码的性能相关
- 探索费米子码的设计

深度评价

优点

1. 数学严谨性

完整的证明结构：从Kac-Rice公式到Laplace方法，逻辑链条清晰
技术细节充分：附录详细展示了Grassmann几何、协方差计算等
引用权威结果：恰当使用BABM22, BBvH23等最新随机矩阵理论成果

2. 创新性

新系综的引入：BHGAE及其复数版本填补了随机矩阵理论的空白
几何洞察：Grassmann流形上的分析揭示了问题的内在几何结构
对偶性的发现：粒子-空穴对称性是物理上深刻的结果，且在数学上自然体现

3. 跨学科价值

连接量子信息、随机矩阵、统计物理、代数几何等多个领域
为费米子系统的纠缠研究提供了坚实的数学基础
推广了自旋玻璃景观复杂度的方法论

4. 实验验证

数值模拟与理论预测高度一致
不同对称类的比较提供了额外洞察
实数/复数情况的对比确认了理论的统一性

不足

1. 上下界的不匹配

目前仅有上界，不能排除真实值更小的可能性
虽然数值结果支持紧性，但缺乏严格证明
相比DM24（后被BS25完成），本文在这方面不够完整

2. 技术假设的必要性

性质3（谱间隙）需要限制在 $u \in K \subset A$ （远离 $E_0(p)$ 的闭集）
这一技术妥协源于结构零，但削弱了结果的普遍性
未来需要更精细的技术处理这一问题

3. 双尺度区域的限制

$\gamma_\alpha(d)$ 的增长速度限制为 $o(\exp(d^{1-\delta}))$
对于增长更快的序列，Laplace方法的适用性未建立
这限制了对某些极端情况的理解

4. 数值实验的范围

仅验证了p≤4的情况
双尺度区域（固定α）的数值验证缺失
高维情况下的计算成本限制了验证的深度

5. 物理解释的深度

对 $\beta(\alpha)$ 函数的物理意义讨论不足
粒子-空穴对称性与物理系统的具体联系可以更详细
与实际费米子系统（如电子气、量子点）的连接未明确

影响力

学术贡献

理论推进：显著扩展了随机张量内射范数的理解，特别是非对称情况
方法论：Grassmann流形上的Kac-Rice公式应用可能启发其他问题
对称类分类：为系统研究所有Young表对应的张量对称类奠定基础

实用价值

量子信息：为费米子纠缠的实验测量提供理论基准
量子计算：可能应用于费米子量子算法的资源估计
凝聚态物理：帮助理解多体费米子系统的纠缠结构

可复现性

高：方法描述详细，数值实验参数明确
代码虽未公开，但FLN22的方法可复现
数学推导可验证（基于已发表的引理）

适用场景

直接适用

费米子量子态的纠缠度量：p个费米子在d维单粒子空间
反对称张量的统计性质：机器学习中的反对称结构
外微分形式的随机理论：微分几何中的应用

潜在扩展

拓扑量子计算：Majorana费米子的纠缠
量子化学：分子轨道的纠缠分析
统计学习：反对称核方法的理论
代数拓扑：随机单纯复形的性质

限制条件

需要高斯假设（i.i.d.标准正态分布）
渐近结果需要d足够大
固定p情况下p不能太大（需 $p \ll d$ ）

参考文献（精选）

核心引用

DM24 Dartois & McKenna: Injective norm of real and complex random tensors（本文直接推广的工作）
BABM22 Ben Arous, Bourgade & McKenna: Exponential growth of random determinants beyond invariance
BBvH23 Bandeira, Boedihardjo & van Handel: Matrix concentration inequalities and free probability

技术基础

AEKN19 Alt, Erdős, Krüger & Nemish: Location of the spectrum of Kronecker random matrices
ABAČ13 Auffinger, Ben Arous & Černý: Random matrices and complexity of spin glasses

量子信息

AS17 Aubrun & Szarek: Alice and Bob meet Banach
FLN22 Fitter, Lancien & Nechita: Estimating the entanglement of random multipartite quantum states

费米子系统

BHK21 Bianchi, Hackl & Kieburg: Page curve for fermionic Gaussian states
GKM12 Grabowski, Kuś & Marmo: Segre maps and entanglement for multipartite systems of indistinguishable particles

总体评价：这是一篇高质量的数学物理论文，在随机张量理论和量子信息的交叉领域做出了实质性贡献。技术上严谨，创新性强，且具有明确的物理意义。尽管存在上下界匹配等未完成的问题，但为后续研究奠定了坚实基础。推荐给对随机矩阵、量子纠缠或统计物理感兴趣的研究者。