If we define
$$ I(a,b) = \int_{-\infty}^{+\infty}\exp\left[-b(x^2-a^2)^2\right]\,dx $$
for $a,b>0$, by setting $c=ba^4$ and $x=az$ we get:
$$ I(a,b) = a \int_{-\infty}^{+\infty}\exp\left[-c(z^2-1)^2\right]\,dz = a\int_{0}^{+\infty}\exp\left[-c(z-1)^2\right]\,\frac{dz}{\sqrt{z}}\stackrel{\text{def}}{=} a\,J(c)$$
and:
$$ J(c) = \int_{-1}^{+\infty}\frac{\exp(-c z^2)}{\sqrt{z+1}}\,dz =\color{blue}{\int_{-1}^{0}\frac{\exp(-cz^2)}{\sqrt{z+1}}\,dz}+\color{red}{\int_{0}^{+\infty}\frac{\exp(-cz^2)}{\sqrt{z+1}}\,dz}$$
where the blue integral can be approximated by expanding the integrand function as a Taylor series and the red integral can be studied by switching to Laplace transforms and getting values of Bessel functions. In any case, the behaviour depends on the magnitude of $\color{green}{ba^4}$.
In terms of modified Bessel functions of the first kind,
$$ I(a,b) = \frac{\pi a}{2 \exp(ba^4/2)}\left[I_{-1/4}(ba^4/2)+I_{1/4}(ba^4/2)\right].$$
It follows that if $ba^4$ is large we have
$$ I(a,b) \approx \frac{\pi a }{\sqrt{\pi b a^4}}=\sqrt{\frac{\pi}{b a^2}}$$
while if $ba^4$ is close to zero we have
$$ I(a,b) \approx \frac{\pi}{2^{3/4}\Gamma(3/4)b^{1/4}}.$$