Bivariate moments of the two-point correlation function for embedded Gaussian unitary ensemble with k-body interactions.

Embedded random matrix ensembles with k-body interactions are well established to be appropriate for many quantum systems. For these ensembles the two point correlation function is not yet derived, though these ensembles are introduced 50 years back. Two-point correlation function in eigenvalues of a random matrix ensemble is the ensemble average of the product of the density of eigenvalues at two eigenvalues, say E and E^{'}. Fluctuation measures such as the number variance and Dyson-Mehta Δ_{3} statistic are defined by the two-point function and so also the variance of the level motion in the ensemble. Recently, it is recognized that for the embedded ensembles with k-body interactions the one-point function (ensemble averaged density of eigenvalues) follows the so called q-normal distribution. With this, the eigenvalue density can be expanded by starting with the q-normal form and using the associated q-Hermite polynomials He_{ζ}(x|q). Covariances S_{ζ}S_{ζ^{'}}[over ¯] (overline representing ensemble average) of the expansion coefficients S_{ζ} with ζ≥1 here determine the two-point function as they are a linear combination of the bivariate moments Σ_{PQ} of the two-point function. Besides describing all these, in this paper formulas are derived for the bivariate moments Σ_{PQ} with P+Q≤8, of the two-point correlation function, for the embedded Gaussian unitary ensembles with k-body interactions [EGUE(k)] as appropriate for systems with m fermions in N single particle states. Used for obtaining the formulas is the SU(N) Wigner-Racah algebra. These formulas with finite N corrections are used to derive formulas for the covariances S_{ζ}S_{ζ^{'}}[over ¯] in the asymptotic limit. These show that the present work extends to all k values, the results known in the past in the two extreme limits with k/m→0 (same as q→1) and k=m (same as q=0).