Solution (source code)

= Solution

The absolutely convergent <matrix exponential> series gives
$$
U^\dagger=(e^{iH})^\dagger=e^{(iH)^\dagger}=e^{-iH}
$$
for Hermitian $H$. The exponents $iH$ and $-iH$ commute, so their exponentials multiply to the exponential of their sum:
$$
\boxed{U^\dagger U=e^{-iH}e^{iH}=I,\qquad UU^\dagger=I.}
$$
Hence $U$ is unitary. Alternatively, an <orthonormal eigenbasis> of $H$ turns the exponential <eigenvalues> into $e^{i\lambda_j}$, all on the <unit circle>.