Is There Still Room for Parameter Free Double Hybrids? Performances of PBE0-DH and B2PLYP over Extended Benchmark Sets

Computed versus experimental energy barriers in solution: Influence of the type of the density functional approximation

Mechanistic investigations at the density functional theory level of organic and organometallic reactions in solution are now broadly accessible and routinely implemented to complement experimental investigations. The selection of an appropriate functional among the plethora of developed ones is the first challenge on the way to reliable energy barrier calculations. To provide guidelines for the choice of an initial and reliable computational level, the performances of commonly used non-empirical (PBE, PBE0, PBE0-DH) and empirical density functionals (BLYP, B3LYP, B2PLYP) were evaluated relative to experimental activation enthalpies. Most reactivity databases to assess density functional performances are primarily based on high level calculations, here a set of experimental activation enthalpies of organic and organometallic reactions performed in solution were selected from the literature. As a general trend, the non-empirical functionals outperform the empirical ones. The most accurate energy barriers are obtained with hybrid PBE0 and double-hybrid PBE0-DH density functionals, both providing similar performance. Regardless of the functional under consideration, the addition of the GD3-BJ empirical dispersion correction does not enhance the accuracy of computed energy barriers.

How Much Electron Donation Is There In Transition Metal Complexes? A Computational Study

The "dative" covalent interactions between metals and ligands in coordination compounds, i.e., metal-to-ligand and ligand-to-metal donation, are manifestations of electron delocalization and subject to errors in approximate calculations. This work addresses the extent of dative bonding/donation in a series of closed-shell transition metal complexes. Several Kohn-Sham density functionals, representing different "rungs" of approximations, along with post-Hartree-Fock methods are assessed in comparison to CCSD(T). Two widely used nonhybrid and global hybrid density functionals (B3LYP, PBE0) tend to produce notably too strong donation. Global hybrids with elevated fractions of exact exchange (40 to 50%) and the range-separated exchange functional CAM-B3LYP tend to perform better for the description of the donation. The performance of a double-hybrid functional is found to be quite satisfactory, correcting errors seen in MP2 calculations. A fast approximate coupled-cluster model (DLPNO-CCSD) also gives a reasonable description of the donation, with a tendency to underestimate its extent.

Electronic Energy and Local Property Errors at QTAIM Critical Points while Climbing Perdew's Ladder of Density-Functional Approximations

We investigate the relationships between electron-density and electronic-energy errors produced by modern exchange-correlation density-functional approximations belonging to all of the rungs of Perdew's ladder. To this aim, a panel of relevant (semi)local properties evaluated at critical points of the electron-density field (as defined within the framework of Bader's atoms-in-molecules theory) are computed on a large selection of molecular systems involved in thermodynamic, kinetic, and noncovalent interaction chemical databases using density functionals developed in a nonempirical and minimally and highly parametrized fashion. The comparison of their density- and energy-based performance, also discussed in terms of density-driven errors, casts light on the strengths and weaknesses of the most recent and efficient density-functional approximations.

Halogen Bonding Involving I2 and d8 Transition-Metal Pincer Complexes

We systematically investigated iodine–metal and iodine–iodine bonding in van Koten’s pincer complex and 19 modifications changing substituents and/or the transition metal with a PBE0–D3(BJ)/aug–cc–pVTZ/PP(M,I) model chemistry. As a novel tool for the quantitative assessment of the iodine–metal and iodine–iodine bond strength in these complexes we used the local mode analysis, originally introduced by Konkoli and Cremer, complemented with NBO and Bader’s QTAIM analyses. Our study reveals the major electronic effects in the catalytic activity of the M–I–I non-classical three-center bond of the pincer complex, which is involved in the oxidative addition of molecular iodine I2 to the metal center. According to our investigations the charge transfer from the metal to the σ* antibonding orbital of the I–I bond changes the 3c–4e character of the M–I–I three-center bond, which leads to weakening of the iodine I–I bond and strengthening of the metal–iodine M–I bond, facilitating in this way the oxidative addition of I2 to the metal. The charge transfer can be systematically modified by substitution at different places of the pincer complex and by different transition metals, changing the strength of both the M–I and the I2 bonds. We also modeled for the original pincer complex how solvents with different polarity influence the 3c–4e character of the M–I–I bond. Our results provide new guidelines for the design of pincer complexes with specific iodine–metal bond strengths and introduce the local vibrational mode analysis as an efficient tool to assess the bond strength in complexes.

Assessing challenging intra- and inter-molecular charge-transfer excitations energies with double-hybrid density functionals

We investigate the performance of a set of recently introduced range-separated double-hybrid functionals, namely ωB2-PLYP, ωB2GP-PLYP, RSX-0DH, and RSX-QIDH models for hard-to-calculate excitation energies. We compare with the parent (B2-PLYP, B2GP-PLYP, PBE0-DH, and PBE-QIDH) and other (DSD-PBEP86) double-hybrid models as well as with some of the most widely employed hybrid functionals (B3LYP, PBE0, M06-2X, and ωB97X). For this purpose, we select a number of medium-sized intra- and inter-molecular charge-transfer excitations, which are known to be challenging to calculate using time-dependent density-functional theory (TD-DFT) and for which accurate reference values are available. We assess whether the high accuracy shown by the newest double-hybrid models is also confirmed for those cases too. We find that asymptotically corrected double-hybrid models yield a superior performance, especially for the inter-molecular charge-transfer excitation energies, as compared to standard double-hybrid models. Overall, the PBE-QIDH and its corresponding range-separated RSX-QIDH functional are recommended for general-purpose TD-DFT applications, depending on whether long-range effects are expected to play a significant role.

Identification and Reactivity of s-cis,s-cis-Dihydroxycarbene, a New [CH2O2] Intermediate

We report the first preparation of the s-cis,s-cis conformer of dihydroxycarbene (1cc) by means of pyrolysis of oxalic acid, isolation of the lower-energy s-trans,s-trans (1tt) and s-cis,s-trans (1ct) product conformers at cryogenic temperatures in a N₂ matrix, and subsequent narrow-band near-infrared (NIR) laser excitation to give 1cc. Carbene 1cc converts quickly to 1ct via quantum-mechanical tunneling with an effective half-life of 22 min at 3 K. The potential energy surface features around 1 were pinpointed by convergent focal point analysis targeting the AE-CCSDT(Q)/CBS level of electronic structure theory. Computations of the tunneling kinetics confirm the time scale of the 1cc → 1ct rotamerization and suggest that direct 1cc → H₂ + CO₂ decomposition may also be a minor pathway. The intriguing latter possibility cannot be confirmed spectroscopically, but hints of it may be present in the measured kinetic profiles.

Range-separated hybrid and double-hybrid density functionals: A quest for the determination of the range-separation parameter

We recently derived a new and simple route to the determination of the range-separation parameter in range-separated exchange hybrid and double-hybrid density functionals by imposing an additional constraint to the exchange-correlation energy to recover the total energy of the hydrogen atom [Brémond et al., J. Chem. Phys. 15, 201102 (2019)]. Here, we thoroughly assess this choice by statistically comparing the derived values of the range-separation parameters to the ones obtained using the optimal tuning (OT) approach. We show that both approaches closely agree, thus, confirming the reliability of ours. We demonstrate that it provides very close performances in the computation of properties particularly prone to the one- and many-electron self-interaction errors (i.e., ionization potentials). Our approach arises as an alternative to the OT procedure, conserving the accuracy and efficiency of a standard Kohn-Sham approach to density-functional theory computation.

A Continuum from Halogen Bonds to Covalent Bonds: Where Do λ3 Iodanes Fit?

The intrinsic bonding nature of λ 3 -iodanes was investigated to determine where its hypervalent bonds fit along the spectrum between halogen bonding and covalent bonding. Density functional theory with an augmented Dunning valence triple zeta basis set ( ω B97X-D/aug-cc-pVTZ) coupled with vibrational spectroscopy was utilized to study a diverse set of 34 hypervalent iodine compounds. This level of theory was rationalized by comparing computational and experimental data for a small set of closely-related and well-studied iodine molecules and by a comparison with CCSD(T)/aug-cc-pVTZ results for a subset of the investigated iodine compounds. Axial bonds in λ 3 -iodanes fit between the three-center four-electron bond, as observed for the trihalide species IF 2 − and the covalent FI molecule. The equatorial bonds in λ 3 -iodanes are of a covalent nature. We explored how the equatorial ligand and axial substituents affect the chemical properties of λ 3 -iodanes by analyzing natural bond orbital charges, local vibrational modes, the covalent/electrostatic character, and the three-center four-electron bonding character. In summary, our results show for the first time that there is a smooth transition between halogen bonding → 3c–4e bonding in trihalides → 3c–4e bonding in hypervalent iodine compounds → covalent bonding, opening a manifold of new avenues for the design of hypervalent iodine compounds with specific properties.

Semi-empirical or non-empirical double-hybrid density functionals: which are more robust?

The development of non-empirical double-hybrid density functionals (DHDFs) is a very active research area with the number of approaches in this field having increased rapidly. At the same time, there is a lack of published work that provides a fair assessment and comparison between non-empirical and semi-empirical DHDFs on an equal footing. Herein, we close this gap and present a thorough analysis of both classes of DHDFs on the large GMTKN55 benchmark database for general main-group thermochemistry, kinetics, and noncovalent interactions [Goerigk et al., Phys. Chem. Chem. Phys., 2017, 19, 32184-32215]. In total, 115 variations of dispersion-corrected and -uncorrected DHDFs are tested, which will be condensed to an in-depth assessment of 31 methods: 19 non-empirical and 12 semi-empirical DHDFs. As such, our study represents the largest DHDF study ever conducted and can serve as an important benchmark informing method developers and users alike. Our results show that semi-empirical DHDFs are the most robust density functional approximations and more reliable and accurate than non-empirical ones. In fact, some non-empirical approaches are even outperformed by hybrid approaches or even dispersion-corrected and -uncorrected MP2 and SCS-MP2. SOS0-PBE0-2-D3(BJ) is the only exception and the only non-empirical DHDF that we can safely recommend for general applicability. However, it is still outperformed by six semi-empirical DHDFs, of which we would like to particularly recommend the following five: ωB97X-2-D3(BJ), DSD-BLYP-D3(BJ), DSD-PBEP86-D3(BJ), B2NC-PLYP-D3(BJ), and B2GPPLYP-D3(BJ). Our findings seriously question current trends in the field and they highlight that novel strategies have to be found in order to outperform the currently best density functional theory methods on the market. We hope that our study can function as an important cornerstone inspiring such a change of direction in the field.

Dissecting the accountability of parameterized and parameter-free single-hybrid and double-hybrid functionals for photophysical properties of TADF-based OLEDs

Organic light emitting diodes (OLEDs) based on thermally activated delayed fluorescence (TADF) emitters are an attractive category of materials that have witnessed a booming development in recent years. In the present contribution, we scrutinize the accountability of parameterized and parameter-free single-hybrid (SH) and double-hybrid (DH) functionals through the two formalisms, full time-dependent density functional theory (TD-DFT) and Tamm-Dancoff approximation (TDA), for the estimation of photophysical properties like absorption energy, emission energy, zero-zero transition energy, and singlet-triplet energy splitting of TADF molecules. According to our detailed analyses on the performance of SHs based on TD-DFT and TDA, the TDA-based parameter-free SH functionals, PBE0 and TPSS0, with one-third of exact-like exchange turned out to be the best performers in comparison to other functionals from various rungs to reproduce the experimental data of the benchmarked set. Such affordable SH approximations can thus be employed to predict and design the TADF molecules with low singlet-triplet energy gaps for OLED applications. From another perspective, considering this point that both the nonlocal exchange and correlation are essential for a more reliable description of large charge-transfer excited states, applicability of the functionals incorporating these terms, namely, parameterized and parameter-free DHs, has also been evaluated. Perusing the role of exact-like exchange, perturbative-like correlation, solvent effects, and other related factors, we find that the parameterized functionals B2π-PLYP and B2GP-PLYP and the parameter-free models PBE-CIDH and PBE-QIDH have respectable performance with respect to others. Lastly, besides the recommendation of reliable computational protocols for the purpose, hopefully this study can pave the way toward further developments of other SHs and DHs for theoretical explorations in the field of OLEDs technology.

Partnering dispersion corrections with modern parameter-free double-hybrid density functionals

The PBE-QIDH and SOS1-PBE-QIDH double-hybrid density functionals are merged with a pair of dispersion corrections, namely the pairwise additive D3(BJ) and the non-local correlation functional VV10, leading to the corresponding dispersion-corrected models.

A look at the density functional theory zoo with the advanced GMTKN55 database for general main group thermochemistry, kinetics and noncovalent interactions

We present the updated and extended GMTKN55 benchmark database for more accurate and extensive energetic evaluation of density functionals and other electronic structure methods with detailed guidelines for method users.

Nonempirical Double-Hybrid Functionals: An Effective Tool for Chemists

Density functional theory (DFT) emerged in the last two decades as the most reliable tool for the description and prediction of properties of molecular systems and extended materials, coupling in an unprecedented way high accuracy and reasonable computational cost. This success rests also on the development of more and more performing density functional approximations (DFAs). Indeed, the Achilles' heel of DFT is represented by the exchange-correlation contribution to the total energy, which, being unknown, must be approximated. Since the beginning of the 1990s, global hybrids (GH) functionals, where an explicit dependence of the exchange-correlation energy on occupied Kohn-Sham orbitals is introduced thanks to a fraction of Hartree-Fock-like exchange, imposed themselves as the most reliable DFAs for chemical applications. However, if these functionals normally provide results of sufficient accuracy for most of the cases analyzed, some properties, such as thermochemistry or dispersive interactions, can still be significantly improved. A possible way out is represented by the inclusion, into the exchange-correlation functional, of an explicit dependence on virtual Kohn-Sham orbitals via perturbation theory. This leads to a new class of functionals, called double-hybrids (DHs). In this Account, we describe our nonempirical approach to DHs, which, following the line traced by the Perdew-Burke-Ernzerhof approach, allows for the definition of a GH (PBE0) and a DH (QIDH) model. In such a way, a whole family of nonempirical functionals, spanning on the highest rungs of the Perdew's quality scale, is now available and competitive with other-more empirical-DFAs. Discussion of selected cases, ranging from thermochemistry and reactions to weak interactions and excitation energies, not only show the large range of applicability of nonempirical DFAs, but also underline how increasing the number of theoretical constraints parallels with an improvement of the DFA's numerical performances. This fact further consolidates the strong theoretical framework of nonempirical DFAs. Finally, even if nonempirical DH approaches are still computationally expensive, relying on the fact that they can benefit of all technical enhancements developed for speeding up post-Hartree-Fock methods, there is substantial hope for their near future routine application to the description and prediction of complex chemical systems and reactions.

Benchmarking Density Functionals on Structural Parameters of Small-/Medium-Sized Organic Molecules

In this Letter we report the error analysis of 59 exchange-correlation functionals in evaluating the structural parameters of small- and medium-sized organic molecules. From this analysis, recently developed double hybrids, such as xDH-PBE0, emerge as the most reliable methods, while global hybrids confirm their robustness in reproducing molecular structures. Notably the M06-L density functional is the only semilocal method reaching an accuracy comparable to hybrids'. A comparison with errors obtained on energetic databases (including thermochemistry, reaction barriers, and interaction energies) indicate that most of the functionals have a coherent behavior, showing low (or high) deviations on both energy and structure data sets. Only a few of them are more prone toward one of these two properties.

Systematic Improvement of Density Functionals through Parameter-Free Hybridization Schemes

It is suggested here that the ultimate accuracy of DFT methods arises from the type of hybridization scheme followed. This idea can be cast into a mathematical formulation utilizing an integrand connecting the noninteracting and the interacting particle system. We consider two previously developed models for it, dubbed as HYB0 and QIDH, and assess a large number of exchange-correlation functionals against the AE6, G2/148, and S22 reference data sets. An interesting consequence of these hybridization schemes is that the error bars, including the standard deviation, are found to markedly decrease with respect to the density-based (nonhybrid) case. This improvement is substantially better than variations due to the underlying density functional used. We thus finally hypothesize about the universal character of the HYB0 and QIDH models.

Assessment of DFT methods for studying acid gas capture by ionic liquids

For the first time, this work reports an analysis of the performance of Density Functional methods for studying acid gas capture (CO₂ and SO₂) by ionic liquids (ILs).

Relative energies of water nanoclusters (H2O)20: comparison of empirical and nonempirical double-hybrids with generalized energy-based fragmentation approach

The applicability of recently developed parameterized and parameter-free double-hybrids for predicting the relative energies of water nanoclusters has been examined.

Double-hybrid density functionals free of dispersion and counterpoise corrections for non-covalent interactions

We have optimized two double-hybrid density functionals (DHDFs) within the frameworks of B2PLYP and mPW2-PLYP against the S22B database. These two functionals are denoted as B2NC-PLYP and mPW2NC-PLYP, where "NC" represents noncovalent interaction. The DHDFs of B2NC-PLYP and mPW2NC-PLYP are optimized free of dispersion and counterpoise corrections with triple-ζ quality basis sets. Combined with the aug-cc-pVTZ basis set, these two functionals are further assessed with the S66 database. According to our computations, both the B2NC-PLYP and mPW2NC-PLYP functionals seem to be competent for investigating noncovalent interactions. Note that the triple-ζ quality basis sets with adequate polarization and diffuse functions should be employed for practical applications. However, different exchange and correlation functionals may be selected and/or modified to reduce the amount of the Fock-exchange in the future.

Computational design of in vivo biomarkers

Fluorescent semiconductor nanocrystals (or quantum dots) are very promising agents for bioimaging applications because their optical properties are superior compared to those of conventional organic dyes. However, not all the properties of these quantum dots suit the stringent criteria of in vivo applications, i.e. their employment in living organisms that might be of importance in therapy and medicine. In our review, we first summarize the properties of an 'ideal' biomarker needed for in vivo applications. Despite recent efforts, no such hand-made fluorescent quantum dot exists that may be considered as 'ideal' in this respect. We propose that ab initio atomistic simulations with predictive power can be used to design 'ideal' in vivo fluorescent semiconductor nanoparticles. We briefly review such ab initio methods that can be applied to calculate the electronic and optical properties of very small nanocrystals, with extra emphasis on density functional theory (DFT) and time-dependent DFT which are the most suitable approaches for the description of these systems. Finally, we present our recent results on this topic where we investigated the applicability of nanodiamonds and silicon carbide nanocrystals for in vivo bioimaging.