Nonadiabatic Exciton and Charge Separation Dynamics at Interfaces of Zinc Phthalocyanine and Fullerene: Orientation Does Matter

Ultrafast Hot Exciton Nonadiabatic Excited-State Dynamics in Green Fluorescent Protein Chromophore Analogue

The ultrafast high-energy nonadiabatic excited-state dynamics of the benzylidenedimethylimidazolinone chromophore dimer has been investigated using an electronic structure method coupled with on-the-fly quantitative wave function analysis to gain insight into the photophysics of hot excitons in biological systems. The dynamical simulation provides a rationalization of the behavior of the exciton in a dimer after the photoabsorption of light to higher-energy states. The results suggest that hot exciton localization within the manifold of excited states is caused by the hindrance of torsional rotation due to imidazolinone (I) or phenolate (P) bonds i.e., ΦI- or ΦP-dihedral rotation, in the monomeric units of a dimer. This hindrance arises due to weak π-π stacking interaction in the dimer, resulting in an energetically uphill excited-state barrier for ΦI- and ΦP-twisted rotation, impeding the isomerization process in the chromophore. Thus, this study highlights the potential impact of the weak π-π interaction in regulating the photodynamics of the green fluorescent protein chromophore derivatives.

Theoretical investigation of distal charge separation in a perylenediimide trimer

An exciton-phonon (ex-ph) model based on our recently developed block interaction product basis framework is introduced to simulate the distal charge separation (CS) process in aggregated perylenediimide (PDI) trimer incorporating the quantum dynamic method, i.e., the time-dependent density matrix renormalization group. The electronic Hamiltonian in the ex-ph model is represented by nine constructed diabatic states, which include three local excited (LE) states and six charge transfer (CT) states from both the neighboring and distal chromophores. These diabatic states are automatically generated from the direct products of the leading localized neutral or ionic states of each chromophore's reduced density matrix, which are obtained from ab initio quantum chemical calculation of the subsystem consisting of the targeted chromophore and its nearest neighbors, thus considering the interaction of the adjacent environment. In order to quantum-dynamically simulate the distal CS process with massive coupled vibrational modes in molecular aggregates, we used our recently proposed hierarchical mapping approach to renormalize these modes and truncate those vibrational modes that are not effectively coupled with electronic states accordingly. The simulation result demonstrates that the formation of the distal CS process undergoes an intermediate state of adjacent CT, i.e., starts from the LE states, passes through an adjacent CT state to generate the intermediates (∼200 fs), and then formalizes the targeted distal CS via further charge transference (∼1 ps). This finding agrees well with the results observed in the experiment, indicating that our scheme is capable of quantitatively investigating the CS process in a realistic aggregated PDI trimer and can also be potentially applied to exploring CS and other photoinduced processes in larger systems.

Exciton Dynamics of TiOPc/WSe2 Heterostructure

The van der Waals (vdW) heterostructures composed of two-dimensional (2D) transition metal dichalcogenides (TMDs) and organic semiconductors demonstrate numerous compelling optoelectronic properties. However, the influence of the vdW epitaxial effect and temperature on the optoelectronic properties and interface exciton dynamics of heterostructures remains unclear. This study systematically investigates the fluorescence properties of TiOPc/WSe2 heterostructure. Comprehensive spectral characterization elucidates that the emission behavior of the TiOPc/WSe2 heterostructure arises from charge/energy transfer at the heterostructure interfaces and the structural ordering of the organic layer on the 2D monolayer WSe2 induced by vdW epitaxy. The interface exciton dynamic features probed by ultrafast transient spectroscopy reveal that the face-to-face molecular stacking configuration of TiOPc exhibits ultrafast exciton dynamics. In particular, we observe picosecond-scale absorption of organic molecular dimer cations, providing direct evidence of interface charge transfer at room temperature. Moreover, energy transfer from the TiOPc to WSe2 may exist based on the tunability in the fluorescence emission of the TiOPc/WSe2 heterostructure as the temperature changes. This study unveils the critical role of vdW epitaxy and temperature in the exciton dynamics of organic/2D TMDs hybrid systems and provides guidance for studying interlayer charge and energy transfer in organic/inorganic heterostructures.

Consistent Construction of the Density Matrix from Surface Hopping Trajectories

Proper construction of the density matrix based on surface hopping trajectories remains a difficult problem. Due to the well-known overcoherence in traditional surface hopping simulations, the electronic wave function cannot be used directly. In this work, we propose a consistent density matrix construction method, which takes the advantage of occupation of active states to rescale the coherence calculated by wave functions and ensures the intrinsic consistency of the density matrix. This new trajectory analysis method can be used for both Tully's fewest switches surface hopping (FSSH) and our recently proposed branching corrected surface hopping (BCSH). As benchmarked in both one- and two-dimensional standard scattering models, the new approach combined with BCSH trajectories achieves highly accurate time-dependent spatial distributions of adiabatic populations and coherence compared to exact quantum results.

Nonadiabatic Dynamics Simulations for Photoinduced Processes in Molecules and Semiconductors: Methodologies and Applications

Nonadiabatic dynamics (NAMD) simulations have become powerful tools for elucidating complicated photoinduced processes in various systems from molecules to semiconductor materials. In this review, we present an overview of our recent research on photophysics of molecular systems and periodic semiconductor materials with the aid of ab initio NAMD simulation methods implemented in the generalized trajectory surface-hopping (GTSH) package. Both theoretical backgrounds and applications of the developed NAMD methods are presented in detail. For molecular systems, the linear-response time-dependent density functional theory (LR-TDDFT) method is primarily used to model electronic structures in NAMD simulations owing to its balanced efficiency and accuracy. Moreover, the efficient algorithms for calculating nonadiabatic coupling terms (NACTs) and spin-orbit couplings (SOCs) have been coded into the package to increase the simulation efficiency. In combination with various analysis techniques, we can explore the mechanistic details of the photoinduced dynamics of a range of molecular systems, including charge separation and energy transfer processes in organic donor-acceptor structures, ultrafast intersystem crossing (ISC) processes in transition metal complexes (TMCs), and exciton dynamics in molecular aggregates. For semiconductor materials, we developed the NAMD methods for simulating the photoinduced carrier dynamics within the framework of the Kohn-Sham density functional theory (KS-DFT), in which SOC effects are explicitly accounted for using the two-component, noncollinear DFT method. Using this method, we have investigated the photoinduced carrier dynamics at the interface of a variety of van der Waals (vdW) heterojunctions, such as two-dimensional transition metal dichalcogenides (TMDs), carbon nanotubes (CNTs), and perovskites-related systems. Recently, we extended the LR-TDDFT-based NAMD method for semiconductor materials, allowing us to study the excitonic effects in the photoinduced energy transfer process. These results demonstrate that the NAMD simulations are powerful tools for exploring the photodynamics of molecular systems and semiconductor materials. In future studies, the NAMD simulation methods can be employed to elucidate experimental phenomena and reveal microscopic details as well as rationally design novel photofunctional materials with desired properties.

Solvent effects on the photoinduced charge separation dynamics of directly linked zinc phthalocyanine-perylenediimide dyads: a nonadiabatic dynamics simulation with an optimally tuned screened range-separated hybrid functional

Herein, we have employed a combination of the optimally tuned screened range-separated hybrid (OT-SRSH) functional, the polarizable continuum model (PCM), and nonadiabatic dynamics (NAMD) simulations to investigate the photoinduced dynamics of directly linked donor-acceptor dyads formed using zinc phthalocyanine (ZnPc) and perylenediimide (PDI), in which ZnPc is the donor while PDI is the acceptor. Our simulations aim to analyze the behavior of these dyads upon local excitation of the ZnPc moiety in the gas phase and in benzonitrile. Our findings indicate that the presence of a solvent can significantly influence the excited state dynamics of ZnPc-PDI dyads. Specifically, the polar solvent benzonitrile effectively lowers the vertical excitation energies of the charge transfer (CT) state from ZnPc to PDI. As a result, the energetic order of the locally excited (LE) states of ZnPc and the CT states is reversed compared to the gas phase. Consequently, the photoinduced electron transfer (PET) dynamics from ZnPc to PDI, which is absent in the gas phase, takes place in benzonitrile with a time constant of 10.4 ps. Importantly, our present work not only qualitatively agrees with experimental results but also provides in-depth insights into the underlying mechanisms responsible for the photoinduced dynamics of ZnPc-PDI. Moreover, this study emphasizes the importance of appropriately considering solvent effects in NAMD simulation of organic donor-acceptor systems, taking into account the distinct excited state dynamics observed in the gas phase and benzonitrile. Furthermore, the combination of the OT-SRSH functional, the PCM solvent model, and nonadiabatic dynamics simulations shows promise as a strategy for investigating the complex excited state dynamics of organic donor-acceptor systems in solvents. These findings will be valuable for the future design of novel organic donor-acceptor structures with improved performance.

Non-negligible roles of charge transfer excitons in ultrafast excitation energy transfer dynamics of a double-walled carbon nanotube

Herein, we employed a developed linear response time dependent density functional theory-based nonadiabatic dynamics simulation method that explicitly takes into account the excitonic effects to investigate photoinduced excitation energy transfer dynamics of a double-walled carbon nanotube (CNT) model with different excitation energies. The E₁₁ excitation of the outer CNT will generate a local excitation (LE) |out*〉 exciton due to its low energy, which does not induce any charge separation. In contrast, the E₁₁ excitation of the inner CNT can generate four kinds of excitons with the LE exciton |in*〉 dominates. In the 500-fs dynamics simulation, the LE exciton |in*〉 and charge transfer (CT) excitons |out^-in⁺〉 and |out⁺in^-〉 are all gradually converted to the |out*〉 exciton, corresponding to a photoinduced excitation energy transfer, which is consistent with experimental studies. Finally, when the excitation energy is close to the E₂₂ state of the outer CNT (∼1.05 eV), a mixed population of different excitons, with the |out*〉 exciton dominated, is generated. Then, photoinduced energy transfer from the outer to inner CNTs occurs in the first 50 fs, which is followed by an inner to outer excitation energy transfer that is completed in 400 fs. The present work not only sheds important light on the mechanistic details of wavelength-dependent excitation energy transfer of a double-walled CNT model but also demonstrates the roles and importance of CT excitons in photoinduced excitation energy transfer. It also emphasized that explicitly including the excitonic effects in electronic structure calculations and nonadiabatic dynamics simulations is significant for correct understanding/rational design of optoelectronic properties of periodically extended systems.

Generalized Ab Initio Nonadiabatic Dynamics Simulation Methods from Molecular to Extended Systems

Nonadiabatic dynamics simulation has become a powerful tool to describe nonadiabatic effects involved in photophysical processes and photochemical reactions. In the past decade, our group has developed generalized trajectory-based ab initio surface-hopping (GTSH) dynamics simulation methods, which can be used to describe a series of nonadiabatic processes, such as internal conversion, intersystem crossing, excitation energy transfer and charge transfer of molecular systems, and photoinduced nonadiabatic carrier dynamics of extended systems with and without spin-orbit couplings. In this contribution, we will first give a brief introduction to our recently developed methods and related numerical implementations at different computational levels. Later, we will present some of our latest applications in realistic systems, which cover organic molecules, biological proteins, organometallic compounds, periodic organic and inorganic materials, etc. Final discussion is given to challenges and outlooks of ab initio nonadiabatic dynamics simulations.

Solvent effects on excited-state relaxation dynamics of paddle-wheel BODIPY-Hexaoxatriphenylene conjugates: Insights from non-adiabatic dynamics simulations

Understanding the excited state dynamics of donor-acceptor (D-A) complexes is of fundamental importance both experimentally and theoretically. Herein, we have first explored the photoinduced dynamics of a recently synthesized paddle-wheel BODIPY-hexaoxatriphenylene (BODIPY is the abbreviation for BF2-chelated dipyrromethenes) conjugates D-A complexes with the combination of both electronic structure calculations and non-adiabatic dynamics simulations. On the basis of computational results, we concluded that the BODIPY-hexaoxatriphenylene (BH) conjugates will be promoted to the local excited (LE) states of the BODIPY fragments upon excitation, which is followed by the ultrafast exciton transfer from LE state to charge transfer (CT). Instead of the photoinduced electron transfer process proposed in previous experimental work, such a exciton transfer process is accompanied with the photoinduced hole transfer from BODIPY to hexaoxatriphenylene. Additionally, solvent effects are found to play an important role in the photoinduced dynamics. Specifically, the hole transfer dynamics is accelerated by the acetonitrile solvent, which can be ascribed to significant influences of the solvents on the charge transfer states, i.e. the energy gaps between LE and CT excitons are reduced greatly and the non-adiabatic couplings are increased in the meantime. Our present work not only provides valuable insights into the underlying photoinduced mechanism of BH, but also can be helpful for the future design of novel donor-acceptor conjugates with better optoelectronic performance.

Ultrafast charge transfer in a nonfullerene all-small-molecule organic solar cell: a nonadiabatic dynamics simulation with optimally tuned range-separated functional

The combination of nonadiabatic dynamics simulation and optimally tuned range-separated functional might be a powerful tool for elucidating the ultrafast charge transfer in nonfullerene all-small-molecule organic solar cells.

Ultrafast Exciton Delocalization and Localization Dynamics of a Perylene Bisimide Quadruple π-­Stack: A Nonadiabatic Dynamics Simulation

Unraveling the photogenerated exciton dynamics of π­stacked molecular aggregates is of great importance for both fundamental studies and industrial applications. Among various π­stacked molecular aggregates, perylene tetracarboxylic acid bisimides (PBI)...

Vibrationally Resolved Absorption Spectra and Exciton Dynamics in Zinc Phthalocyanine Aggregates: Effects of Aggregation Lengths and Remote Exciton Transfer

The vibrationally resolved absorption spectra and exciton dynamics in the α-zinc phthalocyanine aggregates are theoretically investigated by using a non-Markovian stochastic Schrödinger equation. The model Hamiltonian adopted for spectral and dynamic simulations explicitly includes the couplings for both nearest-neighbor and remote exciton transfer, and it is parametrized from first-principles calculations. The results indicate that aggregation lengths and remote exciton transfer significantly influence the relative energy alignment between delocalized Frenkel exciton (FE) and charge transfer (CT) states, which in turn strongly affects the relative intensities of the two absorption peaks in the Q-band region. Analytical formulas are derived to establish quantitative structure-spectra relationships in aggregates, and they offer simple patterns to extract electronic-state properties directly from absorption spectra. The dynamics simulations reveal that the light absorption can directly generate mixed states with both FE and CT features, but it is hard for the photoexcitation from the Q-band region to generate free carriers due to the high energies of charge-separated states.

Chemical Bonding as a New Avenue for Controlling Excited-State Properties and Excitation Energy-Transfer Processes in Zinc Phthalocyanine-Fullerene Dyads

Whether chemical bonding can regulate the excited-state and optoelectronic properties of donor-acceptor dyads has been largely elusive. In this work, we used electronic structure and nonadiabatic dynamics methods to explore the excited-state properties of covalently bonded zinc phthalocyanine (ZnPc)-fullerene (C₆₀ ) dyads with a 6-6 (or 5-6) bonding configuration in which ZnPc is bonded to two carbon atoms shared by the two hexagonal rings (or a pentagonal and a hexagonal ring) in C₆₀ . In both cases, the locally excited (LE) states on ZnPc are spectroscopically bright. However, their different chemical bonding differentiates the electronic interactions between ZnPc and C₆₀ . In the 5-6 bonding configuration, the LE states on ZnPc are much higher in energy than the LE states on C₆₀ . Thus, the excitation energy transfer from ZnPc to C₆₀ is thermodynamically favorable. On the other hand, in the 6-6 bonding configuration, such a process is inhibited because the LE states on ZnPc are the lowest ones. More detailed mechanisms are elucidated from nonadiabatic dynamics simulations. In the 6-6 bonding configuration, no excitation energy transfer was observed. In contrast, in the 5-6 bonding configuration, several LE and charge-transfer (CT) excitons were shown to participate in the energy-transfer process. Further analysis reveals that the photoinduced energy transfer is mediated by a CT exciton, such that electron- and hole-transfer processes take place in a concerted but asynchronous manner in the excitation energy transfer. It is also found that high-level electronic structure methods including exciton effects are indispensable to accurately describe photoinduced energy- and electron-transfer processes. Furthermore, this work opens up new avenues for regulating the excited-state properties of molecular donor-acceptor dyads by means of chemical bonding.

Ultrafast channel I and channel II charge generation processes at a nonfullerene donor-acceptor PTB7:PDI interface is crucial for its excellent photovoltaic performance

Nonfullerene organic solar cells have received much attention in recent years due to their low cost, high absorption coefficient and excellent synthetic flexibility. However, the microscopic photoinduced dynamics at corresponding donor-acceptor interfaces remains unclear. In this work, we have firstly employed state-of-the-art TDDFT-based nonadiabatic dynamics simulations in combination with static electronic structure calculations to explore the ultrafast photoinduced dynamics at a typical nonfullerene donor-acceptor PTB7:PDI interface using a minimal model system (172 atoms). Upon excitation with specific wavelength of light, both PTB7 and PDI can be locally excited to generate |PTB7* and |PDI* excitons due to their high absorption ability and significant overlap in absorption spectrum. After that, these localized excitons gradually convert to charge transfer exciton |PTB7+PDI-, while another |PTB7-PDI+ charge transfer exciton is not involved in the whole process. Along with the exciton conversion, electron transfer from PTB7 to PDI (channel I charge generation) and the hole transfer from PDI to PTB7 (channel II charge generation) occurs simultaneously with time constants of 643 fs and 549 fs respectively. In the same time, D index that measures the centroid distance of electron and hole increases from 1.0 Å to 4.0 Å, which clearly reflects a charge transfer process at the interface. Our present work provides solid evidence that both channel I and channel II charge generation processes play important roles at PTB7:PDI interface, which could be helpful for the design of novel nonfullerene solar cells with better photovoltaic performance.