Locating the peaks of the least energy solutions to an ellyptic system with Neumann boundary conditions

We consider a system of the form -(epsilon)^2 (delta)u+u=g(v), -(epsilon)^2(delta)v+v=f(u) in (omega) with Neumann boundary condition on (delta)(omega), where (omega) is a smooth bounded domain in RN, N>=3 and f,g are power-type nonlinearities having superlinear and subcritical growth at infinity. We prove that the least energy solutions to such a system concentrate, as goes to zero, at a point of the boundary which maximizes the mean curvature of the boundary of (omega).