Thermodynamics of chemical reactions with COSMO-RS: the extreme case of charge separation or recombination.

Many technically relevant chemical processes in the condensed phase involve as elementary reactive steps the formation of ions from neutral species or, as the opposite, recombination of ions. Such reactions that generate or annihilate charge defy the standard gas phase quantum chemical treatment, and also continuum solvation models are only partially able to account for the right amount of stabilization in solution. In this work, for such types of reaction, a solvation treatment involving the COSMO-RS method is assessed, which leads to improved results, i.e., errors of only around 10 kJ/mol for both protic and aprotic solvents. The examples discussed here comprise protolysis reactions and organo halide heterolysis, for both of which a comparison with reliable experimental data is possible. It is observed that for protolysis, the quality of results does not strongly depend on the quantum chemical method used for energy calculation. In contrast, in the case of heterolytic carbon-chlorine bond cleavage, clearly better results are obtained for higher correlated (coupled cluster) methods or the density functional M06-2X, which is well known for its accuracy if applied to organic chemistry. This hints at least that the right answer is obtained for the right reason and not due to a compensation of errors from gas phase thermodynamics with those from the solvation treatment. Problems encountered with certain critical solvents or upon decomposing Gibbs free energies into heats or entropies of reaction are found to relate mostly to the parameterization of the H-bonding term within COSMO-RS.