On the computation of pKa values of organic and organometallic molecules in different solvents