Estimation of electronic coupling in pi-stacked donor-bridge-acceptor systems: correction of the two-state model

Superior End-Group Stacking Promotes Simultaneous Multiple Charge-Transfer Mechanisms in Organic Solar Cells with an All-Fused-Ring Nonfullerene Acceptor

The all-fused-ring acceptor (AFRA) is a success for nonfullerene materials and has attracted considerable attention as its high optical and chemical stability expected to reduce energy loss, and power conversion efficiency (PCE) approaching 15% in constructed all-small-molecule organic solar cells (OSCs). Herein, the intrinsic role of the structure of AFRA F13 and the reason for its high PCE were revealed by comparison with those of typical fused acceptors IDT-IC and Y6. An increased degree of conjugation in F13 leads to broader and red-shifted absorption peaks, facilitating enhancement of the short-circuit current. Multiple charge-transfer mechanisms are mainly attributed to the higher Frenkel exciton (FE) state due to the multiple transition ways for acceptors in the C1-CN:F13 system. The increased number of atoms contributing to the charge-transfer (CT) state facilitated the existence of more superior stacking patterns with easy formation of CT and FE/CT states and a high charge separation rate. It was found using the AFRA is an effective strategy to enhance end-group stacking, enhancing the borrowing of oscillator strength to promote multiple CT mechanisms in the complexes, explaining the high performance of this OSC device. This work is promising to guide designing an efficient AFRA in the future.

A Theoretical Study on the Underlying Factors of the Difference in Performance of Organic Solar Cells Based on ITIC and Its Isomers

Recently, non-fullerene-based organic solar cells (OSCs) have made great breakthroughs, and small structural differences can have dramatic impacts on the power conversion efficiency (PCE). We take ITIC and its isomers as examples to study their effects on the performance of OSCs. ITIC and NFBDT only differed in the side chain position, and they were used as models with the same donor molecule, PBDB-T, to investigate the main reasons for the difference in their performance in terms of theoretical methods. In this work, a detailed comparative analysis of the electronic structure, absorption spectra, open circuit voltage and interfacial parameters of the ITIC and NFBDT systems was performed mainly by combining the density functional theory/time-dependent density functional theory and molecular dynamics simulations. The results showed that the lowest excited state of the ITIC molecule possessed a larger ∆q and more hybrid FE/CT states, and PBDB-T/ITIC had more charge separation paths as well as a larger kCS and smaller kCR. The reason for the performance difference between PBDB-T/ITIC and PBDB-T/NFBDT was elucidated, suggesting that ITIC is a superior acceptor based on a slight modulation of the side chain and providing a guiding direction for the design of superior-performing small molecule acceptor materials.

Modulating the organic photovoltaic properties of non-fullerene acceptors by molecular modification based on Y6: a theoretical study

Developing non-fullerene acceptors (NFAs) by modifying the backbone, side chains and end groups is the most important strategy to improve the power conversion efficiency of organic solar cells (OSCs). Among numerous developed NFAs, Y6 and its derivatives are famous NFAs in the OSC field due to their good performance. Herein, in order to understand the mechanism of tuning the photovoltaic performance by modifying the Y6's center backbone, π-spacer and side-chains, we selected the PM6:Y6 OSC as a reference and systematically studied PM6:AQx-2, PM6:Y6-T, PM6:Y6-2T, PM6:Y6-O, PM6:Y6-1O and PM6:Y6-2O OSC systems based on extensive quantum chemistry calculations. The results indicate that introducing quinoxaline to substitute thiadiazole in the backbone induces a blue-shift of absorption spectra, reduces the charge transfer (CT) distance (Δd) and average electrostatic potential (ESP), and increases the singlet-triplet energy gap (ΔEST), CT excitation energy and the number of CT states in low-lying excitations. Inserting thienyl and dithiophenyl as π spacers generates a red-shift of absorption spectra, enlarges Δd and average ESP, and reduces ΔEST and the number of CT states. Introducing furo[3,2-b]furan for substituting one thieno[3,2-b]thiophene unit in the Y6's backbone causes a red-shift of absorption spectra and increases ΔEST, Δd and average ESP as well as CT excitation energy. Introducing alkoxyl as a side chain results in a blue-shift of absorption spectra, and increases ΔEST, Δd, average ESP, CT excitation energy and the number of CT states. The rate constants calculated using Marcus theory suggest that all the molecular modifications of Y6 reduce the exciton dissociation and charge recombination rates at the heterojunction interface, while introducing furo[3,2-b]furan and alkoxyl enlarges CT rates.

The effects of electronic structures of two non-fullerene systems on their photovoltaic performances

In this contribution, the electronic structures of two polymer donors (PBDB-T and PBDB-T-SF) and two non-fullerene acceptors (ITIC and IT-4F) are researched by density functional theory and time-dependent density functional theory, respectively. The research purpose is to rationalize the relationship between observed experimental performances and structural properties and obtain the effects of structures on their photovoltaic performances. The investigated properties involve in the structure characteristics, absorption spectra, carrier mobilities, and exciton dissociation properties at interfaces to locate the essences of different power conversion efficiency between PBDB-T/ITIC and PBDB-T-SF/IT-4F. The results suggest that both PBDB-T/ITIC and PBDB-T-SF/IT-4F systems have stable structures and relatively high HOMO levels, which benefits to relatively large V_OC values. In addition, the larger PCE of PBDB-T-SF/IT-4F system originates from PBDB-T-SF's large hole transport properties and better exciton dissociation ability. Furthermore, the F and S incorporations enhance hole mobilities and exciton dissociation ability. Consequently, the theoretical results coincide well with the experimental ones.

Molecular tuning of non-fullerene electron acceptors in organic photovoltaics: a theoretical study

On the basis of the famous A–D–A-type non-fullerene acceptor IT-4F, this work investigates the effects of introducing methyl groups and substituting dicyano with O on optoelectronic properties and photovoltaic performances.

The Halogenation Effects of Electron Acceptor ITIC for Organic Photovoltaic Nano-Heterojunctions

Molecular engineering plays a critical role in the development of electron donor and acceptor materials for improving power conversion efficiency (PCE) of organic photovoltaics (OPVs). The halogenated acceptor materials in OPVs have shown high PCE. Here, to investigate the halogenation mechanism and the effects on OPV performances, based on the density functional theory calculations with the optimally tuned screened range-separated hybrid functional and the consideration of solid polarization effects, we addressed the halogenation effects of acceptor ITIC, which were modeled by bis-substituted ITIC with halogen and coded as IT-2X (X = F, Cl, Br), and PBDB-T:ITIC, PBDB-T:IT-2X (X = F, Cl, Br) complexes on their geometries, electronic structures, excitations, electrostatic potentials, and the rate constants of charge transfer, exciton dissociation (ED), and charge recombination processes at the heterojunction interface. The results indicated that halogenation of ITIC slightly affects molecular geometric structures, energy levels, optical absorption spectra, exciton binding energies, and excitation properties. However, the halogenation of ITIC significantly enlarges the electrostatic potential difference between the electron acceptor and donor PBDB-T with the order from fluorination and chlorination to bromination. The halogenation also increases the transferred charges of CT states for the complexes. Meanwhile, the halogenation effects on CT energies and electron process rates depend on different haloid elements. No matter which kinds of haloid elements were introduced in the halogenation of acceptors, the ED is always efficient in these OPV devices. This work provides an understanding of the halogenation mechanism, and is also conducive to the designing of novel materials with the aid of the halogenation strategy.

Assessing Effects of Different π bridges on Properties of Random Benzodithiophene-thienothiophene Donor and Non-fullerene Acceptor Based Active Layer

This report presents the effect of insertion of four different π bridges, furan, thienothiophene, thiophene, and thiazole, into a random benzodithiophene (BDT)-fluorinated-thienothiophene (TT-F) based donor. Starting from a structure of synthesized donor (D)-acceptor (A) random copolymer with 3:1 ratio, we have designed four D-π-A systems with four different π bridges. Structural, optoelectronic, and charge transport/transfer properties of these donors and donor/NDI (NDI = poly[N,N'-bis(2-hexyldecyl)-naphthalene-1,4,5,8-bis(dicarboximide)-2,6-diyl]-alt-5,5'-(2,2'-bithiophene)) blends are investigated using DFT and TD-DFT methodologies. Our results show that the thiazole based TzP1 oligomer has the deepest HOMO value resulting in the highest open circuit voltage among all systems. The maximum absorption wavelengths of π-linked systems are red-shifted compared to the parent molecule. Rates of charge transfer and charge recombination are the highest and smallest in case of the thiazole/NDI blend system. In addition, hole mobilities in thiophene, thienothiophene, and thiazole based systems are larger than in the parent system. The results indicate that the thiazole unit among the four π bridge units is the most suitable for active layer construction.

Tunneling matrix element and tunneling pathways of protein electron transfer calculated with a fragment molecular orbital method

Practical ways to calculate the tunneling matrix elements and analyze the tunneling pathways for protein electron-transfer (ET) reactions with a fragment molecular orbital (FMO) method are presented. The straightforward use of minimal basis sets only for the atoms involved in the covalent bond detachment in FMO can properly describe the ETs through the protein main-chains with the cost-effective two-body corrections (FMO2) without losing the quality of double-zeta basis sets. The current FMO codes have been interfaced with density functional theory, polarizable continuum model, and model core potentials, with which the FMO-based protein ET calculations can consider the effects of electron correlation, solvation, and transition-metal redox centers. The reasonable performance of the FMO-based ET calculations is demonstrated for three different sets of protein-ET model molecules: (1) hole transfer between two tryptophans covalently bridged by a polyalanine linker in the ideal α-helix and β-strand conformations, (2) ET between two plastoquinones covalently bridged by a polyalanine linker in the ideal α-helix and β-strand conformations, and (3) hole transfer between ruthenium (Ru) and copper (Cu) complexes covalently bridged by a stretch of a polyglycine linker as a model for Ru-modified derivatives of azurin.

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.

Photoinduced charge transfer rate of Cy3/C60 blend material

The rates of charge separation and charge recombination of the cyanine dye/C₆₀ heterojunction solar cell in an external electric field were provided using the Marcus and Marcus-Levich-Jortner formalisms. The vibrational mode as another influencing factor was also introduced into the rate expression for the planar heterojunction solar cell. Detailed theoretical analysis of the excited-state of the Cy3/C₆₀ blend was achieved using density functional theory and time-dependent density functional theory. The Gibbs free energy was regulated by an external electric field, while the reorganisation energy presented the opposite conclusion. Frequency analysis was utilised to demonstrate the energy stability of the obtained structures. The rate calculated using the Marcus formalism was greater than that obtained by the Marcus-Levich-Jortner formalism. Consideration of the calculated rates in all vibration modes and at different external electric field strengths indicated that vibrational mode and external electric field played important roles in determining the rates of charge separate and charge recombination, which could provide a more accurate theoretical rate for organic photovoltaic devices.

Modification of NFA-Conjugated Bridges with Symmetric Structures for High-Efficiency Non-Fullerene PSCs

As electron acceptors, non-fullerene molecules can overcome the shortcomings of fullerenes and their derivatives (such as high cost, poor co-solubility, and weak light absorption). The photoelectric properties of two potential non-fullerene polymer solar cells (PSCs) PBDB-T:IF-TN (PB:IF) and PBDB-T:IDT-TN (PB:IDT) are studied by density functional theory (DFT) and time-dependent DFT (TD-DFT). Based on the optimized structure of the ground state, the effects of the electron donor (D) and electron acceptor (A) (D/A) interfaces PBDB-T/IF-TN (PB/IF) and PBDB-T/IDT-TN (PB/IDT) are studied by a quantum-chemical method (QM) and Marcus theory. Firstly, for two non-fullerene acceptors (NFAs) IF-TN and IDT-TN, the NFA IDT-TN has better optical absorption ability and better electron transport ability than IF-TN. Secondly, for the D/A interfaces PB/IF and PB/IDT, they both have high optical absorption and electron transfer abilities, and PB/IDT has better optical absorption and lower exciton binding energy. Finally, some important parameters (open-circuit voltage, voltage loss, fill factor, and power conversion efficiency) are calculated and simulated by establishing the theoretical model. From the above analysis, the results show that the non-fullerene PSC PB:IDT has better photoelectric characteristics than PB:IF.

Photoinduced electron transfer in conjugated molecules containing different donor in an external electric field

Photoinduced electron transfers in different conjugated donor polymer containing with PC₆₀BM as the acceptor (X3:PC₆₀BM and X6:PC₆₀BM) were theoretically investigated. A detailed theoretical analysis of Marcus electron transfer rate in an external electric field was also performed. By discussing the effect of electric field on ΔG, λ, V_DA and charge separation rate, we have come to some conclusions. The results reveal that the free energy change and electronic coupling matrix are significantly influence by the vector properties of external electric field. Thus, the electron transfer process appears critically influenced by external electric fields. Our results will provide meaningful information on enhancing the photoelectric conversion efficiencies of BHJ solar cells.

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.

Photoinduced charge transfer by one and two-photon absorptions: physical mechanisms and applications

We review photoinduced charge transfer in organic solar cells without and with an external electric field and then we introduce the visualization methods of the transition density, charge difference density and transition density matrix for the analysis of the photoinduced charge transfer in a neutral system and a charged system excited by one-photon and two-photon absorption. This review will not only promote a deeper understanding of the available theories and methods of PICT, but also lead to further developments in this field.

Theoretical Study on the Photoelectric Properties of a Class of Copolymers Based on Benzodithiophene for Solar Cells

The structural, electronic, and optical properties of PBDTTBT are comprehensively studied by density functional theory to rationalize the experimentally observed properties. Periodic boundary conditions method is employed to simulate the polymer block and calculate effective charge mass from the band structure calculation to describe charge transport properties. Moreover, both time-dependent density functional theory and a set of multidimensional visualization techniques are used to characterize the exciton dissociation ability in the PCBM: PBDTTBT interface. These theoretical methods and calculation techniques not only promote deep understanding of the connection between chemical structures and optical and electronic properties of the donor-acceptor system but also can be used to rationally design a novel donor-acceptor system. Based on the same calculated methods as PBDTTBT, four copolymers PBDTTTP, PBDTTTO, PBDTTTPD, and PBDTTFPD are designed to study their potentials as donors in polymer BHJ. The results indicate that PBDTTBT’s well conjugation benefits its good stability, and its wide and strong absorption spectra in the range of visible light, appropriate FMO levels, well charge transport, and favorable exciton dissociation lead to its photovoltaic performance. Furthermore, through comparing the four designed polymers with PBDTTBT, we conclude that the four designed polymers have stronger exciton dissociation ability and larger open-circuit voltage and external quantum efficiencies. Consequently, the four designed copolymers are promising candidates for polymer BHJ solar cells.

Non-Fullerene Acceptor-Based Solar Cells: From Structural Design to Interface Charge Separation and Charge Transport

The development of non-fullerene small molecule as electron acceptors is critical for overcoming the shortcomings of fullerene and its derivatives (such as limited absorption of light, poor morphological stability and high cost). We investigated the electronic and optical properties of the two selected promising non-fullerene acceptors (NFAs), IDIC and IDTBR, and five conjugated donor polymers using quantum-chemical method (QM). Based on the optimized structures of the studied NFAs and the polymers, the ten donor/acceptor (D/A) interfaces were constructed and investigated using QM and Marcus semi-classical model. Firstly, for the two NFAs, IDTBR displays better electron transport capability, better optical absorption ability, and much greater electron mobility than IDIC. Secondly, the configurations of D/A yield the more bathochromic-shifted and broader sunlight absorption spectra than the single moiety. Surprisingly, although IDTBR has better optical properties than IDIC, the IDIC-based interfaces possess better electron injection abilities, optical absorption properties, smaller exciton binding energies and more effective electronic separation than the IDTBR-based interfaces. Finally, all the polymer/IDIC interfaces exhibit large charge separation rate (KCS) (up to 1012⁻1014 s-1) and low charge recombination rate (KCR) (<10⁶ s-1), which are more likely to result in high power conversion efficiencies (PCEs). From above analysis, it was found that the polymer/IDIC interfaces should display better performance in the utility of bulk-heterojunction solar cells (BHJ OSC) than polymer/IDTBR interfaces.

A theoretical exploration of the effect of fluorine and cyano substitutions in diketopyrrolopyrrole-based polymer donor for organic solar cells

A series of polymer donor materials 1-5 based on diketopyrrolopyrrole and thiophene unit which have been widely used in organic solar cells (OSCs) were investigated based on quantum chemical calculations. The effect of fluorine and cyano substitutions in polymer donor materials was focused on. Based on the investigation on electronic structures and optical properties of the reported molecules 1 and 2 and the analysis on some parameters relevant to charge dissociation ability at donor/acceptor interface constituted by 1 and 2 with PC₆₁BM such as intermolecular charge transfer and recombination, driving force and Coulombic bound energy, we explained why fluorine substitution can improve OPV efficiency through strengthening eletron-withdrawing ability from a theoretical perspective. Then we designed cyano-substituted polymers 3-5 with the aim of obtaining better photovoltaic donor materials. The results reveal that our attempt to design donor materials which can balance large open-circuit voltage (V_oc) and high short-circuit current (J_sc) in OSCs has worked out. It is worth noting that the substitutions of fluorine and cyano groups synergistically reduce energy gap and HOMO energy level of polymers 3 and 4. Moreover, 3/PC₆₁BM and 4/PC₆₁BM heterojunctions show over 10⁷ and 10⁴ times higher than 1/PC₆₁BM on the ratios of intermolecular charge transfer and recombination rates (k_inter-CT/k_inter-CR). Thus, our work here may provide an efficient strategy to design promising donor materials in OPVs and we hope it could be useful in the future experimental synthesis.

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.

A comparative study of a fluorene-based non-fullerene electron acceptor and PC61BM in an organic solar cell at a quantum chemical level

A comparison of the OSC efficiency between PC61BM and non-fullerene small molecule acceptor-FENIDT was studied by microscopic analysis based on the DFT/TDDFT calculation.

Photoinduced Charge Transport in a BHJ Solar Cell Controlled by an External Electric Field

This study investigated theoretical photoinduced charge transport in a bulk heterojunction (BHJ) solar cell controlled by an external electric field. Our method for visualizing charge difference density identified the excited state properties of photoinduced charge transfer, and the charge transfer excited states were distinguished from local excited states during electronic transitions. Furthermore, the calculated rates for the charge transfer revealed that the charge transfer was strongly influenced by the external electric field. The external electric field accelerated the rate of charge transfer by up to one order when charge recombination was significantly restrained. Our research demonstrated that photoinduced charge transport controlled by an external electric field in a BHJ solar cell is efficient, and the exciton dissociation is not the limiting factor in organic solar cells.Our research should aid in the rational design of a novel conjugated system of organic solar cells.

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.

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.

Guanine holes are prominent targets for mutation in cancer and inherited disease

Single base substitutions constitute the most frequent type of human gene mutation and are a leading cause of cancer and inherited disease. These alterations occur non-randomly in DNA, being strongly influenced by the local nucleotide sequence context. However, the molecular mechanisms underlying such sequence context-dependent mutagenesis are not fully understood. Using bioinformatics, computational and molecular modeling analyses, we have determined the frequencies of mutation at G • C bp in the context of all 64 5'-NGNN-3' motifs that contain the mutation at the second position. Twenty-four datasets were employed, comprising >530,000 somatic single base substitutions from 21 cancer genomes, >77,000 germline single-base substitutions causing or associated with human inherited disease and 16.7 million benign germline single-nucleotide variants. In several cancer types, the number of mutated motifs correlated both with the free energies of base stacking and the energies required for abstracting an electron from the target guanines (ionization potentials). Similar correlations were also evident for the pathological missense and nonsense germline mutations, but only when the target guanines were located on the non-transcribed DNA strand. Likewise, pathogenic splicing mutations predominantly affected positions in which a purine was located on the non-transcribed DNA strand. Novel candidate driver mutations and tissue-specific mutational patterns were also identified in the cancer datasets. We conclude that electron transfer reactions within the DNA molecule contribute to sequence context-dependent mutagenesis, involving both somatic driver and passenger mutations in cancer, as well as germline alterations causing or associated with inherited disease.

Combined density functional theory and Landauer approach for hole transfer in DNA along classical molecular dynamics trajectories

We investigate in detail the charge transport characteristics of DNA wires with various sequences and lengths in the presence of solvent. Our approach combines large-scale quantum/classical molecular dynamics (MD) simulations with transport calculations based on Landauer theory. The quantum mechanical transmission function of the wire is calculated along MD trajectories and thus encodes the influence of dynamical disorder arising from the environment (water, backbone, counterions) and from the internal base dynamics. We show that the correlated fluctuations of the base pair dynamics are crucial in determining the transport properties of the wire and that the effect of fluctuations can be quite different for sequences with low and high static disorders (differences in base ionization potentials). As a result, in structures with high static disorder as is the case of the studied Dickerson dodecamer, the weight of high-transmissive structures increases due to dynamical fluctuations and so does the calculated average transmission. Our analysis further supports the basic intuition of charge-transfer active conformations as proposed by Barton et al. [J. Am. Chem. Soc. 126, 11471 (2004)]. However, not DNA conformations with good stacking contacts leading to large interbase hopping values are necessarily the most important, but rather those where the average fluctuation of ionization potentials along the base stack is small. The reason behind this is that the ensemble of conformations leads to average electronic couplings, which are large enough for sufficient transmission. On the other hand, the alignment of onsite energies is the critical parameter which gates the charge transport.

Charge transport through biomolecular wires in a solvent: bridging molecular dynamics and model Hamiltonian approaches

We present a hybrid method based on a combination of classical molecular dynamics simulations, quantum-chemical calculations, and a model Hamiltonian approach to describe charge transport through biomolecular wires with variable lengths in presence of a solvent. The core of our approach consists in a mapping of the biomolecular electronic structure, as obtained from density-functional based tight-binding calculations of molecular structures along molecular dynamics trajectories, onto a low-dimensional model Hamiltonian including the coupling to a dissipative bosonic environment. The latter encodes fluctuation effects arising from the solvent and from the molecular conformational dynamics. We apply this approach to the case of pG-pC and pA-pT DNA oligomers as paradigmatic cases and show that the DNA conformational fluctuations are essential in determining and supporting charge transport.

CASSCF/CAS-PT2 study of hole transfer in stacked DNA nucleobases

CASSCF and CAS-PT2 calculations are performed for the ground and excited states of radical cations consisting of two and three nucleobases. The generalized Mulliken-Hush approach is employed for estimating electronic couplings for hole transfer in the pi-stacks. We compare the CASSCF results with data obtained within Koopmans' approximation. The calculations show that an excess charge in the ground and excited states in the systems is quite localized on a single base both at the CASSCF level and in Koopmans' picture. However, the CASSCF calculations point to a larger degree of localization and, in line with this, smaller transition dipole moments. The agreement between the CAS-PT2 corrected energy gaps and the values estimated with Koopmans' theorem is better, with the CAS-PT2 calculations giving somewhat smaller gaps. Overall, both factors result in smaller CASSCF/CAS-PT2 couplings, which are reduced by up to 40% of the couplings calculated using Koopmans' approximation. The tabulated data can be used as benchmark values for the electronic couplings of stacked nucleobases. For the base trimers, comparison of the results obtained within two- and three-state models show that the multistate treatment should be applied to derive reliable estimates. Finally, the superexchange approach to estimate the donor acceptor electronic coupling in the stacks GAG and GTG is considered.

Electronic Coupling Mediated by Stacked [Thymine-Hg-Thymine] Base Pairs

Very recently it has been shown that stable metal-mediated base pairs [Thymine-Hg-Thymine] can form in DNA. To estimate the effect of such pairs on the efficiency of charge transfer through DNA, we carry out quantum mechanical calculations of double-stranded pi-stacks GXG, GXXG, and GXXXG, where X = [Thymine-Hg-Thymine] and stacks GT(n)G of canonical base pairs. The charge-transfer efficiency in short duplexes GXG and GTG is found to be similar. However, the donor-acceptor coupling in GXXG and GXXXG is stronger by a factor of 2.5-3.0 than that in GT(n)G (n = 2 and 3), respectively. It is shown that the valence orbitals of Hg atoms do not essentially participate in mediating the electronic coupling for hole transfer; however, they may play an important role in excess electron transfer.