A Gevrey function of order $\rho \geq 1$ is a smooth function $f:\mathbb{R}^n\rightarrow \mathbb{R}$ that satisfies the estimate $$|\partial^\alpha f|\leq AL^{|\alpha|}\alpha!^\rho $$ for some positive constants $A,L$, where $\alpha$ is an arbitrary multi-index.
For $\rho=1$, such functions are easily seen to be analytic and so $\mathcal{G}^1$ does not admit bump functions. On the other hand it is well known that for $\rho>1$ that $\mathcal{G}^\rho$ does admit bump functions.
For a project of mine, I need to construct Gevrey bump functions in arbitrary such classes, and to this end I would like to know:
Q/ How can one find the optimal $\rho$ for the function $\exp(-x^{-\alpha})$?
For $\alpha=1$, we can compute that it suffices to estimate $\sup_{t>0}p_k(t)e^{-t}$, where $p_{k+1}(t)=t^2(p_k(t)-p_k'(t))$ and $p_0=1$. The crude estimate of replacing $p_k$ by its leading term $t^{2k}$ does seem to yield the optimal $\rho=2$ but I have been unable to make this argument rigorous.
Alternatively, I would also welcome either a reference to such a computation or a method of bump function construction where the Gevrey regularity is more transparent.