Integrable matrix theory: Level statistics.

We study level statistics in ensembles of integrable N×N matrices linear in a real parameter x. The matrix H(x) is considered integrable if it has a prescribed number n>1 of linearly independent commuting partners H^{i}(x) (integrals of motion) [H(x),H^{i}(x)]=0, [H^{i}(x),H^{j}(x)]=0, for all x. In a recent work [Phys. Rev. E 93, 052114 (2016)2470-004510.1103/PhysRevE.93.052114], we developed a basis-independent construction of H(x) for any n from which we derived the probability density function, thereby determining how to choose a typical integrable matrix from the ensemble. Here, we find that typical integrable matrices have Poisson statistics in the N→∞ limit provided n scales at least as logN; otherwise, they exhibit level repulsion. Exceptions to the Poisson case occur at isolated coupling values x=x_{0} or when correlations are introduced between typically independent matrix parameters. However, level statistics cross over to Poisson at O(N^{-0.5}) deviations from these exceptions, indicating that non-Poissonian statistics characterize only subsets of measure zero in the parameter space. Furthermore, we present strong numerical evidence that ensembles of integrable matrices are stationary and ergodic with respect to nearest-neighbor level statistics.

A note on Lee-Yang zeros in the negative half-plane.

We obtain lower bounds on the inverse compressibility of systems whose Lee-Yang zeros of the grand-canonical partition function lie in the left half of the complex fugacity plane. This includes in particular systems whose zeros lie on the negative real axis such as the monomer-dimer system on a lattice. We also study the virial expansion of the pressure in powers of the density for such systems. We find no direct connection between the positivity of the virial coefficients and the negativity of the L-Y zeros, and provide examples of either one or both properties holding. An explicit calculation of the partition function of the monomer-dimer system on two rows shows that there are at most a finite number of negative virial coefficients in this case.

Rotationally invariant ensembles of integrable matrices.

We construct ensembles of random integrable matrices with any prescribed number of nontrivial integrals and formulate integrable matrix theory (IMT)-a counterpart of random matrix theory (RMT) for quantum integrable models. A type-M family of integrable matrices consists of exactly N-M independent commuting N×N matrices linear in a real parameter. We first develop a rotationally invariant parametrization of such matrices, previously only constructed in a preferred basis. For example, an arbitrary choice of a vector and two commuting Hermitian matrices defines a type-1 family and vice versa. Higher types similarly involve a random vector and two matrices. The basis-independent formulation allows us to derive the joint probability density for integrable matrices, similar to the construction of Gaussian ensembles in the RMT.