Real Eigenvalues of Asymmetric Wishart Matrices: Expected Number, Global Density and Integrable Structure

We investigate the real eigenvalues of asymmetric Wishart matrices of size $N$, indexed by the rectangular parameter $\nu \in \mathbb{N}$ and the non-Hermiticity parameter $\tau \in [0,1]$. The rectangular parameter $\nu$ is either fixed or proportional to $N$. The non-Hermiticity parameter $\tau$ is either fixed or $\tau = 1 - O(1/N)$, corresponding to the strongly and weakly non-Hermitian regimes, respectively. We establish a decomposition structure for the finite-$N$ correlation kernel of the real eigenvalues, which form Pfaffian point processes. Taking the symmetric limit $\tau = 1$, where the model reduces to the Laguerre orthogonal ensemble, this decomposition structure reduces to the known rank-one perturbation structure established by Adler, Forrester, Nagao, and van Moerbeke, as well as by Widom. Using the decomposition structure, we show that the expected number of real eigenvalues is proportional to $\sqrt{N}$ in the strongly non-Hermitian regime and to $N$ in the weakly non-Hermitian regime, providing explicit leading coefficients in both cases. Furthermore, we derive the limiting real eigenvalue densities, which recovers the Marchenko-Pastur distribution in the symmetric limit.