Let $X$ be a centered Gaussian random variable in a separable Hilbert space ${\mathbb H}$ with covariance operator $Σ.$ We study a problem of estimation of a smooth functional of $Σ$ based on a sample $X_1,\dots ,X_n$ of $n$ independent observations of $X.$ More specifically, we are interested in functionals of the form $\langle f(Σ), B\rangle,$ where $f:{\mathbb R}\mapsto {\mathbb R}$ is a smooth function and $B$ is a nuclear operator in ${\mathbb H}.$ We prove concentration and normal approximation bounds for plug-in estimator $\langle f(\hat Σ),B\rangle,$ $\hat Σ:=n^{-1}\sum_{j=1}^n X_j\otimes X_j$ being the sample covariance based on $X_1,\dots, X_n.$ These bounds show that $\langle f(\hat Σ),B\rangle$ is an asymptotically normal estimator of its expectation ${\mathbb E}_Σ \langle f(\hat Σ),B\rangle$ (rather than of parameter of interest $\langle f(Σ),B\rangle$) with a parametric convergence rate $O(n^{-1/2})$ provided that the effective rank ${\bf r}(Σ):= \frac{{\bf tr}(Σ)}{\|Σ\|}$ (${\rm tr}(Σ)$ being the trace and $\|Σ\|$ being the operator norm of $Σ$) satisfies the assumption ${\bf r}(Σ)=o(n).$ At the same time, we show that the bias of this estimator is typically as large as $\frac{{\bf r}(Σ)}{n}$ (which is larger than $n^{-1/2}$ if ${\bf r}(Σ)\geq n^{1/2}$). In the case when ${\mathbb H}$ is finite-dimensional space of dimension $d=o(n),$ we develop a method of bias reduction and construct an estimator $\langle h(\hat Σ),B\rangle$ of $\langle f(Σ),B\rangle$ that is asymptotically normal with convergence rate $O(n^{-1/2}).$ Moreover, we study asymptotic properties of the risk of this estimator and prove minimax lower bounds for arbitrary estimators showing the asymptotic efficiency of $\langle h(\hat Σ),B\rangle$ in a semi-parametric sense.