The basis vectors for this null-space can be assembled into an (N + 2k) × N matrix U(x) with orthonormalization condition U†U = 1.
For the null-space, the basis vectors can be organized into an (N + 2k) × N matrix U(x), adhering to the orthonormality condition U†U = 1.