Ab initio and density functional theory reinvestigation of gas-phase sulfuric acid monohydrate and ammonium hydrogen sulfate

Theo Kurtén, Markku R Sundberg, Hanna Vehkamäki, Madis Noppel, Johanna Blomqvist, Markku Kulmala
Journal of Physical Chemistry. A 2006 June 8, 110 (22): 7178-88
We have calculated the thermochemical parameters for the reactions H(2)SO(4) + H(2)O <--> H(2)SO(4).H(2)O and H(2)SO(4) + NH(3) <--> H(2)SO(4).NH(3) using the B3LYP and PW91 functionals, MP2 perturbation theory and four different basis sets. Different methods and basis sets yield very different results with respect to, for example, the reaction free energies. A large part, but not all, of these differences are caused by basis set superposition error (BSSE), which is on the order of 1-3 kcal mol(-1) for most method/basis set combinations used in previous studies. Complete basis set extrapolation (CBS) calculations using the cc-pV(X+d)Z and aug-cc-pV(X+d)Z basis sets (with X = D, T, Q) at the B3LYP level indicate that if BSSE errors of less than 0.2 kcal mol(-1) are desired in uncorrected calculations, basis sets of at least aug-cc-pV(T+d)Z quality should be used. The use of additional augmented basis functions is also shown to be important, as the BSSE error is significant for the nonaugmented basis sets even at the quadruple-zeta level. The effect of anharmonic corrections to the zero-point energies and thermal contributions to the free energy are shown to be around 0.4 kcal mol(-1) for the H(2)SO(4).H(2)O cluster at 298 K. Single-point CCSD(T) calculations for the H(2)SO(4).H(2)O cluster also indicate that B3LYP and MP2 calculations reproduce the CCSD(T) energies well, whereas the PW91 results are significantly overbinding. However, basis-set limit extrapolations at the CCSD(T) level indicate that the B3LYP binding energies are too low by ca. 1-2 kcal/mol. This probably explains the difference of about 2 kcal mol(-1) for the free energy of the H(2)SO(4) + H(2)O <--> H(2)SO(4).H(2)O reaction between the counterpoise-corrected B3LYP calculations with large basis sets and the diffusion-based experimental values of S. M. Ball, D. R. Hanson, F. L Eisele and P. H. McMurry (J. Phys. Chem. A. 2000, 104, 1715). Topological analysis of the electronic charge density based on the quantum theory of atoms in molecules (QTAIM) shows that different method/basis set combinations lead to qualitatively different bonding patterns for the H(2)SO(4).NH(3) cluster. Using QTAIM analysis, we have also defined a proton transfer degree parameter which may be useful in further studies.

Full Text Links

Find Full Text Links for this Article


You are not logged in. Sign Up or Log In to join the discussion.

Related Papers

Remove bar
Read by QxMD icon Read

Save your favorite articles in one place with a free QxMD account.


Search Tips

Use Boolean operators: AND/OR

diabetic AND foot
diabetes OR diabetic

Exclude a word using the 'minus' sign

Virchow -triad

Use Parentheses

water AND (cup OR glass)

Add an asterisk (*) at end of a word to include word stems

Neuro* will search for Neurology, Neuroscientist, Neurological, and so on

Use quotes to search for an exact phrase

"primary prevention of cancer"
(heart or cardiac or cardio*) AND arrest -"American Heart Association"