Numerical models of the evolution of interstellar and intergalactic plasmas often assume that the adiabatic parameter {gamma} (the ratio of the specific heats) is constant (5/3 in monoatomic plasmas). However, {gamma} is determined by the total internal energy of the plasma, which depends on the ionic and excitation state of the plasma. Hence, the adiabatic parameter may not be constant across the range of temperatures available in the interstellar medium. We aim to carry out detailed simulations of the thermal evolution of plasmas with Maxwell-Boltzmann and non-thermal ({kappa} and n) electron distributions in order to determine the temperature variability of the total internal energy and of the adiabatic parameter. The plasma, composed of H, He, C, N, O, Ne, Mg, Si, S, and Fe atoms and ions, evolves under collisional ionization equilibrium conditions, from an initial temperature of 10^9^K. The calculations include electron impact ionization, radiative and dielectronic recombinations and line excitation. The ionization structure was calculated solving a system of 112 linear equations using the Gauss elimination method with scaled partial pivoting. Numerical integrations used in the calculation of ionization and excitation rates are carried out using the double-exponential over a semi-finite interval method. In both methods a precision of 10^-15^ is adopted. The total internal energy of the plasma is mainly dominated by the ionization energy for temperatures lower than 8x10^4^K with the excitation energy having a contribution of less than one percent. In thermal and non-thermal plasmas composed of H, He, and metals, the adiabatic parameter evolution is determined by the H and He ionizations leading to a profile in general having three transitions. However, for {kappa} distributed plasmas these three transitions are not observed for {kappa<15} and for {kappa<5} there are no transitions. In general, {gamma} varies from 1.01 to 5/3. Lookup tables of the {gamma} parameter are presented as supplementary material.