Critical analysis of fragment-orbital DFT schemes for the calculation of electronic coupling values

Efficient calculation of electronic coupling integrals with the dimer projection method via a density matrix tight-binding potential

Designing organic semiconductors for practical applications in organic solar cells, organic field-effect transistors, and organic light-emitting diodes requires understanding charge transfer mechanisms across different length and time scales. The underlying electron transfer mechanisms can be efficiently explored using semiempirical quantum mechanical (SQM) methods. The dimer projection (DIPRO) method combined with the recently introduced non-self-consistent density matrix tight-binding potential (PTB) [Grimme et al., J. Chem. Phys. 158, 124111 (2023)] is used in this study to evaluate charge transfer integrals important for understanding charge transport mechanisms. PTB, parameterized for the entire Periodic Table up to Z = 86, incorporates approximate non-local exchange, allowing for efficient and accurate calculations for large hetero-organic compounds. Benchmarking against established databases, such as Blumberger's HAB sets, or our newly introduced JAB69 set and comparing with high-level reference data from ωB97X-D4 calculations confirm that DIPRO@PTB consistently performs well among the tested SQM approaches for calculating coupling integrals. DIPRO@PTB yields reasonably accurate results at low computational cost, making it suitable for screening purposes and applications to large systems, such as metal-organic frameworks and cyanine-based molecular aggregates further discussed in this work.

Impact of Fluoroalkylation on the n-Type Charge Transport of Two Naphthodithiophene Diimide Derivatives

In this work, we investigate two recently synthesized naphthodithiophene diimide (NDTI) derivatives featuring promising n-type charge transport properties. We analyze the charge transport pathways and model charge mobility with the non-adiabatic hopping mechanism using the Marcus-Levich-Jortner rate constant formulation, highlighting the role of fluoroalkylated substitution in α (α-NDTI) and at the imide nitrogen (N-NDTI) position. In contrast with the experimental results, similar charge mobilities are computed for the two derivatives. However, while α-NDTI displays remarkably anisotropic mobilities with an almost one-dimensional directionality, N-NDTI sustains a more isotropic charge percolation pattern. We propose that the strong anisotropic charge transport character of α-NDTI is responsible for the modest measured charge mobility. In addition, when the role of thermally induced transfer integral fluctuations is investigated, the computed electron-phonon couplings for intermolecular sliding modes indicate that dynamic disorder effects are also more detrimental for the charge transport of α-NDTI than N-NDTI. The lower observed mobility of α-NDTI is therefore rationalized in terms of a prominent anisotropic character of the charge percolation pathways, with the additional contribution of dynamic disorder effects.

How the Donor/Acceptor Spin States Affect the Electronic Couplings in Molecular Charge-Transfer Processes?

The electronic coupling matrix element H_AB is an essential ingredient of most electron-transfer theories. H_AB depends on the overlap between donor and acceptor wave functions and is affected by the involved states’ spin. We classify the spin-state effects into three categories: orbital occupation, spin-dependent electron density, and density delocalization. The orbital occupancy reflects the diverse chemical nature and reactivity of the spin states of interest. The effect of spin-dependent density is related to a more compact electron density cloud at lower spin states due to decreased exchange interactions between electrons. Density delocalization is strongly connected with the covalency concept that increases the spatial extent of the diabatic state’s electron density in specific directions. We illustrate these effects with high-level ab initio calculations on model direct donor–acceptor systems relevant to metal oxide materials and biological electron transfer. Obtained results can be used to benchmark existing methods for H_AB calculations in complicated cases such as spin-crossover materials or antiferromagnetically coupled systems.

Excited state diabatization on the cheap using DFT: Photoinduced electron and hole transfer

Excited state electron and hole transfer underpin fundamental steps in processes such as exciton dissociation at photovoltaic heterojunctions, photoinduced charge transfer at electrodes, and electron transfer in photosynthetic reaction centers. Diabatic states corresponding to charge or excitation localized species, such as locally excited and charge transfer states, provide a physically intuitive framework to simulate and understand these processes. However, obtaining accurate diabatic states and their couplings from adiabatic electronic states generally leads to inaccurate results when combined with low-tier electronic structure methods, such as time-dependent density functional theory, and exorbitant computational cost when combined with high-level wavefunction-based methods. Here, we introduce a density functional theory (DFT)-based diabatization scheme that directly constructs the diabatic states using absolutely localized molecular orbitals (ALMOs), which we denote as Δ-ALMO(MSDFT2). We demonstrate that our method, which combines ALMO calculations with the ΔSCF technique to construct electronically excited diabatic states and obtains their couplings with charge-transfer states using our MSDFT2 scheme, gives accurate results for excited state electron and hole transfer in both charged and uncharged systems that underlie DNA repair, charge separation in donor-acceptor dyads, chromophore-to-solvent electron transfer, and singlet fission. This framework for the accurate and efficient construction of excited state diabats and evaluation of their couplings directly from DFT thus offers a route to simulate and elucidate photoinduced electron and hole transfer in large disordered systems, such as those encountered in the condensed phase.

Electron transfer and spin dynamics of the radical-pair in the cryptochrome from Chlamydomonas reinhardtii by computational analysis

In an effort to elucidate the origin of avian magnetoreception, it was postulated that a radical-pair formed in a cryptochrome upon light activation provided the basis for the mechanism that enables an inclination compass sensitive to the geomagnetic field. Photoreduction in this case involves formation of a flavin adenine dinucleotide (FAD)-tryptophan (TRP) radical-pair, following electron transfer within a conserved TRP triad in the cryptochrome. Recently, an animal-like cryptochrome from Chlamydomonas reinhardtii (CraCRY) was analyzed, demonstrating the role of a fourth aromatic residue, which serves as a terminal electron donor in the photoreduction pathway, resulting in the creation of a more distal radical-pair and exhibiting fast electron transfer. In this work, we investigated the electron transfer in CraCRY with a combination of free energy molecular dynamics (MD) simulations, frozen density functional theory, and QM/MM MD simulations, supporting the suggestion of a proton coupled electron transfer mechanism. Spin dynamics simulations discerned details on the dependence of the singlet yield on the direction of the external magnetic field for the [FAD^•- TYRH^•+] and [FAD^•- TYR^•] radical-pairs in CraCRY, in comparison with the previously modeled [FAD^•- TRPH^•+] radical-pair.

An Iterative Fragment Scheme for the ACKS2 Electronic Polarization Model: Application to Molecular Dimers and Chains

The treatment of electrostatic interactions is a key ingredient in the force field-based simulation of condensed phase systems. Most approaches used fixed, site-specific point charges. Yet, it is now clear that many applications of force fields (FFs) demand more sophisticated treatments, prompting the implementation of charge equilibration methods in polarizable FFs to allow the redistribution of charge within the system. One approach allowing both, charge redistribution and site-specific polarization, while at the same time solving methodological shortcomings of earlier methods, is the first-principles-derived atom-condensed Kohn-Sham density functional theory method approximated to the second order (ACKS2). In this work, we present two fragment approaches to ACKS2, termed f-ACKS2 and a self-consistent version, scf-ACKS2, that treat condensed phase systems as a collection of electronically polarizable molecular fragments. The fragmentation approach to ACKS2 not only leads to a more transferable and less system-specific collection of electronic response parameters but also opens up the method to large condensed phase systems. We validate the accuracies of f-ACKS2 and scf-ACKS2 by comparing polarization energies and induced dipole moments for a number of charged hydrocarbon dimers against DFT reference calculations. Finally, we also apply both fragmented ACKS2 variants to calculate the polarization energy for electron-hole pair separation along a chain of anthracene molecules and find excellent agreement with reference DFT calculations.

PyCDFT: A Python package for constrained density functional theory

We present PyCDFT, a Python package to compute diabatic states using constrained density functional theory (CDFT). PyCDFT provides an object-oriented, customizable implementation of CDFT, and allows for both single-point self-consistent-field calculations and geometry optimizations. PyCDFT is designed to interface with existing density functional theory (DFT) codes to perform CDFT calculations where constraint potentials are added to the Kohn-Sham Hamiltonian. Here, we demonstrate the use of PyCDFT by performing calculations with a massively parallel first-principles molecular dynamics code, Qbox, and we benchmark its accuracy by computing the electronic coupling between diabatic states for a set of organic molecules. We show that PyCDFT yields results in agreement with existing implementations and is a robust and flexible package for performing CDFT calculations. The program is available at https://dx.doi.org/10.5281/zenodo.3821097.

Unifying machine learning and quantum chemistry with a deep neural network for molecular wavefunctions

Machine learning advances chemistry and materials science by enabling large-scale exploration of chemical space based on quantum chemical calculations. While these models supply fast and accurate predictions of atomistic chemical properties, they do not explicitly capture the electronic degrees of freedom of a molecule, which limits their applicability for reactive chemistry and chemical analysis. Here we present a deep learning framework for the prediction of the quantum mechanical wavefunction in a local basis of atomic orbitals from which all other ground-state properties can be derived. This approach retains full access to the electronic structure via the wavefunction at force-field-like efficiency and captures quantum mechanics in an analytically differentiable representation. On several examples, we demonstrate that this opens promising avenues to perform inverse design of molecular structures for targeting electronic property optimisation and a clear path towards increased synergy of machine learning and quantum chemistry.

Machine learning models can accurately predict atomistic chemical properties but do not provide access to the molecular electronic structure. Here the authors use a deep learning approach to predict the quantum mechanical wavefunction at high efficiency from which other ground-state properties can be derived.

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

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

Genarris: Random generation of molecular crystal structures and fast screening with a Harris approximation

We present Genarris, a Python package that performs configuration space screening for molecular crystals of rigid molecules by random sampling with physical constraints. For fast energy evaluations, Genarris employs a Harris approximation, whereby the total density of a molecular crystal is constructed via superposition of single molecule densities. Dispersion-inclusive density functional theory is then used for the Harris density without performing a self-consistency cycle. Genarris uses machine learning for clustering, based on a relative coordinate descriptor developed specifically for molecular crystals, which is shown to be robust in identifying packing motif similarity. In addition to random structure generation, Genarris offers three workflows based on different sequences of successive clustering and selection steps: the "Rigorous" workflow is an exhaustive exploration of the potential energy landscape, the "Energy" workflow produces a set of low energy structures, and the "Diverse" workflow produces a maximally diverse set of structures. The latter is recommended for generating initial populations for genetic algorithms. Here, the implementation of Genarris is reported and its application is demonstrated for three test cases.

Fragment-orbital tunneling currents and electronic couplings for analysis of molecular charge-transfer systems

In theoretical charge-transfer research, calculation of the electronic coupling element is crucial for examining the degree of the electronic donor-acceptor interaction. The tunneling current (TC), representing the magnitudes and directions of electron flow, provides a way of evaluating electronic couplings, along with the ability of visualizing how electrons flow in systems. Here, we applied the TC theory to π-conjugated organic dimer systems, in the form of our fragment-orbital tunneling current (FOTC) method, which uses the frontier molecular-orbitals of system fragments as diabatic states. For a comprehensive test of FOTC, we assessed how reasonable the computed electronic couplings and the corresponding TC densities are for the hole- and electron-transfer databases HAB11 and HAB7. FOTC gave 12.5% mean relative unsigned error with regard to the high-level ab initio reference. The shown performance is comparable with that of fragment-orbital density functional theory, which gave the same error by 20.6% or 13.9% depending on the formulation. In the test of a set of nucleobase π stacks, we showed that the original TC expression is also applicable to nondegenerate cases under the condition that the overlap between the charge distributions of diabatic states is small enough to offset the energy difference. Lastly, we carried out visual analysis on the FOTC densities of thiophene dimers with different intermolecular alignments. The result depicts an intimate topological connection between the system geometry and electron flow. Our work provides quantitative and qualitative grounds for FOTC, showing it to be a versatile tool in characterization of molecular charge-transfer systems.

Density Functional Physicality in Electronic Coupling Estimation: Benchmarks and Error Analysis

Electronic coupling estimates from constrained density functional theory configuration interaction (CDFT-CI) depend critically on choice of density functional. In this Letter, the orbital multielectron self-interaction error (OMSIE), vertical electron affinity (VEA), and vertical ionization potential (VIP) are shown to be the key indicators inherited from the density functional that determine the accuracy of electronic coupling estimates. An error metric η is derived to connect the three properties, based on the linear proportionality between electronic coupling and overlap integral, and the hypothesis that the slope of this line is a function of VEA/VIP, η = (1/N_testset)Σ_i^testset|-VE^Ref × OMSIE + ΔVE - ΔVE × OMSIE|_i. Based on η, BH&HLYP and LRC-ωPBEh are suggested as the best functionals for electron and hole transfer, respectively. Error metric η is therefore a useful predictor of errors in CDFT-CI electronic coupling, showing that the physical correctness of the density functional has a direct effect on the accuracy of the electronic coupling.

Charge Transport in Molecular Materials: An Assessment of Computational Methods

The booming field of molecular electronics has fostered a surge of computational research on electronic properties of organic molecular solids. In particular, with respect to a microscopic understanding of transport and loss mechanisms, theoretical studies assume an ever-increasing role. Owing to the tremendous diversity of organic molecular materials, a great number of computational methods have been put forward to suit every possible charge transport regime, material, and need for accuracy. With this review article we aim at providing a compendium of the available methods, their theoretical foundations, and their ranges of validity. We illustrate these through applications found in the literature. The focus is on methods available for organic molecular crystals, but mention is made wherever techniques are suitable for use in other related materials such as disordered or polymeric systems.

Virtual Screening for High Carrier Mobility in Organic Semiconductors

Low carrier mobilities still hamper the use of organic semiconductors in many applications. Using a staged virtual screening approach we compute the electronic couplings and intramolecular reorganization energies as two main descriptors for charge mobility for a set of 95 445 molecular crystals extracted from the Cambridge Structural Database (CSD). As a final step, on the basis of the calculated coupling values, we identify materials with long-range charge percolation pathways. Thus we readily find many acclaimed compounds as well as a number of most promising materials that have not yet been considered for an application in organic electronics. Together with the unique metadata provided in the CSD, the large descriptor database allows us to extract important trends and correlations that will further accelerate the theoretical design and discovery of high mobility organic semiconductors.

Effect of metallation, substituents and inter/intra-molecular polarization on electronic couplings for hole transport in stacked porphyrin dyads

Hole transport properties for stacked porphyrins are analyzed with density functional theory.