Let $a>0$ be a fixed number and consider the Hermite operator (or harmonic oscillator) defined by \begin{equation} Hf(x)=x^2f(x)-f''(x) \end{equation} for any smooth function $f$ compactly supported on the interval $(-a,a)$. Then, we may extend $H$ as an unbounded operator on $L^2(-a,a)$.
Extending from the Wikipedia link : https://en.wikipedia.org/wiki/Dirichlet_eigenvalue I would like to compute the first Dirichlet eigenvalue $\lambda_a$ of the above harmonic operator $H$ on the interval $(-a,a)$.
When working on whole $\mathbb{R}$, it is well-known that the first eigenvalue of $H$ is $1$. What I suspect from connection with the standard Gaussian measure is the following formula: \begin{equation} \lambda_a=\frac{\sqrt{2\pi}}{\int_{-a}^a e^{-\frac{x^2}{2}}dx} \end{equation}
However, I cannot find a way to justify my guess. I tried to look for a relevant reference but all seem to deal with the Laplacian only, let alone the explicit formula for the first Dirichlet eigenvalue..
Could anyone please help me?