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

Basic Information

• Paper ID: 2510.12540
• Title: Empirical bounds for commuting dilations of free unitaries and the universal commuting dilation constant
• Authors: Malte Gerhold (University of Greifswald), Marcel Scherer (Technion), Orr Moshe Shalit (Technion)
• Classification: math.FA cs.NA math.NA math.OA
• Submission Date: October 14, 2025 to arXiv
• Paper Link: https://arxiv.org/abs/2510.12540

Abstract

For a tuple of Hilbert space operators $T$, the "commuting dilation constant" is the smallest number $c$ such that the operators in $T$ are simultaneously compressions of commuting normal operators with norm at most $c$. This paper provides strong numerical evidence that the commuting dilation constant of independent random $N \times N$ unitary matrix pairs converges almost surely to $\sqrt{2}$ (as $N \to \infty$). Under this assumption, it is proven that the commuting dilation constant of any pair of contraction operators is strictly less than 2. The experiments are based on a simple algorithm introduced by the authors for computing dilation constants between matrix tuples.

Research Background and Motivation

1. Problem to be Solved: Determine the exact value of the universal commuting dilation constant $C_2$, which is the smallest constant $c$ such that every pair of contraction operators is simultaneously a compression of commuting normal operators with norm at most $c$.
2. Importance of the Problem: Dilation theory techniques and matrix convex set frameworks have important applications in quantum information, optimal control, mathematical physics, operator systems, and operator algebra theory.
3. Limitations of Existing Methods:
   • Known bounds: $\sqrt{d} \leq C_d \leq \sqrt{2d}$ (for $d$-tuples)
   • For the case $d=2$, the current best upper bound $C_2 \leq \sqrt{2} \cdot 2 = 2\sqrt{2}$ coincides with the trivial bound $C_d \leq d$
   • The lower bound $C_d \geq \sqrt{d}$ is not tight in general
4. Research Motivation: By combining rigorous limit theorems with numerical experiments, provide evidence for $C_2 < 2$ and improve existing upper bounds.

Core Contributions

1. Numerical Evidence: Provides strong numerical evidence that the dilation constant of independent Haar unitary matrix pairs converges to $\sqrt{2}$
2. Theoretical Result: Proves that if this convergence holds, then $C_2 \leq \frac{2\sqrt{2}}{3} < 2$
3. Algorithmic Contribution: Introduces a simple algorithm for computing dilation constants of matrix tuples
4. Rigorous Bounds: Provides strict bounds for finite dimensions in the appendix: $C_2(n) \leq \sqrt{2} + 2\sin(\frac{\pi}{2}(1-\frac{1}{2n})) < 2$

Methodology Details

Task Definition

Given an operator $d$-tuple $u, v$, the dilation constant $c(u,v)$ is defined as: $c(u,v) = \inf\{c : u \prec cv\}$ where $u \prec cv$ means there exist a Hilbert space $H \subseteq K$ and operator tuples $U \in B(H)^d$, $V \in B(K)^d$ such that $u \sim U$, $v \sim V$ and $U = P_H cV|_H$.

Core Theoretical Framework

Theorem 2.1 (Horizontal Convergence and Dilation Constants): If sequences of operator $d$-tuples $(\xi^{(N)})$ and $(\eta^{(N)})$ have matrix ranges that horizontally converge to $\xi^{(\infty)}$ and $\eta^{(\infty)}$, then: $c(\xi^{(\infty)}, \eta^{(\infty)}) \leq \liminf_{N,M \to \infty} c(\xi^{(N)}, \eta^{(M)})$

Corollary 2.4: For collections of random matrices $T^{(N)}$ with strong asymptotic freeness, almost surely: $c(u^f, u_0) \leq \liminf c(T^{(N)}, u_0)$

Algorithm Design

Dilation Constant Computation Algorithm:

1. Construct commuting normal matrices $N = (N_1, N_2)$ whose diagonal values range over all $k^2$ possible pairs in $V_k \times V_k$, where $V_k = \{\exp(\frac{2\pi im}{k}) : m = 0,1,\ldots,k-1\}$
2. Establish a semidefinite programming problem:

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

3. The dilation constant is $c(U,N) = r^{-1}$

Technical Innovations

1. Approximation Strategy: Uses regular $k$-gons to approximate the unit circle with error controlled by the $\cos(\pi/k)$ factor
2. Semidefinite Programming Method: Transforms dilation constant computation into a convex optimization problem
3. Horizontal Convergence Theory: Establishes the relationship between matrix range convergence and dilation constants

Experimental Setup

Dataset

• Independent Haar-distributed $N \times N$ unitary matrix pairs
• Matrix sizes: $N = 10, 15, 20, \ldots, 60$ (main experiments)
• Extended experiments: $N$ up to 300

Evaluation Metrics

• Numerical approximations of the dilation constant $c(U^{(N)}, u_0)$
• Convergence analysis: mean and standard deviation as functions of $N$

Implementation Details

• Programming Language: Julia
• Optimization Package: Convex.jl + SCS solver
• Verification: Some results cross-verified with MOSEK solver
• Parameter Settings: $k = 8, 20, 30$ (number of sides of regular polygons)

Experimental Results

Main Results

1. Convergence Evidence: As $N$ increases, computed dilation constant values cluster around $\sqrt{2} \approx 1.414$
2. Statistical Properties:
   • For $N = 5j$ ($j = 2,\ldots,12$), $k=20$, 100 experimental runs
   • Standard deviation decreases with increasing $N$, indicating convergence
3. Numerical Bounds:
   • All computed values fall below $\frac{\sqrt{2} + \sqrt{3}}{2} \approx 1.57$
   • For $k=8$, $N=125$, all 100 trials cluster around 1.439

Detailed Data Analysis

Figure 1 Analysis: Histograms for different matrix sizes $N$ show that as $N$ increases, dilation constant values become increasingly concentrated around $\sqrt{2}$ (left dashed line) and remain below $\sqrt{2}/\cos(\pi/k)$ (right dashed line).

Figures 4-5 Analysis: For $k=8$, single-sample experiments show the dilation constant stabilizes around 1.439 for $N \geq 40$, with ten random sequences all converging to approximately 1.44.

Theoretical Derivation

Based on experimental results and theoretical relationships: $C_2 = c(u^u, u_0) \leq c(u^u, u^f)c(u^f, u_0) = \frac{2}{\sqrt{3}} \cdot \sqrt{2} = \frac{2\sqrt{2}}{3} < 2$

where the known result $c(u^u, u^f) = \frac{2}{\sqrt{3}}$ is used.

Related Work

1. Historical Development:
   • Ando's theorem: Power dilations of commuting contraction pairs
   • Applications of dilation theory in finite dimensions
2. Evolution of Bounds:
   • $C_d \leq d$ → $C_d \leq \max\{d, 2\sqrt{d}\}$ → $C_d \leq \sqrt{2d}$
   • Lower bound: $C_d \geq \sqrt{d}$
3. Free Probability Theory:
   • Strong asymptotic freeness of Haar unitary matrices
   • Horizontal convergence theory of matrix ranges

Conclusions and Discussion

Main Conclusions

1. Numerical Evidence: Strongly supports $\lim_{N \to \infty} c(U^{(N)}, u_0) = \sqrt{2}$
2. Theoretical Significance: If the conjecture holds, then $C_2 < 2$, improving existing bounds
3. Methodological Contribution: Provides a practical algorithm for computing dilation constants

Limitations

1. Numerical Nature: Results are based on numerical computation, not rigorous proof
2. Computational Constraints: Problem complexity is $O(k^2N^2)$, limiting the range of testable parameters
3. Approximation Error: Systematic error introduced by regular polygon approximation

Future Directions

1. Rigorous Proof: Seek a strict proof of $c(u^f, u_0) = \sqrt{2}$
2. Higher-Dimensional Generalization: Study dilation constants for $d > 2$
3. Algorithm Optimization: Develop more efficient methods for computing dilation constants

In-Depth Evaluation

Strengths

1. Innovative Methodology: Cleverly combines rigorous limit theory with numerical experiments
2. Experimental Sufficiency: Extensive experiments under multiple parameter settings with consistent results
3. Theoretical Rigor: Horizontal convergence theorem provides strict theoretical foundation
4. Practical Value: Algorithm applicable to other dilation constant computation problems

Weaknesses

1. Numerical Dependence: Core conclusions rely on numerical evidence, lacking completely rigorous proof
2. Computational Complexity: Algorithm complexity limits large-scale experiments
3. Error Analysis: Analysis of numerical errors could be more detailed

Impact

1. Theoretical Contribution: Provides significant progress on a long-standing open problem
2. Methodological Value: Demonstrates the application of numerical methods to pure mathematics problems
3. Application Prospects: Potential applications in quantum information and optimization theory

Applicable Scenarios

1. Operator Theory Research: Dilation theory and matrix convex set research
2. Quantum Information: Quantum measurement compatibility problems
3. Optimization Theory: Semidefinite programming and matrix inequality problems

References

The paper cites 63 related references, primarily including:

• Ando (1963): Foundational theory of commuting contraction pairs
• Davidson et al. (2017): Matrix convex sets and dilation theory
• Gerhold & Shalit (2021): Matrix range theory for random matrices
• Helton et al. (2013): Semidefinite programming algorithm foundations

Overall Assessment: This is a high-quality mathematical research paper that provides breakthrough progress on an important open problem through an innovative numerical-theoretical combined approach. Although the core conclusions rely on numerical evidence, its rigorous theoretical framework and sufficient experimental verification make the results highly credible.