$\newcommand{\erf}{\operatorname{erf}}$ $\newcommand{\const}{\mathrm{const}}$
If we have $n$-dimensional periodic function we can decompose it in $n$--dimensional Fourier series:
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_{\mathbb{m}} e^{ i \mathbf{m}\cdot \mathbf{x}} \label{equ-J.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{equ-J.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{equ-J.3} \end{equation} where $\Omega=(0,1)^n$ is $n$-dimensional unit cube.
We need slightly generalize these formulae.
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{equ-J.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$ if \begin{equation} u(\mathbf{x}+\mathbf{y})= u(\mathbf{x})\qquad \forall \mathbf{y}\in \Gamma\ \forall \mathbf{x}. \label{equ-J.5} \end{equation} In the previous section $\Gamma= (2\pi\mathbb{Z})^n$. Let us change coordinate system so that $\Gamma$ becomes $(2\pi\mathbb{Z})^n$, apply (\ref{equ-J.1})--(\ref{equ-J.3}) and then change coordinate system back. We get \begin{equation} u(\mathbf{x})= \sum_{\mathbf{m}\in \Gamma^*} c_{\mathbb{m}} e^{ i \mathbf{m}\cdot \mathbf{x}} \label{equ-J.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{equ-J.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{equ-J.8} \end{equation} where $|\Omega|$ is a volume of $\Omega$ and
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.
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$ (sufficiently fast decaying): \begin{equation} u(\mathbf{x})= \iiint_{\Omega^*} u(\mathbf{k};\mathbf{x})\,d^n\mathbf{k} \label{equ-J.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{n}} u(\mathbf{x}+\mathbf{l}). \label{equ-J.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{equ-J.12} \end{equation} The proof 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}$.