Abstract: Suppose (M, g) is a noncompact Riemannian manifold with nonnegative Ricci curvature, and let hd(M) be the dimension of the space of harmonic functions with polynomial growth of growth order at most d. In this talk, we will review the work in estimating hd(M), and introduce my result on the asymptotic behavior of hd(M) when d is large.