In reply to the Comment by O. Neufeld et al. [Phys. Rev. A 105, 047101 (2022)], we argue that the conclusions of Phys. Rev. A 103, 043106 (2021) remain valid. We disprove the claim that the unphysical even-order harmonics originate from convergence issues related to reflections at the boundary of the simulation box. By additional calculations, we show that such reflections perturb the high-order harmonic generation spectra by oscillations with periods much smaller than the distance between the harmonics. We also demonstrate that the convergence argument of the Comment, in contrast to our multielectron excitations argument, cannot explain why there are no unphysical even-order harmonics in one-electron systems. Moreover, we show that the argument put forward in the Comment to conclude that the time-dependent Kohn-Sham equations are superior to the time-dependent natural Kohn-Sham equations is not valid.