$\renewcommand{\Re}{\operatorname{Re}}$ $\renewcommand{\Im}{\operatorname{Im}}$ $\newcommand{\erf}{\operatorname{erf}}$ $\newcommand{\dag}{\dagger}$ $\newcommand{\const}{\mathrm{const}}$ $\newcommand{\arcsinh}{\operatorname{arcsinh}}$

Appendix 4.B. Multidimensional Fourier series

  1. $2\pi$-periodic case
  2. General case
  3. Special decomposition

$2\pi$-periodic case

Let function $u(\mathbf{x})$, $\mathbf{x}=(x_1,x_2,\ldots,x_n)$ be $2\pi$-periodic with respect to each variable $x_1,x_2,\ldots,x_n$. Then \begin{equation} u(\mathbf{x})= \sum_{\mathbf{m}\in \mathbb{Z}^n} c_{\mathbf{m}} e^{ i \mathbf{m}\cdot \mathbf{x}} \label{eq-4.B.1} \end{equation} with \begin{equation} c_{\mathbf{m}} =(2\pi)^{-n} \iiint_\Omega e^{- i \mathbf{m}\cdot \mathbf{x}} u(\mathbf{x})\,d^nx \label{eq-4.B.2} \end{equation} and \begin{equation} \sum_{\mathbf{m}\in \mathbb{Z}^n} |c_{\mathbf{m}}|^2 =(2\pi)^{-n} \iiint_\Omega |u(\mathbf{x})|^2\,d^n x \label{eq-4.B.3} \end{equation} where $\Omega=(0,1)^n$ is $n$-dimensional unit cube. Here and below we write $n$-dimensional integral as $\iiint$.

We need slightly generalize these formulae.

General case

Definition 1. Let $\Gamma$ be $n$-dimentional lattice. It means that there are $n$ linearly independent vectors $\mathbf{e}_1, \ldots, \mathbf{e}_n$ and \begin{equation} \Gamma = \{ (k_1 \mathbf{e}_1 + k_2 \mathbf{e}_2+\ldots +k_n \mathbf{e}_n:\, k_1,k_2,\ldots,k_n\in \mathbb{Z}\} \label{eq-4.B.4} \end{equation}

Remark 1. The same lattice $\Gamma$ is defined by vectors $\mathbf{e}'_1, \ldots, \mathbf{e}'_n$ with $\mathbf{e}'_j=\sum _k \alpha_{jk}\mathbf{e}_k$ with integer coefficients if and only if the determinant of the matrix $(\alpha_{jk})_{j,k=1,\ldots,n}$ of coefficients is equal $\pm 1$.

Definition 2. Let $\Gamma$ be $n$-dimentional lattice. We call $u(\mathbf{x})$ periodic with respect to $\Gamma$ or simply $\Gamma$-periodic if \begin{equation} u(\mathbf{x}+\mathbf{y})= u(\mathbf{x})\qquad \forall \mathbf{y}\in \Gamma\ \forall \mathbf{x}. \label{eq-4.B.5} \end{equation}

In the previous Subsection $\Gamma= (2\pi\mathbb{Z})^n$. Let us change coordinate system so that $\Gamma$ becomes $(2\pi\mathbb{Z})^n$, apply (\ref{eq-4.B.1})--(\ref{eq-4.B.3}) and then change coordinate system back. We get \begin{equation} u(\mathbf{x})= \sum_{\mathbf{m}\in \Gamma^*} c_{\mathbf{m}} e^{ i \mathbf{m}\cdot \mathbf{x}} \label{eq-4.B.6} \end{equation} with \begin{equation} c_{\mathbf{m}} =|\Omega |^{-1} \iiint_\Omega e^{- i \mathbf{m}\cdot \mathbf{x}} u(\mathbf{x})\,d^n x \label{eq-4.B.7} \end{equation} and \begin{equation} \sum_{\mathbf{m}\in \Gamma^*} |c_{\mathbf{m}}|^2 =|\Omega|^{-1} \iiint_\Omega |u(\mathbf{x})|^2\,d^n x \label{eq-4.B.8} \end{equation} where $|\Omega|$ is a volume of $\Omega$ and

Definition 3.

  1. $\Omega =\{x_1 \mathbf{e}_1+\ldots + x_n \mathbf{e}_n:\, 0< x_1<1,\ldots,0< x_n<1 \}$ is an elementary cell;
  2. $\Gamma^* =\{\mathbf{m}:\, \mathbf{m}\cdot \mathbf{y}\in 2\pi \mathbb{Z}\ \ \forall \mathbf{y}\in \Gamma\}$ is a dual lattice; it could be defined by vectors $\mathbf{e}^*_1,\ldots, \mathbf{e}^*_n$ such that \begin{equation} \mathbf{e}^*_j \cdot \mathbf{e}_k=2\pi \delta_{jk}\quad \forall j,k=1,\ldots,n \label{eq-4.B.9} \end{equation} where $\delta_{jk}$ is a Kronecker symbol;
  3. $\Omega^* =\{k_1 \mathbf{e}_1^*+\ldots + k_n \mathbf{e}_n^*:\, 0< k_1<1,\ldots,0< k_n<1 \}$ is a dual elementary cell.

Remark 2. We prefer to use original coordinate system rather than one with coordinate vectors $(2\pi)^{-1}\mathbf{e}_1,\ldots, (2\pi)^{-1}\mathbf{e}_n$ because the latter is not necessarily orthonormal and in it Laplacian will have a different form.

Special decomposition

These notions are important for studying the band spectrum of the Schrödinger operator $-\Delta +V(\mathbf{x})$ with periodic (with respect to some lattice $\Gamma$) potential in the whole space which has applications to the Physics of crystals. For this the following decomposition is used for functions $u(\mathbf{x})$ in the whole space $\mathbb{R}^n$

Theorem 1. Let $u(\mathbf{x})$ be sufficiently fast decaying function on $\mathbb{R}^n$. Then \begin{equation} u(\mathbf{x})= \iiint_{\Omega^*} u(\mathbf{k};\mathbf{x})\,d^n\mathbf{k} \label{eq-4.B.10} \end{equation} with \begin{equation} u(\mathbf{k};\mathbf{x})= (2\pi)^{-n}|\Omega| \sum_{\mathbf{l}\in \Gamma} e^{-i\mathbf{k}\cdot \mathbf{l}} u(\mathbf{x}+\mathbf{l}). \label{eq-4.B.11} \end{equation} Here $u(\mathbf{k};\mathbf{x})$ is quasiperiodic with quasimomentum $\mathbf{k}$ \begin{equation} u(\mathbf{k};\mathbf{x}+\mathbf{y})= e^{i\mathbf{k}\cdot\mathbf{y}}u(\mathbf{k};\mathbf{x})\qquad \forall \mathbf{y}\in \Gamma\ \forall \mathbf{x}. \label{eq-4.B.12} \end{equation}

Proof. Observe that since $u$ is sufficiently fast decaying series in (\ref{eq-4.B.11}) converges and one can see easily that it defines quasiperiodic with quasimomentum $\mathbf{k}$ function.

The proof of (\ref{eq-4.B.10})is trivial as $\iiint_{\Omega^*} e^{-i\mathbf{k}\cdot \mathbf{l}} \,d^n\mathbf{k}=|\Omega^*|$ as $\mathbf{l}=0$ and $0$ as $0\ne \mathbf{l}\in \Gamma$, and $|\Omega^*|=(2\pi)^n|\Omega|^{-1}$.

$\Leftarrow$  $\Uparrow$  $\Rightarrow$