In this paper we consider an energy operator (a free Hamiltonian), in the second-quantized approach, for the multiparameter quon algebras: aia†j − qija†jai = δij, i, j ∊ I with (qij)i,j∊I any Hermitian matrix of deformation parameters. We obtain an elegant formula for normally ordered (sometimes called Wick-ordered) series expansions of number operators (which determine a free Hamiltonian). As a main result (see theorem 1) we prove that the number operators are given, with respect to a basis formed by 'generalized Lie elements', by certain normally ordered quadratic expressions with coefficients given precisely by the entries of the inverses of Gram matrices of multiparticle weight spaces. (This settles a conjecture by Meljanac S and Perica A (1994 J. Phys. A: Math. Gen.27 4737–44).) These Gram matrices are Hermitian generalizations of the Varchenko matrices, associated with a quantum (symmetric) bilinear form of diagonal arrangements of hyperplanes. The solution of the inversion problem of such matrices in Meljanac S and Svrtan D (1996 Math. Commun.1 1–24 (theorem 2.2.17)), leads to an effective formula for the number operators studied in this paper. The one-parameter case, in the monomial basis, was studied by Zagier, Stanciu and Møller.