Theoretical Insights into Interfacial Electron Transfer between Zinc Phthalocyanine and Molybdenum Disulfide

Ab Initio Modeling of Mixed-Dimensional Heterostructures: A Path Forward

Understanding the electronic structure of mixed-dimensional heterostructures is essential for maximizing their application potential. However, accurately modeling such interfaces is challenging due to the complex interplay between the subsystems. We employ a computational framework integrating first-principles methods, including GW, density functional theory (DFT), and the polarizable continuum model, to elucidate the electronic structure of mixed-dimensional heterojunctions formed by free-base phthalocyanines and monolayer molybdenum disulfide. We assess the impact of dielectric screening across various scenarios, from isolated molecules to organic films on a substrate-supported monolayer. Our findings show that while polarization effects cause significant renormalization of molecular energy levels, band energies and alignments in the most relevant setup can be accurately predicted through DFT simulations of the individual subsystems. Additionally, we analyze orbital hybridization, revealing potential pathways for interfacial charge transfer. This study offers new insights into hybrid inorganic/organic interfaces and provides a practical computational protocol suitable for scaled-up studies.

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.

Optical Enhancement of Indirect Bandgap 2D Transition Metal Dichalcogenides for Multi-Functional Optoelectronic Sensors

The unique electrical and optical properties of transition metal dichalcogenides (TMDs) make them attractive nanomaterials for optoelectronic applications, especially optical sensors. However, the optical characteristics of these materials are dependent on the number of layers. Monolayer TMDs have a direct bandgap that provides higher photoresponsivity compared to multilayer TMDs with an indirect bandgap. Nevertheless, multilayer TMDs are more appropriate for various photodetection applications due to their high carrier density, broad spectral response from UV to near-infrared, and ease of large-scale synthesis. Therefore, this review focuses on the modification of the optical properties of devices based on indirect bandgap TMDs and their emerging applications. Several successful developments in optical devices are examined, including band structure engineering, device structure optimization, and heterostructures. Furthermore, it introduces cutting-edge techniques and future directions for optoelectronic devices based on multilayer TMDs.

Surface hopping modeling of charge and energy transfer in active environments

An active environment is any atomic or molecular system changing a chromophore's nonadiabatic dynamics compared to the isolated molecule. The action of the environment on the chromophore occurs by changing the potential energy landscape and triggering new energy and charge flows unavailable in the vacuum. Surface hopping is a mixed quantum-classical approach whose extreme flexibility has made it the primary platform for implementing novel methodologies to investigate the nonadiabatic dynamics of a chromophore in active environments. This Perspective paper surveys the latest developments in the field, focusing on charge and energy transfer processes.

Enhanced Electron Transfer and Spin Flip through Spin-Orbital Couplings in Organic/Inorganic Heterojunctions: A Nonadiabatic Surface Hopping Simulation

The circumstances of transferred electrons across organic/inorganic interfaces have attracted intensive interest because of the distinctive electronic structure properties of those two components. Leveraging ab initio nonadiabatic molecular dynamics methods in conjunction with spin dynamics induced by spin-orbital couplings (SOCs), this study reports two competitive channels during photoinduced dynamical processes in the prototypical ZnPc/monolayer MoS₂ heterojunction. Interestingly, the electron-transfer and relaxation processes occur simultaneously because of the enhancement of electron-phonon couplings and expansion of dynamical pathways by SOCs, suggesting that the electron-transfer rate and relaxation processes can be tuned by SOCs, hence yielding the performance promotion of photovoltaic and photocatalytic devices. Additionally, approximately half of the transferred electrons flip their spin within 1.6 ps because of strong SOCs in MoS₂, achieving great agreement with experimental measurements. This investigation provides instructive perspectives for designing novel devices and applications based on organic/inorganic heterojunctions, demonstrating the importance of spin dynamics simulations in exploring sophisticated photoinduced processes in materials.

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.

Quasiparticle electronic structure of phthalocyanine:TMD interfaces from first-principles GW

Interfaces formed between monolayer transition metal dichalcogenides and (metallo)phthalocyanine molecules are promising in energy applications and provide a platform for studying mixed-dimensional molecule-semiconductor heterostructures in general. An accurate characterization of the frontier energy level alignment at these interfaces is key in the fundamental understanding of the charge transfer dynamics between the two photon absorbers. Here, we employ the first-principles substrate screening GW approach to quantitatively characterize the quasiparticle electronic structure of a series of interfaces: metal-free phthalocyanine (H2Pc) adsorbed on monolayer MX2 (M = Mo, W; X = S, Se) and zinc phthalocyanine (ZnPc) adsorbed on MoX2 (X = S, Se). Furthermore, we reveal the dielectric screening effect of the commonly used α-quartz (SiO2) substrate on the H2Pc:MoS2 interface using the dielectric embedding GW approach. Our calculations furnish a systematic set of GW results for these interfaces, providing the structure-property relationship across a series of similar systems and benchmarks for future experimental and theoretical studies.

Interfacial photoinduced carrier dynamics tuned by polymerization of coronene molecules encapsulated in carbon nanotubes: bridging type-I and type-II heterojunctions

Carbon nanomaterials play important roles in modern scientific research. Integrating different carbon-based building blocks into nano-hybrid architectures not only takes full advantage of each component, but also brings in novel interfacial properties. Herein, we have employed density functional theory (DFT) calculations to investigate the effects of polymerization degree of coronene molecules encapsulated in single-walled carbon nanotubes (SWNTs) (19,0) on their interfacial properties. The present results reveal that the interfacial properties of the formed heterojunctions are remarkably regulated by the polymerization degree. For example, monomer- and dimer-encapsulated SWNTs are type-I heterojunctions in which interfacial excitation energy transfer is preferred, whereas interfacial charge carrier transfer is favorable in trimer- and polymer-encapsulated SWNTs because they are type-II heterojunctions. On the other hand, we have employed the time-domain nonadiabatic dynamics simulation approach to explore the interfacial carrier dynamics in type-II polymer-encapsulated SWNT heterojunctions. It is found that the electron and hole transfer processes are asymmetric and occur in opposite directions and at different rates. The former takes place from polymers to SWNTs in an ultrafast way (ca. 370 fs), whereas the latter occurs slowly from SWNTs to polymers (ca. 24 ps). A closer analysis uncovers the fact that the different carrier transfer rates mainly originate from the different densities of the acceptor states, energy differences and inter-state couplings between the donor and acceptor states. Finally, the present work demonstrates that the polymerization degree could act as a new regulating strategy to tune the interfacial properties of molecule-encapsulated SWNT heterojunctions.

Rational design of monolayer transition metal dichalcogenide@fullerene van der Waals photovoltaic heterojunctions with time-domain density functional theory simulations

van der Waals heterojunctions formed by transition metal dichalcogenides (TMDs) and fullerenes are promising candidates for novel photovoltaic devices due to the excellent optoelectronic properties of both TMDs and fullerenes. However, relevant experimental and theoretical investigations remain scarce to the best of our knowledge. Herein, we have first employed static density functional theory (DFT) calculations in combination with time-domain density functional theory (TDDFT) based nonadiabatic dynamics simulations to rationally evaluate the photovoltaic performances of four TMD@fullerene heterostructures, i.e. WSe₂@C₆₀, WSe₂@C₇₀, MoTe₂@C₆₀ and MoTe₂@C₇₀, respectively. Our simulation results indicate that the C₇₀-based heterostructures overall have better photoinduced electron transfer efficiencies than their C₆₀-based counterparts, among which the performance of the WSe₂@C₇₀ heterostructure is the best and the electron transfer from WSe₂ to C₇₀ almost accomplishes within 1 ps. In addition, the large build-in potential of about 0.75 eV of WSe₂@C₇₀ is beneficial for the charge separation processes. Our present work not only selects the van der Waals TMD@fullerene heterojunctions that might have excellent photovoltaic properties, but also paves the way for the rational design of novel heterojunctions with better optoelectronic performances with DFT and TDDFT simulations in the future.

Photoinduced electron transfer from carbon nanotubes to fullerenes: C60versus C70

Hybrid carbon materials are found to exhibit novel optoelectronic properties at their interfaces, but the related interfacial carrier dynamics is rarely explored theoretically. In this contribution, we have employed density functional theory (DFT) and DFT-based nonadiabatic dynamics methods to explore photoinduced interfacial electron transfer processes at interfaces between a single-walled carbon nanotube with chiral index (6,5) and C60 or C70 (C60@CNT65 and C70@CNT65). We have found that with low E11 excitation, electron transfer takes place from CNT65 to C60 and C70 in both heterojunctions. This process is ultrafast and completed within about 200 fs, which is consistent with recent experiments. Differently, high E22 excitation does not induce electron injection to C60 in C60@CNT65; instead, "hot" electrons produced within CNT65 will be trapped in its higher conduction band for a while because of slow inter-band relaxation. By contrast, in C70@CNT65, high E22 excitation still can lead to ultrafast electron transfer to C70, but only a comparable amount of electrons are transferred (ca. 30%). Interestingly, electrons either remaining on CNT65 or transferred to C70 are trapped in the higher conduction band for a while, similarly, due to slow inter-band relaxation. The present results could be useful to guide the design of excellent interfaces of mixed-dimensional hybrid carbon materials for various optoelectronic applications.

Emergent Optoelectronic Properties of Mixed-Dimensional Heterojunctions

ConspectusThe electronic dimensionality of a material is defined by the number of spatial degrees of confinement of its electronic wave function. Low-dimensional semiconductor nanomaterials with at least one degree of spatial confinement have optoelectronic properties that are tunable with size and environment (dielectric and chemical) and are of particular interest for optoelectronic applications such as light detection, light harvesting, and photocatalysis. By combining nanomaterials of differing dimensionalities, mixed-dimensional heterojunctions (MDHJs) exploit the desirable characteristics of their components. For example, the strong optical absorption of zero-dimensional (0D) materials combined with the high charge carrier mobilities of two-dimensional (2D) materials widens the spectral response and enhances the responsivity of mixed-dimensional photodetectors, which has implications for ultrathin, flexible optoelectronic devices. MDHJs are highly sensitive to (i) interfacial chemistry because of large surface area-to-volume ratios and (ii) electric fields, which are incompletely screened because of the ultrathin nature of MDHJs. This sensitivity presents opportunities for control of physical phenomena in MDHJs through chemical modification, optical excitation, externally applied electric fields, and other environmental parameters. Since this fast-moving research area is beginning to pose and answer fundamental questions that underlie the fundamental optoelectronic behavior of MDHJs, it is an opportune time to assess progress and suggest future directions in this field.In this Account, we first outline the characteristic properties, advantages, and challenges for low-dimensional materials, many of which arise as a result of quantum confinement effects. The optoelectronic properties and performance of MDHJs are primarily determined by dynamics of excitons and charge carriers at their interfaces, where these particles tunnel, trap, scatter, and/or recombine on the time scales of tens of femtoseconds to hundreds of nanoseconds. We discuss several photophysical phenomena that deviate from those observed in bulk heterojunctions, as well as factors that can be used to vary, probe, and ultimately control the behavior of excitons and charge carriers in MDHJ systems. We then discuss optoelectronic applications of MDHJs, namely, photodetectors, photovoltaics, and photocatalysts, and identify current performance limits compared to state-of-the-art benchmarks. Finally, we suggest strategies to extend the current understanding of dynamics in MDHJs toward the realization of stimuli-driven responses, particularly with respect to exciton delocalization, quantum emission, interfacial morphology, responsivity to external stimuli, spin selectivity, and usage of chemically reactive materials.

Probing Hydrogen-Bonding Properties of a Negatively Charged MoS2 Monolayer by Powder X-ray Diffraction and Density Functional Theory Calculations

The contributions of various noncovalent interactions in stabilization of the assembled and delaminated MoS₂-hexamethylenetetramine (HMTA)-layered compound resulted from the assembly of protonated HMTA molecules and negatively charged 1T-MoS₂ monolayers have been considered on the basis of powder X-ray diffraction pattern modeling, density functional theory calculations, and atoms in molecules quantum theory analysis. The structure with HMTA cations involved in NH···S bonding with MoS₂ layers was concluded to be more advantageous than the alternative one with NH···N bonding between the cations. Delamination was demonstrated to essentially influence the hierarchy of interactions and leads to significant strengthening of the NH···S hydrogen bond established between HMTA and the MoS₂ monolayer surface. The method applied in this study for evaluation of the monolayer MoS₂ properties on the basis of the 3D structure of the MoS₂-organic compound is expected to be helpful to gain insights into the interactions occurring in many MoS₂-based systems.

Modeling nonadiabatic dynamics in condensed matter materials: some recent advances and applications

This review focuses on recent developments in the field of nonadiabatic molecular dynamics (NA-MD), with particular attention given to condensed-matter systems. NA-MD simulations for small molecular systems can be performed using high-level electronic structure (ES) calculations, methods accounting for the quantization of nuclear motion, and using fewer approximations in the dynamical methodology itself. Modeling condensed-matter systems imposes many limitations on various aspects of NA-MD computations, requiring approximations at various levels of theory-from the ES, to the ways in which the coupling of electrons and nuclei are accounted for. Nonetheless, the approximate treatment of NA-MD in condensed-phase materials has gained a spin lately in many applied studies. A number of advancements of the methodology and computational tools have been undertaken, including general-purpose methods, as well as those tailored to nanoscale and condensed matter systems. This review summarizes such methodological and software developments, puts them into the broader context of existing approaches, and highlights some of the challenges that remain to be solved.

Photoinduced Carrier Dynamics at the Interface of Black Phosphorus and Bismuth Vanadate

Two-dimensional (2D) heterostructures of black phosphorus (BP)/bismuth vanadate (BVO) have attracted much attention due to their potential uses in photocatalytic water splitting. However, the interfacial photoinduced electron- and hole-transfer dynamics are not explored computationally. Herein, we have used density functional theory (DFT) calculations and DFT-based fewest-switches surface-hopping dynamics simulations to investigate the light-driven electron and hole dynamics taking place at the interface of BP and the BVO(010) surface. Our results show that the BP monolayer is adsorbed on BVO(010) via van der Waals interaction. Upon irradiation, the electron transfer takes place from BP to BVO(010) within 500 fs but with two distinct processes. In the first phase, the electron transfer proceeds adiabatically and is mainly driven by atomic motions. In the second phase, the electron transfer decays very slowly. The hole-transfer dynamics from BVO(010) to BP exhibits a similar ultrafast decay in the first stage followed by a slow decay; however, there is a comparable amount of hole trapped in a BP state due to a large energy gap from its higher state. These insights may be useful for the design of novel photocatalytic water-splitting materials.

Photoinduced Carrier Dynamics at the Interface of Pentacene and Molybdenum Disulfide

Understanding of photoinduced interfacial carrier dynamics in organic-transition metal dichalcogenides heterostructures is very important for the enhancement of their potential photoelectronic conversion efficiencies. In this work we have used density functional theory (DFT) calculations and DFT-based fewest-switches surface-hopping dynamics simulations to explore the photoinduced hole transfer and subsequent nonadiabatic electron-hole recombination dynamics taking place at the interface of pentacene and MoS₂ in pentacene@MoS₂. Upon photoexcitation the electronic transition mainly occurs on the MoS₂ monolayer, which corresponds to moving an electron to the MoS₂ conduction band. As a result, a hole is left in the valence band. This hole state is energetically lower than certain occupied states of the pentacene molecule; thus, the interfacial hole transfer from MoS₂ to pentacene is favorable in energy. In terms of nonadiabatic dynamics simulations, the hole transfer time to the HOMO-1 state of the pentacene is estimated to be about 600 fs; however, the following hole relaxation process from HOMO-1 to HOMO takes much longer time of ca. 15 ps due to the large energy gap between HOMO-1 and HOMO. Moreover, our results also show that the subsequent radiationless recombination process between the hole transferred to the pentacene molecule and the remaining electron on the MoS₂ CBM needs about 10.2 ns. The computational results shed important mechanistic insights on the interfacial carrier dynamics of mixed-dimensional pentacene@MoS₂. These insights could help to design excellent interfaces for organic-TMDs heterostructures.

Effect of the Interfacial Energy Landscape on Photoinduced Charge Generation at the ZnPc/MoS2 Interface

Monolayer transition-metal dichalcogenide crystals (TMDC) can be combined with other functional materials, such as organic molecules, to form a wide range of heterostructures with tailorable properties. Although a number of works have shown that ultrafast charge transfer (CT) can occur at organic/TMDC interfaces, conditions that would facilitate the separation of interfacial CT excitons into free carriers remain unclear. Here, time-resolved and steady-state photoemission spectroscopy are used to study the potential energy landscape, charge transfer, and exciton dynamics at the zinc phthalocyanine (ZnPc)/monolayer (ML) MoS₂ and ZnPc/bulk MoS₂ interfaces. Surprisingly, although both interfaces have a type-II band alignment and exhibit sub-100 fs CT, the CT excitons formed at the two interfaces show drastically different evolution dynamics. The ZnPc/ML-MoS₂ behaves like typical donor-acceptor interfaces in which CT excitons dissociate into electron-hole pairs. On the contrary, back electron transfer occur at ZnPc/bulk-MoS₂, which results in the formation of triplet excitons in ZnPc. The difference can be explained by the different amount of band bending found in the ZnPc film deposited on ML-MoS₂ and bulk-MoS₂. Our work illustrates that the potential energy landscape near the interface plays an important role in the charge separation behavior. Therefore, considering the energy level alignment at the interface alone is not enough for predicting whether free charges can be generated effectively from an interface.

Nonadiabatic Dynamics Simulations Reveal Distinct Effects of the Thickness of PTB7 on Interfacial Electron and Hole Transfer Dynamics in PTB7@MoS2 Heterostructures

Mixed-dimensional hybrid heterostructures have attracted a lot of experimental attention because they can provide an ideal charge-separated interface for optoelectronic and photonic applications. In this Letter, we have employed first-principles DFT calculations and nonadiabatic dynamics simulations to explore photoinduced interfacial electron and hole transfer processes in two PTB7- nL@MoS₂ models ( n = 1 and 5). The interfacial electron transfer is found to be ultrafast and completes within ca. 10 fs in both PTB7-1L@MoS₂ and PTB7-5L@MoS₂ models, which demonstrates that the electron transfer is not sensitive to the thickness of the PTB7 polymer. Differently, the interfacial hole transfer is sensitive to the thickness of the PTB7 polymer. The transfer time is estimated to be ca. 70 ps in PTB7-1L@MoS₂, while it is significantly accelerated to ca. 1 ps in PTB7-5L@MoS₂. Finally, we have found that the electron transfer is mainly controlled by adiabatic electron evolution, whereas in the hole transfer, nonadiabatic hoppings play a dominant role. These findings are useful for the design of excellent charge-separated interfaces of mixed-dimensional TMD-based heterojunctions for a variety of optoelectronic applications.