Antioxidative Triplet Excitation Energy Transfer in Bacterial Reaction Center Using a Screened Range Separated Hybrid Functional

Excess energy absorbed by photosystems (PSs) can result in photoinduced oxidative damage. Transfer of such energy within the core pigments of the reaction center in the form of triplet excitation is important in regulating and preserving the functionality of PSs. In the bacterial reaction center (BRC), the special pair (P) is understood to act as the electron donor in a photoinduced charge transfer process, triggering the charge separation process through the photoactive branch A pigments that experience a higher polarizing environment. At this work, triplet excitation energy transfer (TEET) in BRC is studied using a computational perspective to gain insights into the roles of the dielectric environment and interpigment orientations. We find in agreement with experimental observations that TEET proceeds through branch B. The TEET process toward branch B pigment is found to be significantly faster than the hypothetical process proceeding through branch A pigments with ps and ms time scales, respectively. Our calculations find that conformational differences play a major role in this branch asymmetry in TEET, where the dielectric environment asymmetry plays only a secondary role in directing the TEET to proceed through branch B. We also address TEET processes asserting the role of carotenoid as the final triplet energy acceptor and in a mutant form, where the branch pigments adjacent to P are replaced by bacteriopheophytins. The necessary electronic excitation energies and electronic state couplings are calculated by the recently developed polarization-consistent framework combining a screened range-separated hybrid functional and a polarizable continuum mode. The polarization-consistent potential energy surfaces are used to parametrize the quantum mechanical approach, implementing Fermi's golden rule expression of the TEET rate calculations.

A Computational Study of the Electronic Energy and Charge Transfer Rates and Pathways in the Tetraphenyldibenzoperiflanthene/Fullerene Interfacial Dyad

The electronic transition rates and pathways underlying interfacial charge separation in tetraphenyldibenzoperiflanthene:fullerene (DBP:C70) blends are investigated computationally. The analysis is based on a polarization-consistent framework employing screened range-separated hybrid functional in a polarizable continuum model to parametrize Fermi's golden rule rate theory. The model considers the possible transitions within the 25 lowest excited states of a DBP:C70 dyad that are accessible by photoexcitation. The different identified pathways contributing to charge carrier generation include electron and hole transfer and backtransfer, exciton transfer, and internal relaxation steps. The larger density of states of C70 appears to explain the previously observed larger efficiency for charge separation through hole transfer mechanism. We also analyze the validity of the high-temperature and short-time semiclassical approximations of the FGR theory, where both overestimated and underestimated Marcus theory based constants can be affected.

Polarization consistent dielectric screening in polarizable continuum model calculations of solvation energies

A polarization consistent framework, where dielectric screening is affected consistently in polarizable continuum model (PCM) calculations, is employed for the study of solvation energies. The computational framework combines a screened range-separated-hybrid functional (SRSH) with PCM calculations, SRSH-PCM, where dielectric screening is imposed in both PCM self-consistent reaction field (SCRF) iterations and the electronic structure Hamiltonian. We begin by demonstrating the impact of modifying the Hamiltonian to include such dielectric screening in SCRF iterations by considering the solutions of electrostatically embedded Hartree-Fock (HF) exact exchange equations. Long-range screened HF-PCM calculations are shown to capture properly the linear dependence of gap energy of frontier orbitals on the inverse of the dielectric constant, whereas unscreened HF-PCM orbital energies are fallaciously semi-constant with respect to the dielectric constant and, therefore, inconsistent with the ionization energy gaps. Similar trends affect density functional theory (DFT) calculations that aim to achieve predictive quality. Importantly, the dielectric screened calculations are shown to significantly affect DFT- and HF PCM-based solvation energies, where screened solvation energies are smaller compared to the unscreened values. Importantly, SRSH-PCM, therefore, appears to reduce the tendency of DFT-PCM to overestimate solvation energies, where we find the effect to increase with the dielectric constant and the polarity of the molecular solute, trends that enhance the quality of DFT-PCM calculations of solvation energy. Understanding the relationship of dielectric screening in the Hamiltonian and DFT-PCM calculations can ultimately benefit on-going efforts for the design of predictive and parameter free descriptions of solvation energies.

Visible light-activatable platinum(IV) prodrugs harnessing CD36 for ovarian cancer therapy

We hereby engineered photoactivatable Pt(IV) metallodrugs that harness CD36 to target ovarian cancer cells. Pt(IV) compounds mimic the structure of fatty acids and take advantage of CD36 as a "Trojan horse" to gain entry into the cells. We confirmed that CD36-dependent entry occurs using graphite furnace atomic absorption spectroscopy with ovarian cancer cells expressing different levels of CD36 and a CD36 inhibitor, SSO. Once the Pt(IV) metallodrugs enter the cancer cells, they can be activated to form Pt(II) with characteristics of cisplatin under visible light (490 nm) irradiation, promoting photoinduced electron transfer from the attached fluorophore to the metal center. This light-induced activation can increase the cytotoxicity of the Pt(IV) metallodrugs by up to 20 times toward ovarian cancer cells, inducing DNA damage and enabling efficient elimination of drug-resistant cancer cells.

Engineering giant excitonic coupling in bioinspired, covalently bridged BODIPY dyads

Strong excitonic coupling in photosynthetic systems is believed to enable efficient light absorption and quantitative charge separation, motivating the development of artificial multi-chromophore arrays with equally strong or even stronger excitonic coupling. However, large excitonic coupling strengths have typically been accompanied by fast non-radiative recombination, limiting the potential of the arrays for solar energy conversion as well as other applications such as fluorescent labeling. Here, we report giant excitonic coupling leading to broad optical absorption in bioinspired BODIPY dyads that have high photostability, excited-state lifetimes at the nanosecond scale, and fluorescence quantum yields of nearly 50%. Through the synthesis, spectroscopic characterization, and computational modeling of a series of dyads with different linking moieties, we show that the strongest coupling is obtained with diethynylmaleimide linkers, for which the coupling occurs through space between BODIPY units with small separations and slipped co-facial orientations. Other linkers allow for broad tuning of both the relative through-bond and through-space coupling contributions and the overall strength of interpigment coupling, with a tradeoff observed in general between the strength of the two coupling mechanisms. These findings open the door to the synthesis of molecular systems that function effectively as light-harvesting antennas and as electron donors or acceptors for solar energy conversion.

Effects of Solvent Dielectric on Thermally Activated Delayed Fluorescence: A Predictive Computational Polarization Consistent Approach

We study computationally thermally activated delayed fluorescence (TADF) in donor-acceptor compounds. The relevant electronic excited states that are strongly affected by the dielectric environment are treated by a polarization consistent framework. The high fidelity potential energy surfaces are used following a quantum-mechanical Fermi's golden rule (FGR) picture to calculate rates of intersystem crossing (ISC) and reverse intersystem crossing (RISC). To demonstrate the potency of the approach, we consider isomers of benzonitrile functionalized tert-butyl-substituted dimethylacridine (DMAC-BN), which were recently found to perform well as TADF emitters. The calculated excited state energies that appear to reproduce well measured spectral trends with respect to the dielectric constant are used to parametrize ISC/RISC FGR rates. The calculated rates reproduce well measured rates, whereas semiclassical based rates are grossly underestimated. In particular, we find in agreement with the recent experimental study [Phys. Rev. Appl.2019, 12, 044021] that the ortho and meta isomers are significantly more effective as TADF emitters. The computational framework provides valuable insight at the molecular level into RISC rates and therefore can contribute to the design of materials of increased TADF efficiency.

Luminescent 1H-1,3-benzazaphospholes

2-R-1H-1,3-Benzazaphospholes (R-BAPs) are an interesting class of σ²P heterocycles containing P[double bond, length as m-dash]C bonds. While closely related 2-R-1,3-benzoxaphospholes (R-BOPs) have been shown to be highly photoluminescent materials depending on specific R substituents, photoluminescence of R-BAPs has been previously limited to an example having a fused carbazole ring system. Here we detail the synthesis and structural characterization of a new R-BAP (3c, R = 2,2'-dithiophene), and compare its photoluminescence against two previously reported R-BAPs (3a, R, R' = Me and 3b, R = 2-thiophene). The significant fluorescence displayed by the thiophene derivatives 3b (φ = 0.53) and 3c (φ = 0.12) stands in contrast to the weakly emissive methyl substituted analogue 3a (φ = 0.08). Comparative computational investigations of 3a-c offer insights into the interplay between structure-function relationships affecting excited state relaxation processes.

Solvent Dependent Nuclear Magnetic Resonance Molecular Parameters Based on a Polarization Consistent Screened Range Separated Hybrid Density Functional Theory Framework

Nuclear magnetic resonance (NMR) properties of solvated molecules are significantly affected by the solvent. We, therefore, employ a polarization consistent framework that efficiently addresses the solvent polarizing environment effects. Toward this goal a dielectric screened range separated hybrid (SRSH) functional is invoked with a polarizable continuum model (PCM) to properly represent the orbital gap in the condensed phase. We build on the success of range separated hybrid (RSH) functionals to address the erroneous tendency of traditional density functional theory (DFT) to collapse the orbital gap. Recently, the impact of RSH that properly opens up the orbital gap in gas-phase calculations on NMR properties has been assessed. Here, we report the use of SRSH-PCM that produces properly solute orbital gaps in calculating isotropic nuclear magnetic shielding and chemical shift parameters of molecular systems in the condensed phase. We show that in contrast to simpler DFT-PCM approaches, SRSH-PCM successfully follows expected dielectric constant trends.

Role of Dielectric Screening in Calculating Excited States of Solvated Azobenzene: A Benchmark Study Comparing Quantum Embedding and Polarizable Continuum Model for Representing the Solvent

The low energy excited states of the conformational isomers of solvated azobenzene are calculated with several DFT methods accounting for the solute-solvent interaction implicitly with the polarizable continuum model or explicitly with subsystem DFT. For the latter, embedding potentials are calculated for 21 sampled snapshots of the solvent molecules. First, we find that accounting for the solvent implicitly or explicitly has little effect on the predicted cis-trans S₁ excitation energy gap. Second, we find that azobenzene's S₁ cis and trans energies are accurate as long as a screened range-separated hybrid exchange-correlation functional is employed. Finally, we also tested a simplified workflow whereby a single, averaged, embedding potential is used. Unfortunately, we find larger deviations against the experiment for the simplified workflow. This highlights a basic flaw in the approach, where the time scale of solvent averaging is much longer than that of the solute's electronic polarization.

Correlating Interfacial Charge Transfer Rates with Interfacial Molecular Structure in the Tetraphenyldibenzoperiflanthene/C70 Organic Photovoltaic System

Organic photovoltaics (OPV) is an emerging solar cell technology that offers vast advantages such as low-cost manufacturing, transparency, and solution processability. However, because the performance of OPV devices is still disappointing compared to their inorganic counterparts, better understanding of how controlling the molecular-level morphology can impact performance is needed. To this end, one has to overcome significant challenges that stem from the complexity and heterogeneity of the underlying electronic structure and molecular morphology. In this Letter, we address this challenge in the context of the DBP/C₇₀ OPV system by employing a modular workflow that combines recent advances in electronic structure, molecular dynamics, and rate theory. We show how the wide range of interfacial pairs can be classified into four types of interfacial donor-acceptor geometries and find that the least populated interfacial geometry gives rise to the fastest charge transfer (CT) rates.

Electronic Spectra of C60 Films Using Screened Range Separated Hybrid Functionals

We study computationally the electronic spectra of C₆₀ thin films using the recently developed density functional theory (DFT) framework combining a screened range separated hybrid (SRSH) functional with a polarizable continuum model (PCM). The SRSH-PCM approach achieves excellent correspondence between the frontier orbital's energy levels and the ionization potential and electron affinity of the molecular system at the condensed phase and consequently leads to high quality electronic excitation energies when used in time-dependent DFT calculations. Our calculated excited states reproduce the experimentally main reported spectral peaks at the 3.6-4.6 eV energy range and when addressing excitonic effects also reproduce the red-shifted spectral feature. Notably, we analyze the low-lying peak at 2.7 eV and associate it to an excitonic state.

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

CTRAMER: An open-source software package for correlating interfacial charge transfer rate constants with donor/acceptor geometries in organic photovoltaic materials

In this paper, we present CTRAMER (Charge-Transfer RAtes from Molecular dynamics, Electronic structure, and Rate theory)-an open-source software package for calculating interfacial charge-transfer (CT) rate constants in organic photovoltaic (OPV) materials based on ab initio calculations and molecular dynamics simulations. The software is based on identifying representative donor/acceptor geometries within interfacial structures obtained from molecular dynamics simulation of donor/acceptor blends and calculating the corresponding Fermi's golden rule CT rate constants within the framework of the linearized-semiclassical approximation. While the methods used are well established, the integration of these state-of-the-art tools originating from different disciplines to study photoinduced CT processes with explicit treatment of the environment, in our opinion, makes this package unique and innovative. The software also provides tools for investigating other observables of interest. After outlining the features and implementation details, the usage and performance of the software are demonstrated with results from an example OPV system.

Simulating energy transfer dynamics in the Fenna-Matthews-Olson complex via the modified generalized quantum master equation

The generalized quantum master equation (GQME) provides a general and formally exact framework for simulating the reduced dynamics of open quantum systems. The recently introduced modified approach to the GQME (M-GQME) corresponds to a specific implementation of the GQME that is geared toward simulating the dynamics of the electronic reduced density matrix in systems governed by an excitonic Hamiltonian. Such a Hamiltonian, which is often used for describing energy and charge transfer dynamics in complex molecular systems, is given in terms of diabatic electronic states that are coupled to each other and correspond to different nuclear Hamiltonians. Within the M-GQME approach, the effect of the nuclear degrees of freedom on the time evolution of the electronic density matrix is fully captured by a memory kernel superoperator, which can be obtained from short-lived (compared to the time scale of energy/charge transfer) projection-free inputs. In this paper, we test the ability of the M-GQME to predict the energy transfer dynamics within a seven-state benchmark model of the Fenna-Matthews-Olson (FMO) complex, with the short-lived projection-free inputs obtained via the Ehrenfest method. The M-GQME with Ehrenfest-based inputs is shown to yield accurate results across a wide parameter range. It is also found to dramatically outperform the direct application of the Ehrenfest method and to provide better-behaved convergence with respect to memory time in comparison to an alternative implementation of the GQME approach previously applied to the same FMO model.

Photoinduced charge transfer in Zn(II) and Au(III)-ligated symmetric and asymmetric bacteriochlorin dyads: A computational study

The excited-state properties and photoinduced charge-transfer (CT) kinetics in a series of symmetrical and asymmetrical Zn- and Au-ligated meso-meso-connected bacteriochlorin (BChl) complexes are studied computationally. BChl derivatives, which are excellent near-IR absorbing chromophores, are found to play a central role in bacterial photosynthetic reaction centers but are rarely used in artificial solar energy harvesting systems. The optical properties of chemically linked BChl complexes can be tuned by varying the linking group and involving different ligated metal ions. We investigate charge transfer in BChl dyads that are either directly linked or through a phenylene ring (1,4-phenylene) and which are ligating Zn or Au ions. The directly linked dyads with a nearly perpendicular arrangement of the BChl units bear markedly different properties than phenylene linked dyads. In addition, we find that the dielectric dependence of the intramolecular CT rate is very strong in neutral Zn-ligated dyads, whereas cationic Au-ligated dyads show negligible dielectric dependence of the CT rate. Rate constants of the photo induced CT process are calculated at the semiclassical Marcus level and are compared to fully quantum mechanical Fermi's golden rule based values. The rates are calculated using a screened range separated hybrid functional that offers a consistent framework for addressing environment polarization. We study solvated systems in two solvents of a low and a high scalar dielectric constant.

Three-state harmonic models for photoinduced charge transfer

A widely used strategy for simulating the charge transfer between donor and acceptor electronic states in an all-atom anharmonic condensed-phase system is based on invoking linear response theory to describe the system in terms of an effective spin-boson model Hamiltonian. Extending this strategy to photoinduced charge transfer processes requires also taking into consideration the ground electronic state in addition to the excited donor and acceptor electronic states. In this paper, we revisit the problem of describing such nonequilibrium processes in terms of an effective three-state harmonic model. We do so within the framework of nonequilibrium Fermi's golden rule (NE-FGR) in the context of photoinduced charge transfer in the carotenoid-porphyrin-C₆₀ (CPC₆₀) molecular triad dissolved in explicit tetrahydrofuran (THF). To this end, we consider different ways for obtaining a three-state harmonic model from the equilibrium autocorrelation functions of the donor-acceptor, donor-ground, and acceptor-ground energy gaps, as obtained from all-atom molecular dynamics simulations of the CPC₆₀/THF system. The quantum-mechanically exact time-dependent NE-FGR rate coefficients for two different charge transfer processes in two different triad conformations are then calculated using the effective three-state model Hamiltonians as well as a hierarchy of more approximate expressions that lead to the instantaneous Marcus theory limit. Our results show that the photoinduced charge transfer in CPC₆₀/THF can be described accurately by the effective harmonic three-state models and that nuclear quantum effects are small in this system.

Cyanide Bridged Platinum-Iron Complexes as Cisplatin Prodrug Systems: Design and Computational Study

The potential role of cyanide-bridged platinum-iron complexes as an anti-cancer Pt(IV) prodrug is studied. We present design principles of a dual-function prodrug that can upon reduction dissociate and release concurrently six cisplatin units and a ferricyanide anion per prodrug unit. The prodrug molecule is a unique complex of hepta metal centers consisting of a ferricyanide core with six Pt(IV) centers each bonded to the Fe(III) core through a cyano ligand. The functionality of the prodrug is addressed through density functional theory (DFT) calculations.

How Well Does a Solvated Octa-acid Capsule Shield the Embedded Chromophore? A Computational Analysis Based on an Anisotropic Dielectric Continuum Model

The optical properties of chromophores embedded in a water-solvated dimer of octa-acid that forms a molecular-shaped capsule are investigated. In particular, we address the anisotropic dielectric environment that appears to blue-shift excitation energies compared to the free aqueous chromophores. Recently we reported that using an effective scalar dielectric constant ε ≈ 3 appears to reproduce the measured spectra of the embedded coumarins, suggesting that the capsule provides a significant, albeit not perfect, screening of the aqueous dielectric environment. Here, we report absorption energies using a theoretical treatment that includes continuum solvation affected by an anisotropic dielectric function reflecting the high-dielectric environment outside of the capsule and the low-dielectric region within. We report time-dependent density functional theory calculations using a range-separated functional with the Poisson boundary conditions that model the anisotropic dielectric environment. Our calculations find that the anisotropic environment due to the water-solvated hydrophobic capsule is equivalent to a homogeneous effective dielectric constant of ≈3. The calculated values also appear to reproduce measured absorption of the embedded coumarin, where we study the effect of the hydrophobic capsule on the excited state.

Efficient Charge Generation via Hole Transfer in Dilute Organic Donor-Fullerene Blends

Efficient organic photovoltaics (OPVs) require broadband charge photogeneration with near-unity quantum yield. This can only be achieved by exploiting all pathways that generate charge. Electron transfer from organic donors to acceptors has been well-studied and is considered the primary path to charge photogeneration in OPVs. In contrast, much less is known about the hole transfer pathway. Here we study charge photogeneration in an archetypal system comprising tetraphenyldibenzoperiflanthene:C₇₀ blends using our recently developed multispectral two-dimensional electronic spectroscopy (M-2DES), supported by time-dependent density functional theory and fully quantum-mechanical Fermi's golden rule rate calculations. Our approach identifies in real time two rapid charge transfer pathways that are confirmed through computational analysis. Surprisingly, we find that both electron and hole transfer occur with comparable rates and efficiencies, facilitated by donor-acceptor electronic interactions. Our results highlight the importance of the hole transfer pathway for optimizing the efficiency of OPV devices employing small-molecule heterojunctions.

On the Role of the Special Pair in Photosystems as a Charge Transfer Rectifier

The special pair, a bacteriochlorophyll a (BChl) dimer found at the core of bacterial reaction centers, is known to play a key role in the functionality of photosystems as a precursor to the photosynthesis process. In this paper, we analyze the inherent affinity of the special pair to rectify the intrapair photo-induced charge transfer (CT). In particular, we show that the molecular environment affects the nuclear geometry, resulting in symmetry breaking between the two possible intrapair CT processes. To this end, we study the relationships of the intrapair CT and the molecular geometry with respect to the effective dielectric constant provided by the molecular environment. We identify the special pair structural feature that breaks the symmetry between the two molecules, leading to CT rectification. Excited state energies, oscillator strengths, and electronic coupling values are obtained via time-dependent density functional theory, employing a recently developed framework based on a screened range-separated hybrid functional within a polarizable continuum model (SRSH-PCM). We analyze the rectification capability of the special pair by calculating the CT rates using a first-principles-based Fermi's golden rule approach.

Vibronic structure of photosynthetic pigments probed by polarized two-dimensional electronic spectroscopy and ab initio calculations

Using polarized 2D spectroscopy and state-of-the-art TDDFT calculations to uncover the vibronic structure of primary photosynthetic pigments and its effect on ultrafast photoexcited dynamics.

Bacteriochlorophyll a (Bchl a) and chlorophyll a (Chl a) play important roles as light absorbers in photosynthetic antennae and participate in the initial charge-separation steps in photosynthetic reaction centers. Despite decades of study, questions remain about the interplay of electronic and vibrational states within the Q-band and its effect on the photoexcited dynamics. Here we report results of polarized two-dimensional electronic spectroscopic measurements, performed on penta-coordinated Bchl a and Chl a and their interpretation based on state-of-the-art time-dependent density functional theory calculations and vibrational mode analysis for spectral shapes. We find that the Q-band of Bchl a is comprised of two independent bands, that are assigned following the Gouterman model to Q_x and Q_y states with orthogonal transition dipole moments. However, we measure the angle to be ∼75°, a finding that is confirmed by ab initio calculations. The internal conversion rate constant from Q_x to Q_y is found to be 11 ps^–1. Unlike Bchl a, the Q-band of Chl a contains three distinct peaks with different polarizations. Ab initio calculations trace these features back to a spectral overlap between two electronic transitions and their vibrational replicas. The smaller energy gap and the mixing of vibronic states result in faster internal conversion rate constants of 38–50 ps^–1. We analyze the spectra of penta-coordinated Bchl a and Chl a to highlight the interplay between low-lying vibronic states and their relationship to photoinduced relaxation. Our findings shed new light on the photoexcited dynamics in photosynthetic systems where these chromophores are primary pigments.

Controlling the Emissive Activity in Heterocyclic Systems Bearing C═P Bonds

The photophysical properties of a series of heteroatom substituted indoles are explored to identify chemical means to control their emissive activity. In particular, we consider impacts of changes in the conjugated backbone, where the C═N bonds of benzoxazoles are replaced by C═P bonds (benzoxaphospholes). The effects of extending the π-conjugation, incorporating various secondary heteroatoms (X-C═P), and enforcing planar rigidity are also examined. Our computational analysis explains the higher fluorescence efficiency observed with extended π-conjugation and highlights the importance of maintaining molecular planarity at both ground- and emissive-state geometries.

Excitonic Interactions in Bacteriochlorin Homo-Dyads Enable Charge Transfer: A New Approach to the Artificial Photosynthetic Special Pair

Excitonically coupled bacteriochlorin (BC) dimers constitute a primary electron donor (special pair) in bacterial photosynthesis and absorbing units in light-harvesting antenna. However, the exact nature of the excited state of these dyads is still not fully understood. Here, we report a detailed spectroscopic and computational investigation of a series of symmetrical bacteriochlorin dimers, where the bacteriochlorins are connected either directly or by a phenylene bridge of variable length. The excited state of these dyads is quenched in high-dielectric solvents, which we attribute to photoinduced charge transfer. The mixing of charge transfer with the excitonic state causes accelerated (within 41 ps) decay of the excited state for the directly linked dyad, which is reduced by orders of magnitude with each additional phenyl ring separating the bacteriochlorins. These results highlight the origins of the excited-state dynamics in symmetric BC dyads and provide a new model for studying the primary processes in photosynthesis and for the development of artificial, biomimetic systems for solar energy conversion.

Enhancing charge mobilities in organic semiconductors by selective fluorination: a design approach based on a quantum mechanical perspective

Selective fluorination of organic semiconducting molecules is proposed as a means to achieving enhanced hole mobility. Naphthalene is examined here as a root molecular system with fluorination performed at various sites. Our quantum chemical calculations show that selective fluorination can enhance attractive intermolecular interactions while reducing charge trapping. Those observations suggest a design principle whereby fluorination is utilized for achieving high charge mobilities in the crystalline form. The utility of this design principle is demonstrated through an application to perylene, which is an important building block of organic semiconducting materials. We also show that a quantum mechanical perspective of nuclear degrees of freedom is crucial for a reliable description of charge transport.

Modification of Molecular Conductance by in Situ Deprotection of Thiol-Based Porphyrin

Acetylthio-protected free base porphyrins are used to form scanning tunneling microscope-molecular break junctions. The porphyrin molecules are deprotected in situ, before the self-assembly. Two types of molecular junctions are formed in the junctions: Au-S-Por-SAc-Au and Au-S-Por-S-Au. Lower conductance values and higher conductance values are observed. Computational modeling attributes the lower conductance to the Au-S-Por-SAc-Au junctions and the higher conductance to the Au-S-Por-S-Au junctions. First-principles calculation suggests that the reduced conductance in the protected porphyrin originates from the presence of the acetyl end groups (-COCH₃), rather than from the elongation of the sulfur-gold (S-Au) bonds at the tip-molecule interface.

Achieving Predictive Description of Molecular Conductance by Using a Range-Separated Hybrid Functional

The conductance of molecular bridges tends to be overestimated by computational studies in comparison to measured values. While this well-established trend may be related to difficulties for achieving robust bridges, the employed computational scheme can also contribute to this tendency. In particular, caveats of the traditional functionals employed in first-principles-based calculations can lead to discrepancies reflected in exaggerated conductance. Here, we show that by employing a range-separated hybrid functional the calculated values are within the same order as the measured conductance for all four considered cases. On the other hand, with B3LYP, which is a widely used functional, the calculated values greatly overestimate the conductance (by about 1-2 orders of magnitude). The improved description of the conductance with a RSH functional builds on achieving a physically meaningful treatment of the quasi particles associated with the frontier orbitals.

Deleterious Effects of Exact Exchange Functionals on Predictions of Molecular Conductance

Kohn-Sham (KS) density functional theory (DFT) describes well the atomistic structure of molecular junctions and their coupling to the semi-infinite metallic electrodes but severely overestimates conductance due to the spuriously large density of charge-carrier states of the KS system. Previous works show that inclusion of appropriate amounts of nonlocal exchange in the functional can fix the problem and provide realistic conductance estimates. Here however we discover that nonlocal exchange can also lead to deleterious effects which artificially overestimate transmittance even beyond the KS-DFT prediction. The effect is a result of exchange coupling between nonoverlapping states of diradical character. We prescribe a practical recipe for eliminating such artifacts.

Unraveling the Mechanism of Photoinduced Charge Transfer in Carotenoid-Porphyrin-C60 Molecular Triad

Photoinduced charge transfer (CT) plays a central role in biologically significant systems and in applications that harvest solar energy. We investigate the relationship of CT kinetics and conformation in a molecular triad. The triad, consisting of carotenoid, porphyrin, and fullerene is structurally flexible and able to acquire significantly varied conformations under ambient conditions. With an integrated approach of quantum calculations and molecular dynamics simulations, we compute the rate of CT at two distinctive conformations. The linearly extended conformation, in which the donor (carotenoid) and the acceptor (fullerene) are separated by nearly 50 Å, enables charge separation through a sequential CT process. A representative bent conformation that is entropically dominant, however, attenuates the CT, although the donor and the acceptor are spatially closer. Our computed rate of CT at the linear conformation is in good agreement with measured values. Our work provides unique fundamental understanding of the photoinduced CT process in the molecular triad.

Ultrafast Charge-Transfer Dynamics at the Boron Subphthalocyanine Chloride/C60 Heterojunction: Comparison between Experiment and Theory

Photoinduced charge-transfer (CT) processes play a key role in many systems, particularly those relevant to organic photovoltaics and photosynthesis. Advancing the understanding of CT processes calls for comparing their rates measured via state-of-the-art time-resolved interface-specific spectroscopic techniques with theoretical predictions based on first-principles molecular models. We measure charge-transfer rates across a boron subphthalocyanine chloride (SubPc)/C60 heterojunction, commonly used in organic photovoltaics, via heterodyne-detected time-resolved second-harmonic generation. We compare these results to theoretical predictions based on a Fermi's golden rule approach, with input parameters obtained using first-principles calculations for two different equilibrium geometries of a molecular donor-acceptor in a dielectric continuum model. The calculated rates (∼2 ps(-1)) overestimate the measured rates (∼0.1 ps(-1)), which is consistent with the expectation that the calculated rates represent an upper bound over the experimental ones. The comparison provides valuable understanding of how the structure of the electron donor-acceptor interface affects the CT kinetics in organic photovoltaic systems.

Molecular structure, spectroscopy, and photoinduced kinetics in trinuclear cyanide bridged complex in solution: a first-principles perspective

We investigate the molecular structure of the solvated complex, [(NC)6Fe-Pt(NH3)4-Fe(CN)6](4-), and related dinuclear and mononuclear model complexes using first-principles calculations. Mixed nuclear complexes in both solution and crystal phases were widely studied as models for charge transfer (CT) reactions using advanced spectroscopical and electrochemical tools. In contrast to earlier interpretations, we find that the most stable gas phase and solvated geometries are substantially different from the crystal phase geometry, mainly due to variance in the underlying oxidation numbers of the metal centers. Specifically, in the crystal phase a Pt(IV) metal center resulting from Fe ← Pt backward electron transfers is stabilized by an octahedral ligand field, whereas in the solution phase a Pt(II) metal complex that prefers a square planar ligand field forms a CT salt by bridging to the iron complexes through long-range electrostatic interactions. The different geometry is shown to be consistent with spectroscopical data and measured CT rates of the solvated complex. Interestingly, we find that the experimentally indicated photoinduced process in the solvated complex is of backward CT (Fe ← Pt).

Donor-to-Donor vs Donor-to-Acceptor Interfacial Charge Transfer States in the Phthalocyanine-Fullerene Organic Photovoltaic System

Charge transfer (CT) states formed at the donor/acceptor heterointerface are key for photocurrent generation in organic photovoltaics (OPV). Our calculations show that interfacial donor-to-donor CT states in the phthalocyanine-fullerene OPV system may be more stable than donor-to-acceptor CT states and that they may rapidly recombine, thereby constituting a potentially critical and thus far overlooked loss mechanism. Our results provide new insight into processes that may compete with charge separation, and suggest that the efficiency for charge separation may be improved by destabilizing donor-to-donor CT states or decoupling them from other states.