We consider the following family of random trigonometric polynomials of the form Sn(θ):=1 n ∑k=1nakcos(kθ)+bksin(kθ), where the variables {ak,bk}k≥1 are i.i.d. and symmetric, defined on a common probability space (Ω,F,P). In the seminal paper [SZ54], Salem and Zygmund proved that for Rademacher coefficients, P−almost surely and for all t∈R limn→+∞ 1 2π∫02πeitSn(θ)dθ=e−t2∕2. To the best of our knowledge, the natural question of the fluctuations in the above almost sure limit has not been tackled so far and is precisely the object of this article. Namely, for general i.i.d. symmetric random coefficients having a finite sixth-moment and for a large class of continuous test functions ϕ we prove that n1 2π∫02πϕ(S n(θ))dθ−∫Rϕ(t)e−t2 2dt 2π→n→∞LawN0,Σϕ2, where the limit variance is given by Σϕ2:=σϕ2+c2(ϕ)22E[a14]−3. Here, the constant σϕ2 is explicit and corresponds to the limit variance in the case of Gaussian coefficients and c2(ϕ) is the coefficient of order 2 in the decomposition of ϕ in the Hermite polynomial basis. It thus turns out that the fluctuations are Gaussian and not universal since they involve the kurtosis of the random coefficients. The method we develop here is robust and allows to use Wiener chaos expansions for proving CLTs for functionals of general, non-necessarily Gaussian, random fields. It combines contributions from the Malliavin–Stein method, noise stability results, and algebraic considerations around Newton sums and the Newton–Girard formula.