Quantum chemical implicit solvent models are used widely to estimate aqueous redox potentials. We compared the accuracy of several popular implicit solvent models (SM8, SMD, C-PCM, IEF-PCM, and COSMO-RS) for the prediction of aqueous single electron oxidation potentials of a diverse test set of neutral organic compounds for which accurate experimental oxidation potential and gas-phase ionization energy data are available. Using a thermodynamic cycle, we decomposed the free energy of oxidation into contributions arising from the gas-phase adiabatic ionization energy, the solvation free energy of the closed-shell neutral species, and the solvation free energy of the radical cation species. For aqueous oxidation potentials, implicit solvent models exhibited mean unsigned errors (MUEs) ranging from 0.27 to 0.50 V, depending on the model. The principal source of error was attributed to the computed solvation free energy of the oxidized radical cation. Based on these results, a recommended implicit solvation approach is the SMD model for the solvation free energy combined with CBS-QB3 for the gas-phase ionization energy. With this approach, the MUE in computed oxidation potentials was 0.27 V, and the MUE in solvation free energy of the charged open-shell species was 0.32 eV. This baseline assessment provides a compiled benchmark test set of vetted experimental data that may be used to judge newly developed solvation models for their ability to produce improved predictions for aqueous oxidation potentials and related properties.

Anilines and phenols are structurally similar compound classes that both are susceptible to oxidation by excited state triplet sensitizers but undergo oxidation by different mechanisms. To gain an understanding of the factors that control the rate of oxidation of anilines and phenols by triplet excited states, a kinetic study was performed on the oxidation of substituted anilines and phenols by methylene blue. The rate constants of one-electron transfer from anilines to triplet state methylene blue and their dependence on the reaction free energy are well fit to a Sandros–Boltzmann model. The observed rate constants are also well modeled when aniline oxidation potentials derived computationally are used. For phenols, the proton-coupled electron transfer rate constants were found to correlate primarily with O–H bond dissociation free energy and secondarily with phenol pK_a. Rate constants for phenols could be modeled using computed bond dissociation free energies. These results provide a basis for predicting aniline and phenol oxidation rates, which could be valuable, for example, in assessing the likely persistence and fate of aniline- and phenol-based aqueous environmental pollutants.

To study the influence of aqueous solvent on the electronic energy levels of dissolved organic molecules, we conducted liquid microjet photoelectron spectroscopy (PES) measurements of the aqueous vertical ionization energies (VIE_aq) of aniline (7.49 eV), veratrole alcohol (7.68 eV), and imidazole (8.51 eV). We also reanalyzed previously reported experimental PES data for phenol, phenolate, thymidine, and protonated imidazolium cation. We then simulated PE spectra by means of QM/MM molecular dynamics and EOM-IP-CCSD calculations with effective fragment potentials, used to describe the aqueous vertical ionization energies for six molecules, including aniline, phenol, veratrole alcohol, imidazole, methoxybenzene, and dimethylsulfide. Experimental and computational data enable us to decompose the VIE_aq into elementary processes. For neutral compounds, the shift in VIE upon solvation, ΔVIE_aq, was found to range from ≈−0.5 to −0.91 eV. The ΔVIE_aq was further explained in terms of the influence of deforming the gas phase solute into its solution phase conformation, the influence of solute hydrogen-bond donor and acceptor interactions with proximate solvent molecules, and the polarization of about 3000 outerlying solvent molecules. Among the neutral compounds, variability in ΔVIE_aq appeared largely controlled by differences in solute–solvent hydrogen-bonding interactions. Detailed computational analysis of the flexible molecule veratrole alcohol reveals that the VIE is strongly dependent on molecular conformation in both gas and aqueous phases. Finally, aqueous reorganization energies of the oxidation half-cell ionization reaction were determined from experimental data or estimated from simulation for the six compounds aniline, phenol, phenolate, veratrole alcohol, dimethylsulfide, and methoxybenzene, revealing a surprising constancy of 2.06 to 2.35 eV.