The Morse property of limit functions appearing in mean field equations on surfaces with boundary

In this paper we study the Morse property for functions related to limit functions of mean field equations on a smooth, compact surface $\mathrm{\Sigma }$ with boundary $\mathrm{\partial }\mathrm{\Sigma }$. Given a Riemannian metric $g$ on $\mathrm{\Sigma }$ we consider functions of the form $$f_g(x):=\sum _{i=1}^m{\sigma }_i^2{R}^g(x_i)+\sum _{i,j=1\text{ı}\mathrm{\ne }j}^m{\sigma }_i{\sigma }_j{G}^g(x_i,x_j)+h(x_1,\dots ,x_m),$$ where ${\sigma }_i\mathrm{\ne }0$ for $i=1,\dots ,m$, $G^g$ is the Green function of the Laplace-Beltrami operator on $(\mathrm{\Sigma },g)$ with Neumann boundary conditions, $R^g$ is the corresponding Robin function, and $h\in {\mathcal{C}}^2({\mathrm{\Sigma }}^m,\mathbb{R})$ is arbitrary. We prove that for any Riemannian metric $g$, there exists a metric $\tilde{g}$ which is arbitrarily close to $g$ and in the conformal class of $g$ such that $f_{\tilde{g}}$ is a Morse function. Furthermore we show that, if all ${\sigma }_i>0$, then the set of Riemannian metrics for which $f_g$ is a Morse function is open and dense in the set of all Riemannian metrics.