Accurate bond dissociation enthalpies by using doubly hybrid XYG3 functional

In this work, we examine the performance of XYG3, a newly developed doubly hybrid density functional (Zhang, Xu, and Goddard III, Proc Natl Acad Sci USA 2009, 106, 4963), to calculate covalent bond dissociation enthalpy (BDE). We use 5 atoms, 32 molecular radicals, and 116 closed-shell molecules to set up 142 bond dissociation reactions. For the total of 148 heats of formation (HOFs) and 142 BDEs, XYG3 leads to mean absolute deviations (MADs) of 1.45 and 1.87 kcal/mol, respectively. In comparison with some other functionals, MADs for HOFs are 2.31 (M06-2X), 2.98 (B2PLYP-D), 3.04 (BMK), 3.96 (B3LYP), 4.47 (B2PLYP), 5.42 (B2GP-PLYP), 6.46 (PBE0), and 29.93 kcal/mol (B3P86), and the corresponding errors for BDEs are 2.06 (M06-2X), 2.25 (BMK), 2.51 (B2PLYP-D), 2.89 (B2GP-PLYP), 3.30 (B3P86), 3.44 (B2PLYP), 3.87 (PBE0), and 6.14 kcal/mol (B3LYP).     

一些密度泛函方法预测电子亲和势

选取了杂化泛函B3LYP, B3PW91, O3LYP, PBE0, 以及与之相对应的GGA泛函BLYP, BPW91, OLYP和PBE, 还选取了能更好地兼顾强相互作用和弱相互作用的X3LYP泛函和在预测NMR的化学位移有较好表现的OPBE泛函, 以及两种meta-GGA泛函VSXC和TPSS, 共12种泛函, 详细地考察了这些泛函在预测EA方面的准确性.     

Ring-opening radical clock reactions: many density functionals have difficulty keeping time

The ability of several density-functional theory methods to describe the kinetics and energetics of a series of ring-opening reactions of cyclopropyl and cyclobutyl-type radicals was explored. PBE, B971 and B3LYP perform quite well in their ability to replicate experiment, based upon the ring opening of cyclopropylcarbinyl, two α-trialkylsilyloxycyclopropylmethyl radicals, pentamethylcyclopropylcarbinyl, cyclobutylcarbinyl and 1-cyclobutylethylcarbinyl. The other functionals tested, which includes BLYP, CAM-B3LYP, BHandHLYP, B2PLYP and B2PLYP-D, as well as functionals designed for kinetics applications, namely MPW1K, BMK and M06-2X, all perform poorly. The latter of these functionals display some integration grid dependencies.     

Assessment of density functionals with long‐range and/or empirical dispersion corrections for conformational energy calculations of peptides

Density functionals with long‐range and/or empirical dispersion corrections, including LC‐ωPBE, B97‐D, ωB97X‐D, M06‐2X, B2PLYP‐D, and mPW2PLYP‐D functionals, are assessed for their ability to describe the conformational preferences of Ac‐Ala‐NHMe (the alanine dipeptide) and Ac‐Pro‐NHMe (the proline dipeptide) in the gas phase and in water, which have been used as prototypes for amino acid residues of peptides. For both dipeptides, the mean absolute deviation (MAD) is estimated to be 0.22–0.40 kcal/mol in conformational energy and 2.0–3.2° in torsion angles ? and ψ using these functionals with the 6‐311++G(d,p) basis set against the reference values calculated at the MP2/aug‐cc‐pVTZ//MP2/aug‐cc‐pVDZ level of theory in the gas phase. The overall performance is obtained in the order B2PLYP‐D ≈ mPW2PLYP‐D > ωB97X‐D ≈ M06‐2X > MP2 > LC‐ωPBE > B3LYP with the 6–311++G(d,p) basis set. The SMD model at the M06‐2X/6‐31+G(d) level of theory well reproduced experimental hydration free energies of the model compounds for backbone and side chains of peptides with MADs of 0.47 and 4.3 kcal/mol for 20 neutral and 5 charged molecules, respectively. The B2PLYP‐D/6‐311++G(d,p)//SMD M06‐2X/6‐31+G(d) level of theory provides the populations of backbone and/or prolyl peptide bond for the alanine and proline dipeptides in water that are consistent with the observed values. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Improved description of nuclear magnetic resonance chemical shielding constants using the M06-L meta-generalized-gradient-approximation density functional

The performance of the M06-L density functional has been tested for four databases of NMR isotropic chemical shielding constants. Comparison with the B3LYP, BLYP, HCTH, KT1, KT2, LSDA, OPBE, OLYP, PBE, TPSS, and VSXC functionals shows that M06-L has improved performance for calculating NMR chemical shielding constants, especially for highly correlated systems. We also found that VSXC and M06-L have encouraging accuracy for calculating (13)C chemical shielding constants, and both functionals perform very well for the chemical shielding constants in the o-benzyne molecule.     

Electronic spectroscopy of UO2(2+), NUO(+) and NUN: an evaluation of time-dependent density functional theory for actinides

The performance of the time-dependent density functional theory (TDDFT) approach has been evaluated for the electronic spectrum of the UO(2)(2+), NUO(+) and NUN molecules. Different exchange-correlation functionals (LDA, PBE, BLYP, B3LYP, PBE0, M06, M06-L, M06-2X, CAM-B3LYP) and the SAOP model potential have been investigated, as has the relative importance of the adiabatic local density approximation (ALDA) to the exchange-correlation kernel. The vertical excitation energies have been compared with reference data obtained using accurate wave-function theory (WFT) methods.     

TD-CI simulation of the electronic optical response of molecules in intense fields II: comparison of DFT functionals and EOM-CCSD

Time-dependent configuration interaction (TD-CI) simulations can be used to simulate molecules in intense laser fields. TD-CI calculations use the excitation energies and transition dipoles calculated in the absence of a field. The EOM-CCSD method provides a good estimate of the field-free excited states but is rather expensive. Linear-response time-dependent density functional theory (TD-DFT) is an inexpensive alternative for computing the field-free excitation energies and transition dipoles needed for TD-CI simulations. Linear-response TD-DFT calculations were carried out with standard functionals (B3LYP, BH&HLYP, HSE2PBE (HSE03), BLYP, PBE, PW91, and TPSS) and long-range corrected functionals (LC-ωPBE, ωB97XD, CAM-B3LYP, LC-BLYP, LC-PBE, LC-PW91, and LC-TPSS). These calculations used the 6-31G(d,p) basis set augmented with three sets of diffuse sp functions on each heavy atom. Butadiene was employed as a test case, and 500 excited states were calculated with each functional. Standard functionals yield average excitation energies that are significantly lower than the EOM-CC, while long-range corrected functionals tend to produce average excitation energies slightly higher. Long-range corrected functionals also yield transition dipoles that are somewhat larger than EOM-CC on average. The TD-CI simulations were carried out with a three-cycle Gaussian pulse (ω = 0.06 au, 760 nm) with intensities up to 1.26 × 10(14) W cm(-2) directed along the vector connecting the end carbons. The nonlinear response as indicated by the residual populations of the excited states after the pulse is far too large with standard functionals, primarily because the excitation energies are too low. The LC-ωPBE, LC-PBE, LC-PW91, and LC-TPSS long-range corrected functionals produce responses comparable to EOM-CC.     

Accurate prediction of vertical electronic transitions of Ni(II) coordination compounds via time dependent density functional theory

Time dependent density functional theory calculations are completed for five Ni(II) complexes formed by polydentate peptides to predict the electronic absorption spectrum. The ligands examined were glycyl‐glycyl‐glycine (GGG), glycyl‐glycyl‐glycyl‐glycine (GGGG), glycyl‐glycyl‐histidine (GGH), glycyl‐glycyl‐cysteine (GGC), and triethylenetetramine (trien). Fifteen functionals and two basis sets were tested. On the basis of the mean absolute percent deviation (MAPD), the ranking among the functionals is: HSE06 ∼ MPW1PW91 ∼ PBE0 > ω‐B97x‐D ∼ B3P86 ∼ B3LYP ∼ CAM‐B3LYP > PBE ∼ BLYP ∼ BP86 > TPSS > TPSSh > BHandHLYP > M06 ≫ M06‐2X. Concerning the basis sets, the triple‐ζ def2‐TZVP performs better than the double‐ζ LANL2DZ. With the functional HSE06 and basis set def2‐TZVP the MAPD with respect to the experimental λ_max is 1.65% with a standard deviation of 1.26%. The absorption electronic spectra were interpreted in terms of vertical excitations between occupied and virtual MOs based on Ni‐d atomic orbitals. The electronic structure of the Ni(II) species is also discussed.     

Assessing alkyl-, silyl-, and halo-substituent effects on the electron affinities of silyl radicals

Neutral anion energy differences for a large class of alpha-substituted silyl radicals have been computed to determine the effect of alkyl, silyl, and halo substituents on their electron affinities. In particular, we report theoretical predictions of the adiabatic electron affinities (AEAs), vertical electron affinities (VEAs), and vertical detachment energies (VDEs) for a series of methyl-, silyl-, and halo-substituted silyl radical compounds. This work utilizes the carefully calibrated DZP++ basis set, in conjunction with the pure BLYP and OLYP functionals, as well as with the hybrid B3LYP, BHLYP, PBE1PBE, MPW1K, and O3LYP functionals. Bromine has the largest effect in stabilizing the anions, and the BLYP/DZP++ AEA for SiBr(3) is 3.29 eV. The other predicted electron affinities are for SiH(3) (1.37 eV), SiH(2)CH(3) (1.09 eV), SiH(2)F (1.54 eV), SiH(2)Cl (1.94 eV), SiH(2)Br (2.05 eV), SiH(2)(SiH(3)) (1.77 eV), SiH(CH(3))(2) (0.92 eV), SiHF(2) (1.86 eV), SiHCl(2) (2.53 eV), SiHBr(2) (2.67 eV), Si(CH(3))(3) (0.86 eV), SiF(3) (2.66 eV), SiCl(3) (3.21 eV), Si(SiH(3))(3) (2.25 eV), and SiFClBr (3.13 eV). For the five silyl radicals where experimental data are available, the BLYP functional gives the most accurate determination of AEAs; the average absolute error is 0.04(1) eV, whereas the corresponding errors for the O3LYP, MPW1K, PBE1PBE, B3LYP, OLYP, and BHLYP functionals are 0.05(8), 0.06(0), 0.06(3), 0.08(5), 0.11(5), and 0.15(3) eV, respectively.     

Vibrational frequency scale factors for density functional theory and the polarization consistent basis sets

Calculated harmonic vibrational frequencies systematically deviate from experimental vibrational frequencies. The observed deviation can be corrected by applying a scale factor. Scale factors for: (i) harmonic vibrational frequencies [categorized into low (<1000 cm^?1) and high (>1000 cm^?1)], (ii) vibrational contributions to enthalpy and entropy, and (iii) zero‐point vibrational energies (ZPVEs) have been determined for widely used density functionals in combination with polarization consistent basis sets (pc‐n, n = 0,1,2,3,4). The density functionals include pure functionals (BP86, BPW91, BLYP, HCTH93, PBEPBE), hybrid functionals with Hartree‐Fock exchange (B3LYP, B3P86, B3PW91, PBE1PBE, mPW1K, BH&HLYP), hybrid meta functionals with the kinetic energy density gradient (M05, M06, M05‐2X, M06‐2X), a double hybrid functional with Møller‐Plesset correlation (B2GP‐PLYP), and a dispersion corrected functional (B97‐D). The experimental frequencies for calibration were from 41 organic molecules and the ZPVEs for comparison were from 24 small molecules (diatomics, triatomics). For this family of basis sets, the scale factors for each property are more dependent on the functional selection than on basis set level, and thus allow for a suggested scale factor for each density functional when employing polarization consistent basis sets (pc‐n, n = 1,2,3,4). A separate scale factor is recommended when the un‐polarized basis set, pc‐0, is used in combination with the density functionals. © 2012 Wiley Periodicals, Inc.     

Performance of Density Functional Theory for Transition Metal Oxygen Bonds

The accuracy of density functional theory (DFT) limits predictions in theoretical catalysis, and strong chemical bonds between transition metals and oxygen pose a particular challenge. We benchmarked 30 diverse density functionals against the bond dissociation enthalpies (BDE) of the 30 MO and 30 MO⁺ diatomic systems of all the 3d, 4d, and 5d metals, to test universality across the d-block as required in comparative studies. Seven functionals, B98, B97-1, B3P86, B2PLYP, TPSSh, B3LYP, and B97-2, display mean absolute errors (MAE) <30 kJ/mol. In contrast, many commonly used functionals such as PBE and RPBE overestimate M−O bonding by +30 kJ/mol and display MAEs from 48–76 kJ/mol. RPBE and OPBE reduce the over-binding of PBE but remain very inaccurate. We identify a linear relationship (p-value 7.6 ⋅ 10⁻⁵) between the precision and accuracy of DFT, i. e. inaccurate functionals tend to produce larger, unpredictable random errors. Some functionals commonly deviate from this relationship: Thus, M06-2X is very precise but not very accurate, whereas B3LYP* and MN15-L are more accurate but less precise than M06-2X. The best-performing hybrids have 10–30 % HF exchange, but this can be relieved by double hybrids (B2PLYP). Most functionals describe trends well, but errors comparing 5d to 4d/3d are ∼10 kJ/mol larger than group-wise errors, due to uncertainties in the spin-orbit coupling corrections for effective core potentials, affecting e. g. Pt/Pd or Au/Ag comparisons.     

DFT,DLPNO-CCSD(T), and NEVPT2 benchmark study of the reaction between ferrocenium and trimethylphosphine

The reaction between ferrocenium and trimethylphosphine was studied using density functional theory (DFT), domain-based local pair natural orbital coupled cluster theory with single-, double-, and perturbative triple excitations (DLPNO-CCSD(T)), and N-electron valence state perturbation theory (NEVPT2). The accuracy of the DFT functionals decreases compared to the DLPNO-CCSD(T) level in the following order: M06-L > TPSS > M06, BLYP > PBE, PBE0, B3LYP > > PWPB95 > > DSD-BLYP. The roles of thermochemical, continuum solvation (SMD), and counterpoise corrections were evaluated. Grimme's D3 empirical dispersion correction is essential for all functionals studied except M06 and M06-L. The reliability of the frequency calculations performed directly within the SMD was confirmed. The systems showed no significant multireference character according to T₁ and T₂ diagnostics and the fractional occupation number (FOD) weighted electron density analysis. The multireference NEVPT2 calculations gave qualitatively valid conclusions about the reaction mechanism. However, a multireference approach is generally not recommended because it requires arbitrary chosen active spaces.     

Accurate heats of formation of polycyclic saturated hydrocarbons predicted by using the XYG3 type of doubly hybrid functionals

Polycyclic saturated hydrocarbons (PSHs) are attractive candidates as hydrocarbon propellants. To assess their potential values, one of the key factors is to determine their energy contents, such as to calculate their heats of formation (HOF). In this work, we have calculated HOFs for a set of 36 PSHs including exo-Tricyclo[5.2.1.0^(2,6)] decane, the principal component of the high-energy density hydrocarbon fuel commonly identified as JP-10. The results from B3LYP, B3LYP-D3BJ, M06-2X, B2PLYP, B2PLYP-D3BJ, and the XYG3 type of doubly hybrid (xDH) functionals are presented. It is demonstrated here that the xDH functionals yield accurate HOFs in good agreement with those from experiments or the G4 theory. In particular, XYGJ-OS, a low scaling xDH functional, is shown to hold the promise for accurate prediction of HOFs for PSHs of larger sizes. © 2018 Wiley Periodicals, Inc.     

Computational studies on ethylene addition to nickel bis(dithiolene)

The density functionals B3LYP, B3PW91, BMK, HSE06, LC-ωPBE, M05, M06, O3LYP, TPSS, ω-B97X, and ω-B97XD are used to optimize key transition states and intermediates for ethylene addition to Ni(edt)(2) (edt = S(2)C(2)H(2)). The efficacy of the basis sets 6-31G**, 6-31++G**, cc-pVDZ, aug-cc-pVDZ, cc-pVTZ, and aug-cc-pVTZ is also examined. The geometric parameters optimized with different basis sets and density functionals are similar and agree well with experimental values. The ω-B97XD functional gives relative energies closest to those from CCSD, while M06 and HSE06 yield results close to those from CCSD(T). CASSCF and CASSCF-PT2 calculation results are also given. Variation of the relative energies from different density functionals appears to arise, in part, from the multireference character of this system, as confirmed by the T1 diagnostic and CASSCF calculations.     

二阶多参考态微扰理论计算电子亲和势

采用二阶多参考微扰理论计算了F, Cl, OH, SH, CN, CH₂和NH₂的电子亲和势.另外,还考察了基函数和完全活性空间大小对电子亲和势精度的影响.通过和CASSCF, CASPT2, CCSD, CCSD(T), B3LYP, X3LYP, M06, HCTH, TPSS, B97D3, mPW2PLYP和B2PLYP的计算结果比较发现,针对目前所用的基函数,二阶多参考态微扰理论的总体计算效果是最好的.     

The shape of gaseous n-butylbenzene: assessment of computational methods and comparison with experiments

Conformational landscape of neutral and ionized n-butylbenzene has been examined. Geometries have been optimized at the B3LYP/6-31G(d), B3LYP/6-31+G(d,p), B3LYP-D/6-31+G(d,p), B2PLYP/6-31+G(d,p), B2PLYP-D/6-31+G(d,p), B97-D/6-31+G(d,p), and M06-2X/6-31+G(d,p) levels. This study is complemented by energy computations using 6-311++G(3df,2p) basis set and CBS-QB3 and G3MP2B3 composite methods to obtain accurate relative enthalpies. Five distinguishable conformers have been identified for both the neutral and ionized systems. Comparison with experimentally determined rotational constants shows that the best geometrical parameters are provided by B3LYP-D and M06-2X functionals, which include an explicit treatment of dispersion effects. Composite G3MP2B3 and CBS-QB3 methods, and B2PLYP-D, B3LYP-D, B97-D, and M06-2X functionals, provide comparable relative energies for the two sets of neutral and ionized conformers of butyl benzene. An exception is noted however for conformer V(+) the stability of which being overestimated by the B3LYP-D and B97-D functionals. The better stability of neutral conformers I, III, and IV, and of cation I(+) , demonstrated by our computations, is in perfect agreement with conclusions based on micro wave, fluorescence, and multiphoton ionization experiments.     

XYG3 and XYGJ-OS performances for noncovalent binding energies relevant to biomolecular structures

XYG3 and XYGJ-OS doubly hybrid density functionals (DHDFs) are applied to calculate intermolecular interaction energies of biological relevance. The results are compared with those of the other types of DHDFs (MC3BB and B2PLYP), as well as their dispersion corrected variants (MC3BB-D B2PLYP-D, and B2PLYP-D3), in addition to those obtained from B3LYP, B3LYP-D and B3LYP-D3. The reference data are taken from the S22, S22x5 and JSCH-2005 benchmark sets consisting mainly of DNA base pairs and amino acid pairs. The overall good agreement with the reference values of the extrapolated CCSD(T) complete basis set limit suggests that the XYG3 and XYGJ-OS functionals are valuable tools for applications in biochemistry.     

用GGA密度泛函及其长程、色散校正方法计算各类氢键的结合能

应用广义梯度近似(GGA) (PW91和PBE)、含动能密度的广义梯度近似(meta-GGA) (M06-L)、杂化泛函(hyper-GGA)(M06-2X、X3LYP和B3LYP)及其长程校正泛函LC-DFT(CAM-B3LYP、LC-ωPBE和ωB97X)和色散校正密度泛函(DFT-D)(ωB97X-D和B97-D),用多种基函数对15种不同强度的传统氢键和非传统氢键体系的结合能进行了系统的计算与分析.并与高精度的CCSD(T)/aug-cc-pVQZ结果比较发现:在上述各类泛函中,对于氢键结合能的计算M06-2X和ωB97X-D泛函较为精确与可靠,且没有必要使用过大的基函数,6-311++G(2d,2p)或aug-cc-pVDZ水平的基组就已足够,各类泛函所计算结合能的基组重叠误差(BSSE)均较小,除ωB97X和ωB97X-D外,其它9种泛函不经BSSE校正也能得到同样甚至更准确的结果.