DSD‐PBEP86‐NL and DOD‐PBEP86‐NL functionals for noncovalent interactions: Basis set effects and tentative applications to large noncovalent systems

Basis set effects on the DSD‐PBEP86‐NL and DOD‐PBEP86‐NL functionals for noncovalent interactions have been extensively studied in this work. The cc‐pVXZ (X = D, T, Q, 5, 6) and augmented aug‐cc‐pVXZ (X = D, T, Q) basis sets are systematically tested without counterpoise (CP) corrections against the well‐known S66 database. Additionally, the basis sets of def2‐TZVPP and def2‐TZVPPD are also examined. Based on our computations, the performances of the aug‐cc‐pVQZ, cc‐pV5Z, and cc‐pV6Z basis sets are very approximate to those obtained with the def2‐QZVP basis set for both the DSD‐PBEP86‐NL and DOD‐PBEP86‐NL functionals. Note that the short‐range attenuation parameters for these two functionals were directly optimized using the def2‐QZVP basis set without CP corrections against the S66 database. Generally speaking, the cc‐pVXZ (X = D, T, Q), aug‐cc‐pVXZ (X = D, T, Q), def2‐TZVPP, and def2‐TZVPPD basis sets favor half CP correction for these two functionals. Nevertheless, the aug‐cc‐pVQZ basis set already performs well without any CP correction, especially for the DOD‐PBEP86‐NL functional. With respect to accuracy and computational cost, the cc‐pVTZ and def2‐TZVPP basis sets with half CP corrections are recommended for these two functionals to evaluate interaction energies of large noncovalent complexes.     

Order of stabilities in water nanoclusters: Insight from some recent double‐hybrid functionals

In the present work, the applicability of some of the recently proposed and modern double‐hybrid (DH) models and other density functional theory (DFT) approximations has been analyzed for a difficult test, the order of stability in low‐energy isomers of water nanoclusters. In particular, we aim to systematically investigate for these functionals the role played by several factors such as dispersion correction, integrand functions upon which the DHs are based, and different spin scaling for the perturbative term in DH calculations of the relative energies for various isomers of water nanoclusters (H₂O)₂₀. From the obtained results, the superior performance of DHs with respect to the functionals from previous rungs is confirmed. It is shown that the dispersion corrected DHs perform better than noncorrected counterparts. Plus, the DH models based on cubic integrand (CI) and quadratic integrand (QI) functions are nearly equivalent in performance. We also find that using only contributions of electron pairs with opposite spin for the perturbative correlation part through scaled opposite spin scheme does not represent a significant improvement on accuracy of the results. Putting all the results together, the dispersion corrected parameterized DHs and parameter‐free DH models involving CI and QI functions outperform other approximations for relative energies of water 20‐mers. Altogether, predicting the correct order of the stability in water nanoclusters may be considered as another Achilles' heel in DFT calculations, although more analyses in this context are still needed. © 2016 Wiley Periodicals, Inc.     

Broken symmetry density functional study of a mixed‐valence unsymmetrical dinuclear iron complex

Density functional theory (DFT) calculations with different exchange‐correlation functionals were performed for a mixed valence Fe(II)/Fe(III) binuclear complex with μ‐methoxo and two μ‐carboxylate bridging ligands, (1) with geometry optimizations being performed for all possible spin multiplicities (M_S = 2, 4, 6, 8, and 10). Within the exchange‐correlation functionals studied, only the hybrid GGA functionals B3P and B3LYP and also the pure GGA functional RPBE, predicts the geometry with high spin (S = 9/2) to be more stable than the geometry with low spin state (S = 1/2) by 20 kcal/mol, in agreement with the experimental findings. These functionals also predict the same stability order for the different spin states, being M_S = 10>8>6>2>4. The meta‐GGA functionals TPSS and TPSSh and also the pure GGA functionals BLYP and BP86 predict different stability orders. The computed average EPR g‐tensor, g_av, of 2.03, at the B3LYP level, is in good agreement with the experimental findings. Heisenberg exchange coupling constants, J, were calculated within the broken‐symmetry formalism, at the B3LYP level, showing that the two iron centers are antiferromagnetic coupling, with a very weak coupling constant of about ?7 cm^?1, in good agreement with the experimental value. Additionally, the effect of using different multiplicities of the reference geometries on the computed J value is discussed. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Highly accurate equilibrium structure of the C2h symmetric N1‐to‐O2 hydrogen‐bonded uracil‐dimer

The highly accurate ab initio equilibrium geometry of the hydrogen‐bonded uracil dimer is derived using a composite geometry extrapolation scheme based on all‐electron, complete basis set extrapolated Møller–Plesset perturbation theory using the jun‐pwCV[T,Q]Z basis sets combined with a valence CCSD(T)/cc‐pVTZ high‐level correction. Geometrical changes on dimerization are discussed and the performance of the several density functional approximations (among others SCAN, ωB97M‐V, DSD‐PBEP86‐D3(BJ), and DSD‐PBEP86‐NL) is evaluated. Orbital‐optimized MP2.5 is discussed as a reduced‐cost alternative to the CCSD(T) gradient in the composite scheme. A new reference interaction energy is calculated with explicitly correlated F12‐CCSD(T).     

Performance of the OP correlation functional in relation to its formulation: Influence of the exchange component and the effect of incorporating same‐spin correlations

In the present study, we have investigated two significant features of the OP correlation functional, namely the incorporation of the exchange functional into itself, and the inclusion of only opposite‐spin (OS) effects. To explore the latter feature, we have compared OP with B95 and a new functional introduced in the present study – the OPB method that combines OP with the same‐spin (SS) component of B95. In general, we find that B95 and OPB perform comparably. Our comparisons of the various DFT procedures suggest that the incorporation of a meta‐GGA (e.g., TPSS) into OP and OPB does not necessarily lead to a chemically more accurate procedure than the use of a related GGA (e.g., PBE). An important finding is the more notable (and somewhat more consistent) improvement in performance with the incorporation of SS correlation, particularly for longer‐range chemical properties. Nonetheless, on average across our test sets of over 800 systems, the difference between the performances of OP versus B95 or OPB is not exceedingly large. By drawing a parallel between these DFT methods and the wavefunction scaled‐MP2‐type methods, we reason that one can further develop the OP functional, and perhaps a wider range of correlation functionals by combining it with the technique of range separation. © 2016 Wiley Periodicals, Inc.     

Parallel implementation of γ‐point pseudopotential plane‐wave DFT with exact exchange

Semi‐local functionals commonly used in density functional theory (DFT) studies of solids usually fail to reproduce localized states such as trapped holes, polarons, excitons, and solitons. This failure is ascribed to self‐interaction which creates a Coulomb barrier to localization. Pragmatic approaches in which the exchange correlation functionals are augmented with small amount of exact exchange (hybrid‐DFT, e.g., B3LYP and PBE0) have shown to promise in rectifying this type of failure, as well as producing more accurate band gaps and reaction barriers. The evaluation of exact exchange is challenging for large, solid state systems with periodic boundary conditions, especially when plane‐wave basis sets are used. We have developed parallel algorithms for implementing exact exchange into pseudopotential plane‐wave DFT program and we have implemented them in the NWChem program package. The technique developed can readily be employed in Γ‐point plane‐wave DFT programs. Furthermore, atomic forces and stresses are straightforward to implement, making it applicable to both confined and extended systems, as well as to Car‐Parrinello ab initio molecular dynamic simulations. This method has been applied to several systems for which conventional DFT methods do not work well, including calculations for band gaps in oxides and the electronic structure of a charge trapped state in the Fe(II) containing mica, annite. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Applications of density functional theory methods in millimeter‐wave spectroscopy

Density functional theory (DFT), using the most common functionals, and ab initio quantum chemistry methods are used to calculate the rotational constants and dipole moments of the astrophysically important molecules HCN, CH₃CN, CH₃CNH⁺, HCCCN, and HCCNC. As far as millimeter‐wave spectroscopy is of interest the DFT methods performed well with most functionals, giving results within ±1% of experiments for rotational constants and ±3% for dipole moments. Analyzing the results obtained with all theoretical models, it may be concluded that the Becke's three‐parameter exchange functional and the gradient‐corrected functional of Lee, Yang, and Paar (B3LYP) and Becke's three‐parameter functional with Perdew–Wang correlational functional [B3PW91/6‐31G(d, p)] give the best performances. A detailed analysis of the electron correlation effects shows that HCCCN is more stable than is HCCNC, by 1.16 eV, with important contribution arising from triple excitations. This result is also compared with those obtained with DFT methods. Despite occasional difficulties, DFT with the currently available functionals are of great utility in quickly assessing spectroscopic parameters of astrophysical interest. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem, 2003     

Determining the role of the underlying orbital‐dependence of PBE0‐DH and PBE‐QIDH double‐hybrid density functionals

We study the orbital‐dependence of three (parameter‐free) double‐hybrid density functionals, namely the PBE0‐DH, the PBE‐QIDH models, and the SOS1‐PBE‐QIDH spin‐opposite‐scaled variant of the latter. To do it, we feed all their energy terms with different sets of orbitals obtained previously from self‐consistent density functional theory calculations using several exchange‐correlation functionals (e.g., PBE, PBE0, PBEH&H), or directly with HF‐PBE orbitals, to see their effect on selected datasets for atomization and reaction energies, the latter proned to marked self‐interaction errors. We find that the PBE‐QIDH double‐hybrid model shows a great consistency, as the best results are always obtained for the set of orbitals corresponding to its hybrid scheme, which prompts us to recommend this model without any other fitting or reparameterization. © 2017 Wiley Periodicals, Inc.     

DFT performance for the hole transfer parameters in DNA π stacks

Recently, we showed that unoccupied Kohn‐Sham (KS) orbitals stemming from DFT calculations of a neutral system can be used to derive accurate estimates of the free energy and electronic couplings for excess electron transfer in DNA (Félix and Voityuk, J Phys Chem A 2008, 112, 9043). In this article, we consider the propagation of radical cation states (hole transfer) through DNA π‐stacks and compare the performance of different exchange‐correlation functionals to estimate the hole transfer (HT) parameters. Two different approaches are used: (1) calculations that use occupied KS orbitals of neutral π stacks of nucleobases, and (2) the time‐dependent DFT method which is applied to the radical cation states of these stacks. Comparison of the calculated parameters with the reference data suggests that the best results are provided by the KS scheme with hybrid functionals (B3LYP, PBE0, and BH&HLYP). The TD DFT approach gives significantly less accurate values of the HT parameters. In agreement with high‐level ab initio results, the KS scheme predicts that the hole in π stacks is confined to a single nucleobase; in contrast, the spin‐unrestricted DFT method considerably overestimates the hole delocalization in the radical cations. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

PW86–PW91 density functional calculation of vertical ionization potentials: Some implications for present‐day functionals

A total of 181 vertical ionization potentials (VIPs) of 41 molecules were calculated by density functional theory (DFT) employing the Perdew–Wang 1986 (PW86) exchange and Perdew–Wang 1991 (PW91) correlation functionals and using the aug‐cc‐pV5Z basis and experimental ground‐state geometries. The overall average absolute deviation (AAD) from experiment was found to be 0.55 eV and only 0.31 eV for linear molecules but 0.86 eV for nonplanar molecules. A number of VIPs were in error by over 2 eV. In particular, DFT performed most poorly when ionization was from an orbital with highly varying density gradients (which arise from the orbital's shape or compactness or through its density being distributed over a number of atoms). Indications are that many or all present‐day functionals suffer from the same failings. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 34–52, 2001     

How reliable is DFT in predicting relative energies of polycyclic aromatic hydrocarbon isomers? comparison of functionals from different rungs of jacob's ladder

Density functional theory (DFT) is the only quantum‐chemical avenue for calculating thermochemical/kinetic properties of large polycyclic aromatic hydrocarbons (PAHs) such as graphene nanoflakes. Using CCSD(T)/CBS PAH isomerization energies, we find that all generalized gradient approximation (GGA) and meta GGA DFT functionals have severe difficulties in describing isomerization energies in PAHs. The poor performance of these functionals is demonstrated by the following root‐mean‐square deviations (RMSDs) obtained for a database of C₁₄H₁₀ and C₁₈H₁₂ isomerization energies. The RMSDs for the GGAs range between 6.0 (BP86‐D3) and 23.0 (SOGGA11) and for the meta GGAs they range between 3.5 (MN12‐L) and 11.9 (τ‐HCTH) kJ mol⁻¹. These functionals (including the dispersion‐corrected methods) systematically and significantly underestimate the isomerization energies. A consequence of this behavior is that they all predict that chrysene (rather than triphenylene) is the most stable C₁₈H₁₂ isomer. A general improvement in performance is observed along the rungs of Jacob's Ladder; however, only a handful of functionals from rung four give good performance for PAH isomerization energies. These include functionals with high percentages (40–50%) of exact Hartree–Fock exchange such as the hybrid GGA SOGGA11‐X (RMSD = 1.7 kJ mol⁻¹) and the hybrid‐meta GGA BMK (RMSD = 1.3 kJ mol⁻¹). Alternatively, the inclusion of lower percentages (20–30%) of exact exchange in conjunction with an empirical dispersion correction results in good performance. For example, the hybrid GGA PBE0‐D3 attains an RMSD of 1.5 kJ mol⁻¹, and the hybrid‐meta GGAs PW6B95‐D3 and B1B95‐D3 result in RMSDs below 1 kJ mol⁻¹. © 2016 Wiley Periodicals, Inc.     

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.     

Photomagnetic and nonlinear optical properties in cis‐trans green fluoroprotein chromophore coupled Bis‐imino nitroxide diradicals

The current study extends an earlier investigation (Bhattacharya, et al., Phys. Chem. Chem. Phys. 2012, 14, 6905) to further explore various photomagnetic and optical properties of bis‐imino nitroxide, that is, (IN)₂‐based green fluorescent protein (GFP) chromophore coupled diradicals revealing new significant features. The conversion mechanisms of selected trans‐isomers into their corresponding cis‐conformers are discussed in detailed using a number of recently‐developed density functional theory (DFT) functionals based on the Minnesota suite of DFT‐models as well as using some other DFT functionals developed earlier. To provide a more in‐depth analysis of variations in magnetic properties as trans‐conformers (singlet ground‐state) convert into their cis‐analogues (triplet ground‐state), the changes in exchange magnetic coupling constants J are compared with the variation of the selected aromaticity indices. The aromaticity indices include the nuclear independent chemical shift [NICS(0)] values calculated at the center of ring structures and the harmonic oscillator model of aromaticity. Furthermore, the investigation of static nonlinear optical response properties in the (IN)₂‐based GFP chromophore coupled diradicals reveal unusually large static first hyperpolarizabilities for these systems which is highly significant for practical applications in optics and optoelectronics. © 2015 Wiley Periodicals, Inc.     

Toward the complete range separation of non‐hybrid exchange–correlation functional

In this study, we use a very simple scheme to achieve range separation of a total exchange–correlation functional. We have utilized this methodology to combine a short‐range pure density functional theory (DFT) functional with a corresponding long‐range pure DFT, leading to a “Range‐separated eXchange–Correlation” (RXC) scheme. By examining the performance of a range of standard exchange–correlation functionals for prototypical short‐ and long‐range properties, we have chosen B‐LYP as the short‐range functional and PBE‐B95 as the long‐range counterpart. The results of our testing using a more diverse range of data sets show that, for properties that we deem to be short‐range in nature, the performance of this prescribed RXC‐DFT protocol does resemble that of B‐LYP in most cases, and vice versa. Thus, this RXC‐DFT protocol already provides meaningful numerical results. Furthermore, we envisage that the general RXC scheme can be easily implemented in computational chemistry software packages. This study paves a way for further refinement of such a range‐separation technique for the development of better performing DFT procedures. © 2015 Wiley Periodicals, Inc.     

Correlation functional in screened‐exchange density functional theory procedures

In the present study, we have explored several prospects for the further development of screened‐exchange density functional theory (SX‐DFT) procedures. Using the performance of HSE06 as our measure, we find that the use of alternative correlation functionals (as oppose to PBEc in HSE06) also yields adequate results for a diverse set of thermochemical properties. We have further examined the performance of new SX‐DFT procedures (termed HSEB‐type methods) that comprise the HSEx exchange and a (near‐optimal) reparametrized B97c (c _OS,0 = c _SS,0 = 1, c _OS,1 = −1.5, c _OS,2 = −0.644, c _SS,1 = −0.5, and c _SS,2 = 1.10) correlation functionals. The different variants of HSEB all perform comparably to or slightly better than the original HSE‐type procedures. These results, together with our fundamental analysis of correlation functionals, point toward various directions for advancing SX‐DFT methods. © 2017 Wiley Periodicals, Inc.     

W4‐17: A diverse and high‐confidence dataset of atomization energies for benchmarking high‐level electronic structure methods

Atomization reactions are among the most challenging tests for electronic structure methods. We use the first‐principles Weizmann‐4 (W4) computational thermochemistry protocol to generate the W4‐17 dataset of 200 total atomization energies (TAEs) with 3σ confidence intervals of 1 kJ mol⁻¹. W4‐17 is an extension of the earlier W4‐11 dataset; it includes first‐ and second‐row molecules and radicals with up to eight non‐hydrogen atoms. These cover a broad spectrum of bonding situations and multireference character, and as such are an excellent benchmark for the parameterization and validation of highly accurate ab initio methods (e.g., CCSD(T) composite procedures) and double‐hybrid density functional theory (DHDFT) methods. The W4‐17 dataset contains two subsets (i) a non‐multireference subset of 183 systems characterized by dynamical or moderate nondynamical correlation effects (denoted W4‐17‐nonMR) and (ii) a highly multireference subset of 17 systems (W4‐17‐MR). We use these databases to evaluate the performance of a wide range of CCSD(T) composite procedures (e.g., G4, G4(MP2), G4(MP2)‐6X, ROG4(MP2)‐6X, CBS‐QB3, ROCBS‐QB3, CBS‐APNO, ccCA‐PS3, W1, W2, W1‐F12, W2‐F12, W1X‐1, and W2X) and DHDFT methods (e.g., B2‐PLYP, B2GP‐PLYP, B2K‐PLYP, DSD‐BLYP, DSD‐PBEP86, PWPB95, ωB97X‐2(LP), and ωB97X‐2(TQZ)). © 2017 Wiley Periodicals, Inc.     

Electronic spectra and intersystem spin‐orbit coupling in 1,2‐ and 1,3‐squaraines

The main photophysical properties of a series of recently synthetized 1,2‐ and 1,3‐squaraines, including absorption electronic spectra, singlet‐triplet energy gaps, and spin‐orbit matrix elements, have been investigated by means of density functional theory (DFT) and time‐dependent DFT approaches. A benchmark of three exchange‐correlation functionals has been performed in six different solvent environments. The investigated 1,2 squaraines have been found to possess two excited triplet states (T₁ and T₂) that lie below the energy of the excited singlet one (S₁). The radiationless intersystem spin crossing efficiency is thus enhanced in both the studied systems and both the transitions could contribute to the excited singlet oxygen production. Moreover, they have a singlet‐triplet energy gap higher than that required to generate the cytotoxic singlet oxygen species. According to our data, these compounds could be used in photodynamic therapy applications that do not require high tissue penetration. © 2014 Wiley Periodicals, Inc.     

Investigation of the use of density functionals in second‐ and third‐row transition metal dimer calculations

We explore the use of density functionals in calculating the equilibrium distances, dissociation energies, and harmonic vibrational frequencies of the homonuclear diatomics of the second‐row transition metals, platinum, and gold. The outermost s–d interconfigurational energies (ICEs) and the outermost s and d ionization potentials (IPs) were also calculated for the second‐ and third‐row transition metal atoms. Compared with the first‐row transition metal dimer calculations (J Chem Phys 2000, 112, 545–553), the binding energies calculated using the combination of the Becke 1988 exchange and the one‐parameter progressive correlation (BOP) functional and Becke's three‐parameter hybrid (B3LYP) functional are in better agreement with the experiment. However, the pure BOP functional still gives the deep and narrow dissociation potential wells for the electron configurations containing high‐angular‐momentum open‐shell orbitals. Analysis of the s–d ICEs and the s and d IPs suggests that the overestimation may be due to the insufficient long‐range interaction between the outermost s and d orbitals in the exchange functional. The hybrid B3LYP functional seems to partly solve this problem for many systems by the incorporation of the Hartree–Fock exchange integral. However, this still leads to an erroneous energy gap between the configurations of fairly different spin multiplicity, probably because of the unbalance of exchange and correlation contributions. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1995–2009, 2001