Singular limits for Liouville-type equations on the flat two-torus

We consider the problem $$\begin{aligned} -\Delta u=\varepsilon ^{2}e^{u}- \frac{1}{|\Omega |}\int _\Omega \varepsilon ^{2} e^{u}+ {4\pi N\over |\Omega |} - 4 \pi N\delta _p, \quad \text{ in} {\Omega }, \quad \int _\Omega u=0 \end{aligned}$$ in a flat two-torus $\Omega $ with periodic boundary conditions, where $\varepsilon >0,\,|\Omega |$ is the area of the $\Omega $ , $N>0$ and $\delta _p$ is a Dirac mass at $p\in \Omega $ . We prove that if $1\le m<N+1$ then there exists a family of solutions $\{u_\varepsilon \}_{\varepsilon }$ such that $\varepsilon ^{2}e^{u_\varepsilon }\rightharpoonup 8\pi \sum _{i=1}^m\delta _{q_i}$ as $\varepsilon \rightarrow 0$ in measure sense for some different points $q_{1}, \ldots , q_{m}$ . Furthermore, points $q_i$ , $i=1,\dots ,m$ are different from $p$ .