Abstract In this article, we are interested in multi-bump solutions of the singularly perturbed problem - ε 2 Δ v + V ( x ) v = f ( v ) in ℝ N . -\varepsilon^{2}\Delta v+V(x)v=f(v)\quad\text{in }\mathbb{R}^{N}. Extending previous results, we prove the existence of multi-bump solutions for an optimal class of nonlinearities f satisfying the Berestycki–Lions conditions and, notably, also for more general classes of potential wells than those previously studied. We devise two novel topological arguments to deal with general classes of potential wells. Our results prove the existence of multi-bump solutions in which the centers of bumps converge toward potential wells as ε → 0 {\varepsilon\rightarrow 0} . Examples of potential wells include the following: the union of two compact smooth submanifolds of ℝ N {\mathbb{R}^{N}} where these two submanifolds meet at the origin and an embedded topological submanifold of ℝ N {\mathbb{R}^{N}} .