distinguished from each other using the La-Cl interionic distance, thereby allowing proper sampling of all the metastable states. In Table , high values of the energy barrier for dissociation can be seen compared to the thermal energy, k B T. This creates a kinetic bottleneck in unbiased simulations, which underlines the motivation to use WMetaD for this system. Further, it can be inferred that PIM predicts lower energy barriers for dissociation compared to those from N-POL MD for all reactions at subcritical as well as supercritical conditions. PIM predicts increasing energy barriers for successive dissociation reactions at both conditions. N-POL MD predicts almost constant energy barriers at 773 K. At 473 K, an increase in energy barrier from reaction 1 to 2 followed by a decrease in reaction 3 is observed.