For $m\geq 2$, let $P_m=4X^{4m}+X^{4m-2}+X^2+4$. I have checked that for $1\leq m \leq 10$,
1) $P_m$ is irreducible over $\mathbb Q$
2) All its roots have modulus $1$
3) Its Galois group has order $4^m(m!)$ (computation made with online MAGMA calculator )
My questions : does anyone know a proof or counterexample for any of those claims ?
What makes this problem interesting in my view is that (according to MAGMA) the Galois group is non-abelian for $m\geq 2$. So we're not in an entirely cyclotomic context here, contrary to what might be expected from the fact that the roots have modulus $1$.
Update : There is a unique polynomial $Q_m$ of degree $m$ such that $P_m(X)=X^{2m}Q_m(X^2+\frac{1}{X^2})$. Some small values of $Q_m$ follow below. As guessed by in user760780's comment below, the Galois group of $Q_m$ is exactly $S_m$.
$$ \begin{array}{lcl} Q_1(y) & = & 4y + 2 \\ Q_2(y) & = & 4y^2 + y - 8 \\ Q_3(y) & = & 4y^3 + y^2 - 12y - 2 \\ Q_4(y) & = & 4y^4 + y^3 - 16y^2 - 3y + 8 \\ Q_5(y) & = & 4y^5 + y^4 - 20y^3 - 4y^2 + 20y + 2 \\ Q_6(y) & = & 4y^6 + y^5 - 24y^4 - 5y^3 + 36y^2 + 5y - 8 \\ Q_7(y) & = & 4y^7 + y^6 - 28y^5 - 6y^4 + 56y^3 + 9y^2 - 28y - 2 \\ Q_8(y) & = & 4y^8 + y^7 - 32y^6 - 7y^5 + 80y^4 + 14y^3 - 64y^2 - 7y + 8 \\ Q_{9}(y) & = & 4y^9 + y^8 - 36y^7 - 8y^6 + 108y^5 + 20y^4 - 120y^3 - 16y^2 + 36y + 2 \\ Q_{10}(y) & = & 4y^{10} + y^9 - 40y^8 - 9y^7 + 140y^6 + 27y^5 - 200y^4 - 30y^3 + 100y^2 + 9y - 8 \\ \end{array} $$