  Operator method and cumulant expansion are used for nonperturbative
calculation of the partition function and the free energy in quantum
statistics. It is shown for Boltzmann diatomic molecular gas with some model
intermolecular potentials that the zeroth order approximation of the proposed
method interpolates the thermodynamic values with rather good accuracy in the
entire range of both the Hamiltonian parameters and temperature. The systematic
procedure for calculation of the corrections to the zeroth order approximation
is also considered.
