I haven't tested this with Pyxaid, but in Pyxaid2 you can definitely do this. Just make sure you use the proper settings (non-relativistic, non-spin-polarized calculations). Then, you don't need to do anything else special other then specifying your (still single!) k-point in the template input file (e.g.
x0.scf.in) in the step2. Mathematically, this should also work with Pyxaid (just ensure do_complete = 0), but this has not been fully explored with that code.