Performance of Frozen Density Embedding for Modeling Hole Transfer Reactions

Evolution of the ionisation energy with the stepwise growth of chiral clusters of [4]helicene

Polycyclic aromatic hydrocarbons (PAHs) are widely established as ubiquitous in the interstellar medium (ISM), but considering their prevalence in harsh vacuum environments, the role of ionisation in the formation of PAH clusters is poorly understood, particularly if a chirality-dependent aggregation route is considered. Here we report on photoelectron spectroscopy experiments on [4]helicene clusters performed with a vacuum ultraviolet synchrotron beamline. Aggregates (up to the heptamer) of [4]helicene, the smallest PAH with helical chirality, were produced and investigated with a combined experimental and theoretical approach using several state-of-the-art quantum-chemical methodologies. The ionisation onsets are extracted for each cluster size from the mass-selected photoelectron spectra and compared with calculations of vertical ionisation energies. We explore the complex aggregation topologies emerging from the multitude of isomers formed through clustering of P and M, the two enantiomers of [4]helicene. The very satisfactory benchmarking between experimental ionisation onsets vs. predicted ionisation energies allows the identification of theoretically predicted potential aggregation motifs and corresponding energetic ordering of chiral clusters. Our structural models suggest that a homochiral aggregation route is energetically favoured over heterochiral arrangements with increasing cluster size, hinting at potential symmetry breaking in PAH cluster formation at the scale of small grains.

Time-Dependent Multilevel Density Functional Theory

We present a novel three-layer approach based on multilevel density functional theory (MLDFT) and polarizable molecular mechanics to simulate the electronic excitations of chemical systems embedded in an external environment within the time-dependent DFT formalism. In our method, the electronic structure of a target system, the chromophore, is determined in the field of an embedded inactive layer, which is treated as frozen. Long-range interactions are described by employing the polarizable fluctuating charge (FQ) force field. The resulting MLDFT/FQ thus accurately describes both electrostatics (and polarization) and non-electrostatic target-environment interactions. The robustness and reliability of the approach are demonstrated by comparing our results with experimental data reported for various organic molecules in solution.

High Electronic Coupling between Cu Complexes and Oxidized Dyes Confirmed by Measurements of Driving Force Dependent Regeneration Kinetics in Minimal Electrolyte System

We confirm fast regeneration kinetics between copper complexes and oxidized organic dyes and the major contribution of electronic coupling (HDA). The highest efficiency of dye-sensitized TiO2 solar cells has been shown by employing Cu complex redox couples. Various groups have reported a fast regeneration rate of oxidized dyes by Cu complexes giving a low driving force attributed to low reorganization energy (λ), but the effect of HDA has not been evaluated. The values of HDA and λ can be derived from driving force dependent transient absorption (TA) measurements. However, analyzing TA decay using Cu complexes is not trivial because accelerated recombination by the presence of Cu2+ complexes and biphasic TA decay often complicates the analysis. Here we employ 16 Cu1+ and Co2+ complexes and two dyes. To simplify the system, i.e., making a minimal electrolyte system, Cu2+ and Co3+ complexes and a common additive of 4-tert-butylpyridine are not used. From the driving force dependent TA decays of oxidized dyes by both Cu1+ and Co2+ complexes, λ for the combination of the Cu complexes and dyes is found to be about 0.15 eV lower than that of Co complexes. Approximately 3 to 5 times higher HDA values of Cu complexes than those of Co complexes are obtained, which is the dominant factor for faster rates. The values vary with the structure of the molecules, showing the possibility of increasing the HDA values further. The higher HDA values of a Cu complex than that of a Co complex are also reproduced by quantum chemical calculations.

Multiscale Frozen Density Embedding/Molecular Mechanics Approach for Simulating Magnetic Response Properties of Solvated Systems

We present a three-layer hybrid quantum mechanical/quantum embedding/molecular mechanics approach for calculating nuclear magnetic resonance (NMR) shieldings and J-couplings of molecular systems in solution. The model is based on the frozen density embedding (FDE) and polarizable fluctuating charges (FQ) and fluctuating dipoles (FQFμ) force fields and permits the accurate ab initio description of short-range nonelectrostatic interactions by means of the FDE shell and cost-effective treatment of long-range electrostatic interactions through the polarizable force field FQ(Fμ). Our approach's accuracy and potential are demonstrated by studying NMR spectra of Brooker's merocyanine in aqueous and nonaqueous solutions.

Orbital-Free Density Functional Theory: An Attractive Electronic Structure Method for Large-Scale First-Principles Simulations

Kohn-Sham Density Functional Theory (KSDFT) is the most widely used electronic structure method in chemistry, physics, and materials science, with thousands of calculations cited annually. This ubiquity is rooted in the favorable accuracy vs cost balance of KSDFT. Nonetheless, the ambitions and expectations of researchers for use of KSDFT in predictive simulations of large, complicated molecular systems are confronted with an intrinsic computational cost-scaling challenge. Particularly evident in the context of first-principles molecular dynamics, the challenge is the high cost-scaling associated with the computation of the Kohn-Sham orbitals. Orbital-free DFT (OFDFT), as the name suggests, circumvents entirely the explicit use of those orbitals. Without them, the structural and algorithmic complexity of KSDFT simplifies dramatically and near-linear scaling with system size irrespective of system state is achievable. Thus, much larger system sizes and longer simulation time scales (compared to conventional KSDFT) become accessible; hence, new chemical phenomena and new materials can be explored. In this review, we introduce the historical contexts of OFDFT, its theoretical basis, and the challenge of realizing its promise via approximate kinetic energy density functionals (KEDFs). We review recent progress on that challenge for an array of KEDFs, such as one-point, two-point, and machine-learnt, as well as some less explored forms. We emphasize use of exact constraints and the inevitability of design choices. Then, we survey the associated numerical techniques and implemented algorithms specific to OFDFT. We conclude with an illustrative sample of applications to showcase the power of OFDFT in materials science, chemistry, and physics.

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.

Towards the description of charge transfer states in solubilised LHCII using subsystem DFT

Light harvesting complex II (LHCII) in plants and green algae have been shown to adapt their absorption properties, depending on the concentration of sunlight, switching between a light harvesting and a non-harvesting or quenched state. In a recent work, combining classical molecular dynamics (MD) simulations with quantum chemical calculations (Liguori et al. in Sci Rep 5:15661, 2015) on LHCII, it was shown that the Chl611-Chl612 cluster of the terminal emitter domain can play an important role in modifying the spectral properties of the complex. In that work the importance of charge transfer (CT) effects was highlighted, in re-shaping the absorption intensity of the chlorophyll dimer. Here in this work, we investigate the combined effect of the local excited (LE) and CT states in shaping the energy landscape of the chlorophyll dimer. Using subsystem Density Functional Theory over the classical [Formula: see text]s MD trajectory we look explicitly into the excitation energies of the LE and the CT states of the dimer and their corresponding couplings. Upon doing so, we observe a drop in the excitation energies of the CT states, accompanied by an increase in the couplings between the LE/LE and the LE/CT states facilitated by a shorter interchromophoric distance upon equilibration. Both these changes in conjunction, effectively produces a red-shift of the low-lying mixed exciton/CT states of the supramolecular chromophore pair.

Integrated Multiscale Multilevel Approach to Open Shell Molecular Systems

We present a novel multiscale approach to study the electronic structure of open shell molecular systems embedded in an external environment. The method is based on the coupling of multilevel Hartree-Fock (MLHF) and Density Functional Theory (MLDFT), suitably extended to the unrestricted formalism, to Molecular Mechanics (MM) force fields (FF). Within the ML region, the system is divided into active and inactive parts, thus describing the most relevant interactions (electrostatic, polarization, and Pauli repulsion) at the quantum level. The surrounding MM part, which is formulated in terms of nonpolarizable or polarizable FFs, permits a physically consistent treatment of long-range electrostatics and polarization effects. The approach is extended to the calculation of hyperfine coupling constants and applied to selected nitroxyl radicals in an aqueous solution.

Subsystem density-functional theory: A reliable tool for spin-density based properties

Subsystem density-functional theory compiles a set of features that allow for efficiently calculating properties of very large open-shell radical systems such as organic radical crystals, proteins, or deoxyribonucleic acid stacks. It is computationally less costly than correlated ab initio wave function approaches and can pragmatically avoid the overdelocalization problem of Kohn-Sham density-functional theory without employing hard constraints on the electron-density. Additionally, subsystem density-functional theory calculations commonly start from isolated fragment electron densities, pragmatically preserving a priori specified subsystem spin-patterns throughout the calculation. Methods based on subsystem density-functional theory have seen a rapid development over the past years and have become important tools for describing open-shell properties. In this Perspective, we address open questions and possible developments toward challenging future applications in connection with subsystem density-functional theory for spin-dependent properties.

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.

Multi-state formulation of the frozen-density embedding quasi-diabatization approach

We present a multi-state implementation of the recently developed frozen-density embedding diabatization (FDE-diab) methodology [D. G. Artiukhin and J. Neugebauer, J. Chem. Phys. 148, 214104 (2018)] in the Serenity program. The new framework extends the original approach such that any number of charge-localized quasi-diabatic states can be coupled, giving an access to calculations of ground and excited state spin-density distributions as well as to excitation energies. We show that it is possible to obtain results similar to those from correlated wave function approaches such as the complete active space self-consistent field method at much lower computational effort. Additionally, we present a series of approximate computational schemes, which further decrease the overall computational cost and systematically converge to the full FDE-diab solution. The proposed methodology enables computational studies on spin-density distributions and related properties for large molecular systems of biochemical interest.

The dynamics of light-induced interfacial charge transfer of different dyes in dye-sensitized solar cells studied by ab initio molecular dynamics

The charge-transport dynamics at the dye-TiO₂ interface plays a vital role for the resulting power conversion efficiency (PCE) of dye sensitized solar cells (DSSCs). In this work, we have investigated the charge-exchange dynamics for a series of organic dyes, of different complexity, and a small model of the semiconductor substrate TiO₂. The dyes studied involve L1, D35 and LEG4, all well-known organic dyes commonly used in DSSCs. The computational studies have been based on ab initio molecular dynamics (aiMD) simulations, from which structural snapshots have been collected. Estimates of the charge-transfer rate constants of the central exchange processes in the systems have been computed. All dyes show similar properties, and differences are mainly of quantitative character. The processes studied were the electron injection from the photoexcited dye, the hole transfer from TiO₂ to the dye and the recombination loss from TiO₂ to the dye. It is notable that the electronic coupling/transfer rates differ significantly between the snapshot configurations harvested from the aiMD simulations. The differences are significant and indicate that a single geometrically optimized conformation normally obtained from static quantum-chemistry calculations may provide arbitrary results. Both protonated and deprotonated dye systems were studied. The differences mainly appear in the rate constant of recombination loss between the protonated and the deprotonated dyes, where recombination losses take place at significantly higher rates. The inclusion of lithium ions close to the deprotonated dye carboxylate anchoring group mitigates recombination in a similar way as when protons are retained at the carboxylate group. This may give insight into the performance-enchancing effects of added salts of polarizing cations to the DSSC electrolyte. In addition, solvent effects can retard charge recombination by about two orders of magnitude, which demonstrates that the presence of a solvent will increase the lifetime of injected electrons and thus contribute to a higher PCE of DSSCs. It is also notable that no simple correlation can be identified between high/low transfer rate constants and specific structural arrangements in terms of atom-atom distances, angles or dihedral arrangements of dye sub-units.

A polarizable three-layer frozen density embedding/molecular mechanics approach

We present a novel multilayer polarizable embedding approach in which the system is divided into three portions, two of which are treated using density functional theory and their interaction is based on frozen density embedding (FDE) theory, and both also mutually interact with a polarizable classical layer described using an atomistic model based on fluctuating charges (FQ). The efficacy of the model is demonstrated by extending the formalism to linear response properties and applying it to the simulation of the excitation energies of organic molecules in aqueous solution, where the solute and the first solvation shell are treated using FDE, while the rest of the solvent is modeled using FQ charges.

Theoretical Assessment of Hinge-Type Models for Electron Donors in Reaction Centers of Photosystems I and II as well as of Purple Bacteria

Hinge-type molecular models for electron donors in reaction centers of Photosystems I and II and purple bacteria were investigated using a two-state computational approach based on frozen-density embedding (FDE). This methodology, dubbed FDE-diab, is known to avoid consequences of the self-interaction error as far as intermolecular phenomena are concerned, which allows a prediction of qualitatively correct spin densities for large biomolecular systems. The calculated spin density distributions are in a good agreement with available experimental results and demonstrated a very high sensitivity to changes in the relative orientation of cofactors and amino acid protonation states. This allows a validation of the previously proposed hinge-type models providing hints on possible protonation states of axial histidine molecules.

Multilevel Density Functional Theory

Following recent developments in multilevel embedding methods, we introduce a novel density matrix-based multilevel approach within the framework of density functional theory (DFT). In this multilevel DFT, the system is partitioned in an active and an inactive fragment, and all interactions are retained between the two parts. The decomposition of the total system is performed upon the density matrix. The orthogonality between the two parts is maintained by solving the Kohn–Sham equations in the MO basis for the active part only, while keeping the inactive density matrix frozen. This results in the reduction of computational cost. We outline the theory and implementation and discuss the differences and similarities with state-of-the-art DFT embedding methods. We present applications to aqueous solutions of methyloxirane and glycidol.

Energy-Based Molecular Orbital Localization in a Specific Spatial Region

We present a novel energy-based localization procedure able to localize molecular orbitals into predefined spatial regions. The method is defined in a multiscale framework based on the multilevel Hartree-Fock approach. In particular, the system is partitioned into active and inactive fragments. The localized molecular orbitals are obtained maximizing the repulsion between the two fragments. The method is applied to several cases including both conjugated and non-conjugated systems. Our multiscale approach is compared with reference values for both ground-state properties, such as dipole moments, and local excitation energies. The proposed approach is useful to extend the application range of high-level electron correlation methods. In fact, the reduced number of molecular orbitals can lead to a large reduction in the computational cost of correlated calculations.

Combining multilevel Hartree–Fock and multilevel coupled cluster approaches with molecular mechanics: a study of electronic excitations in solutions

We present the coupling of different quantum-embedding approaches with a third molecular-mechanics layer, which can be either polarizable or non-polarizable.

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.

Nonlocal Subsystem Density Functional Theory

By invoking a divide-and-conquer strategy, subsystem DFT dramatically reduces the computational cost of large-scale, ab initio electronic structure simulations of molecules and materials. The central ingredient setting subsystem DFT apart from Kohn-Sham DFT is the nonadditive kinetic energy functional (NAKE). Currently employed NAKEs are at most semilocal (i.e., they only depend on the electron density and its gradient), and as a result of this approximation, so far large-scale simulations only included systems composed of weakly interacting subsystems. In this work, we advance the state-of-the-art by introducing fully nonlocal NAKEs in subsystem DFT simulations for the first time. A benchmark analysis based on the S22-5 test set shows that nonlocal NAKEs considerably improve the computed interaction energies and electron densities compared to commonly employed GGA NAKEs, especially when increasing intersubsystem electron density overlap is considered. Most importantly, we resolve the long-standing problem of too attractive interaction energy curves typically resulting from the use of GGA NAKEs.

Detailed balance, internal consistency, and energy conservation in fragment orbital-based surface hopping

We have recently introduced an efficient semi-empirical non-adiabatic molecular dynamics method for the simulation of charge transfer/transport in molecules and molecular materials, denoted fragment orbital-based surface hopping (FOB-SH) [J. Spencer et al., J. Chem. Phys. 145, 064102 (2016)]. In this method, the charge carrier wavefunction is expanded in a set of charge localized, diabatic electronic states and propagated in the time-dependent potential due to classical nuclear motion. Here we derive and implement an exact expression for the non-adiabatic coupling vectors between the adiabatic electronic states in terms of nuclear gradients of the diabatic electronic states. With the non-adiabatic coupling vectors (NACVs) available, we investigate how different flavours of fewest switches surface hopping affect detailed balance, internal consistency, and total energy conservation for electron hole transfer in a molecular dimer with two electronic states. We find that FOB-SH satisfies detailed balance across a wide range of diabatic electronic coupling strengths provided that the velocities are adjusted along the direction of the NACV to satisfy total energy conservation upon a surface hop. This criterion produces the right fraction of energy-forbidden (frustrated) hops, which is essential for correct population of excited states, especially when diabatic couplings are on the order of the thermal energy or larger, as in organic semiconductors and DNA. Furthermore, we find that FOB-SH is internally consistent, that is, the electronic surface population matches the average quantum amplitudes, but only in the limit of small diabatic couplings. For large diabatic couplings, inconsistencies are observed as the decrease in excited state population due to frustrated hops is not matched by a corresponding decrease in quantum amplitudes. The derivation provided here for the NACV should be generally applicable to any electronic structure approach where the electronic Hamiltonian is constructed in a diabatic electronic state basis.

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.

Constrained subsystem density functional theory

Constrained Subsystem Density Fucntional Theory (CSDFT) allows to compute diabatic states for charge transfer reactions using the machinery of the constrained DFT method, and at the same time is able to embed such diabatic states in a molecular environment via a subsystem DFT scheme.

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.