In 2 dimensional models we do not only have the problem of avoiding the singularity but we choose a potential that allows to mimic the average influence that the planet would have on the vertical gas column.
A value often used for the smoothing is $`\epsilon = 0.7H`$ where $`H`$ is the disc scale height.
To select this choice for a simulation in the config file: