In this Wikipedia article on the one dimensional Ising model they use the following factorization
$$ \displaystyle \sum_{\sigma_1,\sigma_2,\cdots,\sigma_N} e^{a\sigma_1} e^{a\sigma_2}\cdots e^{a\sigma_N} = 2\displaystyle\prod_{j=1}^{N}\sum_{\sigma_i}e^{a\sigma_i} = 2 \left[e^{-a}+e^{a}\right]^N \,. $$ The first sum is over all possible sequences ($\sigma_1,\sigma_2,\cdots, \sigma_N)$ where $\sigma_i \in \{-1,1\}$.
Does somebody have a reference for this factorization and/or able to give a proof?
The factor $2$ is not plausible. Taking $N=1$ we have \begin{align*} \sum_{\sigma_1}e^{a\sigma_1}=e^{a}+e^{-a} \end{align*}
Comment:
In (1) we use the product symbol $\prod_{j=1}^N$ to write the product of $N$ factors.
In (2) we reorder the sum by selecting all subsets $S\subseteq [N]=\{1,2,\ldots,N\}$ for which the sign is positive and $[N]\setminus S$ for which the sign is negative.
In (3) we use the fact that there are $\binom{N}{k}$ subsets $S\subseteq [N]$ with size $k$.
in (4) we apply the binomial theorem.