# How to prove consistency and asymptotic normality of the inverse of sample covariance matrix?

Let $X \sim N(0, \Sigma)$ be a $d$-dimensional Gaussian random vector. Suppose I have $n$ samples $X_1, \ldots, X_n$.

Consider the large sample regime where $d$ is fixed and $n$ goes to infinity, the MLE estimator of $\Sigma^{-1}$ is $\hat{\Sigma}^{-1} = (\frac{1}{n}\sum_{i=1}^nX_iX_i^\top)^{-1}$.

My question is how to establish the consistency result that $\hat{\Sigma}^{-1} \rightarrow \Sigma^{-1}$ and the asymptotic distribution of $\hat{\Sigma}^{-1} - \Sigma^{-1}$.

Thank you!

Hint: Look at each entry of $X_i X_i^T$, then, look at their average. Then, use the Law of large numbers. Finally, use the continuous mapping theorem.