We study the existence of solution to the problemwhere Q ≥ 0, κ ∈ (0, ∞) and n ≥ 3. Using ODE techniques Martinazzi for n = 6 and Huang-Ye for n = 4m + 2 proved the existence of solution to the above problem with Q ≡ const > 0 and for every κ ∈ (0, ∞). We extend these results in every dimension n ≥ 5, thus completely answering the problem opened by Martinazzi. Our approach also extends to the case in which Q is non-constant, and under some decay assumptions on Q we can also treat the cases n = 3 and 4.