I am interesting in solutions to Poisson equation

$$\triangle \varphi = 4 \pi \rho \qquad (1)$$

defined on 3-dimensional oriented Riemannian manifold $(M,g)$, where $g$ is metric and $\triangle = \triangle[g]$ is Laplacian, and $\rho$ (density) is a smooth function on $M$ or a distribution. The first simple question: i) is it correct that for compact $M$ equation (1) has a solution if and only if $\int_{M} \rho dVol_g = 0$. If so, the fundamental solution to eq. (1) with $\rho = m \delta_p$, i.e.

$$\triangle \varphi = 4 \pi m \delta_p \qquad (2) $$

where $\delta_p$ is $\delta$-function located at point $p$ and $m \neq 0$ (mass), does not exist. However, in the case when $\rho = m_1 \delta_{p_1} + m_2 \delta_{p_2}$, where $p_1 \neq p_2$ and $m_1 + m_2 =0$ the solution to eq. (1) does exist - is it correct? Another two questions are about the solutions to eq. (2) for the case when : ii) $M = {\mathbb R} \times {\mathbb R} \times S^1$ and ii) $M = {\mathbb R} \times S^1 \times S^1$ with standards metrics (induced from ${\mathbb R}^3$). Could you give me some references, where explicit (analitic) fundamental solutions to Poisson eq. (2) with certain asymptotical conditions at infinity are written for these two cases.