Solution (source code)

= Solution

The phase $f(x)=x(1-\cos x)$ is nonnegative on the interval and vanishes only at its two endpoints. An interval away from both endpoints contributes exponentially little, so both endpoint contributions must be expanded before their orders are ranked.

At the left endpoint,
$$
f(x)=\frac{x^3}{2}-\frac{x^5}{24}+O(x^7).
$$
The <Laplace's method> scale is $x=\sigma^{-1/3}z$. Hence its leading contribution is
$$
\sigma^{-1/3}\int_0^\infty e^{-z^3/2}\,dz
=\frac{2^{1/3}}3\Gamma\left(\frac13\right)\sigma^{-1/3}.
$$
The substitution $v=z^3/2$ evaluates the <integral> by the <Gamma function>. Since there is no quartic term in the phase, the next left-end correction is relatively $O(\sigma^{-2/3})$, and thus absolutely $O(\sigma^{-1})$.

At the right endpoint put $r=2\pi-x$. There
$$
f(2\pi-r)=\pi r^2-\frac{r^3}{2}+O(r^4),
$$
so $r=\sigma^{-1/2}z$ gives
$$
\sigma^{-1/2}\int_0^\infty e^{-\pi z^2}\,dz
=\frac12\sigma^{-1/2}.
$$
The leading right-end integral is a <Gaussian integral>. Its first correction is also absolutely $O(\sigma^{-1})$. The <mixed cubic and quadratic Laplace endpoints> therefore produce the first two global terms
$$
\boxed{\mathcal I(\sigma)=
\frac{2^{1/3}\Gamma(1/3)}{3\sigma^{1/3}}
+\frac1{2\sigma^{1/2}}+O(\sigma^{-1})}.
$$
Keeping only the dominant cubic endpoint would miss the second term, which is larger than that endpoint's own first correction.