If we have a set $X_1,\dots,X_n$ of iid random variables with finite mean $\mu$ and variance $\sigma$, the CLT says that $\sqrt{n}(\overline X_n - \mu) \stackrel{d}{\to} \mathcal{N}(0,\sigma^2)$. If we also have a finite third moment, then the Berry-Essen theorem says that $$ |F_n(x) - \Phi(x)| \le \frac{C E[|X_1|^3]}{\sigma^3 \sqrt{n}}. $$ where $F_n$ is the CDF of $\sqrt{n} \overline X_n /\sigma$.
Do either the CLT or the Berry-Essen theorem allow us to say something about the convergence of the pdf of $\overline X_n$ to the pdf of the normal distribution?
My motivation is that I would like to approximate an expectation of a function of the sample mean for large $n$. So I write $$ \begin{align} E[g(\overline X_n)] & = \int g(x) p_{\overline X_n}(x) dx \\ & = \int g(x) (\varphi(x) - \varphi(x) + p_{\overline X_n}(x)) dx \\ & = \int g(x) \varphi(x)dx + R, \end{align} $$ where the remainder $R$ is $$ R = \int g(x) (p_{\overline X_n}(x) - \varphi(x)) dx. $$ where $\varphi(x)$ is the pdf of the (properly scaled) normal distribution. So I would like to bound the remainder as $n \to \infty$ to show this is a valid approximation. If we were talking about the cdf instead of the pdf, the Berry-Essen inequality could be used to bound the remainder. So is it possible to bound the remainder given that it is in terms of the pdf?
The function $g$ that I am most interested is of the form $g(\overline X_n) = |\overline X_n| \mathbb{1}_{\{|\overline X_n| > c\}}$.