Axial-equatorial equilibrium in substituted cyclohexanes: a DFT perspective on a small but complex problem.

In Chemistry, complexity is not necessarily associated to large systems, as illustrated by the textbook example of axial-equatorial equilibrium in mono-substituted cyclohexanes. The difficulty in modelling such a simple isomerization is related to the need for reproducing the delicate balance between two forces, with opposite effects, namely the attractive London dispersion and the repulsive steric interactions. Such balance is a stimulating challenge for density-functional approximations and it is systematically explored here by considering 20 mono-substituted cyclohexanes. In comparison to highly accurate CCSD(T) reference calculations, their axial-equatorial equilibrium is studied with a large set of 48 exchange-correlation approximations, spanning from semilocal to hybrid to more recent double hybrid functionals. This dataset, called SAV20 (as Steric A-values for 20 molecules), allows to highlight the difficulties encountered by common and more original DFT approaches, including those corrected for dispersion with empirical potentials, the 6-31G*-ACP model, and our cost-effective PBE-QIDH/DH-SVPD protocol, in modeling these challenging interactions. Interestingly, the performance of the approaches considered in this contribution on the SAV20 dataset does not correlate with that obtained with other more standard datasets, such as S66, IDISP or NC15, thus indicating that SAV20 covers physicochemical features not already considered in previous noncovalent interaction benchmarks.

Computational Design of Multiple Resonance-Type BN Molecules for Inverted Singlet and Triplet Excited States.

A computational design of linearly extended multiple resonance (MR)-type BN molecules based on DABNA-1 is proposed herein in the quest to find potential candidates that exhibit a negative singlet-triplet gap (ΔEST) and a large oscillator strength value. The impact of a proper account of the electron correlation in the lowest singlet and triplet excited states is systematically investigated by using double-hybrid functionals within the TD-DFT framework, as well as wavefunction-based methods (EOM-CCSD and SCS-CC2), since this contribution plays an essential role in driving the magnitude of the ΔEST in MR-TADF and inverted singlet-triplet gap compounds. Our results point out a gradual reduction of the ΔEST gap with respect to the increasing sum of the number of B and N atoms, reaching negative ΔEST values for some molecules as a function of their size. The double-hybrid functionals reproduce the gap with only slight deviation compared to available experimental data for DABNA-1, ν-DABNA, and mDBCz and nicely agree with high-level quantum mechanical methods (e.g., EOM-CCSD and SCS-CC2). Larger oscillator strengths are found compared to the azaphenalene-type molecules, also exhibiting the inversion of their singlet and triplet excited states. We hope this study can serve as a motivation for further design of the molecules showing negative ΔEST based on boron- and nitrogen-doped polyaromatic hydrocarbons.

Benchmarking DFT Functionals for Excited-State Calculations of Donor-Acceptor TADF Emitters: Insights on the Key Parameters Determining Reverse Inter-System Crossing.

The importance of intermediate triplet states and the nature of excited states has gained interest in recent years for the thermally activated delayed fluorescence (TADF) mechanism. It is widely accepted that simple conversion between charge transfer (CT) triplet and singlet excited states is too crude, and a more complex route involving higher-lying locally excited triplet excited states has to be invoked to witness the magnitude of the rate of reverse inter-system crossing (RISC) rates. The increased complexity has challenged the reliability of computational methods to accurately predict the relative energy between excited states as well as their nature. Here, we compare the results of widely used density functional theory (DFT) functionals, CAM-B3LYP, LC-ωPBE, LC-ω*PBE, LC-ω*HPBE, B3LYP, PBE0, and M06-2X, against a wavefunction-based reference method, Spin-Component Scaling second-order approximate Coupled Cluster (SCS-CC2), in 14 known TADF emitters possessing a diversity of chemical structures. Overall, the use of the Tamm-Dancoff Approximation (TDA) together with CAM-B3LYP, M06-2X, and the two ω-tuned range-separated functionals LC-ω*PBE and LC-ω*HPBE demonstrated the best agreement with SCS-CC2 calculations in predicting the absolute energy of the singlet S1, and triplet T1 and T2 excited states and their energy differences. However, consistently across the series and irrespective of the functional or the use of TDA, the nature of T1 and T2 is not as accurately captured as compared to S1. We also investigated the impact of the optimization of S1 and T1 excited states on ΔEST and the nature of these states for three different functionals (PBE0, CAM-B3LYP, and M06-2X). We observed large changes in ΔEST using CAM-B3LYP and PBE0 functionals associated with a large stabilization of T1 with CAM-B3LYP and a large stabilization of S1 with PBE0, while ΔEST is much less affected considering the M06-2X functional. The nature of the S1 state barely evolves after geometry optimization essentially because this state is CT by nature for the three functionals tested. However, the prediction of the T1 nature is more problematic since these functionals for some compounds interpret the nature of T1 very differently. SCS-CC2 calculations on top of the TDA-DFT optimized geometries lead to a large variation in terms of ΔEST and the excited-state nature depending on the chosen functionals, further stressing the large dependence of the excited-state features on the excited-state geometries. The presented work highlights that despite good agreement of energies, the description of the exact nature of the triplet states should be undertaken with caution.

SOS1-RSX-QIDH: A spin-opposite-scaled range-separated-exchange quadratic-integrand double-hybrid density functional.

We develop and validate the SOS1-RSX-QIDH density functional, a one-parameter spin-opposite-scaled variant of the range-separated-exchange quadratic-integrand double-hybrid (RSX-QIDH) model. By entering into the family of spin-biased double hybrids, this new density functional benefits from an improved computational scaling that rivals with the one of hybrids, still conserving the accuracy of its RSX-QIDH version. As part of the latter family, this density functional is well-adapted to treat molecular systems that are particularly prone to self-interaction errors in their ground and excited states. In particular, we show that the SOS1-RSX-QIDH model is a good compromise to treat ground-state problems dealing with kinetics and has a real added value when applied to the evaluation of the excited-state properties of equilibrium and out-of-equilibrium molecular complexes. Even if spin-biased double hybrids are recognized to strongly underestimate noncovalent interactions, we notice and recommend coupling SOS1-RSX-QIDH with a nonlocal van der Waals potential, a combination that is here proved to compete with the best density-functional approximations currently in use.

Design Strategies for Diradical Boron/Nitrogen Doped Carbon-Based Materials.

New heterocyclic diradicaloids based on boron and nitrogen-doped polycyclic systems with open-shell ground-states are obtained via concomitant structural and quinoidal extensions, thus allowing to merge the best of both design strategies. A combination of experimental characterization and theoretical calculations have helped disclose their electronic structure, as well as rationalize their associated magnetic and photophysical properties, spanning the chemical space of available molecular templates for cutting-edge applications in organic electronics and spintronics.

Double Hybrids and Noncovalent Interactions: How Far Can We Go?

The accurate evaluation of weak noncovalent interactions in large, that is those containing up to thousand atoms, molecular systems represents a difficult challenge for any quantum chemical method. Indeed, some approximations are often introduced to render affordable these calculations. Here, we consider the PBE-QIDH/DH-SVPD protocol, combining a nonempirical double hybrid functional (PBE-QIDH) with a small basis set (DH-SVPD) tailored for noncovalent interactions with a double aim: (i) explore the robustness and accuracy of this protocol with respect to other Density Functional Approximations; (ii) illustrate how its performances are affected by the computational parameters underlying the calculation of the exact exchange and the Coulomb contribution, as well as the perturbative term. To this end, we consider three data sets, namely S66, L7, and CiM13, incorporating molecules of increasing size. On the bright side, our results suggest that the PBE-QIDH/DH-SVPD protocol is particularly accurate for large systems such as those contained in the CiM13 set (up to more than 1000 atoms and 14 000 basis functions), for which the DLPNO approximation leads to a significant speed-up for the evaluation of the perturbative correlation term. However, our analysis also points out the limit of this statistical exercise, when the quality of the reference data cannot be easily assessed, due to the size of the molecular complexes involved, and when the number of molecules is limited.

Tackling an accurate description of molecular reactivity with double-hybrid density functionals.

In this Communication, we assess a panel of 18 double-hybrid density functionals for the modeling of the thermochemical and kinetic properties of an extended dataset of 449 organic chemistry reactions belonging to the BH9 database. We show that most of DHs provide a statistically robust performance to model barrier height and reaction energies in reaching the "chemical accuracy." In particular, we show that nonempirical DHs, such as PBE0-DH and PBE-QIDH, or minimally parameterized alternatives, such as ωB2PLYP and B2K-PLYP, succeed to accurately model both properties in a balanced fashion. We demonstrate, however, that parameterized approaches, such as ωB97X-2 or DSD-like DHs, are more biased to only one of both properties.

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.

Beyond Chemical Accuracy for Alkane Thermochemistry: The DHthermo Approach.

The so-called protobranching phenomenon, that is the greater stability of branched alkanes with respect to their linear isomers, represents an interesting challenge for approaches based on density functional theory (DFT), since it requires a balanced description of several electronic effects, including (intramolecular) dispersion forces. Here, we investigate this problem using a protocol recently developed based on double-hybrid functionals and a small basis set, DH-SVPD, suited for noncovalent interactions. The energies of bond separation reactions (BSR), defined on the basis of an isodesmic principle, are taken as reference properties for the evaluation of 15 DFT approaches. The obtained results show that error lower than the so-called "chemical accuracy" (<1.0 kcal/mol) can be obtained by the proposed protocol on both relative reaction energies and enthalpies. These results are then verified on the standard BSR36 data set and support the proposition of our computational protocol, named DHthermo, where any DH functional, such as PBE-QIDH or B2PLYP, provides accurate results when coupled to an empirical dispersion correction and the DH-SVPD basis set. This protocol not only gives subchemical accuracy on the thermochemistry of alkanes but it is extremely easy to use with common quantum-chemistry codes.

Investigating the (Poly)Radicaloid Nature of Real-World Organic Compounds with DFT-Based Methods.

Recent advances in the synthesis of stable organic (open-shell) polyradicaloids have opened their application as active compounds for emerging technologies. These systems typically exhibit small energy differences between states with different spin multiplicities, which are intrinsically difficult to calculate by theoretical methods. We thus apply here some DFT-based variants (FT-DFT, SF-DFT, and SF-TDDFT) on a test set of large and real-world molecules, as test systems for which such energy differences are experimentally available, also comparing systematically with RAS-SF results to infer if shortcomings of previous DFT applications are corrected. Additionally, we explore the spin-spin contribution to the ZFS tensor, of high interest for EPR spectroscopy, and derive the spatial extent of the corresponding (photoexcited) triplet state.

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.

Computational Studies of Molecular Materials for Unconventional Energy Conversion: The Challenge of Light Emission by Thermally Activated Delayed Fluorescence.

In this paper we describe the mechanism of light emission through thermally activated delayed fluorescence (TADF)-a process able to ideally achieve 100% quantum efficiencies upon fully harvesting the energy of triplet excitons, and thus minimizing the energy loss of common (i.e., fluorescence and phosphorescence) luminescence processes. If successful, this technology could be exploited for the manufacture of more efficient organic light-emitting diodes (OLEDs) made of only light elements for multiple daily applications, thus contributing to the rise of a sustainable electronic industry and energy savings worldwide. Computational and theoretical studies have fostered the design of these all-organic molecular emitters by disclosing helpful structure-property relationships and/or analyzing the physical origin of this mechanism. However, as the field advances further, some limitations have also appeared, particularly affecting TD-DFT calculations, which have prompted the use of a variety of methods at the molecular scale in recent years. Herein we try to provide a guide for beginners, after summarizing the current state-of-the-art of the most employed theoretical methods focusing on the singlet-triplet energy difference, with the additional aim of motivating complementary studies revealing the stronger and weaker aspects of computational modelling for this cutting-edge technology.

Charge transport parameters for carbon based nanohoops and donor-acceptor derivatives.

The effect of donor-acceptor (D-A) moieties on magnitudes such as reorganization energies and electronic couplings in cycloparaphenylene (CPP) carbon based nanohoops (i.e. conjugated organic molecules with cyclic topology) is highlighted via model computations and analysis of the available crystalline structure of N,N-dimethylaza[8]CPP. For the sake of comparison, intra-molecular and inter-molecular charge transport parameters are concomitantly modelled for the recently determined herringbone polymorph of [6]CPP, along with [8]CPP and [12]CPP. The peculiar contribution of low frequency vibrations to intramolecular reorganization energies is also disclosed by computing the Huang-Rhys factors for the investigated [n]CPPs and the N,N-dimethylaza derivative. In contrast with most planar organic semiconductors where the layer in which molecules are herringbone arranged identifies the high-mobility plane, nanohoops disclose inter-layer electronic couplings larger than the intra-layer counterparts. Charge transfer rate constants modelled with three different approaches (Marcus, Marcus-Levich-Jortner and spectral overlap) suggest that D-A nanohoops, owing to orbital localization, may be more efficient for charge transport than [n]CPPs for suitable solid phase arrangements.

Double-Hybrid Functionals and Tailored Basis Set: Fullerene (C60) Dimer and Isomers as Test Cases.

A computational protocol making use of double hybrid functionals in conjunction with a recently developed basis set tailored to reproduce noncovalent interactions (hereafter named DH-SVPD) is here applied and tested for the evaluation of the properties of C60 fullerenes, namely intermolecular interactions in the weakly bound C60 dimer and relative stabilities of C60 isomers (as described by the C60ISO and iso-C60 data sets). The obtained results suggest that the DH-SVPD performance is very close to that obtained with empirical corrections and larger quadruple-ζ basis for the C60 dimer. In contrast, both approaches (tailored basis set and larger basis with empirical potential) do not reach the envisaged accuracy for the relative stabilities of C60 isomers. Nevertheless, this test well underlined how the DH-SVPD basis set is able to recover the performance obtained by coupling the DH functionals with empirical dispersion corrections and larger basis set, significantly reducing the computational effort for double hybrids and thus enabling expansion of their application domain to larger molecular systems.

sp-hybridized carbon allotrope molecular structures: An ongoing challenge for density-functional approximations.

The recent synthesis of a C18 monocyclic ring constitutes a major breakthrough as a new all-carbon disclosed form. However, modern density functional theory approaches do not lead to the correct experimental polyynic structure and favor the cumulenic one instead. We demonstrate here that this serious drawback can be solved by recently developed range-separated nonempirical schemes, independently of which kind of functional is being applied (i.e., semilocal, hybrid, or double-hybrid).

Range-separated hybrid density functionals made simple.

In this communication, we present a new and simple route to derive range-separated exchange (RSX) hybrid and double hybrid density functionals in a nonempirical fashion. In line with our previous developments [Brémond et al., J. Chem. Theory Comput. 14, 4052 (2018)], we show that by imposing an additional physical constraint to the exchange-correlation energy, i.e., by enforcing to reproduce the total energy of the hydrogen atom, we are able to generalize the nonempirical determination of the range-separation parameter to a family of RSX hybrid density functionals. The success of the resulting models is illustrated by an accurate modeling of several molecular systems and properties, like ionization potentials, particularly prone to the one- and many-electron self-interaction errors.

Quantum-Chemical Insights into the Self-Assembly of Carbon-Based Supramolecular Complexes.

Understanding how molecular systems self-assemble to form well-organized superstructures governed by noncovalent interactions is essential in the field of supramolecular chemistry. In the nanoscience context, the self-assembly of different carbon-based nanoforms (fullerenes, carbon nanotubes and graphene) with, in general, electron-donor molecular systems, has received increasing attention as a means of generating potential candidates for technological applications. In these carbon-based systems, a deep characterization of the supramolecular organization is crucial to establish an intimate relation between supramolecular structure and functionality. Detailed structural information on the self-assembly of these carbon-based nanoforms is however not always accessible from experimental techniques. In this regard, quantum chemistry has demonstrated to be key to gain a deep insight into the supramolecular organization of molecular systems of high interest. In this review, we intend to highlight the fundamental role that quantum-chemical calculations can play to understand the supramolecular self-assembly of carbon-based nanoforms through a limited selection of supramolecular assemblies involving fullerene, fullerene fragments, nanotubes and graphene with several electron-rich π-conjugated systems.

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.

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.