I have a (simple?) question about Fourier transforms.
Consider a 1D Hamiltonian of the form \begin{equation} H = -JS\sum_{j = 1}^{N-1}a_{j+1}^\dagger a_j + a_j^\dagger a_{j+1} - a^\dagger_{j+1}a_{j+1} - a^\dagger_j a_j \end{equation} where $J$ is a coupling between two nearest neighbours, and $S$ is the spin projection along some z-axis, i.e a standard ferromagnetic chain with $N$ lattice sites and lattice spacing $d$.
To diagonalize this one typically introduces the fourier transformed creation/annhilation operators \begin{equation} a_j = \frac{1}{\sqrt{N}} \sum_k e^{ik jd}a_{k}. \end{equation} This is fine as long as we assume periodic boundary conditions such that $a_{j+N}=a_{j}$.
Now consider the case when we let $N\rightarrow \infty$. In this case, it no longer makes sense to use periodic boundary conditions. How then do we define a Fourier transform in order to diagonalize such a problem?
Is it as simple as just writing \begin{equation} a_j = \int dk e^{ikjd}a_{k} \end{equation}?