Discovery and characterization of an acridine radical photoreductant.

Photoinduced electron transfer (PET) is a phenomenon whereby the absorption of light by a chemical species provides an energetic driving force for an electron-transfer reaction1-4. This mechanism is relevant in many areas of chemistry, including the study of natural and artificial photosynthesis, photovoltaics and photosensitive materials. In recent years, research in the area of photoredox catalysis has enabled the use of PET for the catalytic generation of both neutral and charged organic free-radical species. These technologies have enabled previously inaccessible chemical transformations and have been widely used in both academic and industrial settings. Such reactions are often catalysed by visible-light-absorbing organic molecules or transition-metal complexes of ruthenium, iridium, chromium or copper5,6. Although various closed-shell organic molecules have been shown to behave as competent electron-transfer catalysts in photoredox reactions, there are only limited reports of PET reactions involving neutral organic radicals as excited-state donors or acceptors. This is unsurprising because the lifetimes of doublet excited states of neutral organic radicals are typically several orders of magnitude shorter than the singlet lifetimes of known transition-metal photoredox catalysts7-11. Here we document the discovery, characterization and reactivity of a neutral acridine radical with a maximum excited-state oxidation potential of -3.36 volts versus a saturated calomel electrode, which is similarly reducing to elemental lithium, making this radical one of the most potent chemical reductants reported12. Spectroscopic, computational and chemical studies indicate that the formation of a twisted intramolecular charge-transfer species enables the population of higher-energy doublet excited states, leading to the observed potent photoreducing behaviour. We demonstrate that this catalytically generated PET catalyst facilitates several chemical reactions that typically require alkali metal reductants and can be used in other organic transformations that require dissolving metal reductants.

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.

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.

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.

Electronic Spectra of C60 Films Using Screened Range Separated Hybrid Functionals.

We study computationally the electronic spectra of C60 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.

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-C60 (CPC60) 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 CPC60/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 CPC60/THF can be described accurately by the effective harmonic three-state models and that nuclear quantum effects are small in this 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.

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.

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.

Charge transfer rate constants for the carotenoid-porphyrin-C60 molecular triad dissolved in tetrahydrofuran: The spin-boson model vs the linearized semiclassical approximation.

Charge transfer rate constants were calculated for the carotenoid-porphyrin-C60 (CPC60) molecular triad dissolved in explicit tetrahydrofuran. The calculation was based on mapping the all-atom anharmonic Hamiltonian of this system onto the spin-boson Hamiltonian. The mapping was based on discretizing the spectral density from the time correlation function of the donor-acceptor potential energy gap, as obtained from all-atom molecular dynamics simulations. Different spin-boson Hamiltonians were constructed for each of the possible transitions between the three excited electronic states in two different triad conformations. The rate constants of three possible transitions were calculated via the quantum-mechanically exact Fermi's golden rule (FGR), as well as a progression of more approximate expressions that lead to the classical Marcus expression. The advantage of the spin-boson approach is that once the mapping is established, the quantum-mechanically exact FGR and the hierarchy of approximations are known in closed form. The classical Marcus charge transfer rate constants obtained with the spin-boson Hamiltonians were found to reproduce those obtained from all-atom simulations with the linearized semiclassical approximation, thereby confirming the equivalence of the two approaches for this system. Within the spin-boson Hamiltonian, we also found that the quantum-mechanically exact FGR rate constants were significantly enhanced compared to the classical Marcus theory rate constants for two out of three transitions in one of the two conformations under consideration. The results confirm that mapping to the spin-boson model can yield accurate predictions for charge transfer rate constants in a system as complex as CPC60 dissolved in tetrahydrofuran.

A modified approach for simulating electronically nonadiabatic dynamics via the generalized quantum master equation.

We present a modified approach for simulating electronically nonadiabatic dynamics based on the Nakajima-Zwanzig generalized quantum master equation (GQME). The modified approach utilizes the fact that the Nakajima-Zwanzig formalism does not require casting the overall Hamiltonian in system-bath form, which is arguably neither natural nor convenient in the case of the Hamiltonian that governs nonadiabatic dynamics. Within the modified approach, the effect of the nuclear degrees of freedom on the time evolution of the electronic reduced density operator is fully captured by a memory kernel super-operator. A methodology for calculating the memory kernel from projection-free inputs is developed. Simulating the electronic dynamics via the modified approach, with a memory kernel obtained using exact or approximate methods, can be more cost effective and/or lead to more accurate results than direct application of those methods. The modified approach is compared to previously proposed GQME-based approaches, and its robustness and accuracy are demonstrated on a benchmark spin-boson model with a memory kernel which is calculated within the Ehrenfest method.

Combining the mapping Hamiltonian linearized semiclassical approach with the generalized quantum master equation to simulate electronically nonadiabatic molecular dynamics.

The generalized quantum master equation (GQME) provides a powerful framework for simulating electronically nonadiabatic molecular dynamics. Within this framework, the effect of the nuclear degrees of freedom on the time evolution of the electronic reduced density matrix is fully captured by a memory kernel superoperator. In this paper, we consider two different procedures for calculating the memory kernel of the GQME from projection-free inputs obtained via the combination of the mapping Hamiltonian (MH) approach and the linearized semiclassical (LSC) approximation. The accuracy and feasibility of the two procedures are demonstrated on the spin-boson model. We find that although simulating the electronic dynamics by direct application of the two LSC-based procedures leads to qualitatively different results that become increasingly less accurate with increasing time, restricting their use to calculating the memory kernel leads to an accurate description of the electronic dynamics. Comparison with a previously proposed procedure for calculating the memory kernel via the Ehrenfest method reveals that MH/LSC methods produce memory kernels that are better behaved at long times and lead to more accurate electronic dynamics.

A comparative study of different methods for calculating electronic transition rates.

We present a comprehensive comparison of the following mixed quantum-classical methods for calculating electronic transition rates: (1) nonequilibrium Fermi's golden rule, (2) mixed quantum-classical Liouville method, (3) mean-field (Ehrenfest) mixed quantum-classical method, and (4) fewest switches surface-hopping method (in diabatic and adiabatic representations). The comparison is performed on the Garg-Onuchic-Ambegaokar benchmark charge-transfer model, over a broad range of temperatures and electronic coupling strengths, with different nonequilibrium initial states, in the normal and inverted regimes. Under weak to moderate electronic coupling, the nonequilibrium Fermi's golden rule rates are found to be in good agreement with the rates obtained via the mixed quantum-classical Liouville method that coincides with the fully quantum-mechanically exact results for the model system under study. Our results suggest that the nonequilibrium Fermi's golden rule can serve as an inexpensive yet accurate alternative to Ehrenfest and the fewest switches surface-hopping methods.

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.

The Effect of Interfacial Geometry on Charge-Transfer States in the Phthalocyanine/Fullerene Organic Photovoltaic System.

The dependence of charge-transfer states on interfacial geometry at the phthalocyanine/fullerene organic photovoltaic system is investigated. The effect of deviations from the equilibrium geometry of the donor-donor-acceptor trimer on the energies of and electronic coupling between different types of interfacial electronic excited states is calculated from first-principles. Deviations from the equilibrium geometry are found to destabilize the donor-to-donor charge transfer states and to weaken their coupling to the photoexcited donor-localized states, thereby reducing their ability to serve as charge traps. At the same time, we find that the energies of donor-to-acceptor charge transfer states and their coupling to the donor-localized photoexcited states are either less sensitive to the interfacial geometry or become more favorable due to modifications relative to the equilibrium geometry, thereby enhancing their ability to serve as gateway states for charge separation. Through these findings, we eludicate how interfacial geometry modifications can play a key role in achieving charge separation in this widely studied organic photovoltaic system.

Photoinduced Homolytic Bond Cleavage of the Central Si-C Bond in Porphyrin Macrocycles Is a Charge Polarization Driven Process.

Photoinduced cleavage of the bond between the central Si atom in porphyrin macrocycles and the neighboring carbon atom of an axial alkyl ligand is investigated by both experimental and computational tools. Photolysis and electron paramagnetic resonance measurements indicate that the Si-C bond cleavage of Si-phthalocyanine occurs through a homolytic process. The homolytic process follows a low-lying electronic excitation of about 1.8 eV that destabilizes the carbide bond of similar bond dissociation energy. Using electronic structure calculations, we provide insight into the nature of the excited state and the resulting photocleavage mechanism. We explain this process by finding that the electronic excited state is of a charge transfer character from the axial ligand toward the macrocycle in the reverse direction of the ground state polarization. We find that the homolytic process yielding the radical intermediate is energetically the most stable mechanistic route. Furthermore, we demonstrate using our computational approach that changing the phthalocyanine to smaller ring system enhances the homolytic photocleavage of the Si-C bond by reducing the energetic barrier in the relevant excited states.

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).

Communication: Charge-transfer rate constants in zinc-porphyrin-porphyrin-derived dyads: a Fermi golden rule first-principles-based study.

We investigate photoinduced charge transfer (CT) processes within dyads consisting of porphyrin derivatives in which one ring ligates a Zn metal center and where the rings vary by their degree of conjugation. Using a first-principles approach, we show that molecular-scale means can tune CT rates through stabilization affected by the polar environment. Such means of CT tuning are important for achieving high efficiency optoelectronic applications using organic semiconducting materials. Our fully quantum mechanical scheme is necessary for reliably modeling the CT process across different regimes, in contrast to the pervading semi-classical Marcus picture that grossly underestimates transfer in the far-inverted regime.