B3LYP,BLYP and PBE DFT band structures of the nucleotide base stacks

DFT crystal orbital (band structure) calculations have been performed for the nucleotide base stacks of cytosine, thymine, adenine, and guanine arranged in DNA B geometry. The band structures obtained with PBE, BLYP, and B3LYP functionals are presented and compared to other related experimental and theoretical results. The influence of the quality of the basis set on the fundamental gap values was also investigated using Clementi's double ζ, 6‐31G and 6‐31G* basis sets. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Role of nonlocal exchange in molecular crystals: The case of two proton‐ordered phases of ice

We present a periodic density functional theory investigation of twoproton‐ordered phases of ice. Their equilibrium lattice parameters,relative stabilities, formation energies, and densities of states havebeen evaluated. Nine exchange‐correlation functionals, representativeof the generalized gradient approximation (GGA), global hybrids,range‐separated hybrids, meta‐GGA, and hybrid meta‐GGA families havebeen taken into account, considering two oxygen basis sets. Althoughthe hydrogen‐bond network of ice is well reproduced at the B3LYP,M06‐L, or LC‐ wPBE levels, formation energies are only correctlyevaluated with the two former functionals. Band gaps on the other handare only quantitatively reproduced at the B3LYP level. These resultsindicate that this last functional, a de facto reference formolecular calculations, gives in average the most accurate results forthe considered ice properties. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Influence of the exchange‐correlation functional in all‐electron calculations of the vibrational frequencies of corundum (α‐Al2O3)

The equilibrium structural parameters, high‐ and low‐frequency dielectric tensors, Born effective charges, and Γ‐point vibrational frequencies of bulk Al₂O₃ corundum are calculated by using the periodic, ab initio program CRYSTAL, which adopts an all‐electron Gaussian‐type basis set. The effect of basis set and the performance of three different functionals, i.e., LDA, PW91, and B3LYP, are discussed. The mean absolute deviation from the measured frequencies is as small as 7 cm^?1 for both the LDA and B3LYP functionals, indicating that these functionals perform extremely well in this case. The mean absolute deviation increases to 18 cm^?1 when the PW91 functional is used. All three functionals reproduce the equilibrium geometry of corundum to a high level of accuracy, with LDA and B3LYP outperforming PW91 slightly. The comparison of the current all‐electron calculations with previous plane‐wave, pseudo‐potential calculations shows an overall similar performance. The results of isotopic substitution for both Al and O are also presented. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

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.     

The confirmation of accurate combination of functional and basis set for transition‐metal dimers: Fe2, Co2, Ni2, Ru2, Rh2, Pd2, Os2, Ir2, and Pt2

Eight kinds of density functionals named B3LYP, PBE1PBE, B1B95, BLYP, BP86, G96PW91, mPWPW91, and SVWN along with two different valence basis sets (LANL2DZ and CEP‐121g) are employed to study the transition‐metal dimers for the elements of group VIII. By comparing the equilibrium bond distances, vibrational frequencies, and dissociation energies of the ground state of these dimers with the available experimental values and theoretical data, we show that the “pure” DFT methods (G96PW91, BLYP, and BP86) with great‐gradient approximation always give better results relative to the hybrid HF/DFT schemes (B3LYP, PBE1PBE, and B1B95). The striking case found by us is that the G96PW91 functional, which is not tested in previous systemic studies, always predicts the dissociation energy to be well. The Ru₂ and Os₂ dimers are sensitive to not only the functionals employed but also the valence basis sets adopted. The natural bond orbital population is analyzed, and the molecular orbitals of the unpaired electrons are determined. Furthermore, our results indicate that the s and d orbitals of these dimers always hybridize with each other except for Rh₂ and Pt₂ molecules. And by analyzing the electron configuration of the bonding atom, the dissociation limit of the ground state is obtained. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Which density functional is close to CCSD accuracy to describe geometry and interaction energy of small noncovalent dimers? A benchmark study using Gaussian09

A benchmark study on all possible density functional theory (DFT) methods in Gaussian09 is done to locate functionals that agree well with CCSD/aug‐cc‐pVTZ geometry and Ave‐CCSD(T)/(Q‐T) interaction energy (E_int) for small non‐covalently interacting molecular dimers in “dispersion‐dominated” (class 1), “dipole‐induced dipole” (class 2), and “dipole‐dipole” (class 3) classes. A DFT method is recommended acceptable if the geometry showed close agreement to CCSD result (RMSD < 0.045) and E_int was within 80–120% accuracy. Among 382 tested functionals, 1–46% gave good geometry, 13–44% gave good E_int, while 1–33% satisfied geometry and energy criteria. Further screening to locate the best performing functionals for all the three classes was made by counting the acceptable values of energy and geometry given by each functionals. The meta‐generalized gradient approximation (GGA) functional M06L was the best performer with total 14 hits; seven acceptable energies and seven acceptable geometries. This was the only functional “recommended” for at least two dimers in each class. The functionals M05, B2PLYPD, B971, mPW2PLYPD, PBEB95, and CAM‐B3LYP gave 11 hits while PBEhB95, PW91B95, Wb97x, BRxVP86, BRxP86, HSE2PBE, HSEh1PBE, PBE1PBE, PBEh1PBE, and PW91TPSS gave 10 hits. Among these, M05, B971, mPW2PLYPD, Wb97x, and PW91TPSS were among the “recommended” list of at least one dimer from each class. Long‐range correction (LC) of Hirao and coworkers to exchange‐correlation functionals showed massive improvement in geometry and E_int. The best performing LC‐functionals were LC‐G96KCIS and LC‐PKZBPKZB. Our results predict that M06L is the most trustworthy DFT method in Gaussian09 to study small non‐covalently interacting systems. © 2013 Wiley Periodicals, Inc.     

Benchmarking density functionals in conjunction with Grimme's dispersion correction for noble gas dimers (Ne2, Ar2, Kr2, Xe2, Rn2)

Eleven exchange‐correlational functionals of different types corrected for dispersion by Grimme's D3 correction in conjunction with the aug‐cc‐pVTZ basis set were tested on the following noble gas (Ng) dimers: Ne₂, Ar₂, Kr₂, Xe₂, and Rn₂. For comparison, the D2 and D3BJ corrections were probed with the B3LYP functional. From post‐HF wavefunction methods, CCSD(T) theory was also included. The investigated properties involved potential energy curves, equilibrium bond distances, and interaction energies. The B3LYP‐D3, B3LYP‐D3BJ, and PBE0‐D3 functionals performed overall best for bond distances, while B3LYP‐D3 and B97‐D3 performed best for interaction energies. The importance of fortunate error cancellations was seen in the often reduced agreement with reference data upon correction for BSSE. As several functionals performed well selectively for some noble gases (and poorly for others), we also analysed the performance on the Ng₂ dimers individually and recommended DFT‐D3 functionals for the calculation of large clusters of each Ng.     

Benchmarking of density functionals for the kinetics and thermodynamics of the hydrolysis of glycosidic bonds catalyzed by glycosidases

In recent years, there has been an increased interest in understanding the enzymatic mechanism of glycosidases resorting mostly to DFT and DFT/MM calculations. However, the performance of density functionals (DFs) is well known to be system and property dependent. Trends drawn from general studies, despite important to evaluate the quality of the DFs and to pave the way for the development of new DFs, may be misleading when applied to a single specific system/property. To overcome this issue, we carried out a benchmarking study of 40 DFs applied to the geometry optimization and to the electronic barrier height (E _Barrier) and electronic energy of reaction (E _R) of prototypical glycosidase‐catalyzed reactions. Additionally, we report calculations with SCC‐DFTB and four semiempirical MO methods applied to the same problem. We have used a universal molecular model for retaining glycosidases, comprising only a 22‐atoms system that mimics the active site and substrate. High accuracy reference geometries and energies were calculated at the CCSD(T)/CBS//MP2/aug‐cc‐pVTZ level of theory. Most DFs reproduce the reference geometries extremely well, with mean unsigned errors (MUE) smaller than 0.05 Å for bond lengths and 3° for bond angles. Among the DFs, wB97X‐D, CAM‐B3LYP, B3P86, and PBE1PBE have the best performance in geometry optimizations (MUE = 0.02 Å). Conversely, semiempirical MO and SCC‐DFTB methods yielded less accurate geometries (MUE between 0.09 and 0.17 Å). The inclusion of D3 correction has a small, but still relevant, influence in the geometry predicted by some DFs. Regarding E _Barrier, 11 DFs (MPW1B95, CAM‐B3LYP, M06 ‐ 2X, PBE1PBE, wB97X ‐ D, B1B95, BMK, MN12 – SX, M05, M06, and M11) presented errors below 1 kcal.mol^?1, in relation to the reference energy. Most of these functionals belong to the family of hybrid functionals (H‐GGA, HH‐GGA, and HM‐GGA), which shows a positive influence of HF exchange in the determination of E _Barrier. The inclusion of D3 correction has not affected significantly the E _Barrier and E _R. The use of geometries at the accurate but expensive MP2/aug‐cc‐pVTZ level of theory has a small, albeit not insignificant, influence in the E _Barrier when compared with energies calculated with geometries determined with the DFs (usually a few tenths of kcal.mol^?1, with exceptions). In general, semiempirical MO methods and DFTB are associated with larger errors in the determination of E _Barrier, with unsigned errors from 6.9 to 24.7 kcal.mol^?1.     

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.     

One‐bond 29Si‐1H spin‐spin coupling constants in the series of halosilanes: benchmark SOPPA and DFT calculations,relativistic effects,and vibrational corrections

A number of most representative second order polarization propagator approach (SOPPA) based wavefunction methods, SOPPA, SOPPA(CC2) and SOPPA(CCSD), and density functional theory (DFT) based methods, B3LYP, PBE0, KT2, and KT3, have been benchmarked in the calculation of the one‐bond ²⁹Si‐¹H spin‐spin coupling constants in the series of halosilanes SiH_nX_4?n (X = F, Cl, Br, I), both at the non‐relativistic and full four‐parameter Dirac's relativistic levels taking into account vibrational corrections. At the non‐relativistic level, the wavefunction methods showed much better results as compared with those of DFT. At the DFT level, out of four tested functionals, the Perdew, Burke, and Ernzerhof's PBE0 showed best performance. Taking into account, relativistic effects and vibrational corrections noticeably improves wavefunction methods results, but generally worsens DFT results. Copyright © 2013 John Wiley & Sons, Ltd.     

Calculation of the vibration frequencies of alpha-quartz: the effect of Hamiltonian and basis set

The central-zone vibrational spectrum of alpha-quartz (SiO2) is calculated by building the Hessian matrix numerically from the analytical gradients of the energy with respect to the atomic coordinates. The nonanalytical part is obtained with a finite field supercell approach for the high-frequency dielectric constant and a Wannier function scheme for the evaluation of Born charges. The results obtained with four different Hamiltonians, namely Hartree-Fock, DFT in its local (LDA) and nonlocal gradient corrected (PBE) approximation, and hybrid B3LYP, are discussed, showing that B3LYP performs far better than LDA and PBE, which in turn provide better results than HF, as the mean absolute difference from experimental frequencies is 6, 18, 21, and 44 cm(-1), respectively, when a split valence basis set containing two sets of polarization functions is used. For the LDA results, comparison is possible with previous calculations based on the Density Functional Perturbation Theory and usage of a plane-wave basis set. The effects associated with the use of basis sets of increasing size are also investigated. It turns out that a split valence plus a single set of d polarization functions provides frequencies that differ from the ones obtained with a double set of d functions and a set of f functions on all atoms by on average less than 5 cm(-1).     

Analysis of the bond‐valence method for calculating 29Si and 31P magnetic shielding in covalent network solids

²⁹Si and ³¹P magnetic‐shielding tensors in covalent network solids have been evaluated using periodic and cluster‐based calculations. The cluster‐based computational methodology employs pseudoatoms to reduce the net charge (resulting from missing co‐ordination on the terminal atoms) through valence modification of terminal atoms using bond‐valence theory (VMTA/BV). The magnetic‐shielding tensors computed with the VMTA/BV method are compared to magnetic‐shielding tensors determined with the periodic GIPAW approach. The cluster‐based all‐electron calculations agree with experiment better than the GIPAW calculations, particularly for predicting absolute magnetic shielding and for predicting chemical shifts. The performance of the DFT functionals CA‐PZ, PW91, PBE, rPBE, PBEsol, WC, and PBE0 are assessed for the prediction of ²⁹Si and ³¹P magnetic‐shielding constants. Calculations using the hybrid functional PBE0, in combination with the VMTA/BV approach, result in excellent agreement with experiment. © 2016 Wiley Periodicals, Inc.     

Calculations of ionization energies and electron affinities for atoms and molecules: A comparative study with different methods

In the present work, we examined the performance of 36 density functionals, including the newly developed doubly hybrid density functional XYG3 (Y. Zhang, X. Xu, and W. A. Goddard III, Proc. Natl. Acad. Sci, USA, 2009, 106, 4963), to calculate ionization energies (IEs) and electron affinities (EAs). We used the well-established G2-1 set as reference, which contains 14 atoms and 24 molecules for IE, along with 7 atoms and 18 molecules for EA. XYG3 leads to mean absolute deviations (MADs) of 0.057 and 0.080 eV for IEs and EAs, respectively, using the basis set of 6-311 + G (3df,2p). In comparison with some other functionals, MADs for IEs are 0.109 (B2PLYP), 0.119 (M06-2X), 0.159 (X3LYP), 0.161 (PBE), 0.162 (B3LYP), 0.165 (PBE0), 0.173 (TPSS), 0.200 (BLYP), and 0.215 eV (LC-BLYP). MADs for EAs are 0.090 (X3LYP), 0.090 (B2PLYP), 0.102 (PBE), 0.103 (M06-2X), 0.104 (TPSS), 0.105 (BLYP), 0.106 (B3LYP), 0.126 (LC-BLYP), and 0.128 eV (PBE0).