Application of the PKS vibronic coupling model to mixed-valence ferrocene derivatives biferrocene and bis(fulvalene)diiron