We propose a genetic algorithm (GA)-driven density matrix method for calculating the equilibrium geometry and charge storage configurations of hole(bipolaron)-doped Polythiophene (PT) oligomers. A modified version of the Su-Schrieffer-Heeger Hamiltonian is used to describe the PT chain. A population of geometry strings are used to generate the corresponding PT-Hamiltonians which act as generators of corresponding unitary transformations which transform a single trial one-electron density matrix into a population of density matrices for different geometries. As the geometry strings evolve under the action of GA operators, the density matrices also evolve on the fitness landscape. Once the fitness reaches maximum, the optimum geometry string, the corresponding Hamiltonian, and the density matrix lead to equilibrium geometry, energy, charge distribution, band gap, Fermi energy, etc. The bipolaronic defect-induced conducting state is predicted to have a symmetric lattice-like admixture of aromatic and quininoid regions and is characterized by nonzero density of states at Fermi energy and low band gaps. Copyright © Taylor & Francis Group, LLC.
|Journal||Materials and Manufacturing Processes|
|Publisher||TAYLOR & FRANCIS INC|