Mixed interior and boundary bubbling solutions for Neumann problem in R2

Let Ω be a bounded domain in R2 with smooth boundary, we consider the following problem: −Δu+u=λup−1eup, u>0, in Ω, with Neumann boundary condition on ∂Ω, where λ>0 is a small parameter, 0