Random pure states: Quantifying bipartite entanglement beyond the linear statistics

We analyze the properties of entangled random pure states of a quantum system partitioned into two smaller subsystems of dimensions N and M. Framing the problem in terms of random matrices with a fixed-trace constraint, we establish, for arbitrary N <= M, a general relation between the n-point densities and the cross moments of the eigenvalues of the reduced density matrix, i.e., the so-called Schmidt eigenvalues, and the analogous functionals of the eigenvalues of the Wishart-Laguerre ensemble of the random matrix theory. This allows us to derive explicit expressions for two-level densities, and also an exact expression for the variance of von Neumann entropy at finite N, M. Then, we focus on the moments E{K-a} of the Schmidt number K, the reciprocal of the purity. This is a random variable supported on [1, N], which quantifies the number of degrees of freedom effectively contributing to the entanglement. We derive a wealth of analytical results for E{Ka} for N = 2 and 3 and arbitrary M, and also for square N = M systems by spotting for the latter a connection with the probability P(x(min)(GUE) >= root 2N xi) that the smallest eigenvalue x(min)(GUE) of an N x N matrix belonging to the Gaussian unitary ensemble is larger than root 2N xi. As a by-product, we present an exact asymptotic expansion for P(x(min)(GUE) >= root 2N xi) for finite N as xi -> infinity. Our results are corroborated by numerical simulations whenever possible, with excellent agreement.