If you go to 1x1 then at least turn on the Wigner-Seitz scheme, else you are making a Gamma point phonon approximation. I would probably check if some of the machine-learned forcefields like MACE or M3GNet give at least decent phonon dispersions, and if so use them, for a much, much, much shorter calculation time, which hopefully is at least equally qualitative as a very approximate DFT simulation.