The ability of computational methods to predict the structures and energetics that determine the equilibrium of solid state mechanochemical reactions has been assessed. Two previously characterised base-catalysed metathesis reactions between aromatic disulfides are studied using crystal structure prediction methods and lattice energy calculations that combine molecular electronic structure methods with anisotropic atom–atom potentials. We find that lattice energy searches locate three of the six crystal structures as global minima on their respective crystal energy landscapes. The remaining structures are less successfully predicted, due to problems modelling relative conformational energies due to limitations of the density functional theory method for calculating intramolecular energies. Prediction of the overall reaction energies proves challenging for current methods, but the results show promise as a base on which to build more accurate and reliable approaches.