Bound states of a more general exponential screened Coulomb potential
  An alternative approximation scheme has been used in solving the Schrodinger equation to the more general case of exponential screened Coulomb potential, V(r)=-(a/r)\[1+(1+br)e^{-2br}]. The bound state energies of the 1s, $2s, and 3s-states, together with the ground state wave function are obtained analytically upto the second perturbation term. 
Sameer M. Ikhdair and Ramazan Sever