Prediction of electronic couplings for molecular charge transfer using optimally tuned range-separated hybrid functionals

Benchmarking DFT-based excited-state methods for intermolecular charge-transfer excitations

Intermolecular charge-transfer is a highly important process in biology and energy-conversion applications where generated charges need to be transported over several moieties. However, its theoretical description is challenging since the high accuracy required to describe these excited states must be accessible for calculations on large molecular systems. In this benchmark study, we identify reliable low-scaling computational methods for this task. Our reference results were obtained from highly accurate wavefunction calculations that restrict the size of the benchmark systems. However, the density-functional theory based methods that we identify as accurate can be applied to much larger systems. Since targeting charge-transfer states requires the unambiguous classification of an excited state, we first analyze several charge-transfer descriptors for their reliability concerning intermolecular charge-transfer and single out the charge-transfer distance calculated based on the variation of electron density upon excitation (DCT) as an optimal choice for our purposes. In general, best results are obtained for orbital-optimized methods and among those, the maximum overlap method proved to be the most numerically stable variant when using the initial MOs as reference orbitals. Favorable error cancellation with optimally-tuned range-separated hybrid functionals and a rather small basis set can provide an economical yet reasonable wavefunction when using time-dependent density functional theory, which provides relevant information about the excited-state character to be used in the orbital-optimized methods. The qualitative agreement makes these fast calculations attractive for high-throughput screening applications.

Evaluating the interactions between vibrational modes and electronic transitions using frontier orbital energy derivatives

Vibrations affect molecular optoelectronic properties, even at zero kelvin. Accounting for these effects using computational modelling is costly, as it requires many calculations at geometries distorted from equilibrium. Here, we propose a low-cost method for identifying vibrations most strongly coupled to the electronic structure, based on using orbital energy derivatives as a diagnostic.

Electron Localization and Mobility in Monolayer Fullerene Networks

The novel 2D quasi-hexagonal phase of covalently bonded fullerene molecules (qHP C60), the so-called graphullerene, has displayed far superior electron mobilities, if compared to the parent van der Waals three-dimensional crystal (vdW C60). Herein, we present a comparative study of the electronic properties of vdW and qHP C60 using state-of-the-art electronic-structure calculations and a full quantum-mechanical treatment of electron transfer. We show that both materials entail polaronic localization of electrons with similar binding energies (≈0.1 eV) and, therefore, they share the same charge transport via polaron hopping. In fact, we quantitatively reproduce the sizable increment of the electron mobility measured for qHP C60 and identify its origin in the increased electronic coupling between C60 units.

Temperature-transferable tight-binding model using a hybrid-orbital basis

Finite-temperature calculations are relevant for rationalizing material properties, yet they are computationally expensive because large system sizes or long simulation times are typically required. Circumventing the need for performing many explicit first-principles calculations, tight-binding and machine-learning models for the electronic structure emerged as promising alternatives, but transferability of such methods to elevated temperatures in a data-efficient way remains a great challenge. In this work, we suggest a tight-binding model for efficient and accurate calculations of temperature-dependent properties of semiconductors. Our approach utilizes physics-informed modeling of the electronic structure in the form of hybrid-orbital basis functions and numerically integrating atomic orbitals for the distance dependence of matrix elements. We show that these design choices lead to a tight-binding model with a minimal amount of parameters that are straightforwardly optimized using density functional theory or alternative electronic-structure methods. The temperature transferability of our model is tested by applying it to existing molecular-dynamics trajectories without explicitly fitting temperature-dependent data and comparison with density functional theory. We utilize it together with machine-learning molecular dynamics and hybrid density functional theory for the prototypical semiconductor gallium arsenide. We find that including the effects of thermal expansion on the onsite terms of the tight-binding model is important in order to accurately describe electronic properties at elevated temperatures in comparison with experiment.

Ultrafast estimation of electronic couplings for electron transfer between pi-conjugated organic molecules. II

The development of highly efficient methods for the calculation of electronic coupling matrix elements between the electron donor and acceptor is an important goal in theoretical organic semiconductor research. In Paper I [F. Gajdos, S. Valner, F. Hoffmann, J. Spencer, M. Breuer, A. Kubas, M. Dupuis, and J. Blumberger, J. Chem. Theory Comput. 10, 4653 (2014)], we introduced the analytic overlap method (AOM) for this purpose, which is an ultrafast electronic coupling estimator parameterized to and orders of magnitude faster than density functional theory (DFT) calculations at a reasonably small loss in accuracy. In this work, we reparameterize and extend the AOM to molecules containing nitrogen, oxygen, fluorine, and sulfur heteroatoms using 921 dimer configurations from the recently introduced HAB79 dataset. We find again a very good linear correlation between the frontier orbital overlap, calculated ultrafast in an optimized minimum Slater basis, and DFT reference electronic couplings. The new parameterization scheme is shown to be transferable to sulfur-containing polyaromatic hydrocarbons in experimentally resolved dimeric configurations. Our extension of the AOM enables high-throughput screening of very large databases of chemically diverse organic crystal structures and the application of computationally intense non-adiabatic molecular dynamics methods to charge transport in state-of-the-art organic semiconductors, e.g., non-fullerene acceptors.

HAB79: A new molecular dataset for benchmarking DFT and DFTB electronic couplings against high-level ab initio calculations

A new molecular dataset called HAB79 is introduced to provide ab initio reference values for electronic couplings (transfer integrals) and to benchmark density functional theory (DFT) and density functional tight-binding (DFTB) calculations. The HAB79 dataset is composed of 79 planar heterocyclic polyaromatic hydrocarbon molecules frequently encountered in organic (opto)electronics, arranged to 921 structurally diverse dimer configurations. We show that CASSCF/NEVPT2 with a minimal active space provides a robust reference method that can be applied to the relatively large molecules of the dataset. Electronic couplings are largest for cofacial dimers, in particular, sulfur-containing polyaromatic hydrocarbons, with values in excess of 0.5 eV, followed by parallel displaced cofacial dimers. V-shaped dimer motifs, often encountered in the herringbone layers of organic crystals, exhibit medium-sized couplings, whereas T-shaped dimers have the lowest couplings. DFT values obtained from the projector operator-based diabatization (POD) method are initially benchmarked against the smaller databases HAB11 (HAB7-) and found to systematically improve when climbing Jacob's ladder, giving mean relative unsigned errors (MRUEs) of 27.7% (26.3%) for the generalized gradient approximation (GGA) functional BLYP, 20.7% (15.8%) for hybrid functional B3LYP, and 5.2% (7.5%) for the long-range corrected hybrid functional omega-B97X. Cost-effective POD in combination with a GGA functional and very efficient DFTB calculations on the dimers of the HAB79 database give a good linear correlation with the CASSCF/NEVPT2 reference data, which, after scaling with a multiplicative constant, gives reasonably small MRUEs of 17.9% and 40.1%, respectively, bearing in mind that couplings in HAB79 vary over 4 orders of magnitude. The ab initio reference data reported here are expected to be useful for benchmarking other DFT or semi-empirical approaches for electronic coupling calculations.

Stacked Ensemble Machine Learning for Range-Separation Parameters

Density functional theory-based high-throughput materials and drug discovery has achieved tremendous success in recent decades, but its power on organic semiconducting molecules suffered catastrophically from the self-interaction error until the nonempirical but expensive optimally tuned range-separated hybrid (OT-RSH) functionals were developed. An OT-RSH transitions from a short-range (semi)local functional to a long-range Hartree-Fock exchange at a distance characterized by a molecule-specific range-separation parameter (ω). Herein, we propose a stacked ensemble machine learning model that provides an accelerated alternative of OT-RSH based on system-dependent structural and electronic configurations. We trained ML-ωPBE, the first functional in our series, using a database of 1970 molecules with sufficient structural and functional diversity, and assessed its accuracy and efficiency using another 1956 molecules. Compared with nonempirical OT-ωPBE, ML-ωPBE reaches a mean absolute error of 0.00504a₀^-1 for optimal ω's, reduces the computational cost by 2.66 orders of magnitude, and achieves comparable predictive power in optical properties.

Rational design of organic semiconductors with low internal reorganization energies for hole and electron transport: position effect of aza-substitution in phenalenyl derivatives

Amphoteric-redox phenalenyl radical (PLY) is a suitable candidate used to design ambipolar organic materials. Because the singly occupied nonbonding molecular orbital (NBMO) of PLY has a perfect local nonbonding character, its internal reorganization energy (λ) for transporting holes (λ⁺) or electrons (λ^-) is known to be small. Herein, PLY is employed to study the position effect of the aza group on the λ. By adding or extracting an electron from the NBMO, the bond length alterations can be minute. Therefore, the PLY derivatives are also an excellent candidate to study the contributions from the bond angle alterations to the λ. Substituting the aza groups at the β- or α-positions of PLY shows two different trends. When consecutively substituting the aza group at the three β-positions of PLY, the λs are consistently decreased. Contrarily, a series of double functionalization of aza groups at the four α-positions of PLY, the λs are increased. It is because the local bonding or antibonding character in frontier orbitals (FMO) is observed in α2N-PLY and α4N-PLY. As the FMOs of the three β-substituted PLYs and α6N-PLY have perfect local nonbonding character, we found the bond angle alterations are the main contributors of λ. The λs for most aza-PLYs were smaller than 100 meV. Thus, we propose a design rule for substituting aza groups on the parent molecules with strong local nonbonding character in their FMOs. Based on the adiabatic ionization potential and electron affinity, two π-extended PLY derivatives with small λ were recommended for fabricating air-stable ambipolar OFET.

Examination of How Well Long-Range-Corrected Density Functionals Satisfy the Ionization Energy Theorem

We calculated the vertical ionization energies (VIE) of 99 species in two ways to examine the accuracy of several long-range-corrected (LC) hybrid meta functionals in comparison with a gradient approximation (GA), global hybrids, and doubly hybrids. In the category of LC functionals, we examined both those with meta ingredients (i.e., that depend on the kinetic energy density) and those without them. The LC-hybrid meta functionals examined are M11, revM11, M11plus, and ωB97M-V. The reference data used to assess accuracy consist of 95 molecules and 4 atoms in the GW100 set. The two methods studied are the ΔSCF method (involving the difference of neutral and cation self-consistent field (SCF) energies) and the ionization energy theorem (involving the orbital energy of the highest occupied molecular orbital, HOMO). We calculated linear correlation coefficients (r²) and mean absolute deviations (MADs) between each approach and the reference VIE value from the CCSD(T)/def2-TZVPP level of theory. We compared the new LC-hybrid meta calculations to calculations with the 10 functionals in a previous VIE study by Brémond et al. and to the calculations with LC-BLYP (LC-Becke, Lee-Yang-Parr), CAM-B3LYP (Coulomb-attenuating-method Becke-3-parameter Lee-Yang-Parr), LC-ωHPBE, and ωB97X-D. The results show that Minnesota LC-hybrid meta functionals have the smallest mean absolute deviation of ionization energy theorem VIEs with the reference data; the LC-ωHPBE functional also does quite well in this test. This is very encouraging and indicates that LC-hybrid meta functionals would be the best starting points for the tuning strategy that has been shown to be a very good procedure for improving time-dependent density functional calculations, and it also helps explain the good success of LC-hybrid meta functionals for molecular excitation energies.

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.

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.

Electron transfer in extended systems: characterization by periodic density functional theory including the electronic coupling

We describe a new computer implementation of electron transfer (ET) theory in extended systems treated by periodic density functional theory (DFT), including the calculation of the electronic coupling transition element V_AB. In particular, the development opens up the full characterization of electron transfer in the solid state. The approach is valid for any single-determinant wavefunction with localized character representing the electronic structure of the system, from Hartree-Fock (HF) theory, to density functional theory (DFT), hybrid DFT theory, DFT+U theory, and constrained DFT (cDFT) theory. The implementation in CP2K reuses the high-performance functions of the code. The computational cost is equivalent to only one iteration of an HF calculation. We present test calculations for electron transfer in a number of systems, including a 1D-model of ferric oxide, hematite Fe₂O₃, rutile TiO₂, and finally bismuth vanadate BiVO₄.

How To Extract Quantitative Information on Electronic Transitions from the Density Functional Theory "Black Box"

Electronic couplings and vertical excitation energies are crucial determinants of charge and excitation energy transfer rates in a broad variety of processes ranging from biological charge transfer to charge transport through inorganic materials, from molecular sensing to intracellular signaling. Density Functional Theory (DFT) is generally used to calculate these critical parameters, but the quality of the results is unpredictable because of the semiempirical nature of the available DFT approaches. This study identifies a small set of fundamental rules that enables accurate DFT computation of electronic couplings and vertical excitation energies in molecular complexes and materials. These rules are applied to predict efficient DFT approaches to coupling calculations. The result is an easy-to-use guide for reliable DFT descriptions of electronic transitions.