set for the elements in the reactions considered here. SCF calculations and geometry optimizations by Gaussian 16 were conducted with tight convergence criteria. The DFT numerical integration grid had 175 radial shells per atom for H, 250 radial shells per atom for C and O, and 974 angular points per shell with pruning.