In this paper, we are concerned with singularly perturbed Schrödinger equation \(\begin{aligned} -\epsilon ^2\Delta u+V(x)u=f(u),\quad x\in \mathbb {R}^N, \end{aligned}\) where \(\epsilon >0\) is a small parameter, V possesses global minimum points, and f is asymptotically linear at infinity. Using variational methods and construction methods, we reveal the relationship between the number of solutions and the profile of the potential V. In particular, some new tricks and the method of Nehari manifold dependent on a suitable restricted set are introduced to overcome the difficulty resulting from the appearance of asymptotically linear nonlinearity.
