Organic mixed ionic–electronic conductors (OMIECs) are pivotal for bioelectronics; however, the microscopic origins of doping-dependent charge transport remain elusive. In this study, we established a multi-scale computational framework to quantify the distinct intramolecular electronic coupling (Hab) distributions in systems with 25% and 75% doping levels. Our protocol employs automated quantum chemical calculations to exhaustively identify intrinsic local minima, ensuring thermodynamically stable initial conformations. Subsequent Molecular Dynamics (MD) simulations characterize the equilibration timescales and counter-ion dispersion behaviors. The simulation results reveal that the 75% doped system exhibits significantly stronger counter-ion confinement and a distinct electrostatic landscape compared to the 25% system. Finally, hybrid QM/MM calculations integrated with Constrained Density Functional Theory (CDFT) were utilized to evaluate Hab within these specific environments. The computed coupling distributions show a clear correlation with local electrostatic fluctuations induced by differing counter-ion arrangements. These findings indicate that doping-induced environmental disorder is a critical factor modulating intramolecular transport efficiency, providing a theoretical basis for optimizing OMIEC performance through electrostatic engineering.
Gao et al. (Wed,) studied this question.