Confusion relating to discretizing a Lagrangian and the Fourier transform conventions thereof

Question

My question boils down to how the Fourier transform is discretized when we discretize the field $phi(x)$ in a Lagrangian $mathcal{L}$. To put things on a concrete footing, consider the following Lagrangian (from exercise 17.5 of Lancaster and Blundell):
begin{equation}
mathcal{L} = frac{1}{2}left(frac{partialphi(x)}{partial x}right)^2 + frac{m^2}{2}phi(x)^2,
end{equation}
in one dimension. (Please correct any mistakes in this preamble!)
Ordinarially, we would define the Fourier transform relations as
begin{equation}
phi(x) = intfrac{mathrm{d}p}{2pi}tilde{phi}(p)mathrm{e}^{ipcdot x}longleftrightarrow tilde{phi}(p) =intmathrm{d}xphi(x)mathrm{e}^{-ipcdot x},
end{equation}
leading to the delta function relation
begin{equation}
frac{1}{2pi}intmathrm{d}p;mathrm{e}^{ipcdot y}=delta(y).
end{equation}
L&B introduce discretization into the problem by defining the relationships:
begin{equation}
phi_j=frac{1}{sqrt{Na}}sum_ptilde{phi}_pmathrm{e}^{ipja}longleftrightarrow tilde{phi}_p=frac{1}{sqrt{Na}}sum_j{phi}_jmathrm{e}^{-ipja}.
end{equation}
Where $a$ is the lattice spacing, $N$ the number of lattice points and $j$ an index which labels the lattice points. However, this seems to imply that
begin{equation}
frac{1}{Na}sum_kmathrm{e}^{ikja} = delta_{j,0}.
end{equation}
Here is where I am confused. Surely if $j=0$ in this equation, then $sum_k = N$ - so where does the $a$ come from? Generally I am unsure how one derives these discretised Fourier relations. If anyone could point me in the right direction to actually understanding where they come from, in particular relating to the prefactors $1/sqrt{Na}$, then I would be very appreciative.
EDIT: Why $intmathrm{d}xrightarrow Nasum_j$ also confuses me deeply - this time, it seems to me that $mathrm{d}xrightarrow a$, and hence that the $N$ should not be there!

Matt Miguel · Accepted Answer

Intuitively there is a relationship between periodicity and being discrete between dual spaces (time vs frequency or space vs wavenumber).
Periodicity in one space implies being discrete in the dual space and vice versa.
There are four cases:

The continuous Fourier transform is suited when representations in both dual spaces are non-periodic and are non-discrete.

The discrete-time Fourier transform (DTFT) is suited when representation in one space is discrete but not periodic and in its dual space is periodic but not discrete.

The Fourier series is suited when representation in one space is periodic and not discrete and in its dual space is discrete but not periodic.

The discrete Fourier transform (DFT) is suited when representation in a space is both periodic and discrete and it is also that way in the dual space.

Note that 2) and 3) can be used to describe the same situation, just interpreted differently depending on which one you think is the "original/natural" space and which is the transform.
Regarding the $frac{1}{N}$ term - this shows up whenever go from a periodic function to its dual (which must be discrete) - you divide by the period.
In the wiki DTFT inverse transform definition linked to above, there is a $frac{1}{2 pi}$ factor for this factor.
In the wiki Fourier series definition linked to above, there is a $frac{1}{P}$ factor for this factor.
For the DFT, both dual representations are periodic. In some definitions of the DFT, the $frac{1}{N}$ factor is excluded in the forward transform and present in the inverse transform. You can also define a version of the DFT that basically symmetrically shares the square root of the term in both directions $frac{1}{sqrt{N}}$. The only requirements really for the DFT and its inverse are that the exponents have opposite signs and the product of the two scaling factors is $frac{1}{N}$.
Hope that helps with the intuition behind the different Fourier transform variations.

Confusion relating to discretizing a Lagrangian and the Fourier transform conventions thereof

One Answer

Add your own answers!

Ask a Question