A Systematic DFT Study of Two Elementary Steps Involving Hydrogen Peroxide and the Hydroxyl Radical in the Fenton Reaction

2018 ◽  
Author(s):  
Danilo Carmona ◽  
David Contreras ◽  
Oscar A. Douglas-Gallardo ◽  
Stefan Vogt-Geisse ◽  
Pablo Jaque ◽  
...  

The Fenton reaction plays a central role in many chemical and biological processes and has various applications as e.g. water remediation. The reaction consists of the iron-catalyzed homolytic cleavage of the oxygen-oxygen bond in the hydrogen peroxide molecule and the reduction of the hydroxyl radical. Here, we study these two elementary steps with high-level ab-initio calculations at the complete basis set limit and address the performance of different DFT methods following a specific classification based on the Jacob´s ladder in combination with various Pople's basis sets. Ab-initio calculations at the complete basis set limit are in agreement to experimental reference data and identified a significant contribution of the electron correlation energy to the bond dissociation energy (BDE) of the oxygen-oxygen bond in hydrogen peroxide and the electron affinity (EA) of the hydroxyl radical. The studied DFT methods were able to reproduce the ab-initio reference values, although no functional was particularly better for both reactions. The inclusion of HF exchange in the DFT functionals lead in most cases to larger deviations, which might be related to the poor description of the two reactions by the HF method. Considering the computational cost, DFT methods provide better BDE and EA values than HF and post--HF methods with an almost MP2 or CCSD level of accuracy. However, no systematic general prediction of the error based on the employed functional could be established and no systematic improvement with increasing the size in the Pople's basis set was found, although for BDE values certain systematic basis set dependence was observed. Moreover, the quality of the hydrogen peroxide, hydroxyl radical and hydroxyl anion structures obtained from these functionals was compared to experimental reference data. In general, bond lengths were well reproduced and the error in the angles were between one and two degrees with some systematic trend with the basis sets. From our results we conclude that DFT methods present a computationally less expensive alternative to describe the two elementary steps of the Fenton reaction. However, choice of approximated functionals and basis sets must be carefully done and the provided benchmark allows a systematic validation of the electronic structure method to be employed

2018 ◽  
Author(s):  
Danilo Carmona ◽  
David Contreras ◽  
Oscar A. Douglas-Gallardo ◽  
Stefan Vogt-Geisse ◽  
Pablo Jaque ◽  
...  

The Fenton reaction plays a central role in many chemical and biological processes and has various applications as e.g. water remediation. The reaction consists of the iron-catalyzed homolytic cleavage of the oxygen-oxygen bond in the hydrogen peroxide molecule and the reduction of the hydroxyl radical. Here, we study these two elementary steps with high-level ab-initio calculations at the complete basis set limit and address the performance of different DFT methods following a specific classification based on the Jacob´s ladder in combination with various Pople's basis sets. Ab-initio calculations at the complete basis set limit are in agreement to experimental reference data and identified a significant contribution of the electron correlation energy to the bond dissociation energy (BDE) of the oxygen-oxygen bond in hydrogen peroxide and the electron affinity (EA) of the hydroxyl radical. The studied DFT methods were able to reproduce the ab-initio reference values, although no functional was particularly better for both reactions. The inclusion of HF exchange in the DFT functionals lead in most cases to larger deviations, which might be related to the poor description of the two reactions by the HF method. Considering the computational cost, DFT methods provide better BDE and EA values than HF and post--HF methods with an almost MP2 or CCSD level of accuracy. However, no systematic general prediction of the error based on the employed functional could be established and no systematic improvement with increasing the size in the Pople's basis set was found, although for BDE values certain systematic basis set dependence was observed. Moreover, the quality of the hydrogen peroxide, hydroxyl radical and hydroxyl anion structures obtained from these functionals was compared to experimental reference data. In general, bond lengths were well reproduced and the error in the angles were between one and two degrees with some systematic trend with the basis sets. From our results we conclude that DFT methods present a computationally less expensive alternative to describe the two elementary steps of the Fenton reaction. However, choice of approximated functionals and basis sets must be carefully done and the provided benchmark allows a systematic validation of the electronic structure method to be employed


1999 ◽  
Vol 23 (8) ◽  
pp. 502-503
Author(s):  
Branko S. Jursic

High level ab initio and density functional theory studies are performed on highly protonated methane species.


2010 ◽  
Vol 2010 ◽  
pp. 1-6 ◽  
Author(s):  
Prasad Yedlapalli ◽  
Sangyong Lee ◽  
Jae W. Lee

Structure II clathrate hydrates of pure hydrogen and binary hydrates of are studied using ab initio calculations to determine the stable occupancies of small cavities. Ab initio calculations are carried out for a double cavity consisting of one dodecahedron (small cavity) and one hexakaidecahedron (large cavity). These two cavities are attached to each other as in sII hydrates to form a double cavity. One or two molecules are placed in the small cavity and one THF (or 4 molecules) molecule is placed in the large cavity. We have determined the binding energies of the double cavities at the MP2 level using various basis sets (3-21G, 3-21G(2p), 3-21 G(2p), 6-31G, 6-31G(2p), and 6-31 G(2p)). Different basis sets yield different stable occupancies of the small cavity. The results from the highest basis set (6-31 G(2p) with zero point energy corrections) indicate that the single occupancy is slightly more favorable than the double occupancy in both the cases of pure hydrates and THF + double hydrates.


Inorganics ◽  
2019 ◽  
Vol 7 (3) ◽  
pp. 35 ◽  
Author(s):  
Ibon Alkorta ◽  
Anthony Legon

Geometries, equilibrium dissociation energies (De), intermolecular stretching, and quadratic force constants (kσ) determined by ab initio calculations conducted at the CCSD(T)/aug-cc-pVTZ level of theory, with De obtained by using the complete basis set (CBS) extrapolation [CCSD(T)/CBS energy], are presented for the B···BeR2 and B···MgR2 complexes, where B is one of the following Lewis bases: CO, H2S, PH3, HCN, H2O or NH3, and R is H, F or CH3. The BeR2 and MgR2 precursor molecules were shown to be linear and non-dipolar. The non-covalent intermolecular bond in the B···BeR2 complexes is shown to result from the interaction of the electrophilic band around the Be atom of BeR2 (as indicated by the molecular electrostatic potential surface) with non-bonding electron pairs of the base, B, and may be described as a beryllium bond by analogy with complexes such as B···CO2, which contain a tetrel bond. The conclusions for the B···MgR2 series are similar and a magnesium bond can be correspondingly invoked. The geometries established for B···BeR2 and B···MgR2 can be rationalized by a simple rule previously enunciated for tetrel-bonded complexes of the type B···CO2. It is also shown that the dissociation energy, De, is directly proportional to the force constant, kσ, in each B···MR2 series, but with a constant of proportionality different from that established for many hydrogen-bonded B···HX complexes and halogen-bonded B···XY complexes. The values of the electrophilicity, EA, determined from the De for B···BeR2 complexes for the individual Lewis acids, A, reveal the order A = BeF2 > BeH2 > Be(CH3)2—a result that is consistent with the −I and +I effects of F and CH3 relative to H. The conclusions for the MgR2 series are similar but, for a given R, they have smaller electrophilicities than those of the BeR2 series. A definition of alkaline-earth non-covalent bonds is presented.


1977 ◽  
Vol 55 (5) ◽  
pp. 863-868 ◽  
Author(s):  
N. Colin Baird ◽  
Harish B. Kathpal

The important geometrical variables in the structures of the lowest 2A′ and 2A′′ states of the free radicals HCO, CH3CO, NH2CO, HNN, and CH3NN have been determined by ab initio MO calculations using the STO-3G basis set. The energy differences between the states, and the energies of the radicals relative to their decomposition products and relative to their hydrogen atom addition products, are reported using both STO-3G and 4-31G basis sets in the restricted open-shell calculations. The trends in these results and their relation to available experimental data are discussed.


1993 ◽  
Vol 07 (01n03) ◽  
pp. 262-265 ◽  
Author(s):  
M. METHFESSEL ◽  
M. VAN SCHILFGAARDE

A major advance in electronic structure calculations was the combination of local-density techniques with molecular dynamics by Car and Parrinello seven years ago. Unfortunately, application of the Car-Parrinello scheme has been limited essentially to sp materials because only in the plane-wave pseudopotential method forces are trivial to calculate. We present a systematic approach to derive force theorems with desired characteristics within complicated basis sets, which are applicable to all elements of the periodic table equally well. Application to the LMTO basis set yields an accurate force theorem, quite distinct from the Hellman-Feynman form, which is exceptionally insensitive to errors in the trial density. The forces were implemented in a new full-potential LMTO method which is suited to arbitrary geometries. First results for ab-initio molecular dynamics and simulated annealing runs are shown for some random small molecules and small clusters of silver atoms.


2011 ◽  
Vol 76 (4) ◽  
pp. 327-341 ◽  
Author(s):  
Vladimír Špirko ◽  
Xiangzhu Li ◽  
Josef Paldus

Recently generated ground state potential energy curves (PECs) for the nitrogen molecule, as obtained with the reduced multireference (RMR) coupled-cluster (CC) method with singles and doubles (RMR-CCSD), and its version corrected for the secondary triples RMR-CCSD(T), using cc-pVXZ basis sets with X = D, T, and Q, as well as the extrapolated complete basis set (cbs) limit (X. Li and J. Paldus: J. Chem. Phys. 2008, 129, 054104), are compared with both the highly accurate theoretical configuration interaction PEC of Gdanitz (Chem. Phys. Lett. 1998, 283, 253) and analytic PECs obtained by fitting an extensive set of experimental data (R. J. Le Roy et al.: J. Chem. Phys. 2006, 125, 164310). These results are analyzed using a morphing procedure based on the reduced potential curve (RPC) method of Jenč. It is found that an RPC fit of both theoretical potentials can be achieved with only a few parameters. The RMR PECs are found to provide an excellent description of experimentally available vibrational levels, but significantly deviate from those of Gdanitz’s PEC for highly stretched geometries, yet still do provide a qualitatively correct PECs that lie within the region delimited by Le Roy’s analytical PECs.


2013 ◽  
Vol 91 (7) ◽  
pp. 559-572 ◽  
Author(s):  
Jennifer L. Kellie ◽  
Stacey D. Wetmore

When using a hybrid methodology to treat an enzymatic reaction, many factors contribute to selecting the method for the high-level region, which can be complicated by the presence of dispersion-driven interactions such as π–π stacking. In addition, the proper treatment of the reaction center often requires a large number of heavy atoms to be included in the high-level region, precluding the use of ab initio methods such as MP2 as well as large basis sets, in the optimization step. In the present work, popular DFT methods were tested to identify an appropriate functional for treating the high-level region in ONIOM optimizations of reactions catalyzed by nonmetalloenzymes. Eight different DFT methods (B3LYP, B97-2, MPW1K, MPWB1K, BB1K, B1B95, M06-2X, and ωB97X-D) in combination with four double-ζ quality Pople basis sets were tested for their ability to optimize noncovalent interactions (hydrogen bonding and π–π) and characterize reactions (proton transfer, SN2 hydrolysis, and unimolecular cleavage). Although the primary focus of this study is accurate structure determination, energetics were also examined at both the optimization level of theory, and with triple-ζ quality basis set and select (M06-2X or ωB97X-D) methods. If dispersion-driven interactions exist within the active site, then MPWB1K/6-31G(d,p) or M06-2X/6-31+G(d,p) are recommended for the optimization step with subsequent triple-ζ quality single-point energies. However, since dispersion-corrected functionals (M06-2X and ωB97X-D) generally require diffuse functions to yield appropriate geometries, the possible size of the high-level region is greatly limited with these methods. In contrast, if the model is large enough to recover steric constraints on π–π interactions, then B3LYP with a small basis set performs comparatively well for the optimization step and is significantly less computationally expensive. Interestingly, the functionals that afford the best geometries often do not yield the best energetics, which emphasizes the importance of structural benchmark studies.


Sign in / Sign up

Export Citation Format

Share Document