The computational modeling of reactions is simple in theory but can be quite tricky in practice. This article aims at the purpose of providing an assistance to a proper way of describing reactions theoretically and provides rough guidelines to the computational methods involved. Reactions in liquid phase chemical equilibrium can be described theoretically in terms of the Gibbs free energy of reaction. This property can be divided into a sum of three disjunct terms, namely the gas phase reaction energy, the finite temperature contribution to the Gibbs free energy, and the Gibbs free energy of solvation. The three contributions to the Gibbs free energy of reaction can be computed separately, using different theoretico-chemical calculation methods. While some of these terms can be obtained reliably by computationally cheap methods, for others a high level of theory is required to obtain predictions of quantitative quality. To propose workflows which can strike the balance between accuracy and computational cost, a number of benchmarks assessing the precision of different levels of theory is given. As an illustrative example, the low-temperature hydrogenation reaction of acetaldehyde to ethanol in solvent toluene is shown. (C) 2017 American Institute of Chemical Engineers