Ab Initio Molecular Dynamics Simulations Research Articles

The anodic dissolution kinetics of Mg alloys in water based on ab initio molecular dynamics simulations

AbstractThe corrosion susceptibility of magnesium (Mg) alloys presents a significant challenge for their broad application. Although there have been extensive experimental and theoretical investigations, the corrosion mechanisms of Mg alloys are still unclear, especially the anodic dissolution process. Here, a thorough theoretical investigation based on ab initio molecular dynamics and metadynamics simulations has been conducted to clarify the underlying corrosion mechanism of Mg anode and propose effective strategies for enhancing corrosion resistance. Through comprehensive analyses of interfacial structures and equilibrium potentials for Mg(0001)/H2O interface models with different water thicknesses, the Mg(0001)/72 H2O model is identified to be reasonable with −2.17 V vs. standard hydrogen electrode equilibrium potential. In addition, utilizing metadynamics, the free energy barrier for Mg dissolution is calculated to be 0.835 eV, enabling the theoretical determination of anodic polarization curves for pure Mg that aligns well with experimental data. Based on the Mg(0001)/72 H2O model, we further explore the effects of various alloying elements on anodic corrosion resistance, among which Al and Mn alloying elements are found to enhance corrosion resistance of Mg. This study provides valuable atomic‐scale insights into the corrosion mechanism of magnesium alloys, offering theoretical guidance for developing novel corrosion‐resistant Mg alloys.

Mechanistic study of moisture corrosion of FeCr alloys in molten salts by ab-initio molecular dynamics simulations

Molten salts are promising for various energy applications including fuel and solar cells and nuclear energy. These applications face a common challenge: corrosion of structural materials by impurities such as H2O. This work employs ab-initio molecular dynamics simulations to study H2O induced corrosion of FeCr alloys in molten NaF and NaCl salts. H2O is found highly stable in both salts, with infrequent, reversible dissociation into OH− and H+ along with HF or HCl formation. The dissociation tendency correlates positively with the electronegativity and negatively with the size of halogen atoms. Accordingly, H2O reaches the salt/metal interface as a molecule before reacting with metal. Reduction of H+ is found to occur without simultaneous oxidation of specific metal atoms such as Cr, suggesting sequential instead of the commonly proposed concurrent reduction and oxidation. The reduced H atoms prefer to stay at the interface and may re-enter NaF but not NaCl, highlighting the influence of salt chemistry.

Influence of anion structure on Li+ solvation and diffusion in ionic liquid-based electrolyte, and decomposition mechanism of anion

Ionic liquids (ILs) are promising electrolytes for lithium batteries. Particularly, the structure of anions plays a vital role in electrochemical performances. In this work, we chose two asymmetric anions ((fluorosulfonyl)(perfluoroethylsulfonyl)imide (FETI-) and 2,2,2-trifluoro-N-(trifluoromethylsulfonyl)acetamide (TSAC-)) and compared with the symmetric one (bis(trifluoromethanesulfonyl)imide (TFSI-)) to investigate the microcosmic mechanism of the effect of structure on Li+ diffusion through classical molecular dynamics (MD) simulations. We investigated the Li+ solvation structure in electrolytes and found that the asymmetry of anions facilitates the formation of monodentate configuration, but they are not effective in accelerating the Li+ diffusion. On the other hand, we found that the faster Li+ diffusion comes from both the higher diffusion coefficient of anions and the smaller size of Li+ aggregates in electrolytes. In addition, the decomposition mechanism of these anions was investigated by ab initio molecular dynamics (AIMD) simulations. The F directly attached to sulfonyl is unstable, leading to the earlier breaking of the S-F bond compared to S-O and S-N bond, forming early protective components of LiF in solid electrolyte interphase (SEI). This work provides guidance for the design of anions, revealing that the influence of anion structure on Li+ solvation and diffusion in IL-based electrolytes, and decomposition mechanism of anions.

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.

Organic and inorganic sublattice coupling in two-dimensional lead halide perovskites

Two-dimensional layered organic-inorganic halide perovskites have successfully spread to diverse optoelectronic applications. Nevertheless, there remain gaps in our understanding of the interactions between organic and inorganic sublattices that form the foundation of their remarkable properties. Here, we examine these interactions using pump-probe spectroscopy and ab initio molecular dynamics simulations. Unlike off-resonant pumping, resonant excitation of the organic sublattice alters both the electronic and lattice degrees of freedom within the inorganic sublattice, indicating the existence of electronic coupling. Theoretical simulations verify that the reduced bandgap is likely due to the enhanced distortion index of the inorganic octahedra. Further evidence of the mechanical coupling between these two sublattices is revealed through the slow heat transfer process, where the resultant lattice tensile strain launches coherent longitudinal acoustic phonons. Our findings explicate the intimate electronic and mechanical couplings between the organic and inorganic sublattices, crucial for tailoring the optoelectronic properties of two-dimensional halide perovskites.

Ab Initio Study of Novel Phase-Change Heterostructures.

Neuromorphic devices constitute a novel approach to computing that takes inspiration from the brain to unify the processing and storage units. Memories based on phase-change materials (PCMs) are potential candidates for such devices due to their non-volatility and excellent scalability, however their use is hindered by their conductance variability and temporal drift in resistance. Recently, it has been shown that the utilization of phase-change heterostructures consisting of nanolayers of the Sb2Te3 PCM interleaved with a transition-metal dichalcogenide, acting as a confinement material, strongly mitigates these problems. In this work, superlattice heterostructures made of TiTe2 and two prototypical PCMs, respectively GeTe and Ge2Sb2Te5 are considered. By performing ab initio molecular dynamics simulations, it is shown that it is possible to switch the PCMs without destroying the superlattice structure and without diffusion of the atoms of the PCM across the TiTe2 nanolayers. In particular, the model containing Ge2Sb2Te5 shows weak coupling between the two materials during the switching process, which, combined with the high stability of the amorphous state of Ge2Sb2Te5, makes it a very promising candidate for neuromorphic computingapplications.

Theoretical Investigation of Transient Species Following Photodissociation of Ironpentacarbonyl in Ethanol Solution.

Photodissociation of ironpentacarbonyl [1Fe(CO)5] in solution generates transient species in different electronic states, which we studied theoretically. From ab initio molecular dynamics simulations in ethanol solution, the closed-shell parent compound 1Fe(CO)5 is found to interact weakly with the solvent, whereas the irontetracarbonyl [Fe(CO)4] species, formed after photodissociation, has a strongly spin-dependent behavior. It coordinates a solvent molecule tightly in the singlet state [1Fe(CO)4] and weakly in the triplet state [3Fe(CO)4]. From the simulations, we have gained insights into intersystem crossing in solvated irontetracarbonyl based on the distinct structural differences induced by the change in multiplicity. Alternative forms of coordination between 1Fe(CO)4 and functional groups of the ethanol molecule are simulated, and a quantum chemical investigation of the energy landscape for the coordinated irontetracarbonyl gives information about the interconversion of different transient species in solution. Furthermore, insights from the simulations, in which we find evidence of a solvent exchange mechanism, challenge the previously proposed mechanism of chain walking for under-coordinated metal carbonyls in solution.

Computational exploration of two-dimensional vacancy-free boridene sheet and its derivatives: high stabilities and the promise for hydrogen evolution reaction

The recent synthesis of a two-dimensional (2D) MBene sheet, referred to as the boridene sheet (Mo4B6Tz), has ignited considerable interest in exploring 2D transition metal borides. Boridene has an ordered arrangement of metal vacancies, which are pivotal to its stability. Employing first-principles calculations, we explored the stable phases, electronic properties and catalytic abilities of boridene with different vacancy concentrations (Vm). Our results demonstrate that Vm significantly influences the cohesive energies of boridene sheets. Phonon spectrum and ab initio molecular dynamics simulations reveal the high stability of the vacancy-free boridene Mo6B6T6 (T = O, -OH), underscoring their potential for experimental realization. Substituting Mo atoms with Nb, Ta, or W enhances the structural stability of boridene sheets, leading to the identification of four stable variants: Nb6B6F6, Ta6B6F6, Ta6B6O6, and W6B6O6. These boridene sheets exhibit metallic behavior, with five structures displaying near-zero Gibbs free energy for hydrogen atom adsorption, indicating their potential as catalysts for the hydrogen evolution reaction. The uncovering of vacancy-free boridenes and their 2D derivatives greatly broadens the scope of the MBene family.

Unexpected Intermolecular C-H···O Hydrogen Bonds and 1H NMR Chemical Shifts in a Key Linker for Fluorine-18 Labeling of Dimeric Drugs.

The compound 2-{[(trifluoromethyl)sulfonyl]oxy}propane-1,3-diyl bis(4-methylbenzenesulfonate) (TPB) is a crucial intermediate in the synthesis of 18F-radiolabeled cromolyn derivatives. In this work, we combine 1H NMR spectroscopy, X-ray crystallography, ab initio molecular dynamics, and NMR calculations to examine the structure, interactions, and solvation dynamics of the TPB molecule. In CDCl3, the CH2 groups within its glyceryl-derived linker exhibit a single set of proton signals in the 1H NMR measurements. However, when TPB is dissolved in DMSO-d6, distinct splitting patterns emerge despite its seemingly symmetric chemical structure. Crystallographic analysis further unveils the absence of overall symmetry in its three-dimensional arrangement. To elucidate these unique NMR features, we carry out ab initio molecular dynamics simulations and characterize the solvation structures and dynamics of TPB in CHCl3 and DMSO solutions. In contrast to the predominantly nonpolar nature of the CHCl3 solvents, DMSO directly participates in C-H···O hydrogen-bonding interactions with the solute molecule, leading to the splitting of its -CH2 chemical shifts into two distinct distributions. The comprehensive understanding of the structure and solvation interactions of TPB provides essential insights into its application in the radiofluorination reactions of cromolyn derivatives and holds promise for the future development of radiolabeled dimeric drugs.

Ab initio determination of melting and sound velocity of neon up to the deep interior of the Earth.

The thermophysical properties and elemental abundances of the noble gases in terrestrial materials can provide unique insights into the Earth's evolution and mantle dynamics. Here, we perform extensive abinitio molecular dynamics simulations to determine the melting temperature and sound velocity of neon up to 370GPa and 7500K to constrain its physical state and storage capacity, together with to reveal its implications for the deep interior of the Earth. It is found that solid neon can exist stably under the lower mantle and inner core conditions, and the abnormal melting of neon is not observed under the entire temperature (T) and pressure (P) region inside the Earth owing to its peculiar electronic structure, which is substantially distinct from other heavier noble gases. An inspection of the reduction for sound velocity along the Earth's geotherm evidences that neon can be used as a light element to account for the low-velocity anomaly and density deficit in the deep Earth. A comparison of the pair distribution functions and mean square displacements of MgSiO3-Ne and Fe-Ne alloys further reveals that MgSiO3 has a larger neon storage capacity than the liquid iron under the deep Earth condition, indicating that the lower mantle may be a natural deep noble gas storage reservoir. Our results provide valuable information for studying the fundamental behavior and phase transition of neon in a higher T-P regime, and further enhance our understanding for the interior structure and evolution processes inside the Earth.

Computational Insights into Dion-Jacobson Type Oxide Ion Conductors.

Dion-Jacobson type materials have recently emerged as a new structural family of oxide ion conductors, materials important for applications in a variety of electrochemical devices. While some attempts to improve their ionic conductivity have been reported, a detailed understanding of the underlying oxide ion diffusion mechanisms in these materials is still missing. To explore the structure-property relationships leading to the favorable properties, we carried out ab initio molecular dynamics simulations of oxide ion diffusion in CsBi2Ti2NbO10-δ. Our computational study reveals significant out-of-plane dynamics, indicating that the dominant pathway for oxide ion migration is via jumps into and out of the (ab) crystallographic plane. This suggests that further improvement of oxide ion conductivity relative to CsBi2Ti2NbO10-δ could be achieved by enhancing the rotational flexibility of the coordination polyhedra located in the inner perovskite layer, thereby facilitating faster out-of-plane motions.

Hydration of N-Hydroxyurea from Ab Initio Molecular Dynamics Simulations.

N-Hydroxyurea (HU) is an important chemotherapeutic agent used as a first-line treatment in conditions such as sickle cell disease and β-thalassemia, among others. To date, its properties as a hydrated molecule in the blood plasma or cytoplasm are dramatically understudied, although they may be crucial to the binding of HU to the radical catalytic site of ribonucleotide reductase, its molecular target. The purpose of this work is the comprehensive exploration of HU hydration. The topic is studied using ab initio molecular dynamic (AIMD) simulations that apply a first principles representation of the electron density of the system. This allows for the calculation of infrared spectra, which may be decomposed spatially to better capture the spectral signatures of solute-solvent interactions. The studied molecule is found to be strongly hydrated and tightly bound to the first shell water molecules. The analysis of the distance-dependent spectra of HU shows that the E and Z conformers spectrally affect, on average, 3.4 and 2.5 of the closest H2O molecules, respectively, in spheres of radii of 3.7 Å and 3.5 Å, respectively. The distance-dependent spectra corresponding to these cutoff radii show increased absorbance in the red-shifted part of the water OH stretching vibration band, indicating local enhancement of the solvent's hydrogen bond network. The radially resolved IR spectra also demonstrate that HU effortlessly incorporates into the hydrogen bond network of water and has an enhancing effect on this network. Metadynamics simulations based on AIMD methodology provide a picture of the conformational equilibria of HU in solution. Contrary to previous investigations of an isolated HU molecule in the gas phase, the Z conformer of HU is found here to be more stable by 17.4 kJ·mol-1 than the E conformer, pointing at the crucial role that hydration plays in determining the conformational stability of solutes. The potential energy surface for the OH group rotation in HU indicates that there is no intramolecular hydrogen bond in Z-HU in water, in stark contrast to the isolated solute in the gas phase. Instead, the preferred orientation of the hydroxyl group is perpendicular to the molecular plane of the solute. In view of the known chaotropic effect of urea and its N-alkyl-substituted derivatives, N-hydroxyurea emerges as a unique urea derivative that exhibits a kosmotropic ordering of nearby water. This property may be of crucial importance for its binding to the catalytic site of ribonucleotide reductase with a concomitant displacement of a water molecule.

Influence of hydrogen volume/specimen surface area ratio on hydrogen embrittlement sensitivity of X52 pipeline steel

In the transition towards green energy, long distance transportation of gaseous hydrogen by pipeline is considered to be an economically attractive solution. For safe transportation of hydrogen, it is important to conduct relevant mechanical tests in representative hydrogen environments and working conditions of pipelines to demonstrate acceptable hydrogen embrittlement resistance of the steel. One important issue is the testing condition for the ratio of gaseous hydrogen volume over specimen surface area required in mechanical testing since this has not been specified in any testing standards and codes or investigated. This study investigated this issue by performing slow strain rate tensile (SSRT) testing of an X52 pipeline steel in gaseous hydrogen with different ratios of hydrogen volume/specimen surface area, Rv/s, ml/mm2. It was found that, with increasing Rv/s ratio, the area reduction of the steel decreased and hydrogen embrittlement sensitivity increased in the range of the Rv/s ratios investigated. The saturating Rv/s ratio appeared to be ∼2.0 beyond which the effect of Rv/s ratio on hydrogen embrittlement sensitivity was relatively small. The macroscopic and microscopic observations on both the fracture and side surfaces of the SSRT tested specimens also found that the tendency to cleavage fracture increased with increasing Rv/s ratio. The results of the ab initio molecular dynamics (AIMD) simulations performed also predicted a dependence of hydrogen embrittlement on hydrogen volume when Rv/s ratio is relatively small.

Schottky junction-mediated carrier separation in BiOBr/Ti3C2Tx van der Waals heterojunction photocatalyst for increased removal of Cr(VI) under visible-light

The photocatalytic removal of pollutants can integrate efficient solar energy utilization with wastewater treatment. Herein, we report a simple strategy to construct BiOBr/Ti3C2Tx composite photocatalysts consisting of self-assembled van der Waals (vdw) heterojunction structures containing Schottky junctions. The vdW interaction between BiOBr and Ti3C2Tx and the structural stability of the catalysts were validated by density functional theory (DFT) calculations and an ab initio molecular dynamic simulation (AIMD). Compared to BiOBr, the BiOBr/Ti3C2Tx composites exhibited an enhanced light absorption capability and photocurrent response, along with a lower carrier recombination rate and charge-transfer resistance. The BiOBr/Ti3C2Tx activity for the photoreduction of hexavalent chromium (Cr(VI)) under visible-light irradiation considerably surpassed that of the original BiOBr (53.9 %). Using the BTC-0.8 catalyst (with a Ti3C2Tx: BiOBr mass ratio of 0.8 %) optimal reduction efficiency (94.0 %) for Cr(VI) with a reaction rate approximately 3.5 times that of pure BiOBr. The BTC-0.8 catalyst was viable over a wide pH range and maintained an 83.9 % photoreduction rate after four cycling experiments. The exceptional photocatalytic ability of BiOBr/Ti3C2Tx was achieved because Ti3C2Tx served as an electron injection layer that separated photogenerated carriers. The Schottky junction constructed within the BiOBr/Ti3C2Tx vdW heterojunctions suppressed electron backflow, promoting carrier separation. The excellent photocatalytic performance and strong cyclic stability of BiOBr/Ti3C2Tx demonstrated the considerable application potential of this photocatalyst.

Ab initio molecular dynamics study on the local structures and solid/liquid interface in liquid Al-Ti and Al-B-Ti alloys

High-purity aluminum, owing to its excellent physical and chemical properties, finds extensive applications across various fields. A comprehensive understanding of the local structure of liquid Al-based alloys and the characteristics of solute atoms at the forefront of the solid-liquid interface is highly advantageous for the purification of aluminum. Ab initio molecular dynamics (AIMD) simulations were employed to study the local structures, element interactions and electronic structures in Al-Ti and Al-B-Ti melts, together with the behaviors of Ti atom and B-Ti atom pair at the solid/liquid interface forefront of Al. In the Al-B-Ti ternary melt, there are stronger interactions between B and Ti atoms than those between Al and B atoms as well as between Al and Ti atoms. The introduction of B atoms leads to the disruption of Al-Ti bonds and the formation of B-Ti bonds, thereby altering the state of Ti atoms in the melt and then influencing their diffusion properties. The dynamic process of B-Ti bond formation in the melt has been elucidated. The bond angle distribution function shows that the geometric configuration around Ti atoms in the Al melt undergoes significant changes with the addition of B atoms. A strong hybridization occurs between the 3d orbital of Ti atom and the 2p orbital of B atom. Moreover, a notable charge accumulation between B and Ti atoms, indicates the strong interaction between them. The addition of Ti atom can promote the migration of the Al solid/liquid interface. The B-Ti bond is very stable at the forefront of the Al solid/liquid interface, and the formation of B-Ti bond significantly drag the migration of the Al solid/liquid interface.

Theoretically designed M@diaza[2.2.2]cryptand complexes: the role of non-covalent interactions in promoting NLO properties of organic electrides

ABSTRACT Organic excess electron compounds with significant nonlinear optical (NLO) properties are widely employed in optoelectronic applications. Herein, single-alkali metals with diaza[2.2.2] cryptand (M@crypt,M=Li, Na, and K) are investigated for optoelectronic and NLO properties by using the density functional theory. Thermodynamic and kinetic stabilities of present complexes are computed through interaction energy (Eint) and ab-initio molecular dynamic (AIMD) simulations. M@crypt complexes carry excess electrons and mimic molecular electrides. Quantum theory of atoms in molecules (QTAIM) analysis and reduced density gradient (RDG) spectra demonstrate the roles of the weak van der Waals (vdW) interactions between metal and complexant. The remarkable hyperpolarizability (βo) value up to 1.41 × 106 au may be credited to the presence of loosely bound excess electrons. The hyper Rayleigh scattering hyperpolarizability (βHRS) is recorded up to 1.31 × 106 au for the K@crypt. Furthermore, frequency-dependent first-order and second-order hyperpolarizability is more prominent at the applied frequency of ω = 0.042823 au. The electron localizing function (ELF) and localized orbital locator (LOL) analysis further disclose the nature of interaction between alkali metal and complexant. The TD-DFT method is adopted to get excited state parameters and absorbance properties. An electron density difference map (EDDM) is exploited to evaluate the orbital contributions in excited states. Hence, the studied electride may become a promising candidate for NLO materials. We anticipate that the present work will provide insight into further development of molecular electride for optoelectronic applications.

Intermolecular Correlations in Liquid Lactic Acid Based on ab initio Molecular Dynamics Simulations Combined with High-Energy X-ray Diffraction Measurements

Intermolecular Correlations in Liquid Lactic Acid Based on ab initio Molecular Dynamics Simulations Combined with High-Energy X-ray Diffraction Measurements

Novel arsenborane as a sensing material for dichlorofluoromethane and dichlorodifluoromethane vapours – A first-principles study

Presently, two-dimensional (2D) materials are cutting-edge research to detect air pollutants in the environment. In this research study, we used arsenborane as a sensing element towards halomethanes viz. dichlorofluoromethane (Freon 21) and dichlorodifluoromethane (Freon 12). At first, the structural firmness of bare arsenborane is confirmed based on formation energy, phonon band spectrum, and ab initio molecular dynamics simulation. Also, the electronic property of bare stable arsenborane is investigated by plotting bands and projected-density-of-states (PDOS) maps. The energy gap of arsenborane is computed to be 0.927 eV which confirms semiconductor nature. Significantly, the interaction behaviour of Freon 21 and Freon 12 on arsenborane is studied with the influence of relative band gap changes, Mulliken-charge analysis, and adsorption energy. In the present work, the calculated value of adsorption energy is observed in the scope of −0.083 eV to −0.618 eV (physisorption regime). Thus, based on the outcome we suggest that the new 2D material arsenborane is a good base material to detect Freon 21 and Freon 12 vapours in the air environment.

Insight into the effects of gas molecules-adsorbed on 2D-FeS2: A DFT study

Lithium-sulfur (Li-S) batteries are especially competitive in the energy sector due to their excellent performances, like preferable energy density and economic benefits. Studying the adsorption of gas molecules on electrode materials has potential engineering significance for Li-S batteries since they have a highly osmotic potential, which causes unavoidable damage to batteries. In this work, the adsorption phenomenon of common gas molecules (H2O, N2, H2, CO2, and O2) on the two-dimensional pyrite (2D-FeS2) cathode material surface, as well as the effects on the electronic and electrochemical properties, were investigated by the first-principles calculations. The adsorption capabilities were estimated by adsorption energy and Mulliken population analysis. Simulation results demonstrated that whole adsorption energies were less than -1.0 eV and larger than -0.6 eV, which shows a physisorption nature. Among them, the O-S bond of O2/2D-FeS2 has the strongest strength. Electronic structure calculations suggested that 2D-FeS2 maintained good conductivity after gas molecules were adsorbed, achieving efficient transfer between electron, lithium, and sulfur intermediates. Additionally, ab initio molecular dynamics (AIMD) simulations showed that Li+ exhibits excellent diffusion performance and low activation energy at different temperatures. 2D-FeS2 still has a stable electrochemical working window (1.87 ∼ 2.47 V), while the theoretical open current voltage is damaged by gas molecule adsorption. Consequently, this work theoretically reveals the effect of gas molecules on the cathode materials for Li-S batteries, which has guide meaning for engineering.

Chloride Ion-Induced Spatial Separation and Long Recombination Time of Photogenerated Electrons and Holes in Crystalline Carbon Nitride.

Poly(triazine imide)·Li+Cl- (PTI/Li+Cl-) as one of the most reported crystalline carbon nitrides has shown exciting potential for photocatalysis. However, understanding the role of Li+/Cl- in the photoexcited charge transfer in the time and space of PTI is a challenging problem. Here, we have investigated the nonradiative charge recombination of series ion intercalated PTI systems (PTI/Li+X-, where X = F, Cl, Br, and I) by performing the ab initio nonadiabatic molecular dynamics simulations. The results indicate that the intercalated anions in PTI/Li+Cl- and PTI/Li+Br- have the potential to trap holes, separate the electrons and holes, and prolong the nonradiative electron-hole recombination. In particular, ∼70% of holes in PTI/Li+Cl- can transport among interlayers toward the {0001} planes, while most of the electrons are transferred to the {101̅0} planes, exhibiting different transport pathways and directions. Furthermore, PTI/Li+Cl- has an electron-hole recombination time as long as 136 ns, which explains its excellent optoelectronic properties. This work provides a theoretical baseline for the reported facet engineering improvement of crystalline carbon nitride for overall water splitting.