Let's define the harmonic oscillator $H = −∆+x^2$ in a domain $\Omega$ of $\mathbb R^d$. Thus,
we consider the Dirichlet eigenvalue problem
$$ (H - \lambda ^2)u (x) = 0, \ x \in \Omega ;  \  \text{   and    } \ u(x) = 0 \ x \in \partial \Omega.$$  
Denote the spectral projection to the eigenvalue $ k = \lambda^2 $ by $P_k$.  

My question is: Are there any $L^p$ bounds results for eigenfunctions of the Hermite operator $H$ that have the form  
$$ \| P_k u \|_{L^p(\Omega)} \leq k^{\rho(p)} \| u\|_{L^2(\Omega)},$$ where $p$ would be in $[2, 2d/(d-2)]$ and the exponent $\rho$ as function of $1/p$. 

Note that in the case of $\Omega = \mathbb R^d$ there are many results, Karadzhov 1994, Thangavelu 1998, Koch and Tataru 2004.