Let Lqr(Ω) be the usual scale of Sobolev spaces and let ΔN be the Neumann Laplacian in an arbitrary Lipschitz domain Ω. We present an interpolation based approach to the following question: for what range of indices does map isomorphically onto Lqr(Ω)/ℝ?