Multi-peak solutions to coupled Schrödinger systems with Neumann boundary conditions

In this paper, we are concerned with the following two coupled Schrödinger systems in a bounded domain Ω⊂RN(N=2,3)Ω⊂RN(N=2,3) with Neumann boundary conditions {−ε2Δu+u=μ1u3+βuv2,−ε2Δv+v=μ2v3+βu2v,u>0,v>0,∂u/∂n=0,∂v/∂n=0,on  ∂Ω. Suppose the mean curvature H(P)H(P) of the boundary ∂Ω∂Ω has several local minimums or local maximums, we obtain the existence of solutions with multi-peaks to the system with all peaks being on the boundary and all peaks locate either near the local maxima or near the local minima of the mean curvature at the boundary of the domain.