Boundary concentration phenomena for an anisotropic Neumann problem in $\mathbb{R}^2$

Given a smooth bounded domain $\Omega$ in $\mathbb{R}^2$, we study the following anisotropic Neumann problem $$ \begin{cases} -\nabla(a(x)\nabla u)+a(x)u=\lambda a(x) u^{p-1}e^{u^p},\,\,\,\, u>0\,\,\,\,\, \textrm{in}\,\,\,\,\, \Omega,\\[2mm] \frac{\partial u}{\partial\nu}=0\,\, \qquad\quad\qquad\qquad\qquad\qquad\qquad \ \ \ \ \,\qquad\quad\, \textrm{on}\,\,\, \partial\Omega, \end{cases} $$ where $\lambda>0$ is a small parameter, $0More