We provide a general theoretical framework to derive Bernstein-von Mises theorems for matrix functionals. The conditions on functionals and priors are explicit and easy to check. Results are obtained for various functionals including entries of covariance matrix, entries of precision matrix, quadratic forms, log-determinant, eigenvalues in the Bayesian Gaussian covariance/precision matrix estimation setting, as well as for Bayesian linear and quadratic discriminant analysis. are asymptotically the same under the sampling distribution P n θ * . Note that the first one, known as the posterior, is of interest to Bayesians, and the second one is of interest to frequentists in the large sample theory. Applications of Bernstein-von Mises theorem include constructing confidence sets from Bayesian methods with frequentist coverage guarantees.Despite the success of BvM results in the classical parametric setting, little is known about the high-dimensional case, where the unknown parameter is of increasing or even infinite dimensions. The pioneering works of [11] and [13] (see also [17]) showed that generally BvM may not be true in non-classical cases. Despite the negative results, further works on some notions of nonparametric BvM provide some positive answers. See, for example, [22,8,9,24]. In this paper, we consider the question whether it is possible to have BvM results for matrix functionals, such as matrix entries and eigenvalues, when the dimension of the matrix p grows with the sample size n.This paper provides some positive answers to this question. To be specific, we consider a multivariate Gaussian likelihood and put a prior on the covariance matrix. We prove that the posterior distribution has a BvM behavior for various matrix functionals including entries of covariance matrix, entries of precision matrix, quadratic forms, log-determinant, and eigenvalues. All of these conclusions are obtained from a general theoretical framework we provide in Section 2, where we propose explicit easy-to-check conditions on both functionals and priors. We illustrate the theory by both conjugate and non-conjugate priors. A slight extension of the general framework leads to BvM results for discriminant analysis. Both linear discriminant analysis (LDA) and quadratic discriminant analysis (QDA) are considered. This work is inspired by a growing interest in studying the BvM phenomena on a lowdimensional functional of the whole parameter. That is, the asymptotic distribution of √ n(f (θ) −f )|X n , with f being a map from Θ to R d , where d does not grow with n. A special case is the semiparametric setting, where θ = (µ, η) contains both a parametric part µ and a nonparametric part η. The functional f takes the form of f (µ, η) = µ.The works in this field are pioneered by [19] in a right-censoring model and [26] for a general theory in the semiparametric setting. However, the conditions provided by [26] for BvM to hold are hard to check when specific examples are considered. To the best of our knowledge, the first general framework ...