General optimization procedure towards the design of a new family of minimal parameter spin-component-scaled double-hybrid density functional theory

A general optimization procedure towards the development and implementation of a new family of minimal parameter spin-component-scaled double-hybrid (mSD) density functional theory (DFT) is presented. The nature of the proposed exchange-correlation functional establishes a methodology with minimal empiricism. This new family of double-hybrid (DH) density functionals is demonstrated using the PBEPBE functional, illustrating the optimization procedure to the mSD-PBEPBE method, and the performance characteristics shown for a set of non-covalent complexes covering a broad regime of weak interactions. With only two parameters, mSD-PBEPBE and its cost-effective counterpart, RI-mSD-PBEPBE, show a mean absolute error of ca. 0.4 kcal mol^-1 averaged over 66 weak interacting systems. Following a successive 2D-grid refinement for a CBS extrapolation of the coefficients, the optimization procedure can be recommended for the design and implementation of a variety of additional DH methods using any of the plethora of currently available functionals.

Time-Dependent Double-Hybrid Density Functionals with Spin-Component and Spin-Opposite Scaling

For the first time, we combine time-dependent double-hybrid density functional approximations (TD-DHDFAs) for the calculation of electronic excitation energies with the concepts of spin-component and spin-opposite scaling (SCS/SOS) of electron-pair contributions to their nonlocal correlation components. Different flavors of this idea, ranging from standard SCS parameters to fully fitted parameter sets, are presented and tested on six different parent DHDFAs. For cross-validation, we assess those methods on three benchmark sets that cover small- to medium-sized chromophores (up to 78 atoms) and different excitation types. For this purpose, we also introduce new CC3 reference values for the popular Gordon benchmark set that we recommend using in future studies. Our results confirm that already the (unscaled) parent TD-DHDFAs are accurate and outperform some wave function methods. Further introduction of SCS/SOS eliminates extreme outliers, reduces deviation spans from reference values by up to 0.5 eV, aligns the performance of the Tamm-Dancoff approximation (TDA) to that of full TD calculations, and also enables a more balanced description of different excitation types. The best-performing TD-based methods in our cross validation have mean absolute deviations as low as 0.14 eV compared to the time- and resource-intensive CC3 approach. A very important finding is that we also obtained SOS variants with excellent performance, contrary to wave function based methods. This opens a future pathway to highly efficient methods for the optimization of excited-state geometries, particularly when paired with computing strategies such as the Laplace transform. We recommend our SCS- and SOS-based variants for further testing and subsequent applications.

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.

Nonoxido V(IV) Complexes: Prediction of the EPR Spectrum and Electronic Structure of Simple Coordination Compounds and Amavadin

Density functional theory (DFT) calculations of the (51)V hyperfine coupling (HFC) tensor A have been completed for 20 "bare" V(IV) complexes with different donor sets, electric charges, and coordination geometries. Calculations were performed with ORCA and Gaussian software, using functionals BP86, TPSS0, B1LYP, PBE0, B3LYP, B3P, B3PW, O3LYP, BHandHLYP, BHandH, and B2PLYP. Among the basis sets, 6-311g(d,p), 6-311++g(d,p), VTZ, cc-pVTZ, def2-TZVPP, and the "core properties" CP(PPP) were tested. The experimental Aiso and Ai (where i = x or z, depending on the geometry and electronic structure of V(IV) complex) were compared with the values calculated by DFT methods. The results indicated that, based on the mean absolute percentage deviation (MAPD), the best functional to predict Aiso or Ai is the double hybrid B2PLYP. With this functional and the basis set VTZ, it is possible to predict the Aiso and Az of the EPR spectrum of amavadin with deviations of -1.1% and -2.0% from the experimental values. The results allowed us to divide the spectra of nonoxido V(IV) compounds in three types-called "type 1", "type 2", and "type 3", characterized by different composition of the singly occupied molecular orbital (SOMO) and relationship between the values of Ax, Ay, and Az. For "type 1" spectra, Az ≫ Ax ≈ Ay and Az is in the range of (135-155) × 10(-4) cm(-1); for "type 2" spectra, Ax ≈ Ay ≫ Az and Ax ≈ Ay are in the range of (90-120) × 10(-4) cm(-1); and for the intermediate spectra of "type 3", Az > Ay > Ax or Ax > Ay > Az, with Az or Ax values in the range of (120-135) × 10(-4) cm(-1). The electronic structure of the V(IV) species was also discussed, and the results showed that the values of Ax or Az are correlated with the percent contribution of V-dxy orbital in the SOMO. Similarly to V(IV)O species, for amavadin the SOMO is based mainly on the V-dxy orbital, and this accounts for the large experimental value of Az (153 × 10(-4) cm(-1)).

Seeking for Spin-Opposite-Scaled Double-Hybrid Models Free of Fitted Parameters

On the basis of theoretical arguments, a new exchange-correlation energy expression free of any fitted parameter has been proposed for spin-opposite-scaled double-hybrid density functionals (SOS0-DHs). Employing the recently presented DHs, the working expressions for SOS0-DH functionals are obtained and benchmarked numerically against several standard databases. Our test calculations show that for some cases such as interaction energies and barrier heights the SOS0-DHs without dispersion corrections perform better than their non-SOS counterparts. On the other hand, for other properties like atomization energies, the conventional DHs provide reliable results. We hope that the findings of this work can excite further developments of DH functionals in the framework of SOS scheme for a wide variety of applications resolving the failures at a reasonable computational cost. It seems that a bright future lies ahead in this arena.

Theoretical rationalization of the singlet-triplet gap in OLEDs materials: impact of charge-transfer character

New materials for OLED applications with low singlet-triplet energy splitting have been recently synthesized in order to allow for the conversion of triplet into singlet excitons (emitting light) via a Thermally Activated Delayed Fluorescence (TADF) process, which involves excited-states with a non-negligible amount of Charge-Transfer (CT). The accurate modeling of these states with Time-Dependent Density Functional Theory (TD-DFT), the most used method so far because of the favorable trade-off between accuracy and computational cost, is however particularly challenging. We carefully address this issue here by considering materials with small (high) singlet-triplet gap acting as emitter (host) in OLEDs and by comparing the accuracy of TD-DFT and the corresponding Tamm-Dancoff Approximation (TDA), which is found to greatly reduce error bars with respect to experiments thanks to better estimates for the lowest singlet-triplet transition. Finally, we quantitatively correlate the singlet-triplet splitting values with the extent of CT, using for it a simple metric extracted from calculations with double-hybrid functionals, that might be applied in further molecular engineering studies.

Analytical gradients for MP2, double hybrid functionals, and TD-DFT with polarizable embedding described by fluctuating charges

A polarizable quantum mechanics (QM)/ molecular mechanics (MM) approach recently developed for Hartree-Fock (HF) and Kohn-Sham (KS) methods has been extended to energies and analytical gradients for MP2, double hybrid functionals, and TD-DFT models, thus allowing the computation of equilibrium structures for excited electronic states together with more accurate results for ground electronic states. After a detailed presentation of the theoretical background and of some implementation details, a number of test cases are analyzed to show that the polarizable embedding model based on fluctuating charges (FQ) is remarkably more accurate than the corresponding electronic embedding based on a fixed charge (FX) description. In particular, a set of electronegativities and hardnesses has been optimized for interactions between QM and FQ regions together with new repulsion-dispersion parameters. After validation of both the numerical implementation and of the new parameters, absorption electronic spectra have been computed for representative model systems including vibronic effects. The results show remarkable agreement with full QM computations and significant improvement with respect to the corresponding FX results. The last part of the article provides some hints about computation of solvatochromic effects on absorption spectra in aqueous solution as a function of the number of FQ water molecules and on the use of FX external shells to improve the convergence of the results. © 2015 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

Mechanistic Aspects of Aryl-Halide Oxidative Addition, Coordination Chemistry, and Ring-Walking by Palladium

This contribution describes the reactivity of a zero-valent palladium phosphine complex with substrates that contain both an aryl halide moiety and an unsaturated carbon-carbon bond. Although η(2) -coordination of the metal center to a C=C or C≡C unit is kinetically favored, aryl halide bond activation is favored thermodynamically. These quantitative transformations proceed under mild reaction conditions in solution or in the solid state. Kinetic measurements indicate that formation of η(2) -coordination complexes are not nonproductive side-equilibria, but observable (and in several cases even isolated) intermediates en route to aryl halide bond cleavage. At the same time, DFT calculations show that the reaction with palladium may proceed through a dissociation-oxidative addition mechanism rather than through a haptotropic intramolecular process (i.e., ring walking). Furthermore, the transition state involves coordination of a third phosphine to the palladium center, which is lost during the oxidative addition as the C-halide bond is being broken. Interestingly, selective activation of aryl halides has been demonstrated by adding reactive aryl halides to the η(2) -coordination complexes. The product distribution can be controlled by the concentration of the reactants and/or the presence of excess phosphine.

Construction of exchange-correlation functionals through interpolation between the non-interacting and the strong-correlation limit

Drawing on the adiabatic connection of density functional theory, exchange-correlation functionals of Kohn-Sham density functional theory are constructed which interpolate between the extreme limits of the electron-electron interaction strength. The first limit is the non-interacting one, where there is only exchange. The second limit is the strong correlated one, characterized as the minimum of the electron-electron repulsion energy. The exchange-correlation energy in the strong-correlation limit is approximated through a model for the exchange-correlation hole that is referred to as nonlocal-radius model [L. O. Wagner and P. Gori-Giorgi, Phys. Rev. A 90, 052512 (2014)]. Using the non-interacting and strong-correlated extremes, various interpolation schemes are presented that yield new approximations to the adiabatic connection and thus to the exchange-correlation energy. Some of them rely on empiricism while others do not. Several of the proposed approximations yield the exact exchange-correlation energy for one-electron systems where local and semi-local approximations often fail badly. Other proposed approximations generalize existing global hybrids by using a fraction of the exchange-correlation energy in the strong-correlation limit to replace an equal fraction of the semi-local approximation to the exchange-correlation energy in the strong-correlation limit. The performance of the proposed approximations is evaluated for molecular atomization energies, total atomic energies, and ionization potentials.

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.

Astrochemistry of transition metals? The selected cases of [FeN]+, [FeNH]+ and [(CO)2FeN]+: pathways toward CH3NH2 and HNCO

Transition metals are proposed to play an active role in the synthesis of organic compounds containing heteroatoms in astrochemistry.