The pK(a) values of bipyridines in aqueous solution were calculated using a density functional theoretical method in combination with the Poisson-Boltzmann continuum solvation model. Calculated pK(a) values correlates well with experimental results.