The algebra of phonons which we excite by Raman is the Bose-Einstein algebra (commutator of generating vs annihilating operator is 1) which is why we treat them as Bosons and therefore with the Bose-Einstein population. At how low temperatures are you working that it makes a difference from Boltzmann, though?
Dear Zeineb Raddaoui, I guess that the Raman modes are caused by optical phonons, and phonons are bosons, so the correct quantum distribution is the Bose-Einstein statistics.