Theoretical and Computational Chemistry

A Computational Protocol Combining DFT and Cheminformatics for Prediction of pH-Dependent Redox Potentials

Piotr de Silva Technical University of Denmark

Abstract

We present and validate a robust procedure to calculate the redox potentials of organic molecules at any pH value, using widely available quantum chemistry and cheminformatics methods. Using a consistent experimental data set for validation, we explore and compare a few different methods for calculating reaction free energies, the treatment of solvation, and the effect of pH on redox potentials. We find that the B3LYP hybrid functional with COSMO solvation method, in conjunction with thermal contributions evaluated from BLYP gas-phase harmonic frequencies, yields a good prediction of pH=0 redox potentials at a moderate computational cost. To predict how the potentials are affected by pH, we propose an improved version of the Alberty-Legendre transform that allows the construction of a more realistic Pourbaix diagram by taking into account how the protonation state changes with pH.

Content

Thumbnail image of A computational protocol combining DFT and cheminformatics for prediction of pH-dependent redox potentials.pdf
download asset A computational protocol combining DFT and cheminformatics for prediction of pH-dependent redox potentials.pdf 0.82 MB [opens in a new tab]

Supplementary material

Thumbnail image of SI A computational protocol combining DFT and cheminformatics for prediction of pH-dependent redox potentials.pdf
download asset SI A computational protocol combining DFT and cheminformatics for prediction of pH-dependent redox potentials.pdf 0.55 MB [opens in a new tab]
SI A computational protocol combining DFT and cheminformatics for prediction of pH-dependent redox potentials