Poisson equation is given by $\nabla^2V=\frac{\rho}{\epsilon_0}$. Here $\rho$ indicates a volumic charge distribution, which is known in the region $\Omega$ where we solve the Poisson equation.
Is it correct using instead a known surface charge distribution $\sigma$ given in $\Omega$?
I think that the answer is yes, but only in a distributional sense (Dirac delta). In fact the electric field presents a discontinuity near a surface with charge density $\sigma$, then over the surface is not valid the local Gauss law $\vec{\nabla} \cdot \vec{E}=\frac{\rho}{\epsilon_0}$ (unless we interpret $\rho$ as a distribution $\delta$).
Thank you.