Look at the 2-adic integers
$$f(x)= (x^2+x)^{2^n}+1= \prod_{a=1}^{2^n} (x^2+x-\zeta_{2^{n+1}}^{2a+1})\in \Bbb{Q}_2[x]$$
Let $L=\Bbb{Q}_2(\beta)$ with $\beta$ a root of $x^2+x-\zeta_{2^{n+1}}$. It contains $K = \Bbb{Q}_2(\zeta_{2^{n+1}})$ the splitting field of $\Phi_{2^{n+1}}(x)=x^{2^n}+1 \in \Bbb{Q}_2[x]$.
Let $(\pi)$ be $O_K$'s maximal ideal.
From $\zeta_{2^{n+1}}^{2^{n+1}} = 1$ and that $O_K/(\pi)$ is a field with $2^m$ elements we know $\zeta_{2^{n+1}} \equiv 1 \bmod (\pi)$ and $\sum_{l=0}^{2a} \zeta_{2^{n+1}}^l \equiv 2a+1\equiv1 \bmod (\pi)$ which means
$$v_2(\zeta_{2^{n+1}}^{2a+1}-1)=v_2(\zeta_{2^{n+1}}-1)+v_2(\sum_{l=0}^{2a} \zeta_{2^{n+1}}^l)=v_2(\zeta_{2^{n+1}}-1)$$
Together with $$\Phi_{2^{n+1}}(x+1)= (x+1)^{2^n}+1=\prod_{l=1}^{2^n} (x+1-\zeta_{2^{n+1}}^{2a+1}), \qquad \prod_{l=1}^{2^n} (1-\zeta_{2^{n+1}}^{2a+1})=\Phi_{2^{n+1}}(1)=2 $$
we find $v_2(\zeta_{2^{n+1}}-1) =\frac{\sum_{l=1}^{2^n} v_2(1-\zeta_{2^{n+1}}^{2a+1})}{2^n}= \frac{v_2(2)}{2^n}=\frac1{2^n}=\frac1{\deg(\Phi_{2^{n+1}})}$.
Thus $\Phi_{2^{n+1}}$ is irreducible, $K/\Bbb{Q}_2$ is totally ramified and $\pi=1-\zeta_{2^{n+1}}$ is an uniformizer.
We obtain
$$x^2+x-\zeta_{2^{n+1}} \equiv x^2+x-1 \in O_K/(\pi)[x]\cong \Bbb{Z/2Z}[x]$$
Which is irreducible, and since $L$ is $x^2+x-\zeta_{2^{n+1}}\in K[x]$ splitting field it means $$[L:K] = 2, \qquad [L:\Bbb{Q}_2]=[L:K][K:\Bbb{Q}_2] = 2^{n+1}=\deg(f)$$
And hence $f$ is irreducible over $\Bbb{Q}_2$ and it stays irreducible over $\Bbb{Q}$.