Bounds for the gamma function

We improve the upper bound of the following inequalities for the gamma function $\Gamma$ due to H. Alzer and the author. \begin{equation*} \exp\left(-\frac{1}{2}\psi(x+1/3)\right)<\frac{\Gamma(x)}{x^xe^{-x}\sqrt{2\pi}}<\exp\left(-\frac{1}{2}\psi(x)\right). \end{equation*} We also prove the following new inequalities: For $x\geq1$ $$\sqrt{2\pi}x^xe^{-x}\left(x^2+\frac{x}{3}+a_*\right)^{\frac{1}{4}}<\Gamma(x+1)<\sqrt{2\pi}x^xe^{-x}\left(x^2+\frac{x}{3}+a^*\right)^{\frac{1}{4}}$$ with the best possible constants $a_*=\frac{e^4}{4\pi^2}-\frac{4}{3}=0.049653963176...$, and $a^*=1/18=0.055555...$, and for $x\geq0$ \begin{equation*} \exp\left[x\psi\left(\frac{x}{\log (x+1)}\right)\right]\leq\Gamma(x+1)\leq\exp\left[x\psi\left(\frac{x}{2}+1\right)\right], \end{equation*} where $\psi$ is the digamma function.