“…Test calculations estimate the reliability of used methods and the uncertainty related to those energies. The geometrical structures of all reactants, prereactive complexes, transition states, and products were optimized by utilizing the hybrid density functional B3LYP method [39] with cc-pv(T+d)z [40] basis set and MP2 [41] approach with aug-cc-pv(T+d)z basis set, because their reliability and efficiency were confirmed by previous studies on reactions [13,17,18,30,42] and clusters [43][44][45][46][47][48][49] involving sulfuric acid. Moreover, it is also proven [50,51] that the cc-pv(T+d)z and augcc-pv(T+d)z basis sets can produce more accurate results for sulfur-containing compounds.…”