A novel quantum-classical recipe for locating the global minimum on the potential energy surface of a large molecule and simultaneously predicting the associated electronic charge distribution is developed by interfacing the classical particle swarm optimization with a near optimal unitary evolution scheme for the trial one electron density matrix. The unitary transformation is generated by an antihermitian matrix linked to the molecular electronic Hamiltonian at the instantaneous nuclear configurations discovered by the swarm as it flies. The algorithm is used to predict the extensive reorganization of electronic charge distribution and bond lengths in polythiophene oligomers on doping at various levels.antihermitian matrix, neutral and doped polythiophene, one electron density matrix, particle swarm optimization, Su-Schrieffer-Heeger p-electron Hamiltonian † Present address: