2025-11-22T17:31:15.427699

Empirical bounds for commuting dilations of free unitaries and the universal commuting dilation constant

Gerhold, Scherer, Shalit

For a tuple $T$ of Hilbert space operators, the 'commuting dilation constant' is the smallest number $c$ such that the operators of $T$ are a simultaneous compression of commuting normal operators of norm at most $c$. We present numerical experiments giving a strong indication that the commuting dilation constant of a pair of independent random $N{\times}N$ unitary matrices converges to $\sqrt2$ as $N \to \infty$ almost surely. Under the assumption that this is the case, we prove that the commuting dilation constant of an arbitrary pair of contractions is strictly smaller than $2$. Our experiments are based on a simple algorithm that we introduce for the purpose of computing dilation constants between tuples of matrices.

academic

Empirical bounds for commuting dilations of free unitaries and the universal commuting dilation constant

基本信息

论文ID: 2510.12540
标题: Empirical bounds for commuting dilations of free unitaries and the universal commuting dilation constant
作者: Malte Gerhold (University of Greifswald), Marcel Scherer (Technion), Orr Moshe Shalit (Technion)
分类: math.FA cs.NA math.NA math.OA
发表时间: 2025年10月14日提交到arXiv
论文链接: https://arxiv.org/abs/2510.12540

摘要

对于Hilbert空间算子元组 $T$ ，"交换扩张常数"是使得 $T$ 的算子同时为范数至多为 $c$ 的交换正规算子压缩的最小数 $c$ 。本文通过数值实验强有力地表明，独立随机 $N \times N$ 酉矩阵对的交换扩张常数几乎必然收敛到 $\sqrt{2}$ （当 $N \to \infty$ 时）。在此假设下，证明了任意收缩算子对的交换扩张常数严格小于2。实验基于作者引入的用于计算矩阵元组间扩张常数的简单算法。

研究背景与动机

要解决的问题：确定通用交换扩张常数 $C_2$ 的精确值，即使得每对收缩算子都是范数至多为 $c$ 的交换正规算子同时压缩的最小常数 $c$ 。
问题的重要性：扩张理论技术和矩阵凸集框架在量子信息、优化控制、数学物理、算子系统和算子代数理论等领域有重要应用。
现有方法局限性：
- 已知界限： $\sqrt{d} \leq C_d \leq \sqrt{2d}$ （对于 $d$ 元组）
- 对于 $d=2$ 的情况，当前最佳上界 $C_2 \leq \sqrt{2} \cdot 2 = 2\sqrt{2}$ 与简单界限 $C_d \leq d$ 重合
- 下界 $C_d \geq \sqrt{d}$ 在一般情况下不紧致
研究动机：通过结合严格的极限定理与数值实验，提供 $C_2 < 2$ 的证据，改进现有上界。

核心贡献

数值证据：提供强有力的数值证据表明独立Haar酉矩阵对的扩张常数收敛到 $\sqrt{2}$
理论结果：证明了如果该收敛成立，则 $C_2 \leq \frac{2\sqrt{2}}{3} < 2$
算法贡献：引入了计算矩阵元组扩张常数的简单算法
严格界限：在附录中提供了有限维情况的严格界限 $C_2(n) \leq \sqrt{2} + 2\sin(\frac{\pi}{2}(1-\frac{1}{2n})) < 2$

方法详解

任务定义

给定算子 $d$ 元组 $u, v$ ，扩张常数 $c(u,v)$ 定义为： $c(u,v) = \inf\{c : u \prec cv\}$ 其中 $u \prec cv$ 表示存在Hilbert空间 $H \subseteq K$ 和算子元组 $U \in B(H)^d$ ， $V \in B(K)^d$ 使得 $u \sim U$ ， $v \sim V$ 且 $U = P_H cV|_H$ 。

核心理论框架

定理2.1（水平收敛与扩张常数）：如果算子 $d$ 元组序列 $(\xi^{(N)})$ 和 $(\eta^{(N)})$ 的矩阵范围水平收敛到 $\xi^{(\infty)}$ 和 $\eta^{(\infty)}$ ，则： $c(\xi^{(\infty)}, \eta^{(\infty)}) \leq \liminf_{N,M \to \infty} c(\xi^{(N)}, \eta^{(M)})$

推论2.4：对于具有强渐近自由性的随机矩阵集合 $T^{(N)}$ ，几乎必然有： $c(u^f, u_0) \leq \liminf c(T^{(N)}, u_0)$

算法设计

扩张常数计算算法：

构造交换正规矩阵 $N = (N_1, N_2)$ ，其对角线值遍历 $V_k \times V_k$ 的所有 $k^2$ 个可能对，其中 $V_k = \{\exp(\frac{2\pi im}{k}) : m = 0,1,\ldots,k-1\}$

建立半定规划问题：

maximize r
subject to:
C_j ≥ 0, ∀j = 1,...,k²
Σ C_j = I_N
Σ (N_i)_{jj} C_j = rU_i, i = 1,2

扩张常数为 $c(U,N) = r^{-1}$

技术创新点

近似策略：使用正则 $k$ 边形近似单位圆，误差控制在 $\cos(\pi/k)$ 因子内
半定规划方法：将扩张常数计算转化为凸优化问题
水平收敛理论：建立了矩阵范围收敛与扩张常数的关系

实验设置

数据集

独立Haar分布的 $N \times N$ 酉矩阵对
矩阵大小： $N = 10, 15, 20, \ldots, 60$ （主要实验）
扩展实验： $N$ 最大到300

评价指标

扩张常数 $c(U^{(N)}, u_0)$ 的数值近似
收敛性分析：均值和标准差随 $N$ 的变化

实现细节

编程语言：Julia
优化包：Convex.jl + SCS求解器
验证：部分结果用MOSEK求解器交叉验证
参数设置： $k = 8, 20, 30$ （正则多边形边数）

实验结果

主要结果

收敛证据：随着 $N$ 增长，计算得到的扩张常数值聚集在 $\sqrt{2} \approx 1.414$ 附近
统计特性：
- 对于 $N = 5j$ （ $j = 2,\ldots,12$ ）， $k=20$ ，运行100次实验
- 标准差随 $N$ 增长而减小，表明收敛性
数值界限：
- 所有计算值都低于 $\frac{\sqrt{2} + \sqrt{3}}{2} \approx 1.57$
- 对于 $k=8$ ， $N=125$ 的100次试验均聚集在1.439附近