In the present paper, we consider the nonlocal Kirchhoff problem
$$-\left(\epsilon^2a+\epsilon b\int_{{\open R}^{3}}\vert \nabla u \vert^{2}\right)\Delta u+V(x)u=u^{p}, \quad u \gt 0 \quad {\rm in} {\open R}^{3},$$
where a, b>0, 1<p<5 are constants, ϵ>0 is a parameter. Under some mild assumptions on the function V, we obtain multi-peak solutions for ϵ sufficiently small by Lyapunov–Schmidt reduction method. Even though many results on single peak solutions to singularly perturbed Kirchhoff problems have been derived in the literature by various methods, there exist no results on multi-peak solutions before this paper, due to some difficulties caused by the nonlocal term $\left(\int_{{\open R}^{3}} \vert \nabla u \vert^{2}\right)\Delta u$. A remarkable new feature of this problem is that the corresponding unperturbed problem turns out to be a system of partial differential equations, but not a single Kirchhoff equation, which is quite different from most of the elliptic singular perturbation problems.