Let $X=(X_1, X_2, \dots, X_n)$ be a smooth vector field on $\mathbb{R}^n$. The operator $L=(\sum_{i=1}^{m}X_i^2)^p$, where $p$ is an integer, is a degenerated operator. If $X$ satisfies the Hörmander condition, for the case of $p=1$, we have the subelliptic estimates
$$|u|_{s+a}\leq |Lu|_{s} + |u|_{L^2}$$
for $u \in L^2(\Omega)\cap C^{\infty}_0$; here $|\cdot|_{s}$ denotes the Sobolev norm.
Now I want to know what if $p\geq 2$? Moreover, what about general case? i.e. Rockland operator (i.e. higher order operator). I want to get the estimate like the case $p=1$.