Different Density Functional Theory Methods Research Articles

Scalar Relativistic All-Electron and Pseudopotential Ab Initio Study of a Minimal Nitrogenase [Fe(SH)4H]- Model Employing Coupled-Cluster and Auxiliary-Field Quantum Monte Carlo Many-Body Methods.

Nitrogenase is the only enzyme that can cleave the triple bond in N2, making nitrogen available to organisms. The detailed mechanism of this enzyme is currently not known, and computational studies are complicated by the fact that different density functional theory (DFT) methods give very different energetic results for calculations involving nitrogenase models. Recently, we designed a [Fe(SH)4H]- model with the fifth proton binding either to Fe or S to mimic different possible protonation states of the nitrogenase active site. We showed that the energy difference between these two isomers (ΔE) is hard to estimate with quantum-mechanical methods. Based on nonrelativistic single-reference coupled-cluster (CC) calculations, we estimated that the ΔE is 101 kJ/mol. In this study, we demonstrate that scalar relativistic effects play an important role and significantly affect ΔE. Our best revised single-reference CC estimates for ΔE are 85-91 kJ/mol, including energy corrections to account for contributions beyond triples, core-valence correlation, and basis-set incompleteness error. Among coupled-cluster approaches with approximate triples, the canonical CCSD(T) exhibits the largest error for this problem. Complementary to CC, we also used phaseless auxiliary-field quantum Monte Carlo calculations (ph-AFQMC). We show that with a Hartree-Fock (HF) trial wave function, ph-AFQMC reproduces the CC results within 5 ± 1 kJ/mol. With multi-Slater-determinant (MSD) trials, the results are 82-84 ± 2 kJ/mol, indicating that multireference effects may be rather modest. Among the DFT methods tested, τ-HCTH, r2SCAN with 10-13% HF exchange with and without dispersion, and O3LYP/O3LYP-D4, and B3LYP*/B3LYP*-D4 generally perform the best. The r2SCAN12 (with 12% HF exchange) functional mimics both the best reference MSD ph-AFQMC and CC ΔE results within 2 kJ/mol.

Assessment of DFT functionals for a minimal nitrogenase [Fe(SH)4H]- model employing state-of-the-art abinitio methods.

We have designed a [Fe(SH)4H]- model with the fifth proton binding either to Fe or S. We show that the energy difference between these two isomers (∆E) is hard to estimate with quantum-mechanical (QM) methods. For example, different density functional theory (DFT) methods give ∆E estimates that vary by almost 140 kJ/mol, mainly depending on the amount of exact Hartree-Fock included (0%-54%). The model is so small that it can be treated by many high-level QM methods, including coupled-cluster (CC) and multiconfigurational perturbation theory approaches. With extrapolated CC series (up to fully connected coupled-cluster calculations with singles, doubles, and triples) and semistochastic heat-bath configuration interaction methods, we obtain results that seem to be converged to full configuration interaction results within 5 kJ/mol. Our best result for ∆E is 101 kJ/mol. With this reference, we show that M06 and B3LYP-D3 give the best results among 35 DFT methods tested for this system. Brueckner doubles coupled cluster with perturbaitve triples seems to be the most accurate coupled-cluster approach with approximate triples. CCSD(T) with Kohn-Sham orbitals gives results within 4-11 kJ/mol of the extrapolated CC results, depending on the DFT method. Single-reference CC calculations seem to be reasonably accurate (giving an error of ∼5 kJ/mol compared to multireference methods), even if the D1 diagnostic is quite high (0.25) for one of the two isomers.

N2 binding to the E0-E4 states of nitrogenase.

Nitrogenase is the only enzyme that can convert N2 into NH3. The reaction requires the addition of eight electrons and protons to the enzyme and the mechanism is normally described by nine states, E0-E8, differing in the number of added electrons. Experimentally, it is known that three or four electrons need to be added before the enzyme can bind N2. We have used combined quantum mechanical and molecular mechanics methods to study the binding of N2 to the E0-E4 states of nitrogenase, using four different density functional theory (DFT) methods. We test many different structures for the E2-E4 states and study binding both to the Fe2 and Fe6 ions of the active-site FeMo cluster. Unfortunately, the results depend quite strongly on the DFT methods. The TPSS method gives the strongest bonding and prefers N2 binding to Fe6. It is the only method that reproduces the experimental observation of unfavourable binding to the E0-E2 states and favourable binding to E3 and E4. The other three methods give weaker binding, preferably to Fe2. B3LYP strongly favours structures with the central carbide ion triply protonated. The other three methods suggest that states with the S2B ligand dissociated from either Fe2 or Fe6 are competitive for the E2-E4 states. Moreover, such structures with two hydride ions both bridging Fe2 and Fe6 are the best models for E4 and also for the N2-bound E3 and E4 states. However, for E4, other structures are often close in energy, e.g. structures with one of the hydride ions bridging instead Fe3 and Fe7. Finally, we find no support for the suggestion that reductive elimination of H2 from the two bridging hydride ions in the E4 state would enhance the binding of N2.

Investigation of the Structural and Thermochemical Properties of [2.2.2]-Bicyclooctadiene Photoswitches.

Molecular photoswitches can under certain conditions be used to store solar energy in the so-called molecular solar thermal storage systems, which is an interesting technology for renewable energy solutions. The current investigations focus on the performance of seven different density functional theory (DFT) methods (B3LYP, CAM-B3LYP, PBE0, M06-2X, ωB97X-D, B2PLYP, and PBE0DH) when predicting geometries and thermochemical properties of the [2.2.2]-bicyclooctadiene (BOD) photoswitch. We find that all of the investigated DFT methods provide geometries that are in good agreement with those obtained using coupled cluster singles and doubles (CCSD) calculations. The dependence on the employed basis set is not large when predicting geometries. With respect to the thermochemical properties, we find that the M06-2X, CAM-B3LYP, PBE0, and ωB97X-D functionals all predict thermochemical properties that are in good agreement with the results of the CCSD, the CCSD including perturbative triples (CCSD(T)), and the explicitly correlated CCSD-F12 and CCSD(T)-F12 models. Lastly, for energy calculations, we tested the newly developed fourth-order cluster perturbation theory singles and doubles CPS(D-4) model, which in this study provides energy differences that are of CCSD and sometimes also CCSD(T) quality at a relatively low cost. We find that the CPS(D-4) model is an excellent choice for further investigation of BOD derivatives because accurate energies can be obtained routinely using this methodology. From the results, we also note that the predicted storage energies and storage energy densities for the BOD photoswitch are very large compared to other molecular solar thermal storage systems and that these systems could be candidates for such applications.

DFT Study of Photochemical Properties and Radiative Forcing Efficiency Features of the Stereoisomers cis- and trans-CHCl═CH-CF3.

The accurate assessment of radiative forcing efficiency (RFE) of a greenhouse gas is based on the precise knowledge of its structure and infrared absorption spectrum. The present work investigates the UV-vis absorption spectra and IR absorption spectra that are used for the determination of RFE of the short-lived compounds, cis- and trans-CHCl═CHCF3 (CTFP). These investigations were carried out with six different density functional theory (DFT) methods B3LYP, CAM-B3LYP, M06, M06-2X, TPSS0, and ωB97X-D associated to the basis set 6-31G(3df). Therefrom, the relative populations of the two states cis and trans for temperatures over the range 220-370 K at 1 atm and along the atmospheric altitude were assessed. It turns out that trans-CTFP is the abundant component between the two states. This review reveals that B3LYP and M06 reproduce well the experimental results of UV-vis spectra of trans-CTFP. As for cis-CTFP, EOM-CCSD is not well fitted by DFT methods. The cis- to trans-CTFP isomerization leads to the red shift for DFT methods and to the blue shift in regard to EOM-CCSD and experimental results. The IR absorption spectra are well fitted by B3LYP over the range 500-1600 cm-1 and TPSS0 over 1300-2000 cm-1 for both stereoisomers. Moreover, the root-mean-square errors (RMSEs) of frequencies from experimental data are lower for B3LYP and TPSS0 for both systems. The computed IR absorption band strengths over 500-2000 cm-1 for cis- and over 600-1800 cm-1 for trans-CTFP are consistent with the experiment. The relevant descriptor RFEs of the climate effect were calculated using a narrow band model for a constant vertical profile and then corrected with a lifetime factor for different computational methods. The computed values correlate well with the experimental results for both stereoisomers except M06-2X and TPSS0. It is worth noting that, for both systems, the intense radiative forcing spectra are located at frequencies ranging in 1000-1200 cm-1. The lower forcings of trans-CTFP lying in the atmospheric window region 800-1000 cm-1 are greater than those of cis-CTFP. Therefore, RFE(trans-CTFP) = 1.127 RFE(cis-CTFP).

Fragment-Based Approaches for Supramolecular Interaction Energies: Applications to Foldamers and Their Complexes with Anions.

We explore the application of our multilayer Molecules-in-Molecules (MIM) fragment-based method for the study of the energies in supramolecular systems, viz. foldamers and their anion bound complexes. The performance of five different density functional theory (DFT) methods in conjunction with the fragmentation-based method is evaluated against the unfragmented energies for a test set of 5 foldamers (82 to 170 atoms). A systematic protocol has been developed to account for the π···π interactions in such systems in addition to the traditional fragmentation of the system along the backbone comprised of covalently bonded dimer (or trimer) units. We find a significant improvement in the performance of the method on going from a one-layer MIM1 model (errors >10 kcal/mol) to a two-layer MIM2 model (errors 0-2 kcal/mol), due to the incorporation of long-range interactions in the latter approach. Furthermore, we extend the applicability of MIM2 models to determine accurate binding energies of macromolecular receptor-anion complexes. For three different anion bound macrocycles, our MIM2 protocol provides total energies within 1.5 kcal/mol of the unfragmented energies for most of the DFT methods. The corresponding anion binding energies are calculated within 0.5 kcal/mol of the unfragmented binding energies due to systematic error cancellation between the macrocycle and the macrocycle-anion complex. Finally, we have calibrated the absolute accuracy in the calculated binding energies by comparison with unfragmented DLPNO-CCSD(T) calculations on three macromolecule-chloride anion complexes. The most accurate results are obtained using a MIM2 model using DLPNO-CCSD(T) calculations on trimer units as the high level and DFT-D3 (e.g., M06-2X-D3) as the low level of theory, yielding sub kcal/mol errors in the anion binding energies. Our protocol can be an accurate method to calculate anion binding energies for very large supramolecular systems.

A Computational Comparison of Oxygen Atom Transfer Catalyzed by Dimethyl Sulf­oxide Reductase with Mo and W

AbstractA thorough computational study has been performed to investigate oxygen atom transfer (OAT) reactions catalyzed by dimethyl sulfoxide reductase (DMSOR) with a catalytic molybdenum or tungsten ion. Thirteen different density functional theory (DFT) methods have been employed to obtain structural parameters along the reaction pathway, and single‐point energies were computed with local correlation coupled‐cluster methods [LCCSD(T0)]. For both Mo and W, most DFT methods indicate that the enzyme follows a two‐step mechanism with a stable intermediate in which a DMSO molecule coordinates to the metal ion in the +IV oxidation state, and this is also confirmed by the LCCSD(T0) energies. The W‐substituted models have a 26–30 kJ/mol lower activation barrier for the OAT reaction, and the reaction is 63–70 kJ/mol more exothermic than that with Mo. Different DFT methods give widely different activation and reaction energies, which roughly depend on the amount of exact exchange in the method; these differences are also reflected in the structures, especially for the rate‐limiting transition state. Consequently, there is quite a large variation in energies and various energy corrections (thermal, solvation, dispersion, and relativistic; up to 39 kJ/mol) depending on which DFT method is used to obtain the geometries. Therefore, a mechanism predicted by a single method should be viewed with caution.

Determining the absolute configuration of benzopyrenomycin by optical rotation, electronic circular dichroism, and population analysis of different conformations via DFT methods and experiments

In this work, we have studied the absolute configuration of benzopyrenomycin using different density functional theory (DFT) methods, such as optical rotation (OR), electronic circular dichroism (ECD), and conformation distribution analysis for three kinds of Mosher esters at different levels, e.g., B3LYP/6-31G(d), B3LYP/6-311+G(d), and B3LYP/6-311++G(2d,p) in the gas phase and in solution, respectively. Careful investigations for different chiral Mosher esters using DFT theory exhibited the application conditions for certain chiral molecules in absolute configuration determination. Benzopyrenomycin possesses a unique benzo[a]pyrene-type skeleton and shows strong cytotoxicity against various tumor-cell lines.

Theoretical structure and vibrational spectra of ciprofloxacin: Density functional theory study

The molecular structure and vibrational spectra of ciprofloxacin(1-cyclopropyl-6-fluoro-1,4-dihydro-4-oxo-7-(1-piperazinyl)-3-quinolinecarboxylic acid) have been calculated using the different density functional theory (DFT) methods and various basis sets. This paper examines the comparative performance of different DFT methods at various basis sets in predicting molecular and vibrational spectra of ciprofloxacin. The calculation results show that SVWN/LANL2DZ level and SVWN/6-31G level offer the highest certainty to predict the structure and vibrational spectra of ciprofloxacin, respectively.

Density functional theory study on the molecular structure and vibration spectra of fenbufen

This paper examines the comparative performance of different density functional theory (DFT) methods at various basis sets in predicting molecular and vibration spectra of gamma-oxo-[1,1′-biphenyl]-4-butanoic acid (fenbufen). DFT methods including mPW1PW91, HCTH, SVWN, PBEPBE, B3PW91 and B3LYP were investigated. Different basis sets including LANL2DZ, SDD, LANL2MB, CEP-4G, CEP-31G, and CEP-121G were also considered. It is remarkable that the mPW1PW91/LANL2DZ and mPW1PW91/SDD levels are clearly superior to all of the remaining DFT levels in predicting the structure of fenbufen. The calculated results also indicate that SVWN/LANL2DZ level show better performance in the vibration spectra prediction of fenbufen comparing other DFT methods.

Comparison of the performance of different DFT methods in the calculations of the molecular structure and vibration spectra of serotonin (5-hydroxytryptamine, 5-HT)

► We investigated the molecular structure and vibrational spectra of 5-HT by DFT. ► LSDA/6-311++G and 6-311+G * afford the best quality to predict the structure of 5-HT. ► PBEPBE/LANL2DZ show better performance to predict the vibration spectra of 5-HT. Serotonin (5-hydroxytryptamine, 5-HT) is a monoamine neurotransmitter which plays an important role in treating acute or clinical stress. The comparative performance of different density functional theory (DFT) methods at various basis sets in predicting the molecular structure and vibration spectra of serotonin was reported. The calculation results of different methods including mPW1PW91, HCTH, SVWN, PBEPBE, B3PW91 and B3LYP with various basis sets including LANL2DZ, SDD, LANL2MB, 6-31G, 6-311++G and 6-311+G * were compared with the experimental data. It is remarkable that the SVWN/6-311++G and SVWN/6-311+G * levels afford the best quality to predict the structure of serotonin. The results also indicate that PBEPBE/LANL2DZ level show better performance in the vibration spectra prediction of serotonin than other DFT methods.

Comparison of DFT methods for molecular structure and vibration spectra of ofloxacin calculations

Comparison of the performance of different density functional theory (DFT) methods at various basis sets in predicting molecular and vibration spectra of ofloxacin was reported. The methods employed in this study comprise six functionals, namely, mPW1PW91, HCTH, LSDA, PBEPBE, B3PW91 and B3LYP. Different basis sets including LANL2DZ, SDD, LANL2MB, 6-31g, 6-311g and 3-21g were also examined. Comparison between the calculated and experimental data indicates that the mPW1PW91/6-311g level afford the best quality to predict the structure of ofloxacin. The results also indicate that B3LYP/LANL2DZ level show better performance in the vibration spectra prediction of ofloxacin than other DFT methods.

Thermochemistry, Bond Energies, and Internal Rotor Potentials of Dimethyl Tetraoxide

Thermochemical properties of dimethyl tetraoxide (CH(3)OOOOCH(3)), the dimer of the methylperoxy radical, are studied using ab initio and density functional theory methods. Methylperoxy radicals are known to be important intermediates in the tropospheric ozone cycle, and the self-reaction of methylperoxy radicals, which is thought to proceed via dimethyl tetraoxide, leads to significant chain radical termination in this process. Dimethyl tetraoxide has five internal rotors, three of them unique; the potential energy profiles are calculated for these rotors, as well as for those in the CH(3)OO, CH(3)OOO, and CH(3)OOOO radicals. The dimethyl tetraoxide internal rotor profiles show barriers to rotation of 2-8 kcal mol(-1). Using B3LYP/6-31(d) geometries, frequencies, internal rotor potentials, and moments of inertia, we determine entropy and heat capacity values for dimethyl tetraoxide and its radicals. Isodesmic work reactions with the G3B3 and CBS-APNO methods are used; we calculate this enthalpy as -9.8 kcal mol(-1). Bond dissociation energies (BDEs) are calculated for all C-O and O-O bonds in dimethyl tetraoxide, again with the G3B3 and CBS-APNO theoretical methods, and we suggest the following BDEs: 46.0 kcal mol(-1) for CH(3)-OOOOCH(3), 20.0 kcal mol(-1) for CH(3)O-OOOCH(3), and 13.9 kcal mol(-1) for CH(3)OO-OOCH(3). From the BDE calculations and the isodesmic enthalpy of formation for dimethyl tetraoxide, we suggest enthalpies of 2.1, 5.8, and 1.4 kcal mol(-1) for the CH(3)OO, CH(3)OOO, and CH(3)OOOO radicals, respectively. We evaluate the suitability of 10 different density functional theory (DFT) methods for calculating thermochemical properties of dimethyl tetraoxide and its radicals with the 6-31G(d) and 6-311++G(3df,3pd) basis sets, using a variety of work reaction schemes. Overall, the best-performed DFT methods of those tested were TPSSh, BMK, and B1B95. Significant improvements in accuracy were made by moving from atomization to isodesmic work reactions, with most DFT methods yielding errors of less than 2 kcal mol(-1) with the 6-311++G(3df,3pd) basis set for isodesmic calculations on the dimethyl tetraoxide enthalpy. These isodesmic calculations were basis set consistent, with a considerable reduction in error found by using the 6-311++G(3df,3pd) basis set over the 6-31G(d) basis set. This was not the case, however, for atomization and bond dissociation work reactions, where the two basis sets returned similar results. Improved group additivity terms for the O-O-O moiety (O/O2 central atom group) are also determined.

Assessment of Density Functional Theory Methods for the Computation of Heats of Formation and Ionization Potentials of Systems Containing Third Row Transition Metals

The performance of several different density functional theory (DFT) methods, including GGA, hybrid-GGA, meta-GGA, and hybrid-meta-GGA methods, have been assessed in terms of their ability to accurately compute both heats of formation and ionization potentials of systems containing third row transition metals. Two different basis sets were used in this study: 6-31G** and TZVP. It is found that the triple-zeta quality TZVP basis set generally produces the best results for both heats of formation and ionization potentials. One important observation made in this study is that the inclusion of exact exchange terms in DFT methods generally results in more consistently accurate results for both heats of formation and ionization potentials of transition metal systems. In general, DFT methods do not yield good ionization potential results for systems containing titanium or zinc. For heats of formation, it is found that the hybrid-meta-GGA functional, TPSS1KCIS, yields the best overall results when combined with the TZVP basis set, while PBE1PBE (hybrid-GGA) gives the best results for the 6-31G** basis. The hybrid-GGA functional, B3LYP, is found to produce the lowest overall errors for ionization potentials when combined with both 6-31G** and TZVP.

A density functional theory benchmark of the formation enthalpy and first CO dissociation enthalpy of hexacarbonyl complexes of chromium, molybdenum, and tungsten

Reaction energies have been calculated for the hexacarbonyl metal complexes of chromium, molybdenum, and tungsten using six different Density Functional Theory (DFT) methods. Two gas phase reactions have been utilized to benchmark the computed results. The reactions are the heat of formation of the hexacarbonyl complex from its metal dioxide, dicarbon and dioxygen; and the first dissociation enthalpy of the carbonyl ligand. It was found that all DFT methods agreed well with the available experimental data. For the formation reaction, all methods underestimate the formation enthalpy, while non-hybrid methods do better than the methods involving functional hybridization. Including a diffuse function in the basis set decreases the magnitude of the formation enthalpy. For the CO dissociation reaction the LYP methods underestimate the bond enthalpy consistently for all the metals. There was not much difference between using a hybrid and a non-hybrid functionals, except when comparing the PWP91 and B3PW91 methods, but it was found that the level of exchange and correlation affects the result, especially when using the LYP-based methods. A comparison of calculated geometries and frequencies is also made with available experimental and computational data, confirming that the ground state of WO 2 is a singlet.

Binuclear homoleptic copper carbonyls Cu(2)(CO)(x) (x = 1-6): remarkable structures contrasting metal-metal multiple bonding with low-dimensional copper bonding manifolds.

Binuclear homoleptic copper carbonyls Cu(2)(CO)(x) (x = 1-6) have been studied using four different density functional theory methods (DFT) in conjunction with a basis set of extended double-zeta plus polarization quality, labeled as DZP. For each homoleptic binuclear copper carbonyl compound, several stationary point structures are presented, and these structures are characterized in terms of their geometries, thermochemistry, and vibrational frequencies. The optimal unsaturated Cu(2)(CO)(x) (x = 1-6) structures are generated by joining 18-electron tetrahedral, 16-electron trigonal, 14-electron linear copper carbonyl building blocks, and/or bare copper atoms with copper-copper single bonds rather than by joining 18-electron copper carbonyl units with multiple copper-copper bonds. For Cu(2)(CO)(6) the eclipsed and staggered ethane-like structure are virtually degenerate and lie significantly lower in energy than other possible structures. The eclipsed Cu-Cu single bond distance is predicted to be 2.61 A, while that for the staggered structure is 2.65 A. The lowest energy structure for Cu(2)(CO)(5) is the eclipsed ethyl radical-like structure, with r(e)(Cu-Cu) = 2.51 A. The staggered ethyl radical-like structure lies only 0.1 kcal/mol higher in energy, with a Cu-Cu distance shorter by only approximately 0.001 A. For Cu(2)(CO)(4) a methylcarbene-like structure is predicted to lie lowest, with Cu-Cu distance 2.40 A. However, twisted and planar ethylene-like structure lie only 3-5 kcal/mol higher. For Cu(2)(CO)(3) a surprising methylcarbyne-like structure with r(e)(Cu-Cu) = 2.38 A is predicted to lie lowest with all four DFT methods. However, a classical vinyl radical-like lies only 2-4 kcal/mol higher. For Cu(2)(CO)(2) theory predicts a vinylidene-like structure with r(e)(Cu-Cu) = 2.34 A to be essentially degenerate with cis and trans bent acetylene structures with copper-copper distances 2.33 A. Finally, and consistent with earlier theoretical studies, the linear end on Cu-Cu-CO structure with r(e)(Cu-Cu) = 2.27 A is the predicted global minimum for Cu(2)(CO).

An investigation of density functionals: The first-row transition metal dimer calculations

The performance of different density functional theory (DFT) methods was investigated in the calculations of the bond length and the binding energy of the first-low transition metal dimers. The 4s–3d interconfigurational energies and 4s and 3d ionization potentials were also calculated for the first-row transition metal atoms. In general, the hybrid DFT method, B3LYP, yields the bond lengths that are too short compared to the experimental ones. In contrast, the optimized bond lengths by nonhybrid DFT methods such as BOP or PW91 are in good agreement with the experiment. It was also found that nonhybrid DFT methods overestimate the binding energies, because they have a tendency to overstabilize the electron configurations that contain the atomic and molecular orbitals in a higher angular momentum open shell. The hybrid DFT method yields more accurate binding energies, but it estimates rather poor energy gaps between states whose spin multiplicity is quite different.

Properties of phosphorus compounds by density functional theory: CH3P species as a test case

A comparison of different density functional theory (DFT) and molecular orbital (MO) methods for calculating molecular and energetic properties of low-coordinated phosphorus compounds is reported. While DFT methods include both Becke–Lee–Yang–Parr (BLYP and B3LYP) nonlocal functionals, MO methods involve second-order perturbation theory (MP2), quadratic configuration interaction [QCISD(T)], and coupled-cluster theory [CCSD(T)], in conjunction with the 6-31G(d,p), 6-311++G(3df,2p), and 6-311++G(3df,3pd) basis sets. Properties examined include geometrical parameters of the different CH3P equilibrium structures (phosphaethene, phosphinocarbene, methylphosphinidene, and a phosphacarbyne) and relevant transition structures for isomerisations and rearrangements in both the lowest-lying singlet and triplet states, vibrational wave numbers, relative energies, barrier heights, and singlet–triplet energy gaps. In addition, the heat of formation, ionization energy, and proton affinity of phosphaethene are also evaluated. Overall, the B3LYP method, when employed with a large basis set, yields energetic results comparable to the CCSD(T) results. Nevertheless, both DFT methods fail to predict the behavior of the addition/elimination reactions of the hydrogen atom in the triplet state.