We consider a system of N bosons in the limit N → ∞, interacting through singular potentials. For initial data exhibiting Bose-Einstein condensation, the many-body time evolution is well approximated through a quadratic fluctuation dynamics around a cubic non-linear Schrödinger equation of the condensate wave function. We show that these fluctuations satisfy a (multi-variate) central limit theorem.where α (n) ∈ L 2 ⊥ϕ N,t (R 3 ) ⊗sn for all n = 1, · · · , N . Then,This unitary satisfies the following properties proven in [25] U ϕ N,t a * (ϕ N,t )a(ϕ N,t )U * ϕ N,t =N − N + (t) U ϕ N,t a * (ϕ N,t )a(f )U * ϕ N,t = N − N + (t)a(f )