We derive a upper bound on the free energy of a Bose gas system at density ̺ and temperature T . In combination with the lower bound derived previously by Seiringer [15], our result proves that in the low density limit, i.e., when a̺ ≪ 1, where a denotes the scattering length of the pair-interaction potential, the leading term of ∆f the free energy difference per volume between interacting and id...