Assessment of Anharmonicities in Clusters: Developing and Validating a Minimum-Information Partition Function

Precise thermodynamic calculations are essential for understanding the dynamics of cluster systems and new particle formation. However, the widely employed harmonic statistical mechanical approach often falls short in terms of accuracy. In this study, we present an improved statistical model that incorporates vibrational anharmonicity via a novel partition function that requires only one additional system-specific input parameter. In addition to considering vibrational aspects, we also account for anharmonicity related to the configurational space. The role of anharmonicities is thoroughly examined in the case of general clusters, where the complete sets of conformers, mechanically stable spatial arrangements, are known up to clusters composed of 14 monomers. By performing consistent Monte Carlo simulations on these systems, we benchmark the statistical model's efficacy in reproducing key thermodynamic properties (formation free energy and potential energy) in the classical limit. The model exhibits exceptional alignment with simulations, accurately reproducing free energies within a precision of 2kT and reliably capturing cluster melting temperatures. Furthermore, we demonstrate the significance and applicability of the model by reproducing thermodynamic barriers in homogeneous gas-phase nucleation of larger clusters. The transferability of our developed approach extends to more complex molecular systems and bears relevance for atmospheric multicomponent clusters, in particular.