Exploration of density functional methods for one-electron reduction potential of nitrobenzenes

Prediction of the standard potentials for one-electron oxidation of N,N,N',N' tetrasubstituted p-phenylenediamines by calculation

The formal potentials for the reversible one-electron oxidation of N,N,N',N' tetrasubstituted p-phenylenediamines in acetonitrile have been applied as a test set for benchmarking computational methods for a series of compounds with only small structural differences. The aim of the study is to propose a simple method for calculating the standard oxidation potentials, and therefore, the protocol is progressively developed by adding more terms in the energy expression. In addition, the effect of including implicit solvation models (IEFPCM, CPCM, and SMD), larger basis sets, and correlation methods are investigated. The oxidation potentials calculated using the G3MP2B3 approach with IEFPCM resulted in the best fit (R² = 0.9624), but the slope of the correlation line, 0.74, is far from the optimal value, 1.00. B3LYP/6-311++G(d,p) and TPSSh/6-311++G(2d,p) yielded only slightly less consistent data (R² = 0.9388 and R² = 0.9425), but with much better slopes, 1.00 and 0.94, respectively. We conclude that it is important to investigate the basis set size and treatment of electron correlation when calculating oxidation potentials for N,N,N',N' tetrasubstituted p-phenylenediamines.

Directed Design of a AuI Complex with a Reduced Mesoionic Carbene Radical Ligand: Insights from 1,2,3-Triazolylidene Selenium Adducts and Extensive Electrochemical Investigations

Carbene-based radicals are important for both fundamental and applied chemical research. Herein, extensive electrochemical investigations of nine different 1,2,3-triazolylidene selenium adducts are reported. It is found that the half-wave potentials of the first reduction of the selones correlate with their calculated LUMO levels and the LUMO levels of the corresponding triazolylidene-based mesoionic carbenes (MICs). Furthermore, unexpected quasi-reversibility of the reduction of two triazoline selones, exhibiting comparable reduction potentials, was discovered. Through UV/Vis/NIR and EPR spectroelectrochemical investigations supported by DFT calculations, the radical anion was unambiguously assigned to be triazoline centered. This electrochemical behavior was transferred to a triazolylidene-type MIC-gold phenyl complex resulting in a MIC-radical coordinated AuI species. Apart from UV-Vis-NIR and EPR spectroelectrochemical investigations of the reduction, the reduced gold-coordinated MIC radical complex was also formed in situ in the bulk through chemical reduction. This is the first report of a monodentate triazolylidene-based MIC ligand that can be reduced to its anion radical in a metal complex. The results presented here provide design principles for stabilizing radicals based on MICs.

Experimental Validation of Hydrogen Atom Transfer Gibbs Free Energy as a Predictor of Nitroaromatic Reduction Rate Constants

Nitroaromatic compounds (NACs) are a class of prevalent contaminants. Abiotic reduction is an important fate process that initiates NAC degradation in the environment. Many linear free energy relationship (LFER) models have been developed to predict NAC reduction rates. Almost all LFERs to date utilize experimental aqueous-phase one-electron reduction potential ( E_H¹) of NAC as a predictor, and thus, their utility is limited by the availability of E_H¹ data. A promising new approach that utilizes computed hydrogen atom transfer (HAT) Gibbs free energy instead of E_H¹ as a predictor was recently proposed. In this study, we evaluated the feasibility of HAT energy for predicting NAC reduction rate constants. Using dithionite-reduced quinones, we measured the second-order rate constants for the reduction of seven NACs by three hydroquinones of different protonation states. We computed the gas-phase energies for HAT and electron affinity (EA) of NACs and established HAT- and EA-based LFERs for six hydroquinone species. The results suggest that HAT energy is a reliable predictor of NAC reduction rate constants and is superior to EA. This is the first independent, experimental validation of HAT-based LFER, a new approach that enables rate prediction for a broad range of structurally diverse NACs based solely on molecular structures.

Novel 8-nitroquinolin-2(1H)-ones as NTR-bioactivated antikinetoplastid molecules: Synthesis, electrochemical and SAR study

To study the antiparasitic 8-nitroquinolin-2(1H)-one pharmacophore, a series of 31 derivatives was synthesized in 1-5 steps and evaluated in vitro against both Leishmania infantum and Trypanosoma brucei brucei. In parallel, the reduction potential of all molecules was measured by cyclic voltammetry. Structure-activity relationships first indicated that antileishmanial activity depends on an intramolecular hydrogen bond (described by X-ray diffraction) between the lactam function and the nitro group, which is responsible for an important shift of the redox potential (+0.3 V in comparison with 8-nitroquinoline). With the assistance of computational chemistry, a set of derivatives presenting a large range of redox potentials (from -1.1 to -0.45 V) was designed and provided a list of suitable molecules to be synthesized and tested. This approach highlighted that, in this series, only substrates with a redox potential above -0.6 V display activity toward L. infantum. Nevertheless, such relation between redox potentials and in vitro antiparasitic activities was not observed in T. b. brucei. Compound 22 is a new hit compound in the series, displaying both antileishmanial and antitrypanosomal activity along with a low cytotoxicity on the human HepG2 cell line. Compound 22 is selectively bioactivated by the type 1 nitroreductases (NTR1) of L. donovani and T. brucei brucei. Moreover, despite being mutagenic in the Ames test, as most of nitroaromatic derivatives, compound 22 was not genotoxic in the comet assay. Preliminary in vitro pharmacokinetic parameters were finally determined and pointed out a good in vitro microsomal stability (half-life > 40 min) and a 92% binding to human albumin.

Informing Efforts to Develop Nitroreductase for Amine Production

Nitroreductases (NRs) hold promise for converting nitroaromatics to aromatic amines. Nitroaromatic reduction rate increases with Hammett substituent constant for NRs from two different subgroups, confirming substrate identity as a key determinant of reactivity. Amine yields were low, but compounds yielding amines tend to have a large π system and electron withdrawing substituents. Therefore, we also assessed the prospects of varying the enzyme. Several different subgroups of NRs include members able to produce aromatic amines. Comparison of four NR subgroups shows that they provide contrasting substrate binding cavities with distinct constraints on substrate position relative to the flavin. The unique architecture of the NR dimer produces an enormous contact area which we propose provides the stabilization needed to offset the costs of insertion of the active sites between the monomers. Thus, we propose that the functional diversity included in the NR superfamily stems from the chemical versatility of the flavin cofactor in conjunction with a structure that permits tremendous active site variability. These complementary properties make NRs exceptionally promising enzymes for development for biocatalysis in prodrug activation and conversion of nitroaromatics to valuable aromatic amines. We provide a framework for identifying NRs and substrates with the greatest potential to advance.

Predicting reduction potentials of 1,3,6-triphenyl fulvenes using molecular electrostatic potential analysis of substituent effects

The influence of mono- and multiple substituent effect on the reduction potential (E⁰ ) of 1,3,6-triphenyl fulvenes is investigated using B3LYP-SMD/6-311+G(d,p) level density functional theory. The molecular electrostatic potential (MESP) minimum at the fulvene π-system (V_min ) and the change in MESP at any of the fulvene carbon atoms (ΔV_C ) for both neutral and reduced forms are used as excellent measures of substituent effect from the para and meta positions of the 1,3 and 6-phenyl moieties. Substitution at 6-phenyl para position has led to significant change in E⁰ than any other positions. By applying the additivity rule of substituent effects, an equation in ΔV_C is derived to predict E⁰ for multiply substituted fulvenes. Further, E⁰ is predicted for a set of 2000 hexa-substituted fulvene derivatives where the substituents and their positions in the system are chosen in a random way. The calculated E⁰ agreed very well with the experimental E⁰ reported by Godman et al. Predicting E⁰ solely by substituent effect offers a simple and powerful way to select suitable combinations of substituents on fulvene system for light harvesting applications. © 2018 Wiley Periodicals, Inc.

In search of the best DFT functional for dealing with organic anionic species

“And the winner is…” This work assesses the ability of different Density Functional Theory (DFT) functionals for a proper treatment of organic anionic species.

Explicit solvent simulations of the aqueous oxidation potential and reorganization energy for neutral molecules: gas phase, linear solvent response, and non-linear response contributions

First principles simulations were used to predict aqueous one-electron oxidation potentials (Eox) and associated half-cell reorganization energies (λaq) for aniline, phenol, methoxybenzene, imidazole, and dimethylsulfide. We employed quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulations of the oxidized and reduced species in an explicit aqueous solvent, followed by EOM-IP-CCSD computations with effective fragment potentials for diabatic energy gaps of solvated clusters, and finally thermodynamic integration of the non-linear solvent response contribution using classical MD. A priori predicted Eox and λaq values exhibit mean absolute errors of 0.17 V and 0.06 eV, respectively, compared to experiment. We also disaggregate Eox into several well-defined free energy properties, including the gas phase adiabatic free energy of ionization (7.73 to 8.82 eV), the solvent-induced shift in the free energy of ionization due to linear solvent response (-2.01 to -2.73 eV), and the contribution from non-linear solvent response (-0.07 to -0.14 eV). The linear solvent response component is further apportioned into contributions from the solvent-induced shift in vertical ionization energy of the reduced species (ΔVIEaq) and the solvent-induced shift in negative vertical electron affinity of the ionized species (ΔNVEAaq). The simulated ΔVIEaq and ΔNVEAaq are found to contribute the principal sources of uncertainty in computational estimates of Eox and λaq. Trends in the magnitudes of disaggregated solvation properties are found to correlate with trends in structural and electronic features of the solute. Finally, conflicting approaches for evaluating the aqueous reorganization energy are contrasted and discussed, and concluding recommendations are given.

Predicting reduction rates of energetic nitroaromatic compounds using calculated one-electron reduction potentials

The evaluation of new energetic nitroaromatic compounds (NACs) for use in green munitions formulations requires models that can predict their environmental fate. Previously invoked linear free energy relationships (LFER) relating the log of the rate constant for this reaction (log(k)) and one-electron reduction potentials for the NAC (E1NAC) normalized to 0.059 V have been re-evaluated and compared to a new analysis using a (nonlinear) free-energy relationship (FER) based on the Marcus theory of outer-sphere electron transfer. For most reductants, the results are inconsistent with simple rate limitation by an initial, outer-sphere electron transfer, suggesting that the linear correlation between log(k) and E1NAC is best regarded as an empirical model. This correlation was used to calibrate a new quantitative structure-activity relationship (QSAR) using previously reported values of log(k) for nonenergetic NAC reduction by Fe(II) porphyrin and newly reported values of E1NAC determined using density functional theory at the M06-2X/6-311++G(2d,2p) level with the COSMO solvation model. The QSAR was then validated for energetic NACs using newly measured kinetic data for 2,4,6-trinitrotoluene (TNT), 2,4-dinitrotoluene (2,4-DNT), and 2,4-dinitroanisole (DNAN). The data show close agreement with the QSAR, supporting its applicability to other energetic NACs.

Computational electrochemistry: prediction of liquid-phase reduction potentials

This article reviews recent developments and applications in the area of computational electrochemistry. Our focus is on predicting the reduction potentials of electron transfer and other electrochemical reactions and half-reactions in both aqueous and nonaqueous solutions. Topics covered include various computational protocols that combine quantum mechanical electronic structure methods (such as density functional theory) with implicit-solvent models, explicit-solvent protocols that employ Monte Carlo or molecular dynamics simulations (for example, Car-Parrinello molecular dynamics using the grand canonical ensemble formalism), and the Marcus theory of electronic charge transfer. We also review computational approaches based on empirical relationships between molecular and electronic structure and electron transfer reactivity. The scope of the implicit-solvent protocols is emphasized, and the present status of the theory and future directions are outlined.

Oxidation potentials of functionalized sulfone solvents for high-voltage Li-ion batteries: a computational study

New electrolytes with large electrochemical windows are needed to meet the challenge for high-voltage Li-ion batteries. Sulfone as an electrolyte solvent boasts of high oxidation potentials. Here we examine the effect of multiple functionalization on sulfone's oxidation potential. We compute oxidation potentials for a series of sulfone-based molecules functionalized with fluorine, cyano, ester, and carbonate groups by using a quantum chemistry method within a continuum solvation model. We find that multifunctionalization is a key to achieving high oxidation potentials. This can be realized through either a fluorether group on a sulfone molecule or sulfonyl fluoride with a cyano or ester group.

A method to calculate the one-electron reduction potentials for nitroaromatic compounds based on gas-phase quantum mechanics

Nitroaromatic compounds (NACs) are widespread environmental contaminants, and the one-electron reduction potential (E oH) is an important parameter used in modeling their environmental fate. We have identified a method that is both accurate and efficient to predict E oH values for NACs, using gas-phase quantum mechanics (QM) calculations combined with empirical correlations. First, the adiabatic electron affinity (EA) at 0 K is calculated using the B98/MG3S method, and the predictions are scaled by a factor of 0.802 to account for systematic errors in the density functional calculations. Second, the E oH values are predicted from a linear correlation between E oH and EA. Using this method, E oH values were predicted with a mean absolute deviation from measured values of 0.021 V for the 14 NACs used to obtain the correlation and 0.029 V for six additional NACs. This represents a substantial improvement in accuracy over predictions by other QM methods, which are affected by large errors in solvation or aqueous-phase calculations for some compounds.

One-electron standard reduction potentials of nitroaromatic and cyclic nitramine explosives

Extensive studies have been conducted in the past decades to predict the environmental abiotic and biotic redox fate of nitroaromatic and nitramine explosives. However, surprisingly little information is available on one-electron standard reduction potentials (Eo(R-NO2/R-NO2-)). The Eo(R-NO2/R-NO2-) is an essential thermodynamic parameter for predicting the rate and extent of reductive transformation for energetic residues. In this study, experimental (linear free energy relationships) and theoretical (ab initio calculation) approaches were employed to determine Eo(R-NO2/R-NO2-) for nitroaromatic, (caged) cyclic nitramine, and nitroimino explosives that are found in military installations or are emerging contaminants. The results indicate a close agreement between experimental and theoretical Eo(R-NO2/R-NO2-) and suggest a key trend: Eo(R-NO2/R-NO2-) value decreases from di- and tri-nitroaromatic (e.g., 2,4-dinitroanisole) to nitramine (e.g., RDX) to nitroimino compound (e.g., nitroguanidine). The observed trend in Eo(R-NO2/R-NO2-) agrees with reported rate trends for reductive degradation, suggesting a thermodynamic control on the reduction rate under anoxic/suboxic conditions.