I'm trying to find the asymptotic expansion for $$ F(\beta)=\int_{-a}^{a} dx \, \frac{1}{e^{\beta g(x)}+1} \;\;\;\;\;\;\mbox{as}\;\;\;\;\;\; \beta \rightarrow +\infty.$$ where $a>0$ and $g(x)$ is a smooth, bounded and even function.
For example $g(x)=cos(x)$
The expression
$$\frac1{e^{\beta g(x)}+1}$$ corresponds to a sigmoid in the vicinity of the simple roots of $g(x)$, making a transition from $1$ to $0$ or conversely. For large $\beta$, the transition is very sharp and tends to a Heaviside step.
Around a root $r$, using the Taylor development to the first order, $g(x-r)\approx g'_r\,(x-r)$,
$$\int_{-a}^a\frac{dx}{e^{\beta g'_r(x-r)}+1}=2a-\frac1{\beta g'_r}\log\frac{e^{\beta g'_r(a-r)}+1}{e^{\beta g'_r(-a-r)}+1}=a-\frac1{\beta g'_r}\log\frac{e^{\beta g'_ra}e^{-\beta g'_rr}+1}{e^{-\beta g'_rr}+e^{\beta g'_ra}}.$$
When $r=0$, the second term is zero, because the sigmoid is symmetric, and the term $a$ just corresponds to the constant value $1$ in the negatives.
In a more general case of several simple roots, you will get the total length of the intervals where $g$ is negative (a constant), plus a small correction term due to the asymmetric location of the root inside $[-a,a]$.