In this paper, we construct multi-bump solutions for the following coupled system of Schrödinger equations with \(\chi ^{(2)}\) nonlinearity \(\begin{aligned} {\left\{ \begin{array}{ll} -\Delta u_{1}+(1+\varepsilon P(x)) u_{1} = \alpha u_{1} u_{2}, & x \in \mathbb {R}^{N},\\ -\Delta u_{2}+(1+\varepsilon Q(x)) u_{2} = \frac{\alpha }{2}u_{1}^{2}+\beta u_{2}^{2},& x \in \mathbb {R}^{N}, \end{array}\right. } \end{aligned}\) where \(\alpha >0\) , \(\alpha >\beta \) , \(2\le N <6\) and \(\varepsilon >0\) is a small parameter, the positive potentials P(x), Q(x) are continuous and \(\lim \limits _{|x|\rightarrow \infty }P(x)=0, \lim \limits _{|x|\rightarrow \infty }Q(x)=0\) .