Theoretical investigation of the proton affinities of benzazoles in the gas phase and in solution

The ground‐state equilibrium geometries of benzothiazole, benzoxazole, and benzimidazole were optimized at the density functional theory (DFT)/6‐31G** level of theory. Proton affinities on each of the possible sites in the studied series of compounds have been calculated at the DFT/6‐31G**/6‐311++G** level. The results indicate clearly that N‐site protonation is strongly favored over X‐site protonation (X = NH, O, S) for the series studied. Correlation of the computed proton affinities to the energy (E_HOMO) of the highest occupied MO in the gas phase and in solution has been explored and discussed. A comprehensive investigation of the effect of solvent on the process of protonation of the studied compounds has been performed. Different dielectric continuum models (i.e., Onsager, PCM, and IPCM) have been tested; their performance and range of applicability are reported and discussed. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005