Does a Molecule-Specific Density Functional Give an Accurate Electron Density? The Challenging Case of the CuCl Electric Field Gradient

Prediction of 57Fe Mössbauer Nuclear Quadrupole Splittings with Hybrid and Double-Hybrid Density Functionals

As a crucial parameter in Mössbauer spectroscopy, nuclear quadrupole splitting (NQS) exhibits a strong dependence on quantum chemistry methods, which makes accurate theoretical predictions challenging. Meanwhile, the continuous emergence of new density functionals presents opportunities to advance current NQS research. In this study, we evaluate the performance of eleven hybrid density functionals and twelve double-hybrid density functionals, selected from widely used functionals and newly developed functionals, in predicting the NQS values of the 57Fe nuclide for 32 iron-containing molecules within about 70 atoms. The calculations have incorporated scalar relativistic effects using the exact two-component (X2C) Hamiltonian. In general, the double-hybrid functional PBE-0DH demonstrates superior performance compared to the experimental values, achieving a mean absolute error (MAE) of 0.20 mm/s. Meanwhile, rSCAN38 is the best hybrid functional for our database with an MAE = 0.25 mm/s, and it offers a significant advantage in computational efficiency over PBE-0DH. The +/- sign of NQS has also been considered in our error statistics when it has a clear physical meaning; if neglected, the errors of many functionals decrease, but PBE-0DH and rSCAN38 remain unaffected. Notably, when calculating ferrocene [Fe(C5H5)2], which involves strong static correlations, all hybrid functionals that incorporate more than 10% exact exchange fail, while several double-hybrid functionals continue to deliver reliable results. In addition, we encountered two particularly challenging species characterized by strong static correlations: [Fe(H2O)5NO]2+ and FeO2--porphyrin. Unfortunately, none of the density functionals tested in our study yielded satisfactory results for the two cases since the density functional theory (DFT) is a single-determinant approach, and it is imperative to explore large-scale multi-configurational methods for these species. This research offers valuable guidance for selecting density functionals in Mössbauer NQS calculations and serves as a reference point for the future development of new density functionals.

A detailed TD-DFT and intermolecular interaction study of vitamin K in soluble, poorly soluble and insoluble solvents, as well as an ADME and molecular docking study

Vitamin K is one of the most important fat-soluble vitamins and while there are two main types of vitamin K in nature, known as K1 (phylloquinone) and K2 (menaquinones), there is also a synthetic type of vitamin K known as K3 (menadione). Recent studies have shown that it is crucial to know the non-covalent interactions, ADME and molecular docking of molecules in different solvent media. Therefore, we have performed some quantum chemical calculations, ADME and intra-and intermolecular interaction calculations of a number of K1, K2 and K3 such as K1-water (K1 + W), K1-methanol (K1 + M), K1-triacetin (K1 + T), K2-water (K2 + W), K2-methanol (K2 + M), K2-triacetin (K2 + T), K3-water (K3 + W), K3-methanol (K3 + M), K3-triacetin (K3 + T) performed by Density Functional Theory (DFT) and Multiwfn: A multifunctional wavefunction analyzer. Molecular structures, HOMO-LUMO energies, MEP and electronic properties have been calculated and described using DFT at the level of B3LYP/6-311G (d,p) level. The nature of the molecular interactions between vitamin K and solvents such as water, methanol and triacetin were also investigated using topological analyses such as atoms in molecule (AIM), non-covalent interaction index (NCI), reduced density gradient (RDG), Localized orbital locator (LOL) and electron localization function (ELF). In addition, FMO for electronic transitions, MEP for electrophilic and nucleophilic attack, ADME to investigate how a chemical is processed by a living organism, and Fukui functions to determine electron density are explained. Finally, molecular docking was used to determine the biological activity of the vitamin K.

Predicting the Color Polymorphism of ROY from a Time-Dependent Optimally Tuned Screened Range-Separated Hybrid Functional

Polymorphism is a well-known property of molecular crystals, which allows the same molecule to form solids with several crystalline structures that can differ significantly in physical properties. Polymorphs that possess different optical absorption properties in the visible range may exhibit different perceived colors, a phenomenon known as color polymorphism. One striking example of color polymorphism is given by 5-methyl-2-[(2-nitrophenyl)amino]-3-thiophenecarbonitrile, known as ROY for its red-orange-yellow colors. First-principles prediction of color polymorphism may help in polymorph assignment and design but has proven to be challenging. Here, we predict the absorption spectra and simulate the colors of 12 ROY polymorphs using the general, nonempirical method of time-dependent (TD) optimally tuned screened range-separated hybrid (OT-SRSH) functional. For 5 ROY polymorphs with known experimental absorption spectra, we show that the TD-OT-SRSH approach predicts absorption spectra in quantitative agreement with experiment. For all polymorphs, we show that an accurate simulation of the colors is obtained, paving the way to a fully predictive, low-cost calculation of color polymorphism.

Calculation of electric field gradients with the exact two-component (X2C) quasi-relativistic method and its local approximations

When calculating electric field gradients (EFGs), relativistic and electron correlation effects are crucial for obtaining accurate results, and the commonly used density functional methods produce unsatisfactory results, especially for heavy elements and/or strongly correlated systems. In this work, a stand-alone program is presented, which enables calculation of EFGs from the molecular orbitals supplied by an external high accuracy quantum chemical calculation and includes relativistic effects through the exact two-component (X2C) formalism and efficient local approximations to it. Application to BiN and BiP molecules shows that a high precision can be achieved in the calculation of nuclear quadrupole coupling constants of 209Bi by combining advanced ab initio methods with the X2C approach. For seventeen iron compounds, the Mössbauer nuclear quadrupole splittings (NQS) of 57Fe calculated using a double-hybrid functional method are in very good agreement with the experimental values. It is shown that, for strongly correlated molecules, the double-hybrid functionals are much more accurate than the commonly used hybrid functionals. The computer program developed in this study furnishes a useful utility for obtaining EFGs and related nuclear properties with high accuracy.

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.

How does theory compare to experiment for oscillator strengths in electronic spectra? Proposing range-separated hybrids with reliable accountability

As an important quantity in atomic and molecular spectroscopy, oscillator strength should be mentioned. Oscillator strength is linked to the transition dipole moment and consequently to the transition probability between two states, where its magnitude is directly connected to the intensity of the peaks in ultraviolet-visible spectra. However, accurately accounting for oscillator strengths still remains one of the greatest challenges in theory and experiment. Given previous efforts in the context of investigations into oscillator strengths, the related theoretical treatments are relatively limited and have proven to be challenging. In this work, the oscillator strengths in the electronic spectra of organic compounds have thoroughly been investigated with the help of optimally tuned range-separated hybrids (OT-RSHs). In particular, variants of the OT-RSHs combined with the polarizable continuum model (PCM), OT-RSHs-PCM, as well as their screened versions accounting for the screening effects by the electron correlation through the dielectric constant, OT-SRSHs-PCM, are proposed for reliable prediction of the oscillator strengths. The role of the involved ingredients in the proposed methods, namely the underlying density functional approximations, short-range and long-range Hartree-Fock (HF) exchange, as well as the range-separation parameter, has been examined in detail. It is shown that any combination of the parameters in the proposed approximations does not render the reliable oscillator strengths, but a particular compromise among them is needed to describe the experimental data well. Perusing all the results of our developed methods, the best ones are found to be the generalized gradient approximation-based OT-RSHs-PCM, coupled with the linear response theory in the non-equilibrium solvation regime, with the correct asymptotic behavior and incorporating no (low) HF exchange contributions in the short-range part. The best proposed approximations also reveal superior performances not only with respect to their standard counterparts with the default parameters but also as compared to earlier range-separated functionals. Finally, the applicability of the best approximation is also put into broader perspective, where it is used for predicting the oscillator strengths in other sets of compounds not included in the process of developing the approximations. Hopefully, our proposed method can function as an affordable alternative to the expensive wave function-based methods for both theoretical modeling and confirming the experimental observations in the field of electronic spectroscopy.

Understanding Photophysical Properties of Molecules Relevant in Organic Semiconductor Laser Diodes from Electron Localization Function-Tuned and Solvent-Tuned Range-Separated Functionals

Organic semiconductor laser diodes (OSLDs) are prevalent in optoelectronics because of their sustainable energy applications. Organic molecules used in such diodes are usually large; hence, their studies are computationally challenging with high-end benchmark methods. Computational methods with reliable accuracy and efficiency are always indispensable. In the present work, we have applied our computationally inexpensive, nonempirically tuned [electron localization function (ELF*) and solvent (Sol*)] range-separated (RS) functionals to study five molecules used in OSLDs. The emission energies in three different environments [toluene, CBP (4,4'-bis(n-carbazolyl)-1,1'-biphenyl) film, and gas] have been computed with the tuned functionals and compared with the experimental emission energies. ELF* and Sol* functionals can accurately reproduce emission energies in toluene and CBP film environments. On the other hand, both ELF* and IP-tuned functionals with excited-state geometry (IP*) perform better in the gas phase. In addition, a comparative study is performed between time-dependent density functional theory and the Tamm-Dancoff approximation. Along with the emission energy, oscillator strength values have also been reported. Different IP-tuned RS parameters were obtained with the ground- and excited-state geometries. Interestingly, it has been observed that the optimally tuned RS parameter with excited-state geometry (IP*) performs better compared to that with ground-state geometries (IP). Fractional occupation calculations show that the tuned functionals exhibit less localization and delocalization error. The study envisages that ELF* and Sol* functionals can be used to design future candidates for OSLDs.

Koopmans' Theorem-Compliant Long-Range Corrected (KTLC) Density Functional Mediated by Black-Box Optimization and Data-Driven Prediction for Organic Molecules

Density functional theory (DFT) is a significant computational tool that has substantially influenced chemistry, physics, and materials science. DFT necessitates parametrized approximation for determining an expected value. Hence, to predict the properties of a given molecule using DFT, appropriate parameters of the functional should be set for each molecule. Herein, we optimize the parameters of range-separated functionals (LC-BLYP and CAM-B3LYP) via Bayesian optimization (BO) to satisfy Koopmans' theorem. Our results demonstrate the effectiveness of the BO in optimizing functional parameters. Particularly, Koopmans' theorem-compliant LC-BLYP (KTLC-BLYP) shows results comparable to the experimental UV-absorption values. Furthermore, we prepared an optimized parameter dataset of KTLC-BLYP for over 3000 molecules through BO for satisfying Koopmans' theorem. We have developed a machine learning model on this dataset to predict the parameters of the LC-BLYP functional for a given molecule. The prediction model automatically predicts the appropriate parameters for a given molecule and calculates the corresponding values. The approach in this paper would be useful to develop new functionals and to update the previously developed functionals.

Density-Corrected Density Functional Theory for Molecular Properties

Density-corrected (DC) density functional theory (DFT) has been proposed to overcome difficulties related to the self-interaction error. The procedure uses the Hartree-Fock electron density (matrix) non-self-consistently in conjunction with an approximate functional. DC-DFT has so far mainly been tested for total energy differences, whereas other types of molecular properties have not been evaluated systematically. This work focuses on the performance of DC-DFT for molecular properties, namely, dipole moments, static polarizabilities, and electric field gradients (EFGs) at atomic nuclei. Accurate reference data were generated from coupled-cluster theory to assess the performance of DC and self-consistent DFT calculations for twelve molecules, including diatomics with transition metals. DC-DFT does no harm in dipole moment calculations, but it negatively impacts the polarizability in at least one case. DC-DFT performs well for EFGs, even for the difficult case of CuCl.

Insights into the deviation from piecewise linearity in transition metal complexes from supervised machine learning models

Virtual high-throughput screening (VHTS) and machine learning (ML) with density functional theory (DFT) suffer from inaccuracies from the underlying density functional approximation (DFA). Many of these inaccuracies can be traced to the lack of derivative discontinuity that leads to a curvature in the energy with electron addition or removal. Over a dataset of nearly one thousand transition metal complexes typical of VHTS applications, we computed and analyzed the average curvature (i.e., deviation from piecewise linearity) for 23 density functional approximations spanning multiple rungs of "Jacob's ladder". While we observe the expected dependence of the curvatures on Hartree-Fock exchange, we note limited correlation of curvature values between different rungs of "Jacob's ladder". We train ML models (i.e., artificial neural networks or ANNs) to predict the curvature and the associated frontier orbital energies for each of these 23 functionals and then interpret differences in curvature among the different DFAs through analysis of the ML models. Notably, we observe spin to play a much more important role in determining the curvature of range-separated and double hybrids in comparison to semi-local functionals, explaining why curvature values are weakly correlated between these and other families of functionals. Over a space of 187.2k hypothetical compounds, we use our ANNs to pinpoint DFAs for which representative transition metal complexes have near-zero curvature with low uncertainty, demonstrating an approach to accelerate screening of complexes with targeted optical gaps.

Toward highly efficient hyperfluorescence-based emitters through excited-states alignment using novel optimally tuned range-separated models

Hyperfluorescence has recently been introduced as a promising strategy to achieve organic light-emitting diodes (OLEDs) with high color purity and enhanced stability. In this approach, fluorescent emitters (FEs) with strong and narrow band fluorescence are integrated in thin films containing sensitizers exhibiting thermally activated delayed fluorescence (TADF). Toward highly efficient hyperfluorescence-based emitters, the excited-states ordering of the FEs should be well-aligned. Given some recent endeavors in this context, the related theoretical explorations are relatively limited and have proven to be challenging. In this work, alignments of the corresponding excited-states, crucial for both the fast Förster resonance energy transfer and suppression of the Dexter energy transfer from TADF sensitizers to FEs, have theoretically been investigated using optimally tuned range-separated hybrid functionals (OT-RSHs). We have proposed and validated several variants of the models including OT-RSHs, their coupled versions with the polarizable continuum model, OT-RSHs-PCM, as well as the screened versions accounting for the screening effects by the electron correlation through the scalar dielectric constant, OT-SRSHs, for a reliable description of the excited-states ordering in the FEs of the hyperfluorescence-based materials. Particular attention is paid to the influence of the underlying density functional approximations as well as the short- and long-range Hartree-Fock (HF) exchange contributions and the range-separation parameter. Considering a series of experimentally known hyperfluorescence-based emitters as working models, it is unveiled that any combination of the ingredients in the proposed models does not render the correct order of the excited-states of the FEs, but a particular compromise among the involved parameters is needed to more accurately account for the relevant excited-states alignment. Perusing the results of our developed methods, the best ones are found to be the generalized gradient approximation-based OT-RSHs-PCM with the correct asymptotic behavior and incorporating no (low) HF exchange contribution at the short-range regime. The proposed models show superior performances not only with respect to their standard counterparts with the default parameters but also as compared to other range-separated approximations. Accountability of the best-proposed model is also put into broader perspective, where it has been employed for the computational design of several molecules as promising FE candidates prone to be utilized in hyperfluorescence-based materials. Summing up, the proposed models in this study can be recommended for both the theoretical modeling and confirming the experimental observations in the field of hyperfluorescence-based OLEDs.

An optimally tuned range-separated hybrid starting point for ab initio GW plus Bethe–Salpeter equation calculations of molecules

The ab initio GW plus Bethe–Salpeter equation (GW-BSE, where G is the one particle Green's function and W is the screened Coulomb interaction) approach has emerged as a leading method for predicting excitations in both solids and molecules with a predictive power contingent upon several factors. Among these factors are the (1) generalized Kohn–Sham eigensystem used to construct the GW self-energy and to solve the BSE and (2) the efficacy and suitability of the Tamm–Dancoff approximation. Here, we present a detailed benchmark study of low-lying singlet excitations from a generalized Kohn–Sham (gKS) starting point based on an optimally tuned range-separated hybrid (OTRSH) functional. We show that the use of this gKS starting point with one-shot G0W0 and G0W0-BSE leads to the lowest mean absolute errors (MAEs) and mean signed errors (MSEs), with respect to high-accuracy reference values, demonstrated in the literature thus far for the ionization potentials of the GW100 benchmark set and for low-lying neutral excitations of Thiel’s set molecules in the gas phase, without the need for self-consistency. The MSEs and MAEs of one-shot G0W0-BSE@OTRSH excitation energies are comparable to or lower than those obtained with other functional starting points after self-consistency. Additionally, we compare these results with linear-response time-dependent density functional theory (TDDFT) calculations and find GW-BSE to be superior to TDDFT when calculations are based on the same exchange-correlation functional. This work demonstrates tuned range-separated hybrids used in combination with GW and GW-BSE can greatly suppress starting point dependence for molecules, leading to accuracy similar to that for higher-order wavefunction-based theories for molecules without the need for costlier iterations to self-consistency.

Excited-State Properties of Some Thermally Activated Delayed Fluorescence Emitters: Quest for an Accurate and Reliable Computational Method

Thermally activated delayed fluorescence (TADF) finds application in organic light-emitting diodes. The molecules exhibiting TADF are characterized by small singlet-triplet energy gaps that help reverse intersystem crossing. Recently, ionization potential (IP)-tuned range-separated (RS) density functionals have been well accepted for studying excited-state properties. In the present work, two efficient descriptor-based tuning schemes [electron localization function (ELF) and Sol] of RS density functionals have been used to accurately reproduce the excited-state properties of TADF emitters by performing a single self-consistent field calculation. The lowest singlet vertical excitation energies (E_VA(S₁)) and the vertical singlet-triplet energy gaps (ΔE_VST) are computed with ELF-, Sol-, and IP-tuned RS functionals (LC-BLYP, ωB97, ωB97X, and ωB97XD). Encouraging mean absolute deviations from the experimental values with ELF*-, Sol*-, and IP-tuned functionals are observed. Consistent performance of the non-empirical tuned functionals is noted in different solvent dielectrics. In addition to these, fractional occupation calculations have shown that our tuned functionals almost satisfy the energy linearity curve. Thus, ELF*- and Sol*-tuned functionals are promising and reliable alternatives in computing the excited-state properties. Considering the small experimental singlet-triplet gap, we recommend ELF* to calculate E_VA(S₁) and Sol* to calculate ΔE_VST.

Ligand Additivity and Divergent Trends in Two Types of Delocalization Errors from Approximate Density Functional Theory

The predictive accuracy of density functional theory (DFT) is hampered by delocalization errors, especially for correlated systems such as transition-metal complexes. Two complementary strategies have been developed to reduce delocalization error: eliminating the global curvature with change in charge, and applying a linear response Hubbard U as a measure of local curvature at a metal center at fixed charge in a DFT+U framework. We investigate the relationship between the two delocalization error measures as the ligand field strength is varied with the number of strong-field ligands in a series of heteroleptic complexes or by geometrically constraining the metal-ligand bond length in homoleptic octahedral complexes. We show that across these sets of complexes an inverse relationship generally exists between global and local curvatures. We find that effects of ligand substitution on both measures of delocalization are typically additive, but the quantities seldom coincide.

Molecular orbital projectors in non-empirical jmDFT recover exact conditions in transition-metal chemistry

Low-cost, non-empirical corrections to semi-local density functional theory are essential for accurately modeling transition-metal chemistry. Here, we demonstrate the judiciously modified density functional theory (jmDFT) approach with non-empirical U and J parameters obtained directly from frontier orbital energetics on a series of transition-metal complexes. We curate a set of nine representative Ti(III) and V(IV) d1 transition-metal complexes and evaluate their flat-plane errors along the fractional spin and charge lines. We demonstrate that while jmDFT improves upon both DFT+U and semi-local DFT with the standard atomic orbital projectors (AOPs), it does so inefficiently. We rationalize these inefficiencies by quantifying hybridization in the relevant frontier orbitals. To overcome these limitations, we introduce a procedure for computing a molecular orbital projector (MOP) basis for use with jmDFT. We demonstrate this single set of d1 MOPs to be suitable for nearly eliminating all energetic delocalization error and static correlation error. In all cases, MOP jmDFT outperforms AOP jmDFT, and it eliminates most flat-plane errors at non-empirical values. Unlike DFT+U or hybrid functionals, jmDFT nearly eliminates energetic delocalization error and static correlation error within a non-empirical framework.

Optimal Tuning Perspective of Range-Separated Double Hybrid Functionals

We study the optimal tuning of the free parameters in range-separated double hybrid functionals, based on enforcing the exact conditions of piecewise linearity and spin constancy. We find that introducing the range separation in both the exchange and the correlation terms allows for the minimization of both fractional charge and fractional spin errors for singlet atoms. The optimal set of parameters is system specific, underlining the importance of the tuning procedure. We test the performance of the resulting optimally tuned functionals for the dissociation curves of diatomic molecules. We find that they recover the correct dissociation curve for the one-electron system, H2+, and improve the dissociation curves of many-electron molecules such as H2 and Li2, but they also yield a nonphysical maximum and only converge to the correct dissociation limit at very large distances.

Excited-state properties of organic semiconductor dyes as electrically pumped lasing candidates from new optimally tuned range-separated models

Even though many efforts have been devoted to optical lasing in recent years, the realization of lasing by direct electrical excitation of organic semiconductors is hampered mainly due to optical losses from electrical contacts and electrical losses induced by triplets and polarons at high current densities. Hereby, accurately accounting for the electrically pumped organic semiconductor laser diodes (OSLDs) still remains one of the greatest challenges in optoelectronics. In this work, the excited-state characteristics of the organic semiconductor dyes used in the electrically pumped OSLDs have thoroughly been investigated using optimally tuned range-separated hybrids (OT-RSHs). Considering several experimentally known compounds of the electrically pumped OSLDs as working models, several variants of OT-RSHs, their combination forms with the polarizable continuum model (PCM), OT-RSH-PCM, as well as their screened versions accounting for the screening effects by the electron correlation through the scalar dielectric constant, OT-SRSHs, have been proposed for reliable prediction of their emission energies and oscillator strengths in both the gas and solvent phases. The role of involved ingredients in the models, namely, the underlying density functional approximations, short- and long-range exact-like exchange, as well as the range-separation parameter, has been examined in detail. It is shown that the newly designed OT-RSHs with the correct behavior of asymptotic exchange-correlation potential outperform the standard RSHs and other density functionals with both fixed and interelectronic distance-dependent exact-like exchange for describing the excite-state properties of compounds of the electrically pumped OSLDs. Concerning the computational cost of the models, it is unveiled that performing both the optimal tuning procedure and subsequent excited-state computations using OT-RSHs in the gas phase can be considered as a more reliable and affordable framework. Finally, the applicability of the proposed models is also put into a broader perspective for the computational design of several compounds as promising candidates to be used in the OSLD materials. Hopefully, our recommended OT-RSHs can function as efficient models for both the related theoretical modeling and confirming the experimental observations in the field of electrically pumped OSLDs.

A 9.2-GHz clock transition in a Lu(II) molecular spin qubit arising from a 3,467-MHz hyperfine interaction

Spins in molecules are particularly attractive targets for next-generation quantum technologies, enabling chemically programmable qubits and potential for scale-up via self-assembly. Here we report the observation of one of the largest hyperfine interactions for a molecular system, A_iso = 3,467 ± 50 MHz, as well as a very large associated clock transition. This is achieved through chemical control of the degree of s-orbital mixing into the spin-bearing d orbital associated with a series of spin-½ La(II) and Lu(II) complexes. Increased s-orbital character reduces spin-orbit coupling and enhances the electron-nuclear Fermi contact interaction. Both outcomes are advantageous for quantum applications. The former reduces spin-lattice relaxation, and the latter maximizes the hyperfine interaction, which, in turn, generates a 9-GHz clock transition, leading to an increase in phase memory time from 1.0 ± 0.4 to 12 ± 1 μs for one of the Lu(II) complexes. These findings suggest strategies for the development of molecular quantum technologies, akin to trapped ion systems.

Implementation of Relativistic Coupled Cluster Theory for Massively Parallel GPU-Accelerated Computing Architectures

In this paper, we report reimplementation of the core algorithms of relativistic coupled cluster theory aimed at modern heterogeneous high-performance computational infrastructures. The code is designed for parallel execution on many compute nodes with optional GPU coprocessing, accomplished via the new ExaTENSOR back end. The resulting ExaCorr module is primarily intended for calculations of molecules with one or more heavy elements, as relativistic effects on the electronic structure are included from the outset. In the current work, we thereby focus on exact two-component methods and demonstrate the accuracy and performance of the software. The module can be used as a stand-alone program requiring a set of molecular orbital coefficients as the starting point, but it is also interfaced to the DIRAC program that can be used to generate these. We therefore also briefly discuss an improvement of the parallel computing aspects of the relativistic self-consistent field algorithm of the DIRAC program.

Computational Discovery of Transition-metal Complexes: From High-throughput Screening to Machine Learning

Transition-metal complexes are attractive targets for the design of catalysts and functional materials. The behavior of the metal-organic bond, while very tunable for achieving target properties, is challenging to predict and necessitates searching a wide and complex space to identify needles in haystacks for target applications. This review will focus on the techniques that make high-throughput search of transition-metal chemical space feasible for the discovery of complexes with desirable properties. The review will cover the development, promise, and limitations of "traditional" computational chemistry (i.e., force field, semiempirical, and density functional theory methods) as it pertains to data generation for inorganic molecular discovery. The review will also discuss the opportunities and limitations in leveraging experimental data sources. We will focus on how advances in statistical modeling, artificial intelligence, multiobjective optimization, and automation accelerate discovery of lead compounds and design rules. The overall objective of this review is to showcase how bringing together advances from diverse areas of computational chemistry and computer science have enabled the rapid uncovering of structure-property relationships in transition-metal chemistry. We aim to highlight how unique considerations in motifs of metal-organic bonding (e.g., variable spin and oxidation state, and bonding strength/nature) set them and their discovery apart from more commonly considered organic molecules. We will also highlight how uncertainty and relative data scarcity in transition-metal chemistry motivate specific developments in machine learning representations, model training, and in computational chemistry. Finally, we will conclude with an outlook of areas of opportunity for the accelerated discovery of transition-metal complexes.

Molecular DFT+U: A Transferable, Low-Cost Approach to Eliminate Delocalization Error

While density functional theory (DFT) is widely applied for its combination of cost and accuracy, corrections (e.g., DFT+U) that improve it are often needed to tackle correlated transition-metal chemistry. In principle, the functional form of DFT+U, consisting of a set of localized atomic orbitals (AOs) and a quadratic energy penalty for deviation from integer occupations of those AOs, enables the recovery of the exact conditions of piecewise linearity and the derivative discontinuity. Nevertheless, for practical transition-metal complexes, where both atomic states and ligand orbitals participate in bonding, standard DFT+U can fail to eliminate delocalization error (DE). Here, we show that by introducing an alternative valence-state (i.e., molecular orbital or MO) basis to the DFT+U approach, we recover exact conditions in cases for which standard DFT+U corrections have no error-reducing effect. This MO-based DFT+U also eliminates DE where standard AO-based DFT+U is already successful. We demonstrate the transferability of our approach on representative transition-metal complexes with a range of ligand field strengths, electron configurations (i.e., from Sc to Zn), and spin states.

Accurate Molecular Geometries in Complex Excited-State Potential Energy Surfaces from Time-Dependent Density Functional Theory

The interplay of electronic excitations and structural changes in molecules impacts nonradiative decay and charge transfer in the excited state, thus influencing excited-state lifetimes and photocatalytic reaction rates in optoelectronic and energy devices. To capture such effects requires computational methods providing an accurate description of excited-state potential energy surfaces and geometries. We suggest time-dependent density functional theory using optimally tuned range-separated hybrid (OT-RSH) functionals as an accurate approach to obtain excited-state molecular geometries. We show that OT-RSH provides accurate molecular geometries in excited-state potential energy surfaces that are complex and involve an interplay of local and charge-transfer excitations, for which conventional semilocal and hybrid functionals fail. At the same time, the nonempirical OT-RSH approach maintains the high accuracy of parametrized functionals (e.g., B3LYP) for predicting excited-state geometries of small organic molecules showing valence excited states.

Polymorphism Dependent 9-Phosphoanthracene Derivative Exhibiting Thermally Activated Delayed Fluorescence: A Computational Investigation

Polymorphs of anthracene derivatives exhibit diverse photophysical properties that can help to develop efficient organic-based photovoltaic devices. 10-Anthryl-9-phosphoanthracene (10-APA) shows different photophysical behaviors for the solid state due to its variety in crystalline arrangement. Herein, we investigate the ground and excited-state properties of the monomer and two different polymorphs of 10-APA from first-principles. Calculations reveal that strong spin-orbit coupling (SOC) between first excited singlet state (S₁) and triplet manifolds at their S₁-optimized geometries enabling the reverse intersystem crossing (RISC). The electron-vibration coupling (Huang-Rhys factor) in the excited state is the most relevant factor here. For both ISC and RISC, a similarity in Huang-Rhys factors for the molecular vibration along the π···π stacking at low-frequency region makes the rates effective. On the other side, the nonvanishing vibronic relaxation modes provide a relatively slower RISC rate in the red crystal. However, for the red crystal, small reorganization energy (λ) and large Huang-Rhys factor toward S₁ → S₀ conversion reduce nonradiative decay, leading to a prompt fluorescence. As the feasibility of S₁ ↔ T₁ conversion increases in the yellow dimer, it allows a delay in fluorescence emission, leading to thermally activated delayed fluorescence (TADF).

Singlet fission relevant energetics from optimally tuned range-separated hybrids

As a promising idea to design high-efficiency organic photovoltaics, singlet fission (SF) mechanism, i.e., generating two triplet excitons out of a single photon absorption, has recently come into the spotlight. Even though much effort has been devoted to this arena, accurately accounting for the SF process from the theoretical perspective has proven to be challenging. Herein, the SF energetics have thoroughly been investigated with the help of optimally tuned range-separated hybrid functionals (OT-RSHs) in both gas and solvent phases. Taking a series of experimentally known SF chromophores as working models, we have proposed and validated several variants of OT-RSH approximations for the reliable prediction of the energy levels which match the crucial criteria for the SF process, namely, the negative singlet-triplet and triplet-triplet energy gaps. We scrutinize the role of the OT-RSH ingredients, i.e., the underlying density functional approximations, short- and long-range exact-like exchange, as well as the range-separation parameter, for our purpose. The newly designed OT-RSHs outperform the standard RSHs and other related schemes such as screened-exchange approximations as well as other density functionals from different rungs for describing the SF energetics. More importantly, it is unveiled that although the OT-RSH coupled with the polarizable continuum model, OT-RSH-PCM, as well as the screened versions, OT-SRSHs, which account for the screening effect by the electron correlation through the scalar dielectric constant have some advantages over gas-phase computations using OT-RSHs, the energetics criteria of the SF process may not necessarily be satisfied. This in turn corroborates the idea of performing both the optimal tuning procedure and subsequent computations of the SF relevant energetics using OT-RSHs as a more reliable and affordable framework, at least for the present purpose. The applicability of the proposed models is also put into broader perspective, where they are used for the computational design of several chromophores as promising candidates prone to utilization in the SF-based materials. Hopefully, our recommended OT-RSHs can function as efficient models for both the theoretical modeling of SF chromophores and confirming the experimental observations in the field.

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.

Unveiling the role of short-range exact-like exchange in the optimally tuned range-separated hybrids for fluorescence lifetime modeling

We propose and validate several variants of the optimally tuned range-separated hybrid functionals (OT-RSHs) including different density functional approximations for predicting the fluorescence lifetimes of different categories of fluorophores within the time-dependent density functional theory (TD-DFT) framework using both the polarizable continuum and state-specific solvation models. Our main idea originates from performing the optimal tuning in the presence of a contribution of the exact-like exchange at the short-range part, which, in turn, leads to the small values of the range-separation parameter, and computing the fluorescence lifetimes using the models including no or small portions of the short-range exact-like exchange. Particular attention is also paid to the influence of the geometries of emitters on fluorescence lifetime computations. It is shown that our developed OT-RSHs along with the polarizable continuum model can be considered as the promising candidates within the TD-DFT framework for the prediction of fluorescence lifetimes for various fluorophores. We find that the proposed models not only outperform their standard counterparts but also provide reliable data better than or comparable to the conventional hybrid functionals with both the fixed and interelectronic distance-dependent exact-like exchanges. Furthermore, it is also revealed that when the excited state geometries come into play, more accurate descriptions of the fluorescence lifetimes can be achieved. Hopefully, our findings can give impetus for future developments of OT-RSHs for computational modeling of other characteristics in fluorescence spectroscopy as well as for verification of the related experimental observations.

The Sodium Anion Is Strongly Perturbed in the Condensed Phase Even Though It Appears Like a Free Ion in Nuclear Magnetic Resonance Experiments

Solvated sodium anions (Na^-) were thought to behave essentially like isolated gas-phase ions that interact only weakly with their environments. For example, ²³Na NMR signals for solvated Na^- are very sharp, despite the potential for strong quadrupolar broadening. The sharp NMR signals appear to indicate a nearly spherical electron density of the ion. For the present study, ab initio molecular dynamics simulations and quadrupolar relaxation rate calculations were carried out for the Na^-/Na⁺ [2.2.2]cryptand system solvated in methylamine, followed by detailed analyses of the electric field gradient at the sodium nuclei. It is found that Na^- does not behave like a quasi-free ion interacting only weakly with its environment. Rather, the filled 3s shell of Na^- interacts weakly with the ion's own core and the nucleus, causing Na^- to appear in NMR experiments like a free ion.

The effect of N-heterocyclic carbene units on the absorption spectra of Fe(ii) complexes: a challenge for theory

The absorption spectra of five Fe(ii) homoleptic and heteroleptic complexes containing strong sigma-donating N-heterocyclic carbene (NHC) and polypyridyl ligands have been theoretically characterized using a tuned range-separation functional.

Piecewise linearity, freedom from self-interaction, and a Coulomb asymptotic potential: three related yet inequivalent properties of the exact density functional

Three properties of the exact energy functional of DFT are important in general and for spectroscopy in particular, but are not necessarily obeyed by approximate functionals. We explain what they are, why they are important, and how they are related yet inequivalent.

Large-scale comparison of 3d and 4d transition metal complexes illuminates the reduced effect of exchange on second-row spin-state energetics

The origin of distinct 3d vs. 4d transition metal complex sensitivity to exchange is explored over a large data set.

Toward photophysical characteristics of triplet-triplet annihilation photon upconversion: a promising protocol from the perspective of optimally tuned range-separated hybrids

The photon upconversion (UC) process assisted by the triplet-triplet annihilation (TTA) mechanism has recently come into the spotlight. Given the rich collection of efforts in this area, theoretical explorations regarding TTA-UC are relatively limited and have proven to be challenging for its control in devices. In this contribution, the photophysical properties crucial for TTA-UC, such as triplet excited state energies and triplet-triplet energy transfer gaps of the essential ingredients involved in the process, namely sensitizers, annihilators and their pairs, have theoretically been investigated using optimally tuned range-separated hybrid functionals (OT-RSHs) and their screened exchange counterparts, OT-SRSHs. Taking a series of experimentally proven-to-work sensitizer/annihilator pairs as working models, we have constructed and validated several variants of OT-RSHs using both full time-dependent and Tamm-Dancoff formalisms for a reliable description of the TTA-UC photophysics. Given the bimolecular biphotonic nature of the TTA-UC process under study, particular attention is paid to the influence of the factors like the underlying density functional approximations and the tunable parameters such as short- and long-range exact-like exchanges as well as the range-separation parameter for both the sensitizers and annihilators separately. Dissecting all the aspects and relying on the appropriate choices from the tested models, we propose an OT-RSH with the correct asymptotic behavior as a cost-effective yet useful tool for this purpose. Not only against the standard RSHs but also in comparison to the conventional hybrids, the newly developed OT-RSH yields a more reliable description for the TTA-UC energetics in the gas phase and dielectric medium. Accountability of the proposed model has further been confirmed for several theoretically designed sensitizer/annihilator pairs prone to be used in the TTA-UC process. Summing up, in light of this study additional pieces of convincing evidence on the quality of OT-(S)RSHs for computational modeling and experimental verifications of the photophysics of the photon UC based on TTA and other possible technologies are showcased.

Non-empirical, low-cost recovery of exact conditions with model-Hamiltonian inspired expressions in jmDFT

Density functional theory (DFT) is widely applied to both molecules and materials, but well known energetic delocalization and static correlation errors in practical exchange-correlation approximations limit quantitative accuracy. Common methods that correct energetic delocalization errors, such as the Hubbard U correction in DFT+U or Hartree-Fock exchange in global hybrids, do so at the cost of worsening static correlation errors. We recently introduced an alternate approach [Bajaj et al., J. Chem. Phys. 147, 191101 (2017)] known as judiciously modified DFT (jmDFT), wherein the deviation from exact behavior of semilocal functionals over both fractional spin and charge, i.e., the so-called flat plane, was used to motivate functional forms of second order analytic corrections. In this work, we introduce fully nonempirical expressions for all four coefficients in a DFT+U+J-inspired form of jmDFT, where all coefficients are obtained only from energies and eigenvalues of the integer-electron systems. We show good agreement for U and J coefficients obtained nonempirically as compared with the results of numerical fitting in a jmDFT U+J/J' correction. Incorporating the fully nonempirical jmDFT correction reduces and even eliminates the fractional spin error at the same time as eliminating the energetic delocalization error. We show that this approach extends beyond s-electron systems to higher angular momentum cases including p- and d-electrons. Finally, we diagnose some shortcomings of the current jmDFT approach that limit its ability to improve upon DFT results for cases such as weakly bound anions due to poor underlying semilocal functional behavior.

Charge transfer excitation energies from ground state density functional theory calculations

Calculating charge transfer (CT) excitation energies with high accuracy and low computational cost is a challenging task. Kohn-Sham density functional theory (KS-DFT), due to its efficiency and accuracy, has achieved great success in describing ground state problems. To extend to excited state problems, our group recently demonstrated an approach with good numerical results to calculate low-lying and Rydberg excitation energies of an N-electron system from a ground state KS or generalized KS calculations of an (N - 1)-electron system via its orbital energies. In the present work, we explore further the same methodology to describe CT excitations. Numerical results from this work show that performance of conventional density functional approximations (DFAs) is not as good for CT excitations as for other excitations due to the delocalization error. Applying localized orbital scaling correction (LOSC) to conventional DFAs, a recently developed method in our group to effectively reduce the delocalization error, can improve the results. Overall, the performance of this methodology is better than time dependent DFT (TDDFT) with conventional DFAs. In addition, it shows that results from LOSC-DFAs in this method reach similar accuracy to other methods, such as ΔSCF, G₀W₀ with Bethe-Salpeter equations, particle-particle random phase approximation, and even high-level wavefunction methods like CC2. Our analysis shows that the correct 1/R trend for CT excitation can be captured from LOSC-DFA calculations, stressing that the application of DFAs with the minimal delocalization error is essential within this methodology. This work provides an efficient way to calculate CT excitation energies from ground state DFT.

Designing in the Face of Uncertainty: Exploiting Electronic Structure and Machine Learning Models for Discovery in Inorganic Chemistry

Recent transformative advances in computing power and algorithms have made computational chemistry central to the discovery and design of new molecules and materials. First-principles simulations are increasingly accurate and applicable to large systems with the speed needed for high-throughput computational screening. Despite these strides, the combinatorial challenges associated with the vastness of chemical space mean that more than just fast and accurate computational tools are needed for accelerated chemical discovery. In transition-metal chemistry and catalysis, unique challenges arise. The variable spin, oxidation state, and coordination environments favored by elements with well-localized d or f electrons provide great opportunity for tailoring properties in catalytic or functional (e.g., magnetic) materials but also add layers of uncertainty to any design strategy. We outline five key mandates for realizing computationally driven accelerated discovery in inorganic chemistry: (i) fully automated simulation of new compounds, (ii) knowledge of prediction sensitivity or accuracy, (iii) faster-than-fast property prediction methods, (iv) maps for rapid chemical space traversal, and (v) a means to reveal design rules on the kilocompound scale. Through case studies in open-shell transition-metal chemistry, we describe how advances in methodology and software in each of these areas bring about new chemical insights. We conclude with our outlook on the next steps in this process toward realizing fully autonomous discovery in inorganic chemistry using computational chemistry.

Beyond Koopmans' theorem: electron binding energies in disordered materials

The topical review focuses on calculating ionization energies (IE), or electronic polarons in quasi-particle terminology, in large disordered systems, e.g. for a solute dissolved in a molecular solvent. The simplest estimate of the ionization energy is provided by one-electron energies in the Hartree-Fock theory, but the calculated quantities are not accurate. Density functional theory as many-body theory provides a principal opportunity for calculating one-electron energies including correlation and relaxation effects, i.e. the true energies of electronic polarons. We argue that such a principal possibility materializes within the concept of optimally tuned range-separated hybrid functionals (OT-RSH). We describe various schemes for optimal tuning. Importantly, the OT-RSH scheme is investigated for systems capped with dielectric continuum, providing a consistent picture on the QM/dielectric boundary. Finally, some limitations and open issues of the OT-RSH approach are addressed.

Delocalization Errors in Density Functional Theory Are Essentially Quadratic in Fractional Occupation Number

Approximate functionals used in practical density functional theory (DFT) deviate from the piecewise linear behavior of the exact functional for fractional charges. This deviation causes excess charge delocalization, which leads to incorrect densities, molecular properties, barrier heights, band gaps, and excitation energies. We present a simple delocalization function for characterizing this error and find it to be almost perfectly linear vs the fractional electron number for systems spanning in size from the H atom to the C₁₂H₁₄ polyene. This causes the delocalization energy error to be a quadratic polynomial in the fractional electron number, which permits us to assess the comparative performance of 47 popular and recent functionals through the curvature. The quadratic form further suggests that information about a single fractional charge is sufficient to eliminate the principal source of delocalization error. Generalizing traditional two-point information like ionization potentials or electron affinities to account for a third, fractional charge-based data point could therefore permit fitting/tuning of functionals with lower delocalization error.

Dielectric Screening Meets Optimally Tuned Density Functionals

A short overview of recent attempts at merging two independently developed methods is presented. These are the optimal tuning of a range-separated hybrid (OT-RSH) functional, developed to provide an accurate first-principles description of the electronic structure and optical properties of gas-phase molecules, and the polarizable continuum model (PCM), developed to provide an approximate but computationally tractable description of a solvent in terms of an effective dielectric medium. After a brief overview of the OT-RSH approach, its combination with the PCM as a potentially accurate yet low-cost approach to the study of molecular assemblies and solids, particularly in the context of photocatalysis and photovoltaics, is discussed. First, solvated molecules are considered, with an emphasis on the challenge of balancing eigenvalue and total energy trends. Then, it is shown that the same merging of methods can also be used to study the electronic and optical properties of molecular solids, with a similar discussion of the pros and cons. Tuning of the effective scalar dielectric constant as one recent approach that mitigates some of the difficulties in merging the two approaches is considered.

Describing strong correlation with fractional-spin correction in density functional theory

An effective fractional-spin correction is developed to describe static/strong correlation in density functional theory. Combined with the fractional-charge correction from recently developed localized orbital scaling correction (LOSC), a functional, the fractional-spin LOSC (FSLOSC), is proposed. FSLOSC, a correction to commonly used functional approximations, introduces the explicit derivative discontinuity and largely restores the flat-plane behavior of electronic energy at fractional charges and fractional spins. In addition to improving results from conventional functionals for the prediction of ionization potentials, electron affinities, quasiparticle spectra, and reaction barrier heights, FSLOSC properly describes the dissociation of ionic species, single bonds, and multiple bonds without breaking space or spin symmetry and corrects the spurious fractional-charge dissociation of heteroatom molecules of conventional functionals. Thus, FSLOSC demonstrates success in reducing delocalization error and including strong correlation, within low-cost density functional approximation.

Communication: Recovering the flat-plane condition in electronic structure theory at semi-local DFT cost

The flat-plane condition is the union of two exact constraints in electronic structure theory: (i) energetic piecewise linearity with fractional electron removal or addition and (ii) invariant energetics with change in electron spin in a half filled orbital. Semi-local density functional theory (DFT) fails to recover the flat plane, exhibiting convex fractional charge errors (FCE) and concave fractional spin errors (FSE) that are related to delocalization and static correlation errors. We previously showed that DFT+U eliminates FCE but now demonstrate that, like other widely employed corrections (i.e., Hartree-Fock exchange), it worsens FSE. To find an alternative strategy, we examine the shape of semi-local DFT deviations from the exact flat plane and we find this shape to be remarkably consistent across ions and molecules. We introduce the judiciously modified DFT (jmDFT) approach, wherein corrections are constructed from few-parameter, low-order functional forms that fit the shape of semi-local DFT errors. We select one such physically intuitive form and incorporate it self-consistently to correct semi-local DFT. We demonstrate on model systems that jmDFT represents the first easy-to-implement, no-overhead approach to recovering the flat plane from semi-local DFT.

Dipole moments of molecules with multi-reference character from optimally tuned range-separated density functional theory

Dipole moment is the first nonzero moment of the charge density of neutral systems. If a density functional theory (DFT) method is able to yield accurate dipole moments, it should first provide an accurate geometry and then predict a reliable charge distribution for that geometry. In this respect, recent literatures have revealed that most DFT approximations work considerably better for single-reference molecules with respect to multi-reference ones, as may be expected from this fact that DFT utilizes a single configuration state function as reference function to represent the density. Putting together, it seems that as compared to the single-reference systems, situation is slightly more involved in the case of dipole moment calculations of multi-reference molecules. Effort to address this latter issue constitutes the cornerstone of the present investigation. To this end, we rely on a different approach where the new optimally (nonempirically) tuned range-separated hybrid density functionals (OT-RSHs) without invoking any empirical fitting are proposed for predicting the dipole moments of multi-reference molecules containing both main-group elements and transition metals. We have scanned the controlling factors of OT-RSHs like short- and long-range exchange contributions and range-separation parameter with the aim of deriving the best performing models for the purpose. The obtained results unveil that, as compared to the standard range-separated density functionals, our newly developed OT-RSHs not only give an improved description on the dipole moments of the molecules with multi-reference character but also the quality of their predictions is better than other conventional and recently proposed DFT approximations. © 2018 Wiley Periodicals, Inc.