We indicated in our previous work that for QED the role of the scalar potential which appears at the loop level is much smaller than that of the vector potential and is in fact negligible. But the situation is different for QCD, one reason is that the loop effects are more significant because as is much larger than a, and second the non-perturbative QCD effects may induce a sizable scalar potential. In this work, we study phenomenologically the contribution of the scalar potential to the spectra of charmonia, bottomonia and bC（bc） families. Taking into account both vector and scalar potentials, by fitting the well measured charmonia and bottomonia spectra, we re-fix the relevant parameters and test them by calculating other states of not only the eharmonia and bottomonia families, but also the bc family. We also consider the Lamb shift of the spectra.
