I was wondering whether anything is known on the following: Let
$h\_k (x)= (-1)^k e^{x^2/2} \frac {d^k}{dx^k} \, e^{-x^2/2}$, $k \geq 0$, be the classical
Hermite polynomials ($h\_0(x) = 1$, $h\_1(x) = x$,
$h\_2(x) = x^2 -1$, $h\_3(x) = x^3 - 3x$, ...).
Is there anything known on the asymptotic distributional behavior (if any), as $k \to \infty$,
of the sequence
of random variables $X\_k$, $k \geq 0$, such that the law of $X\_k$ has a density
proportional to $h\_k^2 \, e^{-x^2/2}$, $x \in \mathbb{R}$, with respect to the
Lebesgue measure?



Random variables with density distributions given by squared Hermite polynomials
This is related to the probability density of the [quantum harmonic oscillator.](https://en.wikipedia.org/wiki/Quantum_harmonic_oscillator) (The Hermite polynomials are the eigenstates of that problem.) For large $k$ the density
$$P\_k(x)=\frac{1}{\sqrt{2\pi}k!} [h\_k(x)]^2 \, e^{-x^2/2},$$
normalized to unity, has the approximation
$$p\_k(x)=\frac{1}{\pi}(4k-x^2)^{-1/2},\;\;|x|<2\sqrt{k}.$$
The boundaries $\pm 2\sqrt{k}$ are the classical turning points of the harmonic oscillator. In that context $E=4k$ is the energy of the particle and $V(x)=x^2$ is the harmonic potential in which it moves. The large-$k$ approximation is known as the [WKB approximation](https://en.wikipedia.org/wiki/WKB_approximation). More generally, for an arbitrary potential $V(x)$ the density varies as $[E-V(x)]^{-1/2}$.


Alternatively, one can use the asymptotic expansion of the [Hermite polynomials,](https://en.wikipedia.org/wiki/Hermite_polynomials)
$$e^{-x^2/4}h\_k(x) \propto \cos \left(x \sqrt{ k}- \pi k/2 \right)\left(4k-x^2\right)^{-1/4}.$$
Squaring and averaging over the rapidly oscillating $\cos^2$ gives $p\_k(x)$.




---


![](https://i.stack.imgur.com/g0TN6.png)
The plot compares $P\_k$ and $p\_k$ for $k=100$. The function $p\_k$ is a smoothed version of $P\_k$, suitable for averaging quantities that vary little on the scale of $1/\sqrt{k}$. For example, one can check by explicit computation that the average of $x^n$, $n\in\mathbb{N}$, agrees to leading order in $k$,
$$\lim\_{k\rightarrow\infty} k^{-n}\int P\_k(x)x^{2n}\,dx=\frac{2^{2n}\,\Gamma(n+1/2)}{\sqrt{\pi}\,\Gamma(n+1)}=k^{-n}\int p\_k(x) x^{2n}\,dx.$$


