We consider the problem
{epsilon 2 Delta u - u + f(u) = 0 in Omega
u > 0 in Omega, u = 0 on partial derivative Omega,
where Omega is a smooth domain in R-N, not necessarily bounded, epsilon > 0
is a small parameter and f is a superlinear, subcritical nonlinearity. It
is known that this equation possesses a solution that concentrates, as epsi
lon approaches zero, at a maximum of the function d(x) = d(.,partial deriva
tive Omega), the distance to the boundary. We obtain multi-peak solutions o
f the equation given above when the domain Omega presents a distance functi
on to its boundary d with multiple local maxima. We find solutions exhibiti
ng concentration at any prescribed finite set of local maxima, possibly deg
enerate, of d. The proof relies on variational arguments, where a penalizat
ion-type method is used together with sharp estimates of the critical value
s of the appropriate functional. Our main theorem extends earlier results,
including the single peak case. We allow a degenerate distance function and
a more general nonlinearity.