Non-intersecting Squared Bessel Process: Spectral Moments and Dynamical Entanglement Entropy
Youyi Huang, Lu Wei
TL;DR
The paper develops a time-dependent random-matrix-inspired model for entanglement using non-intersecting squared Bessel processes to extend the Hilbert-Schmidt ensemble. It derives recurrence relations for real-order spectral moments $m_k$, enabling systematic computation of entanglement cumulants and linear statistics, including $T_k=\sum_i x_i^k\ln x_i$. Exact finite-$n$ formulae are obtained for the average quantum purity $\mathbb{E}[P]$ and the von Neumann entropy $\mathbb{E}[S]$, with $\mathbb{E}[P]$ given by a closed rational form in $\alpha$, $c$, and $n$, and $\mathbb{E}[S]$ expressed via $\mathbb{E}[T]$ and $m_1$; in the limit $a\to0$ the results recover the Hilbert-Schmidt/Page statistics. The work also develops new biorthogonal polynomial structures, including four-point and structure relations for type I and II multiple orthogonal polynomials, and a Christoffel-Darboux formula for the correlation kernel, enriching the mathematical toolbox for non-intersecting Bessel processes and their spectral properties.
Abstract
Statistical ensembles of reduced density matrices of bipartite quantum systems play a central role in entanglement estimation, but do not capture the non-stationary nature of entanglement relevant to realistic quantum information processing. To address this limitation, we propose a dynamical extension of the Hilbert-Schmidt ensemble, a baseline statistical model for entanglement estimation, arising from non-intersecting squared Bessel processes and perform entanglement estimation via average entanglement entropy and quantum purity. The investigation is enabled by finding spectral moments of the proposed dynamical ensemble, which serves as a new approach for systematic computation of entanglement metrics. Along the way, we also obtain new results for the underlying multiple orthogonal polynomials of modified Bessel weights, including structure and recurrence relations, and a Christoffel-Darboux formula for the correlation kernels.
