The PYXAID Program for Non-Adiabatic Molecular Dynamics in Condensed Matter Systems

Nonadiabatic Dynamics with Exact Factorization: Implementation and Assessment

In this work, we report our implementation of several independent-trajectory mixed-quantum-classical (ITMQC) nonadiabatic dynamics methods based on exact factorization (XF) in the Libra package for nonadiabatic and excited-state dynamics. Namely, the exact factorization surface hopping (SHXF), mixed quantum-classical dynamics (MQCXF), and mean-field (MFXF) are introduced. Performance of these methods is compared to that of several traditional surface hopping schemes, such as the fewest-switches surface hopping (FSSH), branching-corrected surface hopping (BCSH), and the simplified decay of mixing (SDM), as well as conventional Ehrenfest (mean-field, MF) method. Based on a comprehensive set of 1D model Hamiltonians, we find the ranking SHXF ≈ MQCXF > BCSH > SDM > FSSH ≫ MF, with the BCSH sometimes outperforming the XF methods in terms of describing coherences. Although the MFXF method can yield reasonable populations and coherences for some cases, it does not conserve the total energy and is therefore not recommended. We also find that the branching correction for auxiliary trajectories is important for the XF methods to yield accurate populations and coherences. However, the branching correction can worsen the quality of the energy conservation in the MQCXF. Finally, we find that using the time-dependent Gaussian width approximation used in the XF methods for computing decoherence correction can improve the quality of energy conservation in the MQCXF dynamics. The parameter-free scheme of Subotnik for computing the Gaussian widths is found to deliver the best performance in situations where such widths are not known a priori.

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.

Nonadiabatic molecular dynamics with subsystem density functional theory: application to crystalline pentacene

In this work, we report the development and assessment of the nonadiabatic molecular dynamics approach with the electronic structure calculations based on the linearly scaling subsystem density functional method. The approach is implemented in an open-source embedded Quantum Espresso/Libra software specially designed for nonadiabatic dynamics simulations in extended systems. As proof of the applicability of this method to large condensed-matter systems, we examine the dynamics of nonradiative relaxation of excess excitation energy in pentacene crystals with the simulation supercells containing more than 600 atoms. We find that increased structural disorder observed in larger supercell models induces larger nonadiabatic couplings of electronic states and accelerates the relaxation dynamics of excited states. We conduct a comparative analysis of several quantum-classical trajectory surface hopping schemes, including two new methods proposed in this work (revised decoherence-induced surface hopping and instantaneous decoherence at frustrated hops). Most of the tested schemes suggest fast energy relaxation occurring with the timescales in the 0.7-2.0 ps range, but they significantly overestimate the ground state recovery rates. Only the modified simplified decay of mixing approach yields a notably slower relaxation timescales of 8-14 ps, with a significantly inhibited ground state recovery.

Energy-Conserving Surface Hopping for Auger Processes

Auger-type processes are ubiquitous in nanoscale materials because quantum confinement enhances Coulomb interactions, and there exist large densities of states. Modeling Auger processes requires the modification of nonadiabatic (NA) molecular dynamics algorithms to include transitions caused by both NA and Coulomb couplings. The system is split into quantum and classical subsystems, e.g., electrons and vibrations, and as a result, energy conservation becomes nontrivial. In surface hopping, an electronic transition induced by NA coupling is accompanied by a classical velocity readjustment to ensure conservation of the total quantum-classical energy. A different treatment is needed for Auger transitions driven by Coulomb interactions. We develop a nonadiabatic molecular dynamics methodology that meticulously differentiates the energy redistribution accompanying hops induced by the NA coupling and the Coulomb interaction and correctly conserves the total energy at each transition. If the transition is driven by a Coulomb interaction, the hop energy is redistributed within the quantum electronic subsystem only. If the transition is NA, the energy is redistributed between the quantum and classical subsystems. Properly maintaining energy conservation for both types of transitions is crucial to generate a correct order of events, obtain accurate transition times, maintain a proper statistical distribution of state populations, and reach thermodynamic equilibrium. We test the method with biexciton annihilation and Auger-assisted hot electron relaxation in a CdSe quantum dot. The sequence of Auger and phonon-driven processes and the calculated time scales are in excellent agreement with the experimental results. The developed approach can be coupled with any surface-hopping method and provides a crucial practical advance to study charge-carrier dynamics in the nanoscale and condensed matter systems.

Nuclear Quantum Effects Accelerate Charge Separation and Recombination in g-C3N4/TiO2 Heterojunctions

We combined ring-polymer molecular dynamics (MD) and ab initio MD with nonadiabatic MD to study the effects of nuclear quantum effects (NQEs) on interlayer electron transfer and electron-hole recombination at the g-C3N4/TiO2 interface. Our simulations indicate that NQEs significantly affect electron transfer and electron-hole recombination dynamics, accelerating both processes. NQEs deform the g-C3N4 layer and expedite the movement of carbon and nitrogen atoms, thus, enhancing charge delocalization and interlayer coupling. This improved overlap between electronic state wave functions enhances nonadiabatic couplings, facilitating electron transfer and recombination. In addition to the enhanced nonadiabatic couplings accelerating electron transfer, the presence of NQEs narrows the energy gap and delays decoherence by mitigating overall fluctuations, because of restricted TiO2 movements overwhelming enhanced g-C3N4 fluctuations, thereby making the recombination faster. This work provides valuable insights into NQEs in light-element systems and contributes to guiding the development of highly efficient photocatalysts.

Compression of Organic Molecules Coupled with Hydrogen Bonding Extends the Charge Carrier Lifetime in BA2SnI4

Two-dimensional (2D) metal halide perovskites, such as BA2SnI4 (BA═CH3(CH2)3NH3), exhibit an enhanced charge carrier lifetime in experiments under strain. Experiments suggest that significant compression of the BA molecule, rather than of the inorganic lattice, contributes to this enhancement. To elucidate the underlying physical mechanism, we apply a moderate compressive strain to the entire system and subsequently introduce significant compression to the BA molecules. We then perform ab initio nonadiabatic molecular dynamics simulations of nonradiative electron-hole recombination. We observe that the overall lattice compression reduces atomic motions and decreases nonadiabatic coupling, thereby delaying electron-hole recombination. Additionally, compression of the BA molecules enhances hydrogen bonding between the BA molecules and iodine atoms, which lengthens the Sn-I bonds, distorts the [SnI6]4- octahedra, and suppresses atomic motions further, thus reducing nonadiabatic coupling. Also, the elongated Sn-I bonds and weakened antibonding interactions increase the band gap. Altogether, the compression delays the nonradiative electron-hole recombination by more than a factor of 3. Our simulations provide new and valuable physical insights into how compressive strain, accommodated primarily by the organic ligands, positively influences the optoelectronic properties of 2D layered halide perovskites, offering a promising pathway for further performance improvements.

Vertical Dipole Dominates Charge Carrier Lifetime in Monolayer Janus MoSSe

Two-dimensional semiconductor materials with vertical dipoles are promising photocatalysts as vertical dipoles not only promote the electron-hole separation but also enhance the carrier redox ability. However, the influence of vertical dipoles on carrier recombination in such materials, especially the competing relationship between vertical dipoles and band gaps, is not yet clear. Herein, first-principles calculations and nonadiabatic molecular dynamics simulations were combined to clarify the influence of band gap and vertical dipole on the carrier lifetime in Janus MoSSe monolayer. By comparing with the results of MoS2 and MoSe2 as well as exploring the carrier lifetime of MoSSe under strain regulation, it has been demonstrated that the vertical dipole, rather than the band gap, is the dominant factor affecting the carrier lifetime. Strikingly, a linear relationship between the carrier lifetime and vertical dipole is revealed. These findings have important implications for the design of high-performance photocatalysts and optoelectronic devices.

Photoinduced Dynamics of a Single-Walled Carbon Nanotube with a sp3 Defect: The Importance of Excitonic Effects

Herein, we employed linear-response time-dependent functional theory nonadiabatic dynamic simulations to explore the photoinduced exciton dynamics of a chiral single-walled carbon nanotube CNT(6,5) covalently doped with a 4-nitrobenzyl group (CNT65-NO2). The results indicate that the introduction of a sp3 defect leads to the splitting of the degenerate VBM/VBM-1 and CBM/CBM+1 states. Both the VBM upshift and the CBM downshift are responsible for the experimentally observed redshifted E11* trapping state. The simulations reveal that the photoinduced exciton relaxation dynamics completes within 500 fs, which is consistent with the experimental work. On the other hand, we also conducted the nonadiabatic carrier (electron and hole) dynamic simulations, which completely ignore the excitonic effects. The comparison demonstrates that excitonic effects are indispensable. Deep analyses show that such effects induce several dark states, which play an important role in regulating the photoinduced dynamics of CNT65-NO2. The present work demonstrates the importance of including excitonic effects in simulating photoinduced processes of carbon nanotubes. In addition, it not only rationalizes previous experiments but also provides valuable insights that will help in the future rational design of novel covalently doped carbon nanotubes with superior photoluminescent properties.

Prolonged Exciton Lifetime Is Achieved in Porphyrin Nanoring by Template Engineering: A Nonadiabatic Tight Binding Approach

Porphyrin nanoring has been attracting immense attention due to its light harvesting capacity and potential applications in optical, catalysis, sensor, and electronic devices. We demonstrate by nonadiabatic quantum dynamics simulations that the photovoltaic efficiency can be enhanced by template engineering. Altering the hexadentate template (T6) with two tridentate templates (2T3) within the porphyrin ring (P6) cavity accelerated the electron transfer twice and suppressed the electron-hole recombination by nearly three times. The atomistic tight-binding simulation rationalized the dynamics by different localizations of charge of the band edge states, changes in nonadiabatic coupling, alteration in quantum coherence, and involvement of diverse electron-phonon vibrational modes. Further 2T3 templates more strongly hold the P6 ring than T6, reducing the structural fluctuation. As a result, the nonadiabatic coupling becomes weaker and suppresses the carrier recombination. Current atomistic simulation presents a template engineering strategy to enhance the exciton lifetime along with ultrafast charge separation, crucial factors for photovoltaic applications.

Understanding Defects in Perovskite Solar Cells through Computation: Current Knowledge and Future Challenge

Lead halide perovskites with superior optoelectrical properties are emerging as a class of excellent materials for applications in solar cells and light-emitting devices. However, perovskite films often exhibit abundant intrinsic defects, which can limit the efficiency of perovskite-based optoelectronic devices by acting as carrier recombination centers. Thus, an understanding of defect chemistry in lead halide perovskites assumes a prominent role in further advancing the exploitation of perovskites, which, to a large extent, is performed by relying on first-principles calculations. However, the complex defect structure, strong anharmonicity, and soft lattice of lead halide perovskites pose challenges to defect studies. In this perspective, on the basis of briefly reviewing the current knowledge concerning computational studies on defects, this work concentrates on addressing the unsolved problems and proposing possible research directions in future. This perspective particularly emphasizes the indispensability of developing advanced approaches for deeply understanding the nature of defects and conducting data-driven defect research for designing reasonable strategies to further improve the performance of perovskite applications. Finally, this work highlights that theoretical studies should pay more attention to establishing close and clear links with experimental investigations to provide useful insights to the scientific and industrial communities.

Photocatalytic activity of dual defect modified graphitic carbon nitride is robust to tautomerism: machine learning assisted ab initio quantum dynamics

Two-dimensional graphitic carbon nitride (GCN) is a popular metal-free polymer for sustainable energy applications due to its unique structure and semiconductor properties. Dopants and defects are used to tune GCN, and dual defect modified GCN exhibits superior properties and enhanced photocatalytic efficiency in comparison to pristine or single defect GCN. We employ a multistep approach combining time-dependent density functional theory and nonadiabatic molecular dynamics (NAMD) with machine learning (ML) to investigate coupled structural and electronic dynamics in GCN over a nanosecond timescale, comparable to and exceeding the lifetimes of photo-generated charge carriers and photocatalytic events. Although frequent hydrogen hopping transitions occur among four tautomeric structures, the electron-hole separation and recombination processes are only weakly sensitive to the tautomerism. The charge separated state survives for about 10 ps, sufficiently long to enable photocatalysis. The employed ML-NAMD methodology provides insights into rare events that can influence excited state dynamics in the condensed phase and nanoscale materials and extends NAMD simulations from pico- to nanoseconds. The ab initio quantum dynamics simulation provides a detailed atomistic mechanism of photoinduced evolution of charge carriers in GCN and rationalizes how GCN remains photo-catalytically active despite its multiple isomeric and tautomeric forms.

Time-dependent density functional theory with the orthogonal projector augmented wave method

The projector augmented wave (PAW) method of Blöchl linearly maps smooth pseudo wavefunctions to the highly oscillatory all-electron DFT orbitals. Compared to norm-conserving pseudopotentials (NCPP), PAW has the advantage of lower kinetic energy cutoffs and larger grid spacing at the cost of having to solve for non-orthogonal wavefunctions. We earlier developed orthogonal PAW (OPAW) to allow the use of PAW when orthogonal wavefunctions are required. In OPAW, the pseudo wavefunctions are transformed through the efficient application of powers of the PAW overlap operator with essentially no extra cost compared to NCPP methods. Previously, we applied OPAW to DFT. Here, we take the first step to make OPAW viable for post-DFT methods by implementing it in real-time time-dependent (TD) DFT. Using fourth-order Runge-Kutta for the time-propagation, we compare calculations of absorption spectra for various organic and biological molecules and show that very large grid spacings are sufficient, 0.6-0.7 bohr in OPAW-TDDFT rather than the 0.4-0.5 bohr used in traditional NCPP-TDDFT calculations. This reduces the memory and propagation costs by around a factor of 3. Our method would be directly applicable to any post-DFT methods that require time-dependent propagations such as the GW approximation and the Bethe-Salpeter equation.

Real-Time Ab Initio Investigation on Hot Electron Relaxation Dynamics in Silicon

The relaxation of hot electrons in semiconductors is pivotal for both energy harvesting processes and optoelectronics. Utilizing a self-developed non-adiabatic molecular dynamics simulation approach in the momentum space (NAMD_k), we have examined the dynamics of hot electrons in silicon. Whether excited from the Γ or L point, the relaxation dynamics exhibit two distinct stages. Initially, within 100 fs, electrons scatter with phonons throughout the Brillouin zone. Subsequently, over a few picoseconds, they further relax toward the conduction band minimum as a whole. This picture of hot electron relaxation is highly consistent with the quasi-equilibrium hot electron ensemble (HEE) concept. Throughout the hot electron relaxation process, energy transfer to phonons is efficient, leading to time-dependent phonon excitation, which is thoroughly analyzed. This investigation provides a novel perspective on hot electron relaxation in silicon, which holds substantive implications for the realm of photovoltaic and optoelectronic device applications.

Significant Impact of Defect Fluctuation on Charge Dynamics in CsPbI3: A Study Combining Machine Learning with Quantum Dynamics

In this study, we developed a machine-learned force field for CsPbI3 using a neural network potential, enabling molecular dynamics simulations (MD) with ab initio accuracy over nanoseconds. This approach, combined with ab initio MD and nonadiabatic MD, was used to study the charge trapping and recombination dynamics in both pristine and defective CsPbI3. Our simulations revealed key transitions affecting carrier lifetimes, especially in systems with iodine vacancy and interstitial iodine defects. An iodine trimer, formed when iodine replaces cesium, exhibits a high-frequency phonon mode. This mode enhances nonadiabatic coupling, accelerating charge recombination in defective systems compared to pristine ones. In the iodine vacancy system, recombination times varied significantly due to differences in NA coupling and energy gaps. The interplay between nonadiabatic coupling and pure dephasing time is crucial in determining recombination times for interstitial iodine defects. Our findings highlight the role of defect evolution in perovskites, offering insights for enhancing perovskite performance.

The Role of Grain Boundaries on Ion Migration and Charge Recombination in Halide Perovskites

Grain boundaries (GBs) have a significant role in polycrystalline perovskite solar cells (PSCs). However, there is ongoing debate regarding the impact of GBs on the performance and long-term stability of PSCs. Employing the first-principles molecular dynamics for perovskites, the iodine vacancy defect migrations both in bulk and at GBs are investigated. i) The positive iodine vacancy (VI +) is found that have both lower formation energy (1.4 eV) and activation energy (0.18 eV) than those of neutral iodine vacancy (VI), statistically. It indicated the VI + acts as the dominant migrated iodine vacancy rather than VI; ii) the iodine vacancy at GBs has ≈0.48 eV higher activation energy than those in bulk, which leads to the accumulation of iodine vacancy at GBs; iii) the presence of VI + result in a 3-fold increase in charge recombination ratio at GBs, compared to pristine PSCs. Based on quantum molecular dynamics statistical results, which are consistent with experimental measurements, insights into iodine vacancy migration both at GBs and in the bulk are gained. This understanding can be valuable for defects engineering related to ion migration, in order to improve the long-term stability and promote the performance of PSCs.

Impact of large A-site cations on electron-vibrational interactions in 2D halide perovskites: Ab initio quantum dynamics

Using ab initio nonadiabatic molecular dynamics, we study the effect of large A-site cations on nonradiative electron-hole recombination in two-dimensional Ruddlesden-Popper perovskites HA2APb2I7, HA = n-hexylammonium, A = methylammonium (MA), or guanidinium (GA). The steric hindrance created by large GA cations distorts and stiffens the inorganic Pb-I lattice, reduces thermal structural fluctuations, and maintains the delocalization of electrons and holes at ambient and elevated temperatures. The delocalized charges interact more strongly in the GA system than in the MA system, and the charge recombination is accelerated. In contrast, replacement of only some MA cations with GA enhances disorder and increases charge lifetime, as seen in three-dimensional perovskites. This study highlights the key influence of structural fluctuations and disorder on the properties of charge carriers in metal halide perovskites, providing guidance for tuning materials' optoelectronic performance.

Nuclear Quantum Effects Accelerate Charge Recombination but Boost the Stability of Inorganic Perovskites in Mild Humidity

Experiments have demonstrated that mild humidity can enhance the stability of the CsPbBr3 perovskite, though the underlying mechanism remains unclear. Utilizing ab initio molecular dynamics, ring polymer molecular dynamics, and non-adiabatic molecular dynamics, our study reveals that nuclear quantum effects (NQEs) play a crucial role in stabilizing the lattice rigidity of the perovskite while simultaneously shortening the charge carrier lifetime. NQEs reduce the extent of geometric disorder and the number of atomic fluctuations, diminish the extent of hole localization, and thereby improve the electron-hole overlap and non-adiabatic coupling. Concurrently, these effects significantly suppress phonon modes and slow decoherence. As a result, these factors collectively accelerate charge recombination by a factor of 1.42 compared to that in scenarios excluding NQEs. The resulting sub-10 ns recombination time scales align remarkably well with experimental findings. This research offers novel insight into how moisture resistance impacts the stability and charge carrier lifetime in all-inorganic perovskites.

Improved Carrier Separation and Recombination by Ferroelectric Polarization in the CuBiP2Se6/C2N Heterostructure: A Nonadiabatic Molecular Dynamics Study

The rapid recombination of photogenerated carriers heavily restricts the photocatalytic efficiency. Here, we propose a new strategy to improve catalytic efficiency based on the ferroelectric van der Waals heterostructure (CuBiP2Se6/C2N). Combining density functional theory and the nonadiabatic molecular dynamics (NAMD) method, we have systematically analyzed the ground-state properties and carrier dynamics images in the CuBiP2Se6/C2N heterostructure. Our calculations showed that the ferroelectric polarization of CuBiP2Se6 provides the internal driving force for the photogenerated carriers separation. NAMD results demonstrate that the excited-state carrier transfer and recombination processes in the CuBiP2Se6/C2N are consistent with a type II mechanism. Meanwhile, constructing the ferroelectric heterostructure can effectively prolong the carrier lifetime, from ∼65.98 to ∼124.54 ps. Moreover, the high quantum efficiency and tunable band edge positions mean that the CuBiP2Se6/C2N heterostructure is an excellent potential candidate material for photocatalytic water splitting.

Enhanced Charge Separation in Single Atom Cobalt Based Graphitic Carbon Nitride: Time Domain Ab Initio Analysis

In recent years, single atom catalysts have been at the forefront of energy conversion research, particularly in the field of catalysis. Carbon nitrides offer great potential as hosts for stabilizing metal atoms due to their unique electronic structure. We use ab initio nonadiabatic molecular dynamics to study photoexcitation dynamics in single atom cobalt based graphitic carbon nitride. The results elucidate the positive effect of the doped cobalt atom on the electronic structure of GCN. Cobalt doping produces filled midgap states that serve as oxidation centers, advantageous for various redox reactions. The presence of midgap states enables the harvesting of longer wavelength photons, thereby extending the absorption range of solar light. Although doping accelerates charge relaxation overall, charge recombination is significantly slower than charge separation, creating beneficial conditions for catalysis applications. The simulations reveal the detailed microscopic mechanism underlying the improved performance of the doped system due to atomic defects and demonstrate an effective charge separation strategy to construct highly efficient and stable photocatalytic two-dimensional materials.

Efficient photoreduction of carbon dioxide to ethanol using diatomic nitrogen-doped black phosphorus

Successful conversion of CO2 into C2 products requires the development of new catalysts that overcome the difficulties in efficient light harvesting and CO-CO coupling. Herein, density functional theory (DFT) is used to assess the photoreduction properties of nitrogen-doped black phosphorus. The geometric structure, redox potential, first step of hydrogenation activation, CO desorption, and CO-CO coupling are systematically calculated, based on which the diatomic nitrogen-doped black phosphorus (N2@BPV) stands out. The calculated results of the CO2RR pathway demonstrate that N2@BPV has excellent selectivity and high activity for CH3CH2OH production. The results of the time-dependent ab initio nonadiabatic molecular dynamics simulation show that the diatomic N active sites of N2@BPV facilitate charge separation and inhibit electron-hole recombination. In addition, the activation mechanism of CO2 is studied. The main reason for CO2 activation is attributed to the imbalance in electron transfer that destroys the symmetry of CO2. We expect that our study will offer some theoretical guidance in CO2 conversion.

Excitonic Effects of the Excited-State Photocatalytic Reaction at the Molecule/Metal Oxide Interface

Excitonic effects caused by the Coulomb interaction between electrons and holes play a crucial role in photocatalysis at the molecule/metal oxide interface. As an ideal model for investigating the excitonic effect, coadsorption and photodissociation of water and methanol molecules on titanium dioxide involve complex ground-state thermalcatalytic and excited-state photocatalytic reaction processes. Herein, we systemically investigate the excited-state electronic structures of the coadsorption of H2O and CH3OH molecules on a rutile TiO2(110) surface by linear-response time-dependent density functional theory calculations and probe the reaction path for generating HCOOH or CO2, from ground-state and excited-state perspectives. The reaction barriers in excited-state calculations are significantly different from those in ground-state calculations during three processes, with the largest decrease being 0.94 eV for the Ti5c-O-CH2-O-Ti5c formation process.

Estimating Nonradiative Excited-State Lifetimes in Photoactive Semiconducting Nanostructures

The time evolution of the exciton generated by light adsorption in a photocatalyst is an important feature that can be approached from full nonadiabatic molecular dynamics simulations. Here, a crucial parameter is the nonradiative recombination rate between the hole and the electron that form the exciton. In the present work, we explore the performance of a Fermi's golden rule-based approach on predicting the recombination rate in a set of photoactive titania nanostructures, relying solely on the coupling of the ground and first excited state. In this scheme the analysis of the first excited state is carried out by invoking Kasha's rule thus avoiding computationally expensive nonadiabatic molecular dynamics simulations and resulting in an affordable estimate of the recombination rate. Our results show that, compared to previous ones from nonadiabatic molecular dynamics simulations, semiquantitative recombination rates can be predicted for the smaller titania nanostructures, and qualitative values are obtained from the larger ones. The present scheme is expected to be useful in the field of computational heterogeneous photocatalysis whenever a complex and computationally expensive full nonadiabatic molecular dynamics cannot be carried out.

Unveiling the Dual Role of Humidity: The Interplay with Defects Manipulating the Charge Carrier Lifetime in Metal Halide Perovskites

Humidity has exhibited experimentally either beneficial or detrimental effects on the charge carrier lifetime of CH3NH3PbI3 perovskites, leaving the mechanism unresolved. By using ab initio nonadiabatic molecular dynamics simulations, we unveil the dual role of humidity stemming from the complex interplay between water and defects. Beneficially, water passivates iodine vacancies (VI) or grain boundaries (GBs), mitigating electron trapping by reducing nonadiabatic coupling and delaying charge recombination. However, when VI and GBs coexist, water molecules make the two unsaturated lead atoms approach closer and exacerbate electron trapping by deepening the Pb-dimer electron trap that was created by the VI defect, shortening the carrier lifetime to half of pristine CH3NH3PbI3. The study uncovers the origin of the positive and negative effects of humidity on the charge carrier lifetime of perovskites and offers strategies for improving perovskite devices, particularly by avoiding simultaneous point defects and GBs.

Two-dimensional Janus SbTeBr/SbSI heterostructures as multifunctional optoelectronic systems with efficient carrier separation

The stacking of two-dimensional (2D) materials is a highly effective approach in the design of high-performance optoelectronic devices. In this study, we propose a novel Janus monolayer-based 2D/2D van der Waals heterostructure (vdWH) called SbTeBr/SbSI. Starting from its most stable binding configuration, we systematically examined its electronic, optical, mechanical and dynamical properties. The SbTeBr/SbSI vdWH exhibits a type II band arrangement, with an indirect bandgap of 1.28 eV and strong light absorption capabilities in the visible range, achieving an absorption coefficient of 4 × 105 cm-1. These desirable properties suggest that SbTeBr/SbSI holds promise as a material for solar cells, potentially achieving a power conversion efficiency of 8.3%. The dipole-induced electric field in the SbTeBr/SbSI vdWH leads to significant differences in the mobilities of different carriers, which is a critical aspect in suppressing the recombination of photogenerated carriers. Additionally, according to the simulations of nonadiabatic molecular dynamics, a long electron-hole recombination time of 133 ps is predicted. Thus, the SbTeBr/SbSI heterostructure enables efficient charge separation, demonstrating its potential as a high-performance optoelectronic material.

Selective Interfacial Excited-State Carrier Dynamics and Efficient Charge Separation in Borophene-Based Heterostructures

Borophene-based van der Waals heterostructures have demonstrated enormous potential in the realm of optoelectronic and photovoltaic devices, which has sparked a wide range of interest. However, a thorough understanding of the microscopic excited-state electronic dynamics at interfaces is lacking, which is essential for determining the macroscopic optoelectronic and photovoltaic performance of borophene-based devices. In this study, photoexcited carrier dynamics of β12 , χ3 , and α΄ borophene/MoS2 heterostructures are systematically studied based on time-domain nonadiabatic molecular dynamics simulations. Different Schottky contacts are found in borophene/semiconductor heterostructures. The interplay between Schottky barriers, electronic coupling, and the involvement of different phonon modes collectively contribute to the unique carrier dynamics in borophene-based heterostructures. The diverse borophene allotropes within the heterostructures exhibit distinct and selective carrier transfer behaviors on an ultrafast timescale: electrons tunnel into α΄ borophene with an ultrafast transfer rate (≈29 fs) in α΄/MoS2 heterostructures, whereas β12 borophene only allows holes to migrate with a lifetime of 176 fs. The feature enables efficient charge separation and offers promising avenues for applications in optoelectronic and photovoltaic devices. This study provides insight into the interfacial carrier dynamics in borophene-based heterostructures, which is helpful in further design of advanced 2D boron-based optoelectronic and photovoltaic devices.

Rational Design of Full-Color Fluorescent C3N Quantum Dots

Carbon-based quantum dots (QDs) exhibit unique photoluminescence due to size-dependent quantum confinement, giving rise to fascinating full-color emission properties. Accurate emission calculations using time-dependent density functional theory are a time-costing and expensive process. Herein, we employed an artificial neural network (ANN) combined with statistical learning to establish the relationship between geometrical/electronic structures of ground states and emission wavelength for C3N QDs. The emission energy of these QDs can be doubly modulated by size and edge effects, which are governed by the number of C4N2 rings and the CH group, respectively. Moreover, these two structural characteristics also determine the phonon vibration mode of C3N QDs to harmonize the emission intensity and lifetime of hot electrons in the electron-hole recombination process, as indicated by nonadiabatic molecular dynamics simulation. These computational results provide a general approach to atomically precise design the full-color fluorescent carbon-based QDs with targeted functions and high performance.

Ultrafast Charge Transfer and Recombination Dynamics in Monolayer-Multilayer WSe2 Junctions Revealed by Time-Resolved Photoemission Electron Microscopy

The ultrafast carrier dynamics of junctions between two chemically identical, but electronically distinct, transition metal dichalcogenides (TMDs) remains largely unknown. Here, we employ time-resolved photoemission electron microscopy (TR-PEEM) to probe the ultrafast carrier dynamics of a monolayer-to-multilayer (1L-ML) WSe2 junction. The TR-PEEM signals recorded for the individual components of the junction reveal the sub-ps carrier cooling dynamics of 1L- and 7L-WSe2, as well as few-ps exciton-exciton annihilation occurring on 1L-WSe2. We observe ultrafast interfacial hole (h) transfer from 1L- to 7L-WSe2 on an ∼0.2 ps time scale. The resultant excess h density in 7L-WSe2 decays by carrier recombination across the junction interface on an ∼100 ps time scale. Reminiscent of the behavior at a depletion region, the TR-PEEM image reveals the h density accumulation on the 7L-WSe2 interface, with a decay length ∼0.60 ± 0.17 μm. These charge transfer and recombination dynamics are in agreement with ab initio quantum dynamics. The computed orbital densities reveal that charge transfer occurs from the basal plane, which extends over both 1L and ML regions, to the upper plane localized on the ML region. This mode of charge transfer is distinctive to chemically homogeneous junctions of layered materials and constitutes an additional carrier deactivation pathway that should be considered in studies of 1L-TMDs found alongside their ML, a common occurrence in exfoliated samples.

Resolving the Puzzle of Charge Carrier Lifetime in ZnO by Revisiting the Role of Oxygen Vacancy

Zinc oxide (ZnO) is a wide bandgap prototypical n-type semiconductor due to the presence of intrinsic oxygen vacancies (VO). The VO can readily transfer to the most energetically favorable +2 charged VO (VO2+) by losing two electrons mediated by the metastable VO1+ defect. Nevertheless, the influence of charged VO on the charge dynamics in ZnO and the underlying mechanisms remain elusive. By performing nonadiabatic molecular dynamics simulations of the charge trapping and recombination processes, we show that both VO1+ and VO2+ slow down the nonradiative electron-hole recombination via assisted defect states and, thus, extending charge carrier lifetime compared to pristine ZnO. Our study contributes to identifying the different recombination pathways that take place in VO1+ and VO2+ of n-type ZnO systems, providing useful guidance for designing high-performance ZnO-based devices.

Nonadiabatic Molecular Dynamics with Non-Condon Effect of Charge Carrier Dynamics

Nonradiative multiphonon transitions play a crucial role in understanding charge carrier dynamics. To capture the non-Condon effect in nonadiabatic molecular dynamics (NA-MD), we develop a simple and accurate method to calculate noncrossing and crossing k-point NA coupling in momentum space on an equal footing and implement it with a trajectory surface hopping algorithm. Multiple k-point MD trajectories can provide sufficient nonzero momentum multiphonons coupled to electrons, and the momentum conservation is maintained during nonvertical electron transition. The simulations of indirect bandgap transition in silicon and intra- and intervalley transitions in graphene show that incorporation of the non-Condon effect is needed to correctly depict these types of charge dynamics. In particular, a hidden process is responsible for the delayed nonradiative electron-hole recombination in silicon: the thermal-assisted rapid trapping of an excited electron at the conduction band minimum by a long-lived higher energy state through a nonvertical transition extends charge carrier lifetime, approaching 1 ns, which is about 1.5 times slower than the direct bandgap recombination. For graphene, intervalley scattering takes place within about 225 fs, which can occur only when the intravalley relaxation proceeds to about 50 fs to gain enough phonon momentum. The intra- and intervalley scattering constitute energy relaxation, which completes within sub-500 fs. All the simulated time scales are in excellent agreement with experiments. The study establishes the underlying mechanisms for a long-lived charge carrier in silicon and valley scattering in graphene and underscores the robustness of the non-Condon approximation NA-MD method, which is suitable for rigid, soft, and large defective systems.

Activating two-dimensional semiconductors for photocatalysis: a cross-dimensional strategy

The emerging two-dimensional (2D) semiconductors substantially extend materials bases for versatile applications such as semiconductor photocatalysis demanding semiconductive matrices and large surface areas. The dimensionality, while endowing 2D semiconductors the unique properties to host photocatalytic functionality of pollutant removal and hydrogen evolution, hurdles the activation paths to form heterogenous photocatalysts where the photochemical processes are normally superior over these on the mono-compositional counterparts. In this perspective, we present a cross-dimensional strategy to employ thenD (n= 0-2) clusters or nanomaterials as activation partners to boost the photocatalytic activities of the 2D semiconductors. The formation principles of heterogenous photocatalysts are illustrated specifically for the 2D matrices, followed by selection criteria of them among the vast 2D database. The computer investigations are illustrated in the density functional theory route and machine learning benefitted from the vast samples in the 2D library. Synthetic realizations and characterizations of the 2D heterogenous systems are introduced with an emphasis on chemical methods and advanced techniques to understand materials and mechanistic studies. The perspective outlooks cross-dimensional activation strategies of the 2D materials for other applications such as CO2removal, and materials matrices in other dimensions which may inspire incoming research within these fields.

Controlling Charge Carrier Dynamics in Porphyrin Nanorings by Optically Active Templates

Understanding the dynamics of photogenerated charge carriers is essential for enhancing the performance of solar and optoelectronic devices. Using atomistic quantum dynamics simulations, we demonstrate that a short π-conjugated optically active template can be used to control hot carrier relaxation, charge carrier separation, and carrier recombination in light-harvesting porphyrin nanorings. Relaxation of hot holes is slowed by 60% with an optically active template compared to that with an analogous optically inactive template. Both systems exhibit subpicosecond electron transfer from the photoactive core to the templates. Notably, charge recombination is suppressed 6-fold by the optically active template. The atomistic time-domain simulations rationalize these effects by the extent of electron and hole localization, modification of the density of states, participation of distinct vibrational motions, and changes in quantum coherence. Extension of the hot carrier lifetime and reduction of charge carrier recombination, without hampering charge separation, demonstrate a strategy for enhancing efficiencies of energy materials with optically active templates.

Blocking recombination centers by controlling the charge density of a sulfur vacancy in antimony trisulfide

By performing nonadiabatic molecular dynamics combined with ab initio time-domain density functional theory, we have explored the effects of the charge density of a sulfur vacancy on charge trapping and recombination in antimony trisulfide (Sb2S3). The simulations demonstrate that, compared to an antimony vacancy, the sulfur vacancy generates a high charge density trap state within the band gap. This state acts as the recombination center and provides new channels for charge carrier relaxation. Filling the sulfur vacancy with electron donors elevates the defect state to the Fermi level due to the introduced extra electrons. In contrast, the electron acceptor lowers the charge density of the sulfur vacancy by capturing its local electrons, eliminating the charge recombination center and extending the photo-generated charge carrier lifetime. Additionally, compared with electron injection, hole injection can also decrease the charge density of the trap state via neutralizing its local electronic states, eliminate the trap state within the band gap, and suppress nonradiative electron-hole recombination. This study is expected to shed new light on the blocking recombination centers and provide valuable insights into the design of high-performance solar cells.

Exploration of the origin of the excellent charge-carrier dynamics in Ruddlesden-Popper oxysulfide perovskite Y2Ti2O5S2

Although the efficient separation of electron-hole (e-h) pairs is one of the most sought-after electronic characteristics of materials, due to thermally induced atomic motion and other factors, they do not remain separated during the carrier transport process, potentially leading to rapid carrier recombination. Here, we utilized real-time time-dependent density functional theory in combination with nonadiabatic molecular dynamics (NAMD) to explore the separated dynamic transport path within Ruddlesden-Popper oxysulfide perovskite Y2Ti2O5S2 caused by the dielectric layer and phonon frequency difference. The underlying origin of the efficient overall water splitting in Y2Ti2O5S2 is systematically explored. We report the existence of the bi-directional e-h separate-path transport, in which, the electrons transport in the Ti2O5 layer and the holes diffuse in the rock-salt layer. This is in contrast to the conventional e-h separated distribution with a crowded transport channel, as observed in SrTiO3 and hybrid perovskites. Such a unique feature finally results in a long carrier lifetime of 321 ns, larger than that in the SrTiO3 perovskite (160 ns) with only one carrier transport channel. This work provides insights into the carrier transport in lead-free perovskites and yields a novel design strategy for next-generation functionalized optoelectronic devices.

Time-Domain View of Polaron Dynamics in Metal Oxide Photocatalysts

The polaron is a fundamental physical phenomenon in transition metal oxides (TMOs), and it has been studied extensively for decades. However, the implication of a polaron on photochemistry is still ambiguous. As such, understanding the fundamental properties and controlling the dynamics of polarons at the atomistic level is desired. In this Perspective, we seek to highlight the recent advances in studying small polarons in TMOs, with a particular focus on nonadiabatic molecular dynamics at the ab initio level, and discuss the implications for photocatalysis from the aspects of the structure, intrinsic physical properties, formation, migration, and recombination of small polarons. Finally, various methods were proposed to advance our understanding of manipulating the small-polaron dynamics, and strategies to design high-performance TMO-based photoelectrodes were discussed.

Lattice Distortion and Low-Frequency Anharmonic Phonons Suppress Charge Recombination in Lead Halide Perovskites upon Pseudohalide Doping: Time-Domain Ab Initio Analysis

Perovskite solar cells have witnessed a surge in interest as a promising technology for low-cost, high-efficiency photovoltaics with certified power conversion efficiencies beyond 25%. However, their commercial development is hindered by poor stability and nonradiative losses that restrict their approach to the theoretical efficiency limit. Using ab initio nonadiabatic molecular dynamics, we demonstrate that nonradiative charge recombination is suppressed when the iodide in formamidinium lead iodide (FAPbI3) is partially replaced with pseudohalide anions (SCN-, BF4-, and PF6-). The replacement breaks the symmetry of the system and creates local structural distortion and dynamic disorder, decreasing electron-hole overlap and nonadiabatic electron-vibrational coupling. The charge carrier lifetime is found to increase with increased structural distortion and is the longest for PF6-. This work is fundamentally relevant to the design of high-performance perovskite materials for optoelectronic applications.

Electron- versus Spin-Phonon Coupling Governs the Temperature-Dependent Carrier Dynamics in the Topological Insulator Bi2Te3

Ultrafast charge and spin dynamics have immense effects on the applications of topological insulators (TIs). By performing spin-adiabatic nonadiabatic molecular dynamics simulations in the presence of electron-phonon (e-ph) and spin-phonon couplings, we investigate temperature-dependent intra- and interband charge and spin relaxation dynamics via the bulk and surface paths in the three-dimensional TI Bi2Te3. The e-ph coupling dominates charge relaxation in the bulk path, and the relaxation rate is positively correlated with temperature due to the large energy gaps and weak spin polarization. Conversely, the relaxation dynamics exhibits an opposite temperature dependence in the surface path because of electron re-excitation and spin mismatching induced by spin-phonon coupling, which arises from small energy gaps and strong spin polarization. The two mechanisms rationalize the charge carriers being long-lived in the bulk and surface phases at low and room temperature, respectively. Additionally, strong thermal fluctuations of the topological states' magnetic moments destroy the spin-momentum locking and trigger backscattering at room temperature.

Ultrafast many-body bright-dark exciton transition in anatase TiO2

The momentum-forbidden dark excitons can have a pivotal role in quantum information processing, Bose-Einstein condensation, and light-energy harvesting. Anatase TiO2 with an indirect band gap is a prototypical platform to study bright to momentum-forbidden dark exciton transition. Here, we examine, by GW plus the real-time Bethe-Salpeter equation combined with the nonadiabatic molecular dynamics (GW + rtBSE-NAMD), the many-body transition that occurs within 100 fs from the optically excited bright to the strongly bound momentum-forbidden dark excitons in anatase TiO2. Comparing with the single-particle picture in which the exciton transition is considered to occur through electron-phonon scattering, within the GW + rtBSE-NAMD framework, the many-body electron-hole Coulomb interaction activates additional exciton relaxation channels to notably accelerate the exciton transition in competition with other radiative and nonradiative processes. The existence of dark excitons and ultrafast bright-dark exciton transitions sheds insights into applications of anatase TiO2 in optoelectronic devices and light-energy harvesting as well as the formation process of dark excitons in semiconductors.

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.

Twisting Aromaticity and Photoinduced Dynamics in Hexapole Helicenes

Curved aromatic molecules are attractive electronic materials, where an additional internal strain uniquely modifies their structure, aromaticity, dynamics, and optical properties. Helicenes are examples of such twisted conjugated systems. Herein, we analyze the photoinduced dynamics in different stereoisomers of a hexapole helicene by using nonadiabatic excited-state molecular dynamics simulations. We explore how changes in symmetry and structural distortion modulate the intramolecular energy redistribution. We find that distinct helical assembly leads to different rigid distorted structures that in turn impact the nonradiative energy relaxation and ultimately formation of the self-trapped exciton. Subsequently, the value of the twisting angles relative to the central triphenylene core structure controls the global molecular aromaticity and electronic localization during the internal conversion process. Our work sheds light on how the future synthesis of novel curved aromatic compounds can be directed to attain specific desired electronic properties through the modulation of their twisted aromaticity.

The sp3 Defect Decreases Charge Carrier Lifetime in (8,3) Single-Walled Carbon Nanotubes

A recent experimental approach introduces sp3 defects into single-walled carbon nanotubes (SWNTs) through controlled functionalization with guanine, resulting in a decrease in charge carrier lifetime. However, the physical mechanism behind this phenomenon remains unclear. We employ nonadiabatic molecular dynamics to systematically model the nonradiative recombination process of electron-hole pairs in SWNTs with sp3 defects generated by a guanine molecule. We demonstrate that the introduction of sp3 defects creates an overlapping channel between the highest occupied (HOMO) and lowest unoccupied molecular orbital (LUMO), significantly enhancing the nonadiabatic (NA) coupling and leading to a 4.7-fold acceleration in charge carrier recombination compared to defect-free SWNTs. The charge carrier recombination slows significantly at a lower temperature (50 K) due to the weakening of the NA coupling. Our results rationalize the accelerated recombination of charge carriers in SWNTs with sp3 defects in experiments and contribute to a deeper understanding of the carrier dynamics in SWNTs.

Lessons Learned from Catalysis to Qubits: General Strategies to Build Accessible and Accurate First-Principles Models of Point Defects

Defects and dopants play critical roles in defining the properties of a material. Achieving a mechanistic understanding of how such properties arise is challenging with current experimental methods, and computational approaches suffer from significant modeling limitations that frequently require a posteriori fitting. Consequently, the pace of dopant discovery as a means of tuning material properties for a particular application has been slow. However, recent advances in computation have enabled researchers to move away from semiempirical schemes to reposition density functional theory as a predictive tool and improve the accessibility of highly accurate first-principles methods to all researchers. This Perspective discusses some of these recent achievements that provide more accurate first-principles geometric, thermodynamic, optical, and electronic properties simultaneously. Advancements related to supercells, basis sets, functionals, and optimization protocols, as well as suggestions for evaluating the quality of a computational model through comparison to experimental data, are discussed. Moreover, recent computational results in the fields of energy materials, heterogeneous catalysis, and quantum informatics are reviewed along with an evaluation of current frontiers and opportunities in the field of computational materials chemistry.

Convergence of Time-Derivative Nonadiabatic Couplings in Plane-Wave DFT Calculations

Accurate prediction of charge carrier relaxation rates is essential to design molecules and materials with the desired photochemical properties for applications like photocatalysis and solar energy conversion. Nonadiabatic molecular dynamics allows one to simulate the relaxation process of excited charge carriers. Plane-wave density functional theory (DFT) calculations make the time-derivative nonadiabatic couplings (TNACs) simple to compute because the basis is independent of the atomic positions. However, the effect of the kinetic energy cutoff for the plane-wave basis on the accuracy of the dynamics has not been studied. Here, we examine the effect of the kinetic energy cutoff on the TNACs and decay time scales for the prototypical model system of tetracene. These calculations show that the choice of kinetic energy cutoff can change the relaxation time by up to 30%. The relaxation times of states that have small TNACs to other states or are far from degenerate are more sensitive to the kinetic energy cutoff than those of states with large TNACs or near degeneracies. A kinetic energy cutoff of 60 Ry is sufficient for all states to reach semiquantitative agreement (absolute error <10%) with the decay times of our 110 Ry reference data, and a cutoff of 80 Ry is required for all states to reach quantitative agreement (absolute error <2%).

Substrate Ferroelectric Proximity Effects Have a Strong Influence on Charge Carrier Lifetime in Black Phosphorus

By stacking monolayer black phosphorus (MBP) with nonpolarized and ferroelectric polarized bilayer hexagonal boron nitride (h-BN), we demonstrate that ferroelectric proximity effects have a strong influence on the charge carrier lifetime of MBP using nonadiabatic (NA) molecular dynamics simulations. Through enhancing the motion of phosphorus atoms, ferroelectric polarization enhances the overlap of electron-hole wave functions that improves NA coupling and decreases the bandgap, resulting in a rapid electron-hole recombination completing within a quarter of nanoseconds, which is two times shorter than that in nonpolarized stackings. In addition to the dominant in-plane Ag2 mode in free-standing MBP, the out-of-plane high-frequency Ag1 and low-frequency interlayer breathing modes presented in the heterojunctions drive the recombination. Notably, the resonance between the breathing mode within bilayer h-BN and the B1u mode of MBP provides an additional nonradiative channel in ferroelectric stackings, further accelerating charge recombination. These findings are crucial for charge dynamics manipulation in two-dimensional materials via substrate ferroelectric proximity effects.

The impact of spatially heterogeneous chemical doping on the electronic properties of CdSe quantum dots: insights from ab initio computation

The introduction of copper (Cu) impurity in semiconductor CdSe quantum dots (QDs) gives rise to unique photoluminescence (PL) bands exhibiting distinctive characteristics, like broad line width, significant Stokes shift, and complex temporal decay. The atomistic origins of these spectral features are yet to be understood comprehensively. We employed multiple computational techniques to systematically study the impact of the spatial heterogeneity of Cu atoms on the stability and photophysical properties, including the emission linewidth of doped QDs under ambient conditions. The Cu substitution introduces a spin-polarized intragap state, the energetic position of which is strongly dependent on the dopant location and causes spectral broadening in QD ensembles. Furthermore, the dopant dynamics under ambient conditions are significantly influenced by the specific arrangement of Cu within the QDs. The dynamic electronic structures of surface-doped CdSe illustrate more pronounced perturbations and vary the mid-gap state position more drastically than those of the core-doped QDs. Vibronic coupling broadens the photoluminescence peaks associated with the conduction band-to-defect level transition for individual QDs. These insights into the dynamic structure-photophysical property relationship suggest viable approaches, such as tuning the operational temperature and selective co-doping, to enhance the functional performances of doped CdSe QDs strategically.

Direct Observation of Photon Induced Giant Band Renormalization in 2D PdSe2 Dichalcogenide by Transient Absorption Spectroscopy

Insight into fundamental light-matter interaction as well as underlying photo-physical processes is crucial for the development of novel optoelectronic devices. Palladium diselenide (PdSe2 ), an important representative of emerging 2D noble metal dichalcogenides, has gain considerable attention owing to its unique optical, physical, and chemical properties. In this study, 2D PdSe2 nanosheets (NSs) are prepared using the liquid-phase exfoliation method. A broadband carrier relaxation dynamics from visible to near-infrared bands are revealed using a time-resolved transient absorption spectrometer, giving results that indicate band filling and bandgap renormalization (BGR) effects in the 2D PdSe2 NSs. The observed blue-shift of the transient absorption spectra at the primary stage and the subsequent red-shift can be ascribed to this BGR effect. These findings reveal the many-body character of the 2D TMDs material and may hold keys for applications in the field of optoelectronics and ultrafast photonics.

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.

Data-Driven Discovery of a Covalent Organic Framework Heterojunction as Efficient Photocatalysts for Overall Solar Water Splitting

Searching for highly efficient visible-light photocatalysts is a high-cost and time-consuming process in the water splitting field. The integration of data-driven screening based on the database and density functional theory calculations represents a promising approach. In this study, we first present a topologically assembled single-layer covalent organic framework (COF) that is used to build a COF heterojunction database via AA stacking. Then we propose a systematic search procedure for COF heterojunctions as overall solar water splitting photocatalysts, including suitable band gap (screen 1), appropriate band edge position (screen 2), spontaneous catalytic reactions for water splitting (screen 3), and efficient separation of photogenerated electrons and holes (proof). Finally, we successfully identify 1 heterojunction from the pool of 222 items as an efficient photocatalyst for overall solar water splitting. Clearly, this kind of data-driven screening procedure, based on a COF heterojunction database, opens up new avenues and inspires the development of high-performance photocatalysts.

Twist Angle-Dependent Intervalley Charge Carrier Transfer and Recombination in Bilayer WS2

A twist angle at a van der Waals junction provides a handle to tune its optoelectronic properties for a variety of applications, and a comprehensive understanding of how the twist modulates electronic structure, interlayer coupling, and carrier dynamics is needed. We employ time-dependent density functional theory and nonadiabatic molecular dynamics to elucidate angle-dependent intervalley carrier transfer and recombination in bilayer WS2. Repulsion between S atoms in twisted configurations weakens interlayer coupling, increases the interlayer distance, and softens layer breathing modes. Twisting has a minor influence on K valleys while it lowers Γ valleys and raises Q valleys because their wave functions are delocalized between layers. Consequently, the reduced energy gaps between the K and Γ valleys accelerate the hole transfer in the twisted structures. Intervalley electron transfer proceeds nearly an order of magnitude faster than hole transfer. The more localized wave functions at K than Q values and larger bandgaps result in smaller nonadiabatic couplings for intervalley recombination, making it 3-4 times slower in twisted than high-symmetry structures. B2g breathing, E2g in-plane, and A1g out-of-plane modes are most active during intervalley carrier transfer and recombination. The faster intervalley transfer and extended carrier lifetimes in twisted junctions are favorable for optoelectronic device performance.

Surface Hydroxylation during Water Splitting Promotes the Photoactivity of BiVO4(010) Surface by Suppressing Polaron-Mediated Charge Recombination

Polaron-based electron transport restricts the photoelectrochemical (PEC) water splitting efficiency of BiVO4. However, the location and dynamics of polarons are significantly dependent on the surface hydroxylation. By performing ab initio nonadiabatic molecular dynamics simulations, we demonstrated that hydroxylation of BiVO4(010) surface greatly alleviates the detrimental effect of oxygen-vacancy-induced electron polaron (EP). Surface hydroxylation stabilizes the EP at the surface to facilitate water splitting, makes the polaron a shallow localized state, and reduces the intensity of high-frequency V-O bond stretching vibrations. By decreasing the nonadiabatic coupling and decoherence time, the charge carrier lifetimes are extended by 1-3 orders of magnitude depending on the hydroxylation coverage. Our study not only reveals that the surface hydroxylation mitigated detrimental impacts of polarons in metal oxides but also provided valuable insights into the benign effect of intermediate species on the photocatalytic reactivity.

Synergetic Role of Thermal Catalysis and Photocatalysis in CO2 Reduction on Cu2/MoS2

Effective activation of CO2 is a primarily challenging issue in CO2 reduction to value-added hydrocarbon chemicals, due to the large energy gap between the highest-occupied and lowest-unoccupied molecular orbitals (HOMO-LUMO). Here, we employ state-of-the-art first-principles calculations to explore the synergetic role of thermal catalysis and photocatalysis in CO2 reduction, on typical single-atom scale catalyst, i.e., Cu2 magic cluster on a semiconducting two-dimensional MoS2 substrate. It is identified that only about 1% of the hot electrons excited from the MoS2 substrate by at least 6.3 eV photons may be trapped by the inert CO2 molecule at the expense of 400 fs. Moreover, the physisorption-to-chemisorption transition of CO2 can be observed within 500 fs upon overcoming an about 0.05 eV energy barrier. Contrastingly, upon chemisorption, the activated CO2δ- species may trap about 7% of the hot electron excited from the MoS2 substrate by about 2.5 eV visible photons, with a cost of 140 fs.