Theoretical Studies of Photoinduced Electron Transfer in Dye-Sensitized TiO2

Data-Driven Design of Electroactive Spacer Molecules to Tune Charge Carrier Dynamics in Layered Halide Perovskite Heterostructures

Crafting rational heterojunctions with nanostructured materials is instrumental in fostering effective interfacial charge separation and transport for optoelectronics. Layered halide perovskites (LHPs) that form heterojunctions between organic spacer molecules and inorganic metal halide layers exhibit tunable photophysics owing to their customizable band alignment. However, controlling photogenerated carrier dynamics by strategically designing layered perovskite heterojunctions remains largely unexplored. We combine a data-driven approach with time-domain density functional theory (TD-DFT) and non-adiabatic molecular dynamics (NAMD) to screen and select electronically active spacer dications (A') that introduce a type-II heterojunction in the lead iodide-based Dion-Jacobson phase LHPs. The composition-structure-electronic property correlations reveal that the number of nitrogens in aromatic heterocycles is the key factor in designing electron-accepting spacers in these perovskites. The detailed atomistic simulations validate the design strategy further by modeling (A')PbI4 perovskites, which incorporate three different screened electroactive A' spacers. The computed excited charge carrier dynamics illustrate the phonon-mediated ultrafast interfacial electron transfer from the inorganic conduction band edge to the lower-lying unoccupied orbitals of spacers, exhibiting photoluminescence quenching in these (A')PbI4 perovskites. The spatially separated electrons and holes at the type-II heterojunction interface prolong the excited charge carrier lifetime, boosting the carrier transport and exciton dynamics. Our work illustrates a robust in silico approach for designing LHPs with exciting optoelectronic properties originating from their fine-tuned heterojunctions.

Decoherence ensures convergence of non-adiabatic molecular dynamics with number of states

Non-adiabatic (NA) molecular dynamics (MD) is a powerful approach for studying far-from-equilibrium quantum dynamics in photophysical and photochemical systems. Most NA-MD methods are developed and tested with few-state models, and their validity with complex systems involving many states is not well studied. By modeling intraband equilibration and interband recombination of charge carriers in MoS2, we investigate the convergence of three popular NA-MD algorithms, fewest switches surface hopping (FSSH), global flux surface hopping (GFSH), and decoherence induced surface hopping (DISH) with the number of states. Only the standard DISH algorithm converges with the number of states and produces Boltzmann equilibrium. Unitary propagation of the wave function in FSSH and GFSH violates the Boltzmann distribution, leads to internal inconsistency between time-dependent Schrödinger equation state populations and trajectory counts, and produces non-convergent results. Introducing decoherence in FSSH and GFSH by collapsing the wave function fixes these problems. The simplified version of DISH that omits projecting out the occupied state and is applicable to few-state systems also causes problems when the number of states is increased. We discuss the algorithmic application of wave function collapse and Boltzmann detailed balance and provide detailed FSSH, GFSH, and DISH flow charts. The use of convergent NA-MD methods is highly important for modeling complicated quantum processes involving multiple states. Our findings provide the basis for investigating quantum dynamics in realistic complex systems.

Wavelength-Dependent Activity of Oxygen Species in Propane Conversion on Rutile TiO2(110)

Photocatalytic oxidative dehydrogenation of propane (C3H8) into propene (C3H6) under mild conditions holds great potential in the chemical industry, but understanding how active species participate in C3H8 conversion remains a significant challenge. Here, the wavelength-dependent activities of bridging oxygen (Ob2-) and the Ti5c-bound oxygen adatom (OTi2-) of model rutile (R) TiO2(110) in C3H8 conversion have been investigated. Under 257 and 343 nm irradiation, hole-trapped OTi- and Ob- can abstract the hydrogen atom of C3H8, forming the CH3CH•CH3 radical and C3H6. However, the rate of C3H8 conversion with hole-trapped Ob- is strongly dependent on the wavelength, primarily producing the C3H7• radical. In the case of hole-trapped OTi-, C3H6 is the main product, which is nearly independent of wavelength. The differences in the wavelength-dependent activity and product selectivity are likely due to dynamic control rather than thermodynamic control. The result provides a deeper understanding of the dynamic processes involved in the conversion of light alkanes in TiO2 photocatalysis.

Emergence of Z-Scheme Photocatalysis for Total Water Splitting: An Improvised Route to High Efficiency

Photocatalytic water splitting to spontaneously produce H2 and O2 is a long-standing goal in solar energy conversion, presenting a significant challenge without using sacrificial electron donors or external biases. Inspired by natural photosynthesis, the design of artificial Z-scheme photocatalytic systems is at the forefront of this field. These systems achieve higher redox potential by separating photogenerated electrons and holes through a fast interlayer recombination process between valence and conduction band edges. Z-scheme photocatalysis involves using two different semiconductors with distinct bandgap energies. Here, we explore potential systems based on two-dimensional (2D) heterostructures composed of carbon, nitrogen, or similar main group elements. The advantages and disadvantages of these systems are discussed, with a focus on enhancing their efficiency through strategic design. Special emphasis is placed on the dynamics of excited charge carrier transfer and recombination processes, which are crucial for developing efficient photocatalytic systems for overall water splitting.

Exploration of 4-substituted thiophene-based azo dyes for dye-sensitized solar cells and non-linear optical materials: synthesis and an in silico approach

This work examines the effects of changing the secondary donors' donating strengths at the thiophene ring's fourth position in a range of dyes designated SR1 to SR9. The DFT results indicate that the molecular planarity is greatly affected by the placement of the secondary donor at position four, which changes the charge transfer (CT) characteristics in the thiophene-azo-salicylic acid backbone. These results are corroborated by TD-DFT analysis, which indicates that as the secondary donor's donating strength increases, so does the vertical absorption maximum. Based on the computed photovoltaic characteristics, these dyes perform better from SR1 to SR5. Research on dye@TiO2 clusters suggests a potential for binding with TiO2, which might cause the dye@TiO2 clusters' absorbance to shift red. Additionally, computed linear and non-linear optical (NLO) properties exhibit similar trends to those observed for DSSC performance. The experimental results, which include HOMO and LUMO energies as well as initial absorption in dimethylformamide (DMF), match very well with the patterns seen in DFT and TD-DFT calculations. The high thermal stability of SR1 to SR4 is indicated by thermogravimetric analysis (TGA), indicating their practical applicability in non-linear optics (NLO) and DSSC applications. The efficiencies of the produced DSSCs vary; SR4 has the highest efficiency (4.50 ± 0.1), and SR1 has the lowest (0.37 ± 0.1). The combined effects of theoretical and experimental DSSC results demonstrate the vital role of secondary donors in influencing molecular characteristics and NLO and DSSC performance.

Energy-Conserving and Thermally Corrected Neglect of Back-Reaction Approximation Method for Nonadiabatic Molecular Dynamics

In this work, the energy-conserving and thermally corrected neglect of the back-reaction approximation approach for nonadiabatic molecular dynamics in extended atomistic systems is developed. The new approach introduces three key corrections to the original method: (1) it enforces the total energy conservation, (2) it introduces an explicit coupling of the system to its environment, and (3) it introduces a renormalization of nonadiabatic couplings to account for a difference between the instantaneous nuclear kinetic energy and the kinetic energy of guiding trajectories. In the new approach, an auxiliary kinetic energy variable is introduced as an independent dynamical variable. The new approach produces nonzero equilibrium populations, whereas the original neglect of the back-reaction approximation method does not. It yields population relaxation time scales that are favorably comparable to the reference values, and it introduces an explicit and controllable way of dissipating energy into a bath without an assumption of the bath being at equilibrium.

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.

A promising small-sized near-infrared absorbing zwitterionic dye for DSSC and NLO applications: DFT and TD-DFT approaches

Herein we investigate three quinoid zwitterionic dye sensitizers having donor-donor (4-dimethylaniline; ZIDM), donor-acceptor (4-dimethylaniline and 4-benzoic acid; ZIMCA), and acceptor-acceptor (4-benzoic acid; ZIDCA) that can be used in dye sensitized-solar cells and non-linear optical (NLO) application through density-functional theory (DFT) and time-dependent-DFT computations. ZIDM showed better charge transfer than ZIMCA and ZIDCA, which showed similar trends in chemical potential, electrophilicity index, hardness, and hyperhardness. The higher values of open circuit voltage, light harvesting efficiency, lower binding, and adsorption energy values for the dye to bind with the TiO2 cluster were observed for ZIDM. The results suggest that these dyes can easily hold with the TiO2 cluster through the monodentate binding mode possible between Ti and oxygen of the zwitterionic backbone. The examination of the linear and NLO properties of these dyes revealed that ZIDM has a higher α0 = 80.64 × 10-24 esu, β0 = 448.54 × 10-30 esu, and γ = 2219.23 × 10-36 esu in DCM. Similarly, higher values of molecular hyperpolarizability of 1335.0 × 10-48 esu and 8818.3 × 10-48 esu were observed in gas and DCM for ZIDM than ZIMCA and ZIDCA.

Enhancing photocatalytic overall water-splitting performance on dual-active-sites of the Co-P@MoS2 catalysts: a DFT study

The rational construction of photocatalysts possesses tremendous potential to solve the energy crisis and environmental pollution; however, designing a catalyst for solar-driven overall water-splitting remains a great challenge. Herein, we propose a new MoS2-based photocatalyst (Co-P@MoS2), which skillfully uses the cobalt (Co) atom to stimulate in-plane S atoms and employs the phosphorus (P) atom to stabilize the basal plane by forming the Co-P bands. Using density functional theory (DFT), it was found that oxygen evolution reaction (OER) and hydrogen evolution reaction (HER) can occur at the P site and S2 site of the Co-P@MoS2, respectively, and the dual-active sites successfully makes a delicate balance between the adsorption and dissociation of hydrogen. Furthermore, the improved overall water-splitting performance of Co-P@MoS2 was verified by analyzing the results of the electron structure and the dynamics of photogenerated carries. It was found that the imbalance of electron transfer caused by the introduction of the Co atom was the main contributor to the catalytic activity of Co-P@MoS2. Our study broadens the idea of developing photocatalysts for the overall water-splitting.

Characterization of Type I/II g-C3N4/MoS2 van der Waals Heterostructures: A New Theoretical Insight

The charge transfer mechanism of the g-C3N4/MoS2 heterojunction is still disputed. Some regard it as a type I pathway, some regard it as a type II pathway, and still some regard it as a Z-scheme pathway. Especially, the results obtained by density functional theory (DFT) calculations are not totally in agreement. Here, we constructed four g-C3N4/MoS2 heterojunctions on the basis of the aperture alignment modes of g-C3N4 and MoS2. Their morphology and photocatalytic activity were investigated via first-principles and excited state dynamics simulations. By systemically comparing the interfacial binding energy and electronic structure (e.g., band structure, electrostatic potential, and band edge positions) of g-C3N4/MoS2 heterojunctions, we found that both type I and type II band alignment structures could be obtained. Moreover, the calculated lifetimes of interlayer photogenerated electrons and holes show that type II g-C3N4/MoS2 tends to favor a general type II pathway rather than a Z-scheme pathway. This study could provide a deep understanding of the photocatalytic mechanism of g-C3N4/MoS2 van der Waals heterostructures, which will be of great use for applications in photocatalysis.

Oxide- and Silicate-Water Interfaces and Their Roles in Technology and the Environment

Interfacial reactions drive all elemental cycling on Earth and play pivotal roles in human activities such as agriculture, water purification, energy production and storage, environmental contaminant remediation, and nuclear waste repository management. The onset of the 21st century marked the beginning of a more detailed understanding of mineral aqueous interfaces enabled by advances in techniques that use tunable high-flux focused ultrafast laser and X-ray sources to provide near-atomic measurement resolution, as well as by nanofabrication approaches that enable transmission electron microscopy in a liquid cell. This leap into atomic- and nanometer-scale measurements has uncovered scale-dependent phenomena whose reaction thermodynamics, kinetics, and pathways deviate from previous observations made on larger systems. A second key advance is new experimental evidence for what scientists hypothesized but could not test previously, namely, interfacial chemical reactions are frequently driven by "anomalies" or "non-idealities" such as defects, nanoconfinement, and other nontypical chemical structures. Third, progress in computational chemistry has yielded new insights that allow a move beyond simple schematics, leading to a molecular model of these complex interfaces. In combination with surface-sensitive measurements, we have gained knowledge of the interfacial structure and dynamics, including the underlying solid surface and the immediately adjacent water and aqueous ions, enabling a better definition of what constitutes the oxide- and silicate-water interfaces. This critical review discusses how science progresses from understanding ideal solid-water interfaces to more realistic systems, focusing on accomplishments in the last 20 years and identifying challenges and future opportunities for the community to address. We anticipate that the next 20 years will focus on understanding and predicting dynamic transient and reactive structures over greater spatial and temporal ranges as well as systems of greater structural and chemical complexity. Closer collaborations of theoretical and experimental experts across disciplines will continue to be critical to achieving this great aspiration.

Titanium(IV) Surface Complexes Bearing Chelating Catecholato Ligands for Enhanced Band-Gap Reduction

Protonolysis reactions between dimethylamido titanium(IV) catecholate [Ti(CAT)(NMe₂)₂]₂ and neopentanol or tris(tert-butoxy)silanol gave catecholato-bridged dimers [(Ti(CAT)(OCH₂tBu)₂)(HNMe₂)]₂ and [Ti(CAT){OSi(OtBu)₃}₂(HNMe₂)₂]₂, respectively. Analogous reactions using the dimeric dimethylamido titanium(IV) (3,6-di-tert-butyl)catecholate [Ti(CATtBu₂-3,6)(NMe₂)₂]₂ yielded the monomeric Ti(CATtBu₂-3,6)(OCH₂tBu)₂(HNMe₂)₂ and Ti(CATtBu₂-3,6)[OSi(OtBu)₃]₂(HNMe₂)₂. The neopentoxide complex Ti(CATtBu₂-3,6)(OCH₂tBu)₂(HNMe₂)₂ engaged in further protonolysis reactions with Si-OH groups and was consequentially used for grafting onto mesoporous silica KIT-6. Upon immobilization, the surface complex [Ti(CATtBu₂-3,6)(OCH₂tBu)₂(HNMe₂)₂]@[KIT-6] retained the bidentate chelating geometry of the catecholato ligand. This convergent grafting strategy was compared with a sequential and an aqueous approach, which gave either a mixture of bidentate chelating species with a bipodally anchored Ti(IV) center along with other physisorbed surface species or not clearly identifiable surface species. Extension of the convergent and aqueous approaches to anatase mesoporous titania (m-TiO₂) enabled optical and electronic investigations of the corresponding surface species, revealing that the band-gap reduction is more pronounced for the bidentate chelating species (convergent approach) than for that obtained via the aqueous approach. The applied methods include X-ray photoelectron spectroscopy, ultraviolet photoelectron spectroscopy, and solid-state UV/vis spectroscopy. The energy-level alignment for the surface species from the aqueous approach, calculated from experimental data, accounts for the well-known type II excitation mechanism, whereas the findings indicate a distinct excitation mechanism for the bidentate chelating surface species of the material [Ti(CATtBu₂-3,6)(OCH₂tBu)₂(HNMe₂)₂]@[m-TiO₂].

Effect of regio-specific arylamine substitution on novel π-extended zinc salophen complexes: density functional and time-dependent density functional study on DSSC applications

A series of π-extended salophen-type Schiff-base zinc(ii) complexes, e.g., zinc-salophen complexes (ZSC), were investigated toward potential applications for dye-sensitized solar cells. The ZSC dyes adopt linear-, X-, or π-shaped geometries either with the functionalization of 1 donor/1 acceptor or 2 donors/2 acceptors to achieve a push-pull type molecular structure. The frontier molecular orbitals, light-harvesting properties as well as charge transfer characters against regio-specific substitution of donor/acceptor groups were studied by using density functional theory (DFT) and time-dependent density functional theory (TDDFT). The results reveal that all ZSC dyes of D-ZnS-π-A geometry (where D, S, and A denote to donor, salophen ligand, and acceptor, respectively) exhibit relatively lower HOMO energy compared to the structurally resembled porphyrin dye YD2-o-C8. Natural transition orbital (NTO) and electron-hole separation (EHS) approaches clearly differentiate the linear type YD-series dyes from CL-, AJ1-, and AJ2-series dyes because of poor charge transfer (CT) properties. In contrast, the π-shaped AJ2-series and X-shaped AJ1-series dyes outperform the others in a manner of stronger CT characteristics, broadened UV-vis absorption as well as tunable bandgap simply via substitution of p-ethynylbenzoic acids (EBAs) and arylamine donors at salophen 7,8- and 2,3,12,13-positions, respectively. Both EHS and calculated exciton binding energies suggest the strength of CT character for ZSC dyes with an amino donor in the trend TPA > AN > DPA. This work has provided clear illustration toward molecular design of efficient dyes featuring a zinc-salophen backbone.

Dynamics of Electron Transfers in Photosensitization Reactions of Zinc Porphyrin Derivatives

Photocatalytic systems for CO₂ reduction operate via complicated multi-electron transfer (ET) processes. A complete understanding of these ET dynamics can be challenging but is key to improving the efficiency of CO₂ conversion. Here, we report the ET dynamics of a series of zinc porphyrin derivatives (ZnPs) in the photosensitization reactions where sequential ET reactions of ZnPs occur with a sacrificial electron donor (SED) and then with TiO₂. We employed picosecond time-resolved fluorescence spectroscopy and femtosecond transient absorption (TA) measurement to investigate the fast ET dynamics concealed in the steady-state or slow time-resolved measurements. As a result, Stern-Volmer analysis of fluorescence lifetimes evidenced that the reaction of photoexcited ZnPs with SED involves static and dynamic quenching. The global fits to the TA spectra identified much faster ET dynamics on a few nanosecond-time scales in the reactions of one-electron reduced species (ZnPs^•-) with TiO₂ compared to previously measured minute-scale quenching dynamics and even diffusion rates. We propose that these dynamics report the ET dynamics of ZnPs^•- formed at adjacent TiO₂ without involving diffusion. This study highlights the importance of ultrafast time-resolved spectroscopy for elucidating the detailed ET dynamics in photosensitization reactions.

A quantitative model of charge injection by ruthenium chromophores connecting femtosecond to continuous irradiance conditions

A kinetic framework for the ultrafast photophysics of tris(2,2-bipyridine)ruthenium(II) phosphonated and methyl-phosphonated derivatives is used as a basis for modeling charge injection by ruthenium dyes into a semiconductor substrate. By including the effects of light scattering, dye diffusion, and adsorption kinetics during sample preparation and the optical response of oxidized dyes, quantitative agreement with multiple transient absorption datasets is achieved on timescales spanning femtoseconds to nanoseconds. In particular, quantitative agreement with important spectroscopic handles-the decay of an excited state absorption signal component associated with charge injection in the UV region of the spectrum and the dynamical redshift of a ∼500 nm isosbestic point-validates our kinetic model. Pseudo-first-order rate coefficients for charge injection are estimated in this work, with an order of magnitude ranging from 10¹¹ to 10¹² s^-1. The model makes the minimalist assumption that all excited states of a particular dye have the same charge injection coefficient, an assumption that would benefit from additional theoretical and experimental exploration. We have adapted this kinetic model to predict charge injection under continuous solar irradiation and find that as many as 68 electron transfer events per dye per second take place, significantly more than prior estimates in the literature.

How Good Is the Vibronic Hamiltonian Repetition Approach for Long-Time Nonadiabatic Molecular Dynamics?

Multiple applied studies of slow nonadiabatic processes in nanoscale and condensed matter systems have adopted the "repetition" approximation in which long trajectories for such simulations are obtained by concatenating shorter trajectories, directly available from ab initio calculations, many times. Here, we comprehensively assess this approximation using model Hamiltonians with parameters covering a wide range of regimes. We find that state transition time scales may strongly depend on the length of the repeated data, although the convergence is not monotonic and may be slow. The repetition approach may under- or overestimate the time scales by a factor of ≤7-8, does not directly depend on the dispersion of energy gap and nonadiabatic coupling (NAC) frequencies, but may depend on the magnitude of the NACs. We suggest that the repetition-based nonadiabatic dynamics may be inaccurate in simulations with very small NACs, where intrinsic transition times are on the order of ≥100 ps.

Synthesis of innovative triphenylamine-functionalized organic photosensitizers outperformed the benchmark dye N719 for high-efficiency dye-sensitized solar cells

Herein, we present a thorough photovoltaic investigation of four triphenylamine organic sensitizers with D-π-A configurations and compare their photovoltaic performances to the conventional ruthenium-based sensitizer N719. SFA-5-8 are synthesized and utilized as sensitizers for dye-sensitized solar cell (DSSC) applications. The effects of the donor unit (triphenylamine), π-conjugation bridge (thiophene ring), and various acceptors (phenylacetonitrile and 2-cyanoacetamide derivatives) were investigated. Moreover, this was asserted by profound calculations of HOMO (highest occupied molecular orbital) and LUMO (lowest unoccupied molecular orbital) energy levels, the molecular electrostatic potential (MEP), and natural bond orbital (NBO) that had been studied for the TPA-sensitizers. Theoretical density functional theory (DFT) was performed to study the distribution of electron density between donor and acceptor moieties. The sensitization by the absorption of sensitizers SFA-5-8 leads to an obvious enhancement in the visible light absorption (300-750 nm) as well as a higher photovoltaic efficiency in the range of (5.53-7.56%). Under optimized conditions, SFA-7 showed outstanding sensitization of nanocrystalline TiO2, resulting in enhancing the visible light absorption and upgrading the power conversion efficiency (PCE) to approximately 7.56% over that reported for the N719 (7.29%). Remarkably, SFA-7 outperformed N719 by 4% in the total conversion efficiency. Significantly, the superior performance of SFA-7 could be mainly ascribed to the higher short-circuit photocurrents (Jsc) in parallel with larger open-circuit voltages (Voc) and more importantly, the presence of different anchoring moieties that could enhance the ability to fill the gaps on the surface of the TiO2 semiconductor. That could be largely reflected in the overall enhancement in the device efficiency. Moreover, the theoretical electronic and photovoltaic properties of all studied sensitizers have been compared with experimental results. All the 2-cyanoacrylamide derivative sensitizers demonstrated robust photovoltaic performance.

Ultrafast Dynamics of Photoinduced Electron Transfer in Bay-Aryl-Substituted Perylene Diimide Derivatives

Blends of donors and acceptors have been widely used in bulk-heterojunction solar cells to have exciton formation and charge separation by photoinduced electron transfer (PET). In this work, we have synthesized perylene diimide (PDI)-based materials having different aryl substituents at the bay positions (4-Anisyl-PDI, CBZ-N-Ph-PDI, and 4-Pyridyl-PDI) to understand the excited-state dynamics of electron transfer. The detailed photophysics was studied using steady-state as well as ultrafast dynamics of the excited states in different solvents. CBZ-N-Ph-PDI showed tremendous effects of the solvent on the electronic properties compared with the other two derivatives. The emission quantum yield of CBZ-N-Ph-PDI decreases drastically in dichloromethane and other polar solvents, indicating strong electron transfer. DFT calculations showed that in CBZ-N-Ph-PDI the HOMO is centered mostly on the N-phenylcarbazole and the LUMO is on the electron-poor PDI moieties. In addition, the energy levels of the HOMO and HOMO-1 in CBZ-N-Ph-PDI are estimated to be identical. The free energy change for charge separation (ΔG_CS) was calculated using electrochemical and photophysical data and found to be negative for CBZ-N-Ph-PDI. The ground- and excited-state dipole moment ratios suggest that the excited state of 4-Pyridyl-PDI (1.90) is less polar than that of 4-Anisyl-PDI (3.67), which provides an idea of the lower possibility of charge separation in 4-Anisyl-PDI and 4-Pyridyl-PDI. Ultrafast photodynamics studies of 4-Anisyl-PDI, CBZ-N-Ph-PDI, and 4-Pyridyl-PDI showed fast electron transfer only in CBZ-N-Ph-PDI and not in the other PDI derivatives. It was also observed that electron transfer is faster in DCM and THF than in toluene. Ultrafast dynamics studies showed the presence of an equilibrium between electron transfer and decay from the singlet excited state. Ultrafast studies also showed the features of the N-phenylcarbazole cation and PDI anion, further confirming the intramolecular electron transfer in CBZ-N-Ph-PDI.

Electron Transfer at Quantum Dot–Metal Oxide Interfaces for Solar Energy Conversion

Electron transfer at a donor–acceptor quantum dot–metal oxide interface is a process fundamentally relevant to solar energy conversion architectures as, e.g., sensitized solar cells and solar fuels schemes. As kinetic competition at these technologically relevant interfaces largely determines device performance, this Review surveys several aspects linking electron transfer dynamics and device efficiency; this correlation is done for systems aiming for efficiencies up to and above the ∼33% efficiency limit set by Shockley and Queisser for single gap devices. Furthermore, we critically comment on common pitfalls associated with the interpretation of kinetic data obtained from current methodologies and experimental approaches, and finally, we highlight works that, to our judgment, have contributed to a better understanding of the fundamentals governing electron transfer at quantum dot–metal oxide interfaces.

Promoting Water Activation by Photogenerated Holes in Monolayer C2N

In photocatalytic reactions, the activation of H₂O is very important for achieving high energy conversion efficiency. However, its activation mechanism under photoirradiation is still not fully understood. Here, on the basis of first-principles calculations, the role of photogenerated holes on the activation of H₂O is investigated in a typical photocatalytic material C₂N. The H₂O molecule adsorbs at the six-membered N pore of C₂N with a dual H-bonding configuration. Due to the electrostatic repulsion between the O atom of H₂O and six N atoms of C₂N, the energy level of the H₂O molecule's highest occupied molecular orbital is raised significantly to exceed the valence band maximum of C₂N, so that the photogenerated holes in C₂N can be quickly captured by the H₂O molecule. The captured photogenerated holes boost the activation of H₂O and reduce the dissociation energy barrier from 1.61 to 0.69 eV. Besides, p-type defects of C₂N have similar effects as photogenerated holes.

Ultrafast Interlayer Charge Separation, Enhanced Visible-Light Absorption, and Tunable Overpotential in Twisted Graphitic Carbon Nitride Bilayers for Water Splitting

Moiré pattern superlattice formed by 2D van der Waals layered structures have attracted great attention for diverse applications. In experiments, the enhancement of catalytic performance in twisted bilayer systems is reported while its mechanism remains unclear. From high-accuracy first-principles and time-dependent ab initio nonadiabatic molecular dynamics calculations, ultrafast interlayer charge transfer within 120 fs, excellent charge separation, improved visible-light absorption, and satisfactory overpotentials for the hydrogen evolution and oxygen evolution reactions in twisted graphitic carbon nitride (g-C₃ N₄ ) bilayers are found, which are beneficial to photocatalytic, photo-electrocatalytic, or electrocatalytic water splitting. This work provides insightful guidance to advanced nanocatalysis based on twisted layered materials.

Real Space-Real Time Evolution of Excitonic States Based on the Bethe-Salpeter Equation Method

We introduce a method for constructing localized excitations and simulating the real time dynamics of excitons at the Many-Body Perturbation Theory Bethe-Salpeter Equation level. We track, on the femto-seconds scale, electron injection from a photoexcited dye into a semiconducting slab. From the time-dependent many-body wave function we compute the spatial evolution of the electron and of the hole; full electron injection is attained within 5 fs. Time-resolved analysis of the electron density and electron-hole interaction energy hints at a two-step charge transfer mechanism through an intermediary partially injected state. We adopt the Von-Neumann entropy for analyzing how the electron and hole entangle. We find that the excitation of the dye-semiconductor model may be represented by a four-level system and register a decrease in entanglement upon electron injection. At full injection, the electron and the hole exhibit only a small degree of entanglement indicative of pure electron and hole states.

Structure and Energetics of Dye-Sensitized NiO Interfaces in Water from Ab Initio MD and Large-Scale GW Calculations

The energy-level alignment across solvated molecule/semiconductor interfaces is a crucial property for the correct functioning of dye-sensitized photoelectrodes, where, following the absorption of solar light, a cascade of interfacial hole/electron transfer processes has to efficiently take place. In light of the difficulty of performing X-ray photoelectron spectroscopy measurements at the molecule/solvent/metal-oxide interface, being able to accurately predict the level alignment by first-principles calculations on realistic structural models would represent an important step toward the optimization of the device. In this respect, dye/NiO surfaces, employed in p-type dye-sensitized solar cells, are undoubtedly challenging for ab initio methods and, also for this reason, much less investigated than the n-type dye/TiO₂ counterpart. Here, we consider the C343-sensitized NiO surface in water and combine ab initio molecular dynamics (AIMD) simulations with GW (G₀W₀) calculations, performed along the MD trajectory to reliably describe the structure and energetics of the interface when explicit solvation and finite temperature effects are accounted for. We show that the differential perturbative correction on the NiO and molecule states obtained at the GW level is mandatory to recover the correct (physical) interfacial energetics, allowing hole transfer from the semiconductor valence band to the highest occupied molecular orbital (HOMO) of the dye. Moreover, the calculated average driving force quantitatively agrees with the experimental estimate.

Description of Sudden Polarization in the Excited Electronic States with an Ensemble Density Functional Theory Method

Sudden polarization (SP) is one of the manifestations of electron transfer in the electronically excited states of molecules. Proposed initially to explain the unusual reactivity of photoexcited olefins, SP often occurs in the excited states of molecules possessing strongly correlated diradical ground state. Theoretical description of SP involves mixing between the singly excited and the doubly excited zwitterionic states, which makes it inaccessible with the use of the popular linear-response time-dependent density functional theory methods. In this work, an extended variant of the state-interaction state-averaged spin-restricted ensemble-referenced Kohn-Sham (SI-SA-REKS, or SSR) method is applied to study SP in a number of organic diradical systems. To this end, the analytical derivative formalism is derived and implemented for the SSR(3,2) method (see the main text for explanation of the acronym), which enables the automatic geometry optimization and obtains the relaxed density matrices as well as the electron binding energies and respective Dyson's orbitals. Application of the new method to SP in the lowest singlet excited state of ethylene agrees with the results obtained previously with the use of multireference methods of wavefunction theory. A number of interesting manifestations of SP are observed, such as the charge transfer in photoexcited tetramethyleneethene (TME) diradical mediated by the vibrational motion and conductivity switching in the excited state of a donor-acceptor dyad placed in an external electric field.

Real-time GW-BSE investigations on spin-valley exciton dynamics in monolayer transition metal dichalcogenide

We develop an ab initio nonadiabatic molecular dynamics (NAMD) method based on GW plus real-time Bethe-Salpeter equation (GW + rtBSE-NAMD) for the spin-resolved exciton dynamics. From investigations on MoS₂, we provide a comprehensive picture of spin-valley exciton dynamics where the electron-phonon (e-ph) scattering, spin-orbit interaction (SOI), and electron-hole (e-h) interactions come into play collectively. In particular, we provide a direct evidence that e-h exchange interaction plays a dominant role in the fast valley depolarization within a few picoseconds, which is in excellent agreement with experiments. Moreover, there are bright-to-dark exciton transitions induced by e-ph scattering and SOI. Our study proves that e-h many-body effects are essential to understand the spin-valley exciton dynamics in transition metal dichalcogenides and the newly developed GW + rtBSE-NAMD method provides a powerful tool for exciton dynamics in extended systems with time, space, momentum, energy, and spin resolution.

Effect of Ligand Functionalization on the Rate of Charge Carrier Recombination in Metal-Organic Frameworks: A Case Study of MIL-125

Ligand functionalization is a powerful approach for modifying the electronic structure of metal-organic frameworks when targeting the optimal electronic properties for photocatalysis and photovoltaics. However, its effect on the charge carrier lifetimes and recombination pathways remains unexplored. In this work, first-principles simulations, including nonadiabatic molecular dynamics, are performed for the representative TiO₂-based metal-organic framework systems MIL-125-X to unravel the impact of ligand functionalization on the nonradiative electron-hole recombination process, decoherence rates, and phonon modes giving the largest contribution to the nonradiative decay. Nonradiative recombination rates, simulated using the PBE0 density functional, are in excellent agreement with experiment. The ligand functionalization in MIL-125-X influences the recombination rates, unraveling the trend opposite to the evolution of the band gap and affecting the nonadiabatic coupling coefficients. Ligand modification impacts the phonon modes, which contribute most to the recombination process, altering the distribution between soft phonon modes and vibrational modes associated with specific structural motifs.

Symmetry Effects in Photoinduced Electron Transfer in Chlorin‐Quinone Dyads: Adiabatic Suppression in the Marcus Inverted Region

In donor–acceptor dyads undergoing photoinduced electron transfer (PET), a direction or pathway for electron movement is usually dictated by the redox properties and the separation distance between the donor and acceptor subunits, while the effect of symmetry is less recognized. We have designed and synthesized two isomeric donor–acceptor assemblies in which electronic coupling between donor and acceptor is altered by the orbital symmetry control with the reorganization energy and charge transfer exothermicity being kept unchanged. Analysis of the optical absorption and luminescence spectra, supported by the DFT and TD‐DFT calculations, showed that PET in these assemblies corresponds to the Marcus inverted region (MIR) and has larger rate for isomer with weaker electronic coupling. This surprising observation provides the first experimental evidence for theoretically predicted adiabatic suppression of PET in MIR, which unambiguously controlled solely by symmetry.

Accurate Computation of Nonadiabatic Coupling with Projector Augmented-Wave Pseudopotentials

Synergy of nonadiabatic molecular dynamics with real-time time-dependent density functional theory has led to significant progress in modeling excited-state dynamics in nanoscale and condensed matter systems over the past decade. Nonadiabatic coupling (NAC) is the central quantity in such simulations, and its accurate and efficient evaluation is an enduring challenge in time-dependent Kohn-Sham theory, particularly in conjunction with planewave basis sets and projector augmented-wave (PAW) pseudopotentials because of the complexity of the PAW "all-electron" wave function. We report a method for rigorous evaluation of the NAC with PAW wave functions and demonstrate an efficient approximation to the rigorous NAC that gives comparable accuracy. As a validation, we intensely examine the NAC matrix elements calculated using both pseudo- and all-electron wave functions under the PAW formalism in six representative systems. The approximate NAC obtained with pseudowave functions is close to the exact all-electron NAC, with the largest deviations observed when subshell d-electrons are involved in the transitions. The developed approach provides a rigorous and convenient methodology for the numerical computation of NAC in the Kohn-Sham theory framework.

Interlayer coupling prolonged the photogenerated carrier lifetime of few layered Bi2OS2 semiconductors

Layered semiconductors with broad photoabsorption, a long carrier lifetime and high carrier mobility are of crucial importance for high-performance optoelectronic and photovoltaic devices; however it is hard to satisfy these requirements simultaneously in a system due to the opposite dependence on the layer thickness. Herein, by means of ab initio time-domain nonadiabatic molecular dynamic simulations, we find a new mechanism in Bi2OS2 nanosheets inducing an anomalous layer-dependent property of carrier lifetimes, which makes the few layered Bi2OS2 a possible system for fulfilling the above requirements concurrently. It is revealed that the interlayer dipole-dipole interaction in few layered Bi2OS2 effectively breaks the two-fold degenerate orbitals of [BiS2] layers, which not only cuts down the overlap of the electron and hole wave functions, but also accelerates the electron decoherence process. This significantly suppresses the electron-hole recombination and prolongs the photogenerated carrier lifetime of few layered Bi2OS2. The mechanism unveiled here paves a possible way for developing advanced optoelectronic and photovoltaic devices through engineering interlayer dipole-dipole coupling.

Non-adiabatic Excited-State Molecular Dynamics: Theory and Applications for Modeling Photophysics in Extended Molecular Materials

Optically active molecular materials, such as organic conjugated polymers and biological systems, are characterized by strong coupling between electronic and vibrational degrees of freedom. Typically, simulations must go beyond the Born-Oppenheimer approximation to account for non-adiabatic coupling between excited states. Indeed, non-adiabatic dynamics is commonly associated with exciton dynamics and photophysics involving charge and energy transfer, as well as exciton dissociation and charge recombination. Understanding the photoinduced dynamics in such materials is vital to providing an accurate description of exciton formation, evolution, and decay. This interdisciplinary field has matured significantly over the past decades. Formulation of new theoretical frameworks, development of more efficient and accurate computational algorithms, and evolution of high-performance computer hardware has extended these simulations to very large molecular systems with hundreds of atoms, including numerous studies of organic semiconductors and biomolecules. In this Review, we will describe recent theoretical advances including treatment of electronic decoherence in surface-hopping methods, the role of solvent effects, trivial unavoided crossings, analysis of data based on transition densities, and efficient computational implementations of these numerical methods. We also emphasize newly developed semiclassical approaches, based on the Gaussian approximation, which retain phase and width information to account for significant decoherence and interference effects while maintaining the high efficiency of surface-hopping approaches. The above developments have been employed to successfully describe photophysics in a variety of molecular materials.

Fundamentals of TiO2 Photocatalysis: Concepts, Mechanisms, and Challenges

Photocatalysis has been widely applied in various areas, such as solar cells, water splitting, and pollutant degradation. Therefore, the photochemical mechanisms and basic principles of photocatalysis, especially TiO₂ photocatalysis, have been extensively investigated by various surface science methods in the last decade, aiming to provide important information for TiO₂ photocatalysis under real environmental conditions. Recent progress that provides fundamental insights into TiO₂ photocatalysis at a molecular level is highlighted. Insights into the structures of TiO₂ and the basic principles of TiO₂ photocatalysis are discussed first, which provides the basic concepts of TiO₂ photocatalysis. Following this, details of the photochemistry of three important molecules (oxygen, water, methanol) on the model TiO₂ surfaces are presented, in an attempt to unravel the relationship between charge/energy transfer and bond breaking/forming in TiO₂ photocatalysis. Lastly, challenges and opportunities of the mechanistic studies of TiO₂ photocatalysis at the molecular level are discussed briefly, as well as possible photocatalysis models.

Suppression of Electron-Hole Recombination by Intrinsic Defects in 2D Monoelemental Material

The Shockley-Read-Hall (SRH) model, in which the deep trap defect states in the band gap are proposed as nonradiative electron-hole (e-h) recombination centers, has been widely used to describe the nonradiative e-h recombination through the defects in semiconductor. By using the ab initio nonadiabatic molecular dynamics method, we find that the SRH model fails to describe the e-h recombination behavior for defects in 2D monoelemental material such as monolayer black phosphorus (BP). Through the investigation of three intrinsic defects with shallow and deep defect states in monolayer BP, it is found that, surprisingly, none of these defects significantly accelerates the e-h recombination. Further analysis shows that because monolayer BP is a monoelemental material, the distinct impurity phonon, which often induces fast e-h recombination, is not formed. Moreover, because of the flexibility of 2D material, the defects scatter the phonons present in pristine BP, generating multiple modes with lower frequencies compared with the pristine BP, which further suppresses the e-h recombination. We propose that the conclusion can be extended to other monoelemental 2D materials, which is important guidance for the future design of functional semiconductors.