Multi-bump solutions for a Kirchhoff problem type

In this paper, we are going to study the existence of solution for the following Kirchhoff problem $$\{\begin{array}{l}M({\displaystyle {\int }_{{\mathbb{R}}^3}\mathrm{\mid }\mathrm{\nabla }u{\mathrm{\mid }}^{2}dx+{\displaystyle {\int }_{{\mathbb{R}}^3}\lambda a(x)+1)u^{2}dx)(-\mathrm{\Delta }u+(\lambda a(x)+1)u)=f(u)in{\mathbb{R}}^3,}}\\ \text{mbox}\\ u\in H^1({\mathbb{R}}^3)\mathrm{.}\end{array}$$ Assuming that the nonnegative function $a(x)$ has a potential well with $int(a^{-1}(\{0\}))$ consisting of $k$ disjoint components ${\mathrm{\Omega }}_1,{\mathrm{\Omega }}_2,\mathrm{.}\mathrm{.}\mathrm{.}\mathrm{.}\mathrm{.},{\mathrm{\Omega }}_k$ and the nonlinearity $f(t)$ has a subcritical growth, we are able to establish the existence of positive multi-bump solutions by variational methods.