We study concentration phenomena for the system

\[\varepsilon^2 \Delta v - v - \delta \phi v + \gamma v^{p} =0,\quad \Delta\phi+ \delta v^2=0 \]

in the unit ball $B_1$ of $\mathbb{R}^3$ with Dirichlet boundary conditions. Here $\varepsilon,\ \delta,\ \gamma >0$ and $p>1$. We prove the existence of positive radial solutions $(v_{\varepsilon}, \phi_{\varepsilon})$ such that $v_{\varepsilon}$ concentrates at a distance $({\varepsilon}/{2}) |{\rm log}\, {\varepsilon} |$ away from the boundary $\partial B_1$ as the parameter $\varepsilon$ tends to 0. The approach is based on a combination of Lyapunov–Schmidt reduction procedure together with a variational method.