Benchmarking of DFT methods using experimental free energies and volumes of activation for the cycloaddition of alkynes to cuboidal Mo3S4 clusters

Here, the kinetics of the concerted [3 + 2] cycloaddition reaction between the [Mo3(μ3‐S)(μ‐S)3Cl3(dmen)3]+ (dmen = N,N′‐dimethyl‐ethylenediamine) ([1]+) cluster and various alkynes to form dithiolene derivatives is thoroughly studied, with measurements at different temperatures and pressures allowing the determination of the free energies and volumes of activation. These parameters, together with the available single‐crystal X‐ray diffraction structures, are used to test a number of commonly used density functional theory (DFT) methods from Jacob's ladder, as well as the effects associated with the size of the basis sets, the way in which solvent effects are taken into account, or the inclusion of dispersion effects. Overall, a protocol that leads to average deviations between experimental and computed ΔV‡ and ΔG‡ values similar to the uncertainty of the experimental measurements is obtained. Despite the pressure dependence of the rate of reactions allows obtaining its activation volume, this magnitude is rarely studied from a computational viewpoint. Along with their free energy barriers, here, the activation volume for the reaction of Mo3S4 clusters with alkynes is determined and subsequently used to obtain a computational protocol featuring average errors of 0.7 kcal/mol and 2 cm3/mol, respectively.