Nonadiabatic Molecular Dynamics Simulations Research Articles

Ultrafast hydrogen production in boron/oxygen-codoped graphitic carbon nitride revealed by nonadiabatic dynamics simulations.

Graphitic carbon nitride (g-C3N4 or GCN) shows promise in photocatalytic water splitting, despite facing the challenge of rapid electron-hole recombination. In this study, we investigated the influence of boron/oxygen codoping on the photocatalytic performance of GCN systems for hydrogen generation. First-principles calculations and nonadiabatic molecular dynamics (NAMD) simulations were employed to reveal that the recombination time of photogenerated carriers could be increased by 16% to 64% in the codoped systems compared to the pristine GCN. The time-dependent density functional theory (TDDFT) scheme was utilized to select energy windows and initiate dynamics in cluster models of B/O co-doped heptazine with water molecules. Notably, we observed efficient direct photodissociation of hydrogen atoms from water molecules within 60 fs and proton hops within the hydrogen-bonded network within 80 fs in the co-doped system, diverging from the previously proposed mechanism for pristine heptazine in NAMD simulations. This discovery underscores the significant role of faster proton-coupled electron transfer (PCET) reactions and rapid radiationless relaxation in achieving high photocatalytic efficiency in water splitting. Our work enhances the understanding of the internal mechanism of highly efficient photocatalysts for water splitting and provides a new design strategy for doped GCN.

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.

Exploring hot carrier relaxation dynamics in MoSi2N4 monolayer for photocatalytic applications: Insights into electron-phonon coupling mechanisms

The relaxation process of hot carriers in MoSi2N4 monolayer remains unclear due to the need to sample electron–phonon interaction at different momenta. Utilizing density functional theory (DFT) combined with nonadiabatic molecular dynamics (NAMD), we determined that the material possesses an indirect bandgap of approximately 2.33 eV, primarily originating from the d orbitals of inner Mo atoms. NAMD simulations unveil rapid relaxation of hot electrons and holes, exhibiting lifetimes of approximately 166 fs and 146 fs, respectively. The rapid relaxation of hot carriers couples with the vibrations of outer Si and N atoms, with no inter-valley scattering observed between the K and K′ valleys. Notably, the coupling between the d orbitals of inner Mo atoms and the vibrations of outer Si and N atoms plays a pivotal role in determining the relaxation dynamics of hot carriers, shedding light on the intricate interplay between electronic and vibrational modes in this material. Our research reveals the potential application value of MoSi2N4 monolayer material as a photocatalyst and elucidates a relaxation mechanism governing the dynamics of excited carriers.

Spontaneous polarization dynamics in V-doped monolayer MoS2

Valleytronics has drawn attention both for its fundamental aspects and potential applications in device technology. Lifting the valley degeneracy is an attractive route to achieve stable valley polarization and spontaneous valley polarization. Our first-principles calculations on V-doped monolayer MoS2 reveal a substantial permanent valley splitting of 174 meV, driven by the combined effects of spin-orbit coupling (SOC) and magnetic coupling with the doped atom. Non-adiabatic molecular dynamics simulations show that the spontaneous valley polarization dynamics depends on the SOC, electron-phonon (e-ph) coupling and the scattering with the impurity, and it varies notably with different initial spin orientations of excited holes. For spin-up holes, the process involves e-ph scattering and spin-flip induced by spinor-phonon scattering over several picoseconds. On the other hand, spin-down hole excitation leads to spontaneous polarization through e-ph and impurity scattering on a timescale shorter by an order of magnitude. This study provides insights into the distinct valley polarization dynamics of carriers with different spins, contributing to the design of ultrafast valleytronic and spintronic devices.

Hot Carrier Controlled Nitrogen Fixation Reaction in Metal-Free Boron-Anchored Aza-COF: Insight from Nonadiabatic Molecular Dynamics Simulation.

Designing highly efficient photocatalysts for the production of renewable energy is a challenging task that necessitates simultaneous control of chemical activity and photocarrier dynamics for a particular reaction. To this end, we have investigated the catalytic mechanism and real-time photocarrier dynamics of the nitrogen reduction reaction (NRR) at the metal-free boron-functionalized 2D aza-COF (B-aza-COF), an inexpensive and environmentally friendly semiconductor. By employing density functional theory (DFT) and time-dependent ab initio nonadiabatic molecular dynamics simulation, we have investigated the electronic structure, light harvesting ability, free energy change, and dynamics of photoexcited carriers. Our calculated results reveal that the gas phase N2 molecule can be effectively reduced into NH3 on B-aza-COF under UV-visible light. Therefore, our investigation on the design of efficient photocatalysts for the nitrogen reduction reaction (NRR) provides a cost-effective opportunity for the sustainable production of NH3.

Effects of Surface Defects on Performance and Dynamics of CsPbI2Br Perovskite: First-Principles Nonadiabatic Molecular Dynamics Simulations.

Inorganic mixed-halogen perovskites exhibit excellent photovoltaic properties and stability; yet, their photoelectric conversion efficiency is limited by inherent surface defects. In this work, we study the impact of defects on properties of CsPbI2Br slabs using first-principles calculations, focusing on specific defects such as I vacancy (VI), I interposition (Ii), and I substitution by Pb (PbI). Our findings reveal that these defects affect the geometric and optoelectronic properties as well as dynamics of charge carriers of slabs. We employ two theoretical frameworks (surface hopping and Redfield theory) of nonadiabatic molecular dynamics simulations to comprehensively study relaxation processes and obtain consistent results. The presence of VI reduces carrier lifetimes, while the influence of PbI on carrier lifetimes is negligible. In contrast, Ii defects lead to prolonged carrier lifetimes. These insights provide valuable guidance for the rational design of perovskite photovoltaic devices, aiming to enhance their efficiency and stability.

Hot carrier relaxation dynamics of an aza-covalent organic framework during photoexcitation: An insight from abinitio quantum dynamics.

In order to develop an efficient metal-free solar energy harvester, we herein performed the electronic structure calculation, followed by the hot carrier relaxation dynamics of two dimensional (2D) aza-covalent organic framework by time domain density functional calculations in conjunction with non-adiabatic molecular dynamics (NAMD) simulation. The electronic structure calculation shows that the aza-covalent organic framework (COF) is a direct bandgap semiconductor with acute charge separation and effective optical absorption in the UV-visible region. Our study of non-adiabatic molecular dynamics simulation predicts the sufficiently prolonged electron-hole recombination process (6.8 nanoseconds) and the comparatively faster electron (22.48ps) and hole relaxation (0.51ps) dynamics in this two-dimensional aza-COF. According to our theoretical analysis, strong electron-phonon coupling is responsible for the rapid charge relaxation, whereas the electron-hole recombination process is slowed down by relatively weak electron-phonon coupling, relatively lower non-adiabatic coupling, and quick decoherence time. We do hope that our results of NAMD simulation on exciton relaxation dynamics will be helpful for designing photovoltaic devices based on this two dimensional aza-COF.

Interlayer Carrier Dynamics in Two-Dimensional Perovskites Determined by the Length of Conjugated Organic Cations.

Unlocking the restricted interlayer carrier transfer in a two-dimensional perovskite is a crucial means to achieve the harmonization of efficiency and stability in perovskite solar cells. In this work, the effects of conjugated organic molecules on the interlayer carrier dynamics of 2D perovskites were investigated through nonadiabatic molecular dynamics simulations. We found that elongated conjugated organic cations contributed significantly to the accelerated interlayer carrier dynamics, originating from lowered transport barrier and boosted π-p coupling between organic and inorganic layers. Utilizing conjugated molecules of moderate length as spacer cations can yield both superior efficiency and exceptional stability simultaneously. However, conjugated chains that are too long lead to structural instability and stronger carrier recombination. The potential of conjugated chain-like molecules as spacer cations in 2D perovskites has been demonstrated in our work, offering valuable insights for the development of high-performance perovskite solar cells.

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.

Reveal long-lived hot electrons in 2D indium selenide and ferroelectric-regulated carrier dynamics of InSe/α-In2Se3/InSe heterostructure.

As typical representatives of group III chalcogenides, InSe, α-In2Se3, and β'-In2Se3 have drawn considerable interest in the domain of photoelectrochemistry. However, the microscopic mechanisms of carrier dynamics in these systems remain largely unexplored. In this work, we first reveal that hot electrons in the three systems have different cooling rate stages and long-lived hot electrons, through the utilization of density functional theory calculations and nonadiabatic molecular dynamics simulations. Furthermore, the ferroelectric polarization of α-In2Se3 weakens the nonadiabatic coupling of the nonradioactive recombination, successfully competing with the narrow bandgap and slow dephasing process, and achieving both high optical absorption efficiency and long carrier lifetime. In addition, we demonstrate that the ferroelectric polarization of α-In2Se3 not only enables the formation of the double type-II band alignment in the InSe/α-In2Se3/InSe heterostructure, with the top and bottom InSe sublayers acting as acceptors and donors, respectively, but also eliminates the hindrance of the built-in electric field at the interface, facilitating an ultrafast interlayer carrier transfer in the heterojunction. This work establishes an atomic mechanism of carrier dynamics in InSe, α-In2Se3, and β'-In2Se3 and the regulatory role of the ferroelectric polarization on the charge carrier dynamics, providing a guideline for the design of photoelectronic materials.

Hydrogen generation from direct Z-scheme for photocatalytic overall water splitting with the SiSe/SnSe2 and SiSe/SnSSe heterostructures

The direct Z-schemes of the photocatalytic overall water splitting with the SiSe/SnSe2, SiSe/SnSSe-I, and SiSe/SnSSe-II heterostructures are constructed based on the density functional theory calculations and nonadiabatic molecular dynamics (NAMD) simulations. The maximum solar-to-hydrogen efficiency (ηʹSTH) reaches 19.18% and can be promoted to 28.71% under tensile biaxial strains. NAMD simulations indicate the transfer of the electron for hydrogen evolution reaction (HER) and hole for oxygen evolution reaction (OER) for SiSe/SnSe2 are apparently slower than those for the other two heterostructures, implicating the reduction and oxidation activities of this heterostructure are well-protected. Moreover, the shortest electron-hole recombination time is attributed to SiSe/SnSSe-I, indicating that it holds the best photocatalytic performance. Remarkably, the Gibbs free energies indicate that HER and OER with SiSe/SnSSe-I can spontaneously proceed, while OERs can but HERs cannot spontaneously proceed with the other two heterostructures. Therefore, the newfound heterostructures, especially SiSe/SnSSe-I, are promising candidates in photocatalytic overall water splitting.

Deciphering the Internal Conversion Processes Involved in the Photochemical Ring-Opening of 1,3-Cyclohexadiene by Symmetric sp2-Carbon Substitutions.

Chemical substituents hold the potential to markedly influence the photochemical behavior in molecular systems and assist in gaining a comprehensive understanding of nonadiabatic phenomena. In this study, we have conducted a comparative analysis of the influence of chemical substituents on the photochemical ring-opening of 1,3-cyclohexadiene (CHD), considering four systems: CHD, 2,3-dimethylcyclohexadiene (CHD-Me2-1), 1,4-dimethylcyclohexadiene (CHD-Me2-2), and 1,2,3,4-tetramethylcyclohexadiene (CHD-Me4), using electronic structure theory calculations and nonadiabatic molecular dynamics simulations. Employing extended multistate complete active space second-order perturbation (XMS-CASPT2) theory, we optimized reactants, S1 states, conical intersections (CIns), and products, revealing structural and energetic variations consistent with prior research. Nonadiabatic molecular dynamics simulation was used to gain insights into photochemical dynamics at state-averaged complete active space self-consistent field (SA-CASSCF) theory. CHD-Me4 exhibited reduced carbon-carbon single bond rupture rates, responsible for ring-opening, due to substituent proximity. Further, CHD-Me2-2 and CHD-Me4 displayed prolonged excited-state relaxation times, highlighting notable substituents' impact. Analysis of kinetic energy profiles of specific carbon atoms also revealed restrained atomic displacements, particularly in CHD-Me2-2 and CHD-Me4. These findings advance our understanding of how substituents modulate photochemical reactions in cyclohexadiene derivatives, guiding new molecular design and future research.

Ultrafast electronic relaxation pathways of the molecular photoswitch quadricyclane

The light-induced ultrafast switching between molecular isomers norbornadiene and quadricyclane can reversibly store and release a substantial amount of chemical energy. Prior work observed signatures of ultrafast molecular dynamics in both isomers upon ultraviolet excitation but could not follow the electronic relaxation all the way back to the ground state experimentally. Here we study the electronic relaxation of quadricyclane after exciting in the ultraviolet (201 nanometres) using time-resolved gas-phase extreme ultraviolet photoelectron spectroscopy combined with non-adiabatic molecular dynamics simulations. We identify two competing pathways by which electronically excited quadricyclane molecules relax to the electronic ground state. The fast pathway (<100 femtoseconds) is distinguished by effective coupling to valence electronic states, while the slow pathway involves initial motions across Rydberg states and takes several hundred femtoseconds. Both pathways facilitate interconversion between the two isomers, albeit on different timescales, and we predict that the branching ratio of norbornadiene/quadricyclane products immediately after returning to the electronic ground state is approximately 3:2.

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.

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.

Investigating the photo-relaxation mechanism of 6-azauracil through ab initio nonadiabatic molecular dynamics simulations

On-the-fly ab initio molecular dynamics simulations were performed to study the decay dynamics of photo-excited 6-azauracil. The Zhu–Nakamura trajectory surface hopping method was used to determine the nonadiabatic transitions across different electronic states. The lifetime of the photoexcited S2 state was estimated to be 201 fs, which is in line with the experimentally measured lifetime (Hua et al., 2015). The results of the molecular dynamics simulations revealed that the internal conversion from S2 to S1 is facilitated by the destabilization of the ground state caused by molecular ring puckering.

Unraveling the effect of aromaticity for the dynamics of excited states of single benzene fluorophores.

The photophysical properties of a series of recently synthesized single benzene fluorophores were investigated using ensemble density functional theory calculations. The energetic stability of the ground and excited state species were counterposed against the aromaticity index derived from local vibrational modes. It was found that the large Stokes shift of the fluorophores (up to ca. 5800 cm ) originates from the effect of electron donating and electron withdrawing substituents rather than -delocalization and related (anti-)aromaticity. On the basis of nonadiabatic molecular dynamics simulations, the absence of fluorescence from one of the regioisomers was explained by the occurrence of easily accessible S /S conical intersections below the vertical excitation energy level. It is demonstrated in the manuscript that the analysis of local mode force constants and the related aromaticity index represent a useful tool for the characterization of -delocalization effects in -conjugated compounds.

Modeling pH-Dependent Biomolecular Photochemistry.

The tuning mechanism of pH can be extremely challenging to model computationally in complex biological systems, especially with respect to the photochemical properties. This article reports a protocol aimed at modeling pH-dependent photodynamics using a combination of constant-pH molecular dynamics and semiclassical nonadiabatic molecular dynamics simulations. With retinal photoisomerization in Anabaena sensory rhodopsin (ASR) as a testbed, we show that our protocol produces pH-dependent photochemical properties, such as the isomerization quantum yield or decay rates. We decompose our results into single-titrated residue contributions, identifying some key tuning amino acids. Additionally, we assess the validity of the single protonation state picture to represent the system at a given pH and propose the most populated protein charge state as a compromise between cost and accuracy.

Hydrogen production of overall water splitting with direct Z-scheme driven by antimonene and arsenide nanoribbon heterostructures: Insight from electronic properties and carrier nonadiabatic dynamics

Two nanoribbon heterostructures are suitable for driving efficient direct Z-schemes of the photocatalytic overall water splitting for hydrogen production and are screened out of nine heterostructures constructed from edge-passivated antimonene (SbNR-X) and arsenide (AsNR-X, X = F, Cl, CN) nanoribbons. The solar-to-hydrogen (STH) efficiency and the Gibbs free energies (ΔGs) of the redox reactions at different sites are calculated to identify the preferable Z-schemes. Meanwhile, the transfer and recombination of the photogenerated carriers are explored by using nonadiabatic molecular dynamics (NAMD) simulations. The obtained band alignments and built-in electric fields of the two considered heterostructures can match the requirements of the photocatalytic Z-scheme, and the corresponding STH efficiencies can reach 15.04 % and 30.53 %, respectively. The oxidation evolution reactions can proceed spontaneously, and hydrogen evolution reactions are feasible in thermodynamics. NAMD simulations show that the SbNR-CN/AsNR-CN nanoribbon heterostructure has faster interlayer recombination but slower transfer of the photogenerated electron from the SbNR-CN to AsNR-CN nanoribbons and of the photogenerated hole along the reversible direction, indicating that this heterostructure has a higher carrier utilization rate. These findings confirm that the SbNR-CN/AsNR-CN nanoribbon heterostructure is a promising candidate for developing high STH efficiency materials with direct Z-schemes.

Generalization of One-Center Nonorthogonal Configuration Interaction Singles to Open-Shell Singlet Reference States: Theory and Application to Valence-Core Pump-Probe States in Acetylacetone.

We formulate a one-center nonorthogonal configuration interaction singles (1C-NOCIS) theory for the computation of core excited states of an initial singlet state with two unpaired electrons. This model, which we refer to as 1C-NOCIS two-electron open-shell (2eOS), is appropriate for computing the K-edge near-edge X-ray absorption spectra (NEXAS) of the valence excited states of closed-shell molecules relevant to pump-probe time-resolved (TR) NEXAS experiments. With the inclusion of core-hole relaxation effects and explicit spin adaptation, 1C-NOCIS 2eOS requires mild shifts to match experiment, is free of artifacts due to spin contamination, and can capture the high-energy region of the spectrum beyond the transitions into the singly occupied molecular orbitals (SOMOs). Calculations on water and thymine illustrate the different key features of excited-state NEXAS, namely, the core-to-SOMO transitions as well as shifts and spin-splittings in the transitions analogous to those of the ground state. Simulations of the TR-NEXAS of acetylacetone after excitation to its π → π* singlet excited state at the carbon K-edge, an experiment carried out recently, showcase the ability of 1C-NOCIS 2eOS to efficiently simulate NEXAS based on nonadiabatic molecular dynamics simulations.