A general strategy for the accurate computation of structural
and
spectroscopic properties of biomolecule building blocks in the gas
phase has been further improved and validated with a special reference
to tautomeric equilibria. The main improvements concern the use of
the cc-pVTZ-F12 basis set in both DFT and CCSD(T)-F12 computations,
the inclusion of core-valence correlation in geometry optimizations
by double hybrid functionals, and the use of the cc-pVQZ-F12 basis
set for complete basis set extrapolation at the MP2-F12 level. The
resulting model chemistry is applied to the challenging problem of
cytosine tautomers in the gas phase. The results are in remarkable
agreement with experiment concerning both rotational and vibrational
spectroscopic parameters and permit their unbiased interpretation
in terms of structural and thermochemical features. Together with
the intrinsic interest of the studied molecule, the accuracy of the
results obtained at reasonable cost without any empirical parameter
suggests that the proposed composite method can be profitably employed
for accurate investigations of other molecular bricks of life.