A path integral methodology for obtaining thermodynamic properties of nonadiabatic systems using Gaussian mixture distributions