Numerical method for N electrons bound to a polar quantum dot with a Coulomb impurity

A numerical method is proposed to calculate the Frohlich Hamiltonian containing N electrons bound to polar quantum dot with a Coulomb impurity without transformation to the coordination frame of the center of mass and by direct diagonalization. As an example to demonstrate the formalism of this method, the low-lying spectra of three interacting electrons bound to an on-center Coulomb impurity, both for accepter and donor, are calculated and analyzed in a polar quantum dot under a perpendicular magnetic field. Taking polaron effect into account, the physical meaning of the phonon-induced terms, both self-square terms and cross terms of the Hamiltonian are discussed. The calculation can also be applied to systems containing particles with opposite charges, such as excitons.