Excited state diabatization on the cheap using DFT: Photoinduced electron and hole transfer

Electronic Couplings and Conversion Dynamics between Localized and Charge Transfer Excitations from Many-Body Green's Functions Theory

We investigate the determination of electronic coupling between localized excitations (LEs) and charge-transfer (CT) excitations based on many-body Green's functions theory in the GW approximation with the Bethe-Salpeter equation (GW-BSE). Using a small molecule dimer system, we first study the influence of different diabatization methods, as well as different model choices within GW-BSE, such as the self-energy models or different levels of self-consistency, and find that these choices affect the LE-CT couplings only minimally. We then consider a large-scale low-donor morphology formed from rubrene and fullerene and evaluate the LE-CT couplings based on coupled GW-BSE-molecular mechanics calculations. For these disordered systems of bulky molecules, we observe differences in the couplings based on the Edmiston-Ruedenberg diabatization compared to the more approximate Generalize Mulliken-Hush and fragment charge difference diabatization formalisms. In a kinetic model for the conversion between LE and CT states, these differences affect the details of state populations in an intermediate time scale but not the final populations.

Eclipsed and Twisted Excimers of Pyrene and 2-Azapyrene: How Nitrogen Substitution Impacts Excimer Emission

Due to their unique photophysical and electronic properties, pyrene and its analogues have been the subject of extensive research in recent decades. The propensity of pyrene and its derivatives to form excimers has found wide application in various fields. Nitrogen-substituted pyrene derivatives display similar photophysical properties, but for them, excimer emission has not been reported to date. Here, we use time-dependent density functional theory (TD-DFT) calculations to investigate the low-lying exciton states of dimers of pyrene and 2-azapyrene. The excimer equilibrium structures are determined and the contribution of charge transfer (CT) excitations and intermolecular interactions to the exciton states is disclosed using a diabatization procedure. The study reveals that the dimers formed by the two molecules have quite similar exciton-state patterns, in which the relevant CT contributions govern the formation of excimer states, along with the La/Lb state inversion. In contrast with pyrene, the dipole-dipole interactions in 2-azapyrene stabilize the dark eclipsed excimer structure and increase the barrier for conversion into a bright twisted excimer. It is suggested that these differences in the nitrogen-substituted derivative might influence the excimer emission properties.

Tryptophan to Tryptophan Hole Hopping in an Azurin Construct

Electron transfer (ET) between neutral and cationic tryptophan residues in the azurin construct [ReI(H126)(CO)3(dmp)](W124)(W122)CuI (dmp = 4,7-Me2-1,10-phenanthroline) was investigated by Born-Oppenheimer quantum-mechanics/molecular mechanics/molecular dynamics (QM/MM/MD) simulations. We focused on W124•+ ← W122 ET, which is the middle step of the photochemical hole-hopping process *ReII(CO)3(dmp•-) ← W124 ← W122 ← CuI, where sequential hopping amounts to nearly 10,000-fold acceleration over single-step tunneling (ACS Cent. Sci. 2019, 5, 192-200). In accordance with experiments, UKS-DFT QM/MM/MD simulations identified forward and reverse steps of W124•+ ↔ W122 ET equilibrium, as well as back ET ReI(CO)3(dmp•-) → W124•+ that restores *ReII(CO)3(dmp•-). Strong electronic coupling between the two indoles (≥40 meV in the crossing region) makes the productive W124•+ ← W122 ET adiabatic. Energies of the two redox states are driven to degeneracy by fluctuations of the electrostatic potential at the two indoles, mainly caused by water solvation, with contributions from the protein dynamics in the W122 vicinity. ET probability depends on the orientation of Re(CO)3(dmp) relative to W124 and its rotation diminishes the hopping yield. Comparison with hole hopping in natural systems reveals structural and dynamics factors that are important for designing efficient hole-hopping processes.

Intermolecular Interactions and Charge Resonance Contributions to Triplet and Singlet Exciton States of Oligoacene Aggregates

Intermolecular interactions modulate the electro-optical properties of molecular materials and the nature of low-lying exciton states. Molecular materials composed by oligoacenes are extensively investigated for their semiconducting and optoelectronic properties. Here, we analyze the exciton states derived from time-dependent density functional theory (TDDFT) calculations for two oligoacene model aggregates: naphthalene and anthracene dimers. To unravel the role of inter-molecular interactions, a set of diabatic states is selected, chosen to coincide with local (LE) and charge-transfer (CT) excitations within a restricted orbital space including two occupied and two unoccupied orbitals for each molecular monomer. We study energy profiles and disentangle inter-state couplings to disclose the (CT) character of singlet and triplet exciton states and assess the influence of inter-molecular orientation by displacing one molecule with respect to the other along the longitudinal translation coordinate. The analysis shows that (CT) contributions are relevant, although comparably less effective for triplet excitons, and induce a non-negligible mixed character to the low-lying exciton states for eclipsed monomers and for small translational displacements. Such (CT) contributions govern the L_a/L_b state inversion occurring for the low-lying singlet exciton states of naphthalene dimer and contribute to the switch from H- to J-aggregate type of the strongly allowed B_b transition of both oligoacene aggregates.

Photochemical Method for Laser Absorption

During the laser application process, laser energy is usually converted into heat energy, causing high temperature, which affects the (high-speed) aircraft in routine flight. A completely novel photochemical method was investigated to potentially minimize the energy effect of the laser beam. Ag nanoparticles/C3N4 were synthesized by an ultra-low temperature reduced deposit method with Ag mean diameters of 5-25 nm for photofixation of N2. The absorption performance of laser can be improved by using appropriate charge density and small size Ag metal particles. The energy absorption rate was 7.1% over Ag/C3N4 (-40) at 5 mJ/cm2 of laser energy.

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

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

Diabatic States of Molecules

Quantitative simulations of electronically nonadiabatic molecular processes require both accurate dynamics algorithms and accurate electronic structure information. Direct semiclassical nonadiabatic dynamics is expensive due to the high cost of electronic structure calculations, and hence it is limited to small systems, limited ensemble averaging, ultrafast processes, and/or electronic structure methods that are only semiquantitatively accurate. The cost of dynamics calculations can be made manageable if analytic fits are made to the electronic structure data, and such fits are most conveniently carried out in a diabatic representation because the surfaces are smooth and the couplings between states are smooth scalar functions. Diabatic representations, unlike the adiabatic ones produced by most electronic structure methods, are not unique, and finding suitable diabatic representations often involves time-consuming nonsystematic diabatization steps. The biggest drawback of using diabatic bases is that it can require large amounts of effort to perform a globally consistent diabatization, and one of our goals has been to develop methods to do this efficiently and automatically. In this Feature Article, we introduce the mathematical framework of diabatic representations, and we discuss diabatization methods, including adiabatic-to-diabatic transformations and recent progress toward the goal of automatization.

Dynamic signatures of electronically nonadiabatic coupling in sodium hydride: a rigorous test for the symmetric quasi-classical model applied to realistic, ab initio electronic states in the adiabatic representation

Sodium hydride (NaH) in the gas phase presents a seemingly simple electronic structure making it a potentially tractable system for the detailed investigation of nonadiabatic molecular dynamics from both computational and experimental standpoints. The single vibrational degree of freedom, as well as the strong nonadiabatic coupling that arises from the excited electronic states taking on considerable ionic character, provides a realistic chemical system to test the accuracy of quasi-classical methods to model population dynamics where the results are directly comparable against quantum mechanical benchmarks. Using a simulated pump-probe type experiment, this work presents computational predictions of population transfer through the avoided crossings of NaH via symmetric quasi-classical Meyer-Miller (SQC/MM), Ehrenfest, and exact quantum dynamics on realistic, ab initio potential energy surfaces. The main driving force for population transfer arises from the ground vibrational level of the D¹Σ⁺ adiabatic state that is embedded in the manifold of near-dissociation C¹Σ⁺ vibrational states. When coupled through a sharply localized first-order derivative coupling most of the population transfers between t = 15 and t = 30 fs depending on the initially excited vibronic wavepacket. While quantum mechanical effects are expected due to the reduced mass of NaH, predictions of the population dynamics from both the SQC/MM and Ehrenfest models perform remarkably well against the quantum dynamics benchmark. Additionally, an analysis of the vibronic structure in the nonadiabatically coupled regime is presented using a variational eigensolver methodology.

Path Integrals for Nonadiabatic Dynamics: Multistate Ring Polymer Molecular Dynamics

This review focuses on a recent class of path-integral-based methods that simulate nonadiabatic dynamics in the condensed phase using only classical molecular dynamics trajectories in an extended phase space. Specifically, a semiclassical mapping protocol is used to derive an exact, continuous, Cartesian variable path-integral representation for the canonical partition function of a system in which multiple electronic states are coupled to nuclear degrees of freedom. Building on this exact statistical foundation, multistate ring polymer molecular dynamics methods are developed for the approximate calculation of real-time thermal correlation functions. The remarkable promise of these multistate ring polymer methods, their successful applications, and their limitations are discussed in detail.Expected final online publication date for the Annual Review of Physical Chemistry, Volume 73 is April 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Impact of Charge-Resonance Excitations on CT-Mediated J-Type Aggregation in Singlet and Triplet Exciton States of Perylene Di-Imide Aggregates: A TDDFT Investigation

The modulation of intermolecular interactions upon aggregation induces changes in excited state properties of organic molecules that can be detrimental for some optoelectronic applications but can be exploited for others. The time-dependent density functional theory (TDDFT) is a cost-effective approach to determining the exciton states of molecular aggregates, and it has been shown to provide reliable results when coupled with the appropriate choice of the functional. Here we apply a general procedure to analyze the aggregates’ exciton states derived from TDDFT calculations in terms of diabatic states chosen to coincide with local (LE) and charge-transfer (CT) excitations within a restricted orbital space. We apply the approach to study energy profiles, interstate couplings, and the charge-transfer character of singlet and triplet exciton states of perylene di-imide aggregates (PDI). We focus on the intermolecular displacement along the longitudinal translation coordinate, which mimics different amounts of slip-stacking observed in PDI crystals. The analysis, in terms of symmetry-adapted Frenkel excitations (FE) and charge-resonance (CR) states and their interactions, discloses how the interchange of the H/J character for small longitudinal shifts, previously reported for singlet exciton states, also occurs for triplet excitons.

Photoredox Chemistry with Organic Catalysts: Role of Computational Methods

Organic catalysts have the potential to carry out a wide range of otherwise thermally inaccessible reactions via photoredox routes. Early demonstrated successes of organic photoredox catalysts include one-electron CO₂ reduction and H₂ generation via water splitting. Photoredox systems are challenging to study and design owing to the sheer number and diversity of phenomena involved, including light absorption, emission, intersystem crossing, partial or complete charge transfer, and bond breaking or formation. Designing a viable photoredox route therefore requires consideration of a host of factors such as absorption wavelength, solvent, choice of electron donor or acceptor, and so on. Quantum chemistry methods can play a critical role in demystifying photoredox phenomena. Using one-electron CO₂ reduction with phenylene-based chromophores as an illustrative example, this perspective highlights recent developments in quantum chemistry that can advance our understanding of photoredox processes and proposes a way forward for driving the design and discovery of organic catalysts.

Addressing the Frenkel and charge transfer character of exciton states with a model Hamiltonian based on dimer calculations: Application to large aggregates of perylene bisimide

To understand the influence of interchromophoric arrangements on photo-induced processes and optical properties of aggregates, it is fundamental to assess the contribution of local excitations [charge transfer (CT) and Frenkel (FE)] to exciton states. Here, we apply a general procedure to analyze the adiabatic exciton states derived from time-dependent density functional theory calculations, in terms of diabatic states chosen to coincide with local excitations within a restricted orbital space. In parallel, motivated by the need of cost-effective approaches to afford the study of larger aggregates, we propose to build a model Hamiltonian based on calculations carried out on dimers composing the aggregate. Both approaches are applied to study excitation energy profiles and CT character modulation induced by interchromophore rearrangements in perylene bisimide aggregates up to a tetramer. The dimer-based approach closely reproduces the results of full-aggregate calculations, and an analysis in terms of symmetry-adapted diabatic states discloses the effects of CT/FE interactions on the interchange of the H-/J-character for small longitudinal shifts of the chromophores.