On the distribution of perturbations of propagated Schrödinger eigenfunctions

Let $(M,g_0)$ be a compact Riemmanian manifold of dimension $n$. Let $P_0 (\h) := -\h^2\Delta_{g}+V$ be the semiclassical Schr\"{o}dinger operator for $\h \in (0,\h_0]$, and let $E$ be a regular value of its principal symbol $p_0(x,\xi)=|\xi|^2_{g_0(x)} +V(x)$. Write $\varphi_\h$ for an $L^2$-normalized eigenfunction of $P(\h)$, $P_0(\h)\varphi_\h =E(\h)\varphi_\h$ and $E(\h) \in [E-o(1),E+ o(1)]$. Consider a smooth family of perturbations $g_u$ of $g_0$ with $u$ in the ball $\mathcal B^k(\varepsilon) \subset \mathbb R^k$ of radius $\varepsilon>0$. For $P_{u}(\h) := -\h^2 \Delta_{g_u} +V$ and small $|t|$, we define the propagated perturbed eigenfunctions $$\varphi_\h^{(u)}:=e^{-\frac{i}{\h}t P_u(\h)} \varphi_\h.$$ We study the distribution of the real part of the perturbed eigenfunctions regarded as random variables $$\Re (\varphi^{(\cdot)}_\h(x)):\mathcal B^{k}(\varepsilon) \to \mathbb R \quad \quad \text{for}\;\, x\in M.$$ In particular, when $(M,g)$ is ergodic, we compute the $h \to 0^+$ asymptotics of the variance $\text{Var} [\Re (\varphi^{(\cdot)}_\h(x))] $ and show that all odd moments vanish as $h \to 0^+.$