Machine learning-accelerated screening of hydroquinone analogs for proton-coupled electron transfer

Proton-coupled electron transfer (PCET) mediated by hydroquinone and related molecules is key to natural and artificial energy conversion. The reactivity of these molecules depends on their bond dissociation free energy (BDFE), but studying the relationship between structure and thermochemistry across this chemical space has been limited by challenging experimental setup and high computational expense. Here, we present the first use of the AIMNet2 neural network potential to calculate average BDFE (BDFE_avg) values for the 2H⁺/2e^- dehydrogenation of about 200 000 hydroquinone-like compounds, including vicinal diamines, diols, and dithiols. Benchmarking against DFT calculations for 168 substituted ortho-phenylenediamines (opda) shows good agreement (R ² ∼ 0.84). Our analysis finds that the BDFE_avg of diamines ranges from 50 to 80 kcal mol^-1 and can be systematically tuned by modifying the backbone and N-substitution: electron-withdrawing groups raise BDFE_avg by up to 15 kcal mol^-1, while lower aromaticity in furan and thiophene backbones decreases BDFE_avg by approximately 10 kcal mol^-1 compared to the phenyl systems (∼65 kcal mol^-1). Validation through cyclic voltammetry and reactivity studies with quinone oxidants for selected compounds supports the computational results. This extensive thermochemical database and a web-based prediction tool developed as a result of this work will offer valuable resources for designing PCET reagents for catalysis, energy storage, and biomedical uses.