Isodesmic reaction modelling