The thermodynamic integration/path integral Monte Carlo (TI/PIMC) method of calculating the temperature dependence of the equilibrium constant quantum mechanically is applied to O + HCl ⇌ OH + Cl reaction. The method is based upon PIMC simulations for energies of the reactants and the products and subsequently on thermodynamic integration for the ratios of partition functions. PIMC calculations are performed with the primitive approximation (PA) and the Takahashi–Imada approximation (TIA).