Multistate Effects in Calculations of the Electronic Coupling Element for Electron Transfer Using the Generalized Mulliken−Hush Method

From Many-Body Ab Initio to Effective Excitonic Models: A Versatile Mapping Approach Including Environmental Embedding Effects

We present an original multistate projective diabatization scheme based on Green's function formalisms that allows the systematic mapping of many-body ab initio calculations onto effective excitonic models. This method inherits the ability of the Bethe-Salpeter equation to describe Frenkel molecular excitons and intermolecular charge-transfer states equally well, as well as the possibility for an effective description of environmental effects in a QM/MM framework. The latter is found to be a crucial element in order to obtain accurate model parameters for condensed phases and to ensure their transferability to excitonic models for extended systems. The method is presented through a series of examples illustrating its quality, robustness, and internal consistency.

Rational design of anthocyanidins-directed near-infrared two-photon fluorescent probes

Recently, two-photon fluorescent probes based on anthocyanidin molecules have attracted extensive attention due to their outstanding photophysical properties. However, there are only a few two-photon excited fluorescent probes that really meet the requirements of relatively long emission wavelengths (>600 nm), large two-photon absorption (TPA) cross-sections (300 GM), significant Stokes shift (>80 nm), and high fluorescence intensity. Herein, the photophysical properties of a series of anthocyanidins with the same substituents but different fluorophore skeletons are investigated in detail. Compared with b-series molecules, a-series molecules with a six-membered ring in the backbone have a slightly higher reorganization energy. This results in more energy loss upon light excitation, enabling the reaction products to detect NTR through a larger Stokes shift. More importantly, there is very little decrease in fluorescence intensity as the Stokes shift increases. These features are extremely valuable for high-resolution NTR detection. In light of this, novel 2a-n (n = 1-5) compounds are designed, which are accomplished by inhibiting the twisted intramolecular charge transfer (TICT) effect through alkyl cyclization, azetidine ring and extending π conjugation. Among them, 2a-3 gains a long emission spectrum (λem = 691.4 nm), noticeable TPA cross-section (957 GM), and large Stokes shift (110 nm), indicating that it serves as a promising candidate for two-photon fluorescent dyes. It is hoped that this work will offer some insightful theoretical direction for the development of novel high performance anthocyanin fluorescent materials.

The three kingdoms-Photoinduced electron transfer cascades controlled by electronic couplings

Excited states are the key species in photocatalysis, while the critical parameters that govern their applications are (i) excitation energy, (ii) accessibility, and (iii) lifetime. However, in molecular transition metal-based photosensitizers, there is a design tension between the creation of long-lived excited (triplet), e.g., metal-to-ligand charge transfer (3MLCT) states and the population of such states. Long-lived triplet states have low spin-orbit coupling (SOC) and hence their population is low. Thus, a long-lived triplet state can be populated but inefficiently. If the SOC is increased, the triplet state population efficiency is improved-coming at the cost of decreasing the lifetime. A promising strategy to isolate the triplet excited state away from the metal after intersystem crossing (ISC) involves the combination of transition metal complex and an organic donor/acceptor group. Here, we elucidate the excited state branching processes in a series of Ru(II)-terpyridyl push-pull triads by quantum chemical simulations. Scalar-relativistic time-dependent density theory simulations reveal that efficient ISC takes place along 1/3MLCT gateway states. Subsequently, competitive electron transfer (ET) pathways involving the organic chromophore, i.e., 10-methylphenothiazinyl and the terpyridyl ligands are available. The kinetics of the underlying ET processes were investigated within the semiclassical Marcus picture and along efficient internal reaction coordinates that connect the respective photoredox intermediates. The key parameter that governs the population transfer away from the metal toward the organic chromophore either by means of ligand-to-ligand (3LLCT; weakly coupled) or intra-ligand charge transfer (3ILCT; strongly coupled) states was determined to be the magnitude of the involved electronic coupling.

Perturbative Expansion of Nonorthogonal Product Approach for Charge Transfer States

Modeling of the excited states of multichromophoric systems is crucial for the understanding of photosynthesis functioning. The excitonic Hamiltonian method is widely used for such calculations. Excited states of the combined system are constructed from the wave functions of individual chromophores while their interactions are described by excitonic couplings. In the current study, we enhance a previously proposed nonorthogonal product approach to incorporate dynamic correlation effects accounted for by the multireference perturbation theory. We discuss the problems of constructing the excitonic Hamiltonian including charge transfer states for the molecular systems where the overlap contribution to the excitonic couplings is non-negligible. The benchmark calculations were performed for a model system. It was shown that the overlap component of the excitonic coupling is of great importance. The enhanced method provides an accurate description of the excited state energies and other properties.

A photoinduced mixed valence photoswitch

The ground state and photoinduced mixed valence states (GSMV and PIMV, respectively) of a dinuclear (Dp⁴⁺) ruthenium(II) complex bearing 2,2'-bipyridine ancillary ligands and a 2,2':4',4'':2'',2'''-quaterpyridine (Lp) bridging ligand were investigated using femtosecond and nanosecond transient absorption spectroscopy, electrochemistry and density functional theory. It was shown that the electronic coupling between the transiently light-generated Ru(II) and Ru(III) centers is H_DA ∼ 450 cm^-1 in the PIMV state, whereas the electrochemically generated GSMV state showed H_DA ∼ 0 cm^-1, despite virtually identical Ru-Ru distances. This stemmed from the changes in dihedral angles between the two bpy moieties of Lp, estimated at 30° and 4° for the GSMV and PIMV states, respectively, consistent with a through-bond rather than a through-space mechanism. Electronic coupling can be turned on by using visible light excitation, making Dp⁴⁺ a competitive candidate for photoswitching applications. A novel strategy to design photoinduced charge transfer molecular switches is proposed.

Charge-Separation and Charge-Recombination Rate Constants in a Donor-Acceptor Buckybowl-Based Supramolecular Complex: Multistate and Solvent Effects

The kinetics of the nonradiative photoinduced processes (charge-separation and charge-recombination) experimented in solution by a supramolecular complex formed by an electron-donating bowl-shaped truxene-tetrathiafulvalene (truxTTF) derivative and an electron-accepting fullerene fragment (hemifullerene, C30H12) has been theoretically investigated. The truxTTF·C30H12 heterodimer shows a complex decay mechanism after photoexcitation with the participation of several low-lying excited states of different nature (local and charge-transfer excitations) all close in energy. In this scenario, the absolute rate constants for all of the plausible charge-separation (CS) and charge-recombination (CR) channels have been successfully estimated using the Marcus-Levich-Jortner (MLJ) rate expression, electronic structure calculations, and a multistate diabatization method. The outcomes suggest that for a reasonable estimate of the CS and CR rate constants, it is necessary to include the following: (i) optimally tuned long-range (LC) corrected density functionals, to predict a correct energy ordering of the low-lying excited states; (ii) multistate effects, to account for the electronic couplings; and (iii) environmental solvent effects, to provide a proper stabilization of the charge-transfer excited states and accurate external reorganization energies. The predicted rate constants have been incorporated in a simple but insightful kinetic model that allows estimating global CS and CR rate constants in line with the most generalized three-state model used for the CS and CR processes. The values computed for the global CS and CR rates of the donor-acceptor truxTTF·C30H12 supramolecular complex are found to be in good agreement with the experimental values.

Extended Mulliken-Hush Method with Applications to the Theoretical Study of Electron Transfer

A novel adiabatic-to-diabatic (ATD) transformation strategy, namely, the extended Mulliken-Hush (XMH) method, is proposed to evaluate diabatic properties including electronic couplings, potential energy surfaces, and their crossings. The XMH method is developed by adopting our recently proposed ATD transformation formula of a general vectorial physical observable, in which a useful ATD transformation is further determined by using an auxiliary dipole between localized frontier orbitals as a simple approximation of the diabatic transition dipole. The XMH method is simple and practical that provides a flexible way to construct diabatic states. To some extent, it can be regarded as an extension of the generalized Mulliken-Hush (GMH) method since the latter takes a stronger approximation, in which the diabatic transition dipole is assumed to be vanishing. Test calculations on the HeH₂⁺ system show that the electronic couplings predicted by the XMH method are closer to the ones calculated by the valence bond block-diagonalization approach than the GMH ones since the XMH method takes into account both the magnitude and direction of the diabatic transition dipole, which is consistent with the properties of this molecule. In the study of electron transfer in the two kinds of donor-bridge-acceptor systems, the XMH method maintains the simplicity of the GMH method and gives reasonable results even when the latter fails, wherein the diabatic transition dipole is nearly perpendicular to the difference of the initial and final adiabatic dipoles. More importantly, the XMH method can be easily combined with high-level electronic structure methods, in which the properties of the ground and excited states may be more accurately calculated, and hence, one may expect that further development of the XMH method would result in a general computational model for studying electron transfer reactions.

Molecular Dynamics Simulations of Bimolecular Electron Transfer: the Distance-Dependent Electronic Coupling

Understanding the distance dependence of the parameters underpinning Marcus theory is imperative when interpreting the results of experiments on electron transfer (ET). Unfortunately, most of these parameters are difficult or impossible to access directly with experiments, necessitating the use of computer simulations to model them. In this work, we use molecular dynamics simulations in conjunction with constrained density functional theory calculations to study the distance dependence of the electronic coupling matrix element, |H_RP|, for bimolecular ET. Contrary to what is typically assumed for such intermolecular reactions, we find that the magnitude of |H_RP| does not decay exponentially with the center-of-mass separation of the reactants, r_COM. The addition of other simple measures of donor/acceptor (D/A) orientation did not improve the correlation of |H_RP| with r_COM. Using the minimum distance separation, r_min, of the reactants as the structural descriptor allowed the system to be partitioned into high-coupling/close-contact and low-coupling/non-contact regimes, but large fluctuations of |H_RP| were still found for the close-contact reactant pairs. Despite the persistent large fluctuations of |H_RP|, its mean value was found to decay piecewise exponentially with increasing r_min, which was attributed to significant changes in the average D/A pair structure. The results herein advise one to use caution when interpreting the experimental results derived from spherical reactant models of bimolecular ET.

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.

Structure Dependence of Kinetic and Thermodynamic Parameters in Singlet Fission Processes

Singlet fission-whereby one absorbed photon generates two coupled triplet excitons-is a key process for increasing the efficiency of optoelectronic devices by overcoming the Shockley-Queisser limit. A crucial parameter is the rate of dissociation of the coupled triplets, as this limits the number of free triplets subsequently available for harvesting and ultimately the overall efficiency of the device. Here we present an analysis of the thermodynamic and kinetic parameters for this process in parallel and herringbone dimers measured by electron paramagnetic resonance spectroscopy in coevaporated films of pentacene in p-terphenyl. The rate of dissociation is higher for parallel dimers than for their herringbone counterparts, as is the rate of recombination to the ground state. DFT calculations, which provide the magnitude of the electronic coupling as well as the distribution of molecular orbitals for each geometry, suggest that weaker triplet coupling in the parallel dimer is the driving force for faster dissociation. Conversely, localization of the molecular orbitals and a stronger triplet-triplet interaction result in slower dissociation and recombination. The identification and understanding of how the intermolecular geometry promotes efficient triplet dissociation provide the basis for control of triplet coupling and thereby the optimization of one important parameter of device performance.

The Influence of Single, Tandem, and Clustered DNA Damage on the Electronic Properties of the Double Helix: A Theoretical Study

Oxidatively generated damage to DNA frequently appears in the human genome as the effect of aerobic metabolism or as the result of exposure to exogenous oxidizing agents, such as ionization radiation. In this paper, the electronic properties of single, tandem, and clustered DNA damage in comparison with native ds-DNA are discussed as a comparative analysis for the first time. A single lesion—8-oxo-7,8-dihydro-2′-deoxyguanosine (G^oxo), a tandem lesion—(5′S) and (5′R) 5′,8-cyclo-2′-deoxyadenosine (cdA), and the presence of both of them in one helix turn as clustered DNA damage were chosen and taken into consideration. The lowest vertical and adiabatic potential (VIP ~ 5.9 and AIP ~ 5.5 eV, respectively) were found for G^oxo, independently of the discussed DNA lesion type and their distribution within the double helix. Moreover, the VIP and AIP were assigned for ds-trimers, ds- dimers and single base pairs isolated from parental ds-hexamers in their neutral and cationic forms. The above results were confirmed by the charge and spin density population, which revealed that G^oxo can be considered as a cation radical point of destination independently of the DNA damage type (single, tandem, or clustered). Additionally, the different influences of cdA on the charge transfer rate were found and discussed in the context of tandem and clustered lesions. Because oligonucleotide lesions are effectively produced as a result of ionization factors, the presented data in this article might be valuable in developing a new scheme of anticancer radiotherapy efficiency.

Uniform potential difference scheme to evaluate effective electronic couplings for superexchange electron transfer in donor-bridge-acceptor systems

This article proposes an ab initio quantum chemical method to evaluate the effective electronic coupling that determines the rate of superexchange electron transfer in donor-bridge-acceptor (D-B-A) systems. The method utilizes the fragment charge difference to define electronic diabatic states and to apply an electrostatic potential in a form of a uniform potential difference that mimics solvation effects on the relative energies of the electronic states. The two-state generalized Mulliken-Hush method is used to obtain the effective electronic coupling as the nondiagonal element of the effective Hamiltonian that is derived based on the Green's function approach and the quasi-degenerate perturbation theory. A theoretical basis is provided for the dependence of the calculated effective electronic coupling on the applied potential and for how to find the optimal potential to give the desired effective electronic coupling that coincides with the result of the minimum energy splitting method. The method is applied to typical D-B-A molecules and gives the effective electronic couplings in reasonable agreement with the experimental estimates.

Clustered DNA Damage: Electronic Properties and Their Influence on Charge Transfer. 7,8-Dihydro-8-Oxo-2'-Deoxyguaosine Versus 5',8-Cyclo-2'-Deoxyadenosines: A Theoretical Approach

Approximately 3 × 10¹⁷ DNA damage events take place per hour in the human body. Within clustered DNA lesions, they pose a serious problem for repair proteins, especially for iron–sulfur glycosylases (MutyH), which can recognize them by the electron-transfer process. It has been found that the presence of both 5′,8-cyclo-2′-deoxyadenosine (cdA) diastereomers in the ds-DNA structure, as part of a clustered lesion, can influence vertical radical cation distribution within the proximal part of the double helix, i.e., d[~oxoGcAoxoG~] (7,8-dihydro-8-oxo-2′-deoxyguaosine - ^oxodG). Here, the influence of cdA, “the simplest tandem lesion”, on the charge transfer through ds-DNA was taken into theoretical consideration at the M062x/6-31+G** level of theory in the aqueous phase. It was shown that the presence of (5′S)- or (5′R)-cdA leads to a slowdown in the hole transfer by one order of magnitude between the neighboring dG→^oxodG in comparison to “native” ds-DNA. Therefore, it can be concluded that such clustered lesions can lead to defective damage recognition with a subsequent slowing down of the DNA repair process, giving rise to an increase in mutations. As a result, the unrepaired, ^oxodG: dA base pair prior to genetic information replication can finally result in GC → TA or AT→CG transversion. This type of mutation is commonly observed in human cancer cells. Moreover, because local multiple damage sites (LMSD) are effectively produced as a result of ionization factors, the presented data in this article might be useful in developing a new scheme of radiotherapy treatment against the background of DNA repair efficiency.

Electronic couplings and rates of excited state charge transfer processes at poly(thiophene-co-quinoxaline)–PC71BM interfaces: two- versus multi-state treatments

Multi-state effects should be considered when calculating electronic couplings at local polymer–fullerene interfaces with the non-tuned and optimally tuned long-range corrected functionals.

Kinetics teach that electronic coupling lowers the free-energy change that accompanies electron transfer

Electron-transfer theories predict that an increase in the quantum-mechanical mixing (H_DA) of electron donor and acceptor wavefunctions at the instant of electron transfer drives equilibrium constants toward unity. Kinetic and equilibrium studies of four acceptor-bridge-donor (A-B-D) compounds reported herein provide experimental validation of this prediction. The compounds have two redox-active groups that differ only by the orientation of the aromatic bridge: a phenyl-thiophene bridge (p) that supports strong electronic coupling of H_DA > 1,000 cm^-1; and a xylyl-thiophene bridge (x) that prevents planarization and decreases H_DA < 100 cm^-1 without a significant change in distance. Pulsed-light excitation allowed kinetic determination of the equilibrium constant, K_eq In agreement with theory, K_eq(p) were closer to unity compared to K_eq(x). A van't Hoff analysis provided clear evidence of an adiabatic electron-transfer pathway for p-series and a nonadiabatic pathway for x-series. Collectively, the data show that the absolute magnitude of the thermodynamic driving force for electron transfers are decreased when adiabatic pathways are operative, a finding that should be taken into account in the design of hybrid materials for solar energy conversion.

Approximate DFT-based methods for generating diabatic states and calculating electronic couplings: models of two and more states

Four types of density functional theory (DFT)-based approaches are assessed in this work for the approximate construction of diabatic states and the evaluation of electronic couplings between these states. These approaches include the constrained DFT (CDFT) method, the constrained noninteracting electron (CNE) model to post-process Kohn-Sham operators, the approximate block-diagonalization (BD) of the Kohn-Sham operators, and the generalized Mulliken-Hush method. It is shown that the first three approaches provide a good description for long-distance intermolecular electron transfer (ET) reactions. On the other hand, inconsistent results were found when applying these approaches to intramolecular ET in strongly coupled, mixed-valence systems. Model analysis shows that this discrepancy is caused by the inappropriate use of the two-state model rather than the defects of the approaches themselves. The situation is much improved when more states are included in the model electronic Hamiltonian. The CNE and BD approaches can thus serve as efficient and robust alternatives for building ET models based on DFT calculations.

Investigation of NLO Properties of Fluorescent BORICO Dyes: a Comprehensive Experimental and Theoretical Approach

BORICO dyes with N, N-diethyl as a strong donor and BF₂ complexed iminocoumarin six member core as strong acceptor are investigated as an efficient non linear optical chromophores. Extended π-conjugation over iminocoumarin moiety is useful to make ICT character of BORICO dyes more significant and is established on the scale of Generalised Mulliken Hush analysis scale. Bond length alternation and bond order alternation values for three BORICO chromophores estimates the cyanine like framework for optimal non linear optical response. The frontier molecular orbital diagrams obtained from density functional theory calculations shows that there is charge transfer from donor to accepter as well as effective overlap between them making the basis for optimal NLO response of BORICO chromophores. The theoretical values of linear and non linear optical responses for three BORICO NLOphores obtained by using three different functionals B3LYP, CAMB3LYP and BHandHLYP with 6-311+g(d,p) basis set are quite consistent for the values of static dipole moment (μ), linear polarizability (α) and first hyperpolarizability (β). However in case of the γ values calculation, compare to the similar values obtained by CAMB3LYP and BHandHLYP functionals, B3LYP overestimates the same. The vibrational motions play decisive role in the overall non linear optical properties of BORICO chromophores.

Nonlinear Optical Properties of Pyrene Based Fluorescent Hemicurcuminoid and their BF2 Complexes -Spectroscopic and DFT Studies

The acetyl acetone and benzoyl acetone pyrene based difluoroboron scaffolds and their respective decomplexed congeners were investigated for their nonlinear optical properties in various microenvironment. The geometries of the chromophore were optimized and electronic excitation properties were estimated using time dependant density functional theory. The vertical excitation properties were found to be complimentary to the experimental results. Nolinear optical properties were evaluated using solvatochromism and Density Functional Theory. The experimental nonlinear optical properties well correlated using Mulliken Hush parameters and also compared with calculated values. CAM-B3LYP results were found to be superior over the B3LYP and BHandHLYP methods. The positive solvatochromism in emission were well explained by multilinear regression and bond length alteration.

Photoinduced Electron Transfer in Organic Solar Cells

Electron transfer (ET) is the key process in light-driven charge separation reactions in organic solar cells. The current review summarizes the progress in theoretical modelling of ET in these materials. First we give an account of ET, with a description originating from Marcus theory. We systematically go through all the relevant parameters and show how they depend on different material properties, and discuss the consequences such dependencies have for the performance of the devices. Finally, we present a set of visualization methods which have proven to be very useful in analyzing the elementary processes in absorption and charge separation events. Such visualization tools help us to understand the properties of the photochemical and photobiological systems in solar cells.

DFT-based Green's function pathways model for prediction of bridge-mediated electronic coupling

New LMO-GFM methodology enables intuitive understanding of electron tunneling in terms of through-bond and through-space interactions.

A multi-state fragment charge difference approach for diabatic states in electron transfer: extension and automation

The electron transfer (ET) rate prediction requires the electronic coupling values. The Generalized Mulliken-Hush (GMH) and Fragment Charge Difference (FCD) schemes have been useful approaches to calculate ET coupling from an excited state calculation. In their typical form, both methods use two eigenstates in forming the target charge-localized diabatic states. For problems involve three or four states, a direct generalization is possible, but it is necessary to pick and assign the locally excited or charge-transfer states involved. In this work, we generalize the 3-state scheme for a multi-state FCD without the need of manual pick or assignment for the states. In this scheme, the diabatic states are obtained separately in the charge-transfer or neutral excited subspaces, defined by their eigenvalues in the fragment charge-difference matrix. In each subspace, the Hamiltonians are diagonalized, and there exist off-diagonal Hamiltonian matrix elements between different subspaces, particularly the charge-transfer and neutral excited diabatic states. The ET coupling values are obtained as the corresponding off-diagonal Hamiltonian matrix elements. A similar multi-state GMH scheme can also be developed. We test the new multi-state schemes for the performance in systems that have been studied using more than two states with FCD or GMH. We found that the multi-state approach yields much better charge-localized states in these systems. We further test for the dependence on the number of state included in the calculation of ET couplings. The final coupling values are converged when the number of state included is increased. In one system where experimental value is available, the multi-state FCD coupling value agrees better with the previous experimental result. We found that the multi-state GMH and FCD are useful when the original two-state approach fails.

On the performance of the Kohn–Sham orbital approach in the calculation of electron transfer parameters. The three state model

The electron transfer parameters for the 3-state G⋯Ind system can be obtained using an efficient Kohn–Sham orbital based scheme.

Multistate treatments of the electronic coupling in donor-bridge-acceptor systems: insights and caveats from a simple model

We use a simple one-dimensional delta function electronic structure model (dfm) to investigate the results of a pair of multistate diabatization techniques (i.e., based on n states, with n ≥ 2) for linear DBA and DBBA (donor-bridge-acceptor) electron-transfer systems. In particular, we focus on the physical meaning of the couplings obtained from multistate methods and their relationship to two-state (n = 2) coupling elements. On the basis of the simple dfm approach, which allows exact as well as finite basis set treatment and has no many-electron effects, we conclude that for orthogonal diabatic states, it is difficult to assign clear physical significance to multistate matrix elements for coupling beyond nearest-neighbor contacts. The implications of these results for more complex multistate many-electron treatments are discussed. It is emphasized that physically meaningful coupling elements must involve states that are orthogonal, either explicitly or implicitly.

Electronic coupling for charge transfer in donor-bridge-acceptor systems. Performance of the two-state FCD model

Electronic coupling is a key parameter that determines the rate of electron transfer reactions and electrical conductivity of molecular wires. To examine the performance of a two-state approach based on the orthogonal transformation of adiabatic states to diabatic states, we compare the effective donor-acceptor coupling V(DA) computed with three different approaches in model donor-bridge-acceptor (D-B-A) systems. It is found that V(DA) derived with the two-state method accounts properly for both the direct and superexchange interactions. The approach becomes, however, less accurate with the increasing energy difference of the donor and acceptor states. We suggest a simple diagnostic to identify the situation when the estimated coupling might be inaccurate and consider how to improve the performance of the two-state scheme in such a case.

Electronic couplings in DNA pi-stacks: multistate effects

In this study, we employ a multistate generalized Mulliken-Hush approach for calculating electronic couplings V(da) for charge transfer (CT) in DNA pi-stacks consisting of three, four, and five base pairs. In these systems the guanine donor and acceptor sites are separated by several (AT) pairs. The Hartree-Fock calculations of the stacks are carried out with the standard 6-31G basis sets. All possible superexchange pathways are accounted for. We examine electronic couplings estimated using the two-state and multistate models. Although for some systems the two-state scheme provides reasonable estimates of V(da), in general this simple model fails to reproduce the electronic couplings calculated with the multistate approach. The two-state treatment of pi-stacks with a tunneling gap less than 0.3 eV, for instance, GAAG and GAAAG, may lead to invalid estimates of V(da). We consider the dependence of V(da) on the length and composition of the bridge. The calculations show that V(da) is less sensitive to the arrangement of nucleobases in the bridge, as can be predicted on the basis of electronic couplings between adjacent base pairs.