Compute configuration integral

Dear LAMMPS community,
I mean the configuration integral Q_N in the partition function:
изображение
Is it possible to compute it at runtime for given atoms (atom IDs or group) and a specified temperature? What is the way to do it?
Or probably it is possible to compute the whole partition function?

Maybe start here: Extraction of configurational entropy from molecular simulations via an expansion approximation | The Journal of Chemical Physics | AIP Publishing

Please note that the quantity you are thinking about is not trivial to calculate, and you are likely to get stuck without the help of a local, experienced mentor or supervisor.