In this work, we explore the use of constrained density functional theory for the calculation of the charge transfer parameters of 2,1,3-benzothiadiazole (BTZ), a promising redoxmer, in acetonitrile (MeCN) solvent. The BTZ molecule has been studied as an anolyte in redox flow batteries, where charge transfer is a crucial process. It is highly desirable to simulate ab initio charge-transfer parameters, given their accuracy in predicting electron-transfer rates and structure-activity relationships. This work explores the state of the art in charge-transfer simulation for this process. Constrained density functional theory (DFT) calculations are used to predict charge-transfer free energies and electronic couplings, which are crucial for evaluating charge transfer within the Marcus theory. Based on the simulations, we find that electronic coupling fluctuates rapidly with time and also depends on the difference between the donor and acceptor state energies (reaction gap energy). Based on our evaluation of Marcus theory, BTZ has a predicted self-exchange reaction rate constant on the order of 0.5 M-1 s-1 at 1 M concentration in MeCN. Our work demonstrates the utility of constrained DFT for providing physical insight into a charge-transfer process, while also highlighting current limitations in computational and algorithmic capacity in achieving desirable system sizes and levels of ergodicity in molecular dynamics simulations. A significant conclusion of this work is that time-dependent sampling of electronic coupling as a function of the reaction gap energy, as described herein, is essential for future predictions of charge and electron transfer.
Howard et al. (Wed,) studied this question.