Asymptotic infinitesimal freeness of covariance matrices

We consider $n\times n$ covariance matrices $M=\frac{1}{n}X{X}^{\ast }$ where $X=(x_{i,j})$ is a matrix whose entries are independent complex random variables with $\mathbb{E}(x_{i,j})=0$ and $\mathbb{E}(\mathrm{\mid }{x}_{i,j}{\mathrm{\mid }}^2)=1$. We derive a $\frac{1}{n}$ expansion of the mixed moments, $\frac{1}{n}\mathbb{E}(\text{Tr}(M^{(r_1)}\cdots M^{(r_q)}))$, of the form $a_0+a_1\frac{1}{n}+O(\frac{1}{n^2})$. This permits us to find explicit formulas for the moments and infinitesimal moments of several covariance matrices where we allow repetition. As an application of our formulas, we derive asymptotic freeness and infinitesimal freeness of independent covariance matrices under a fourth-moment condition. This generalizes previous results for the Wishart ensemble in which $x_{i,j}$ is complex Gaussian.