Let's say I have a single particle hamiltonian in a periodic potential, for example a 1D lattice such that: $$H = -\frac{\partial_x^2}{2m} + V(x) $$ with $ V(x+a) = V(x)$ where $a$ is the lattice spacing between the atoms or sites. It is known by Bloch's theorem that a solution to such a system will have the form $$\psi_{k}(x)=e^{ikx}u_k(x)$$ where $u_k(x+a)=u_k(x)$.
My questions is about the boundary conditions. If we take $$\psi(x+Na) = \psi(x)$$ we get, if $N$ is large enough, a lot of different values for $k$ in the first Brillouin zone: $$k=\frac{2\pi n}{N} \text{ with }-\frac{\pi}{a}<k<\frac{\pi}{a},$$ so we get a band of possible states.
In this case we can define Wannier functions which using Fourier over the wave-functions: $$\phi(x-R) = \sum_k e^{-ik R} \psi_k(x)$$ where the summation is over all the $k$'s in the first Brillouin zone.
But if I take the B.C $$\psi(x+a) = \psi(x)$$ I get a single value for the momentum in each Brillouin zone $$k = 0, \pm 2\pi , \pm 4\pi,...$$ Is it still possible to define a Wannier function for such a state? I mean what will the Fourier be like if we have a single possible value of $k$??