To parameterize the electronic Hamiltonian of Eq. ( ), we adopt the FPHD method to obtain the corresponding Hamiltonian parameters including the transition dipole moments, states energies and couplings between the various different states in the quasi-diabatic representation. The central idea of FPHD is to search for orthonormal electronic states that maximally localize electron and hole densities in terms of predefined molecular fragments.