skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on September 3, 2025

Title: Breaking the size limitation of nonadiabatic molecular dynamics in condensed matter systems with local descriptor machine learning
Nonadiabatic molecular dynamics (NA-MD) is a powerful tool to model far-from-equilibrium processes, such as photochemical reactions and charge transport. NA-MD application to condensed phase has drawn tremendous attention recently for development of next-generation energy and optoelectronic materials. Studies of condensed matter allow one to employ efficient computational tools, such as density functional theory (DFT) and classical path approximation (CPA). Still, system size and simulation timescale are strongly limited by costly ab initio calculations of electronic energies, forces, and NA couplings. We resolve the limitations by developing a fully machine learning (ML) approach in which all the above properties are obtained using neural networks based on local descriptors. The ML models correlate the target properties for NA-MD, implemented with DFT and CPA, directly to the system structure. Trained on small systems, the neural networks are applied to large systems and long timescales, extending NA-MD capabilities by orders of magnitude. We demonstrate the approach with dependence of charge trapping and recombination on defect concentration in MoS2. Defects provide the main mechanism of charge losses, resulting in performance degradation. Charge trapping slows with decreasing defect concentration; however, recombination exhibits complex dependence, conditional on whether it occurs between free or trapped charges, and relative concentrations of carriers and defects. Delocalized shallow traps can become localized with increasing temperature, changing trapping and recombination behavior. Completely based on ML, the approach bridges the gap between theoretical models and realistic experimental conditions and enables NA-MD on thousand-atom systems and many nanoseconds.  more » « less
Award ID(s):
2154367
PAR ID:
10583509
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
PNAS
Date Published:
Journal Name:
Proceedings of the National Academy of Sciences
Volume:
121
Issue:
36
ISSN:
0027-8424
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Defects, such as halide interstitials, act as charge recombination centers, induce degradation of halide perovskites, and create major obstacles to applications of these materials. Alkali metal dopants greatly improve perovskite performance. Using ab initio nonadiabatic molecular dynamics, it is demonstrated that alkalis bring favorable effects. The formation energy of halide interstitials increases by up to a factor of four in the presence of alkali dopants, and therefore, defect concentration decreases. When defects are present, alkali metals strongly bind to them. Halide interstitials introduce mid‐gap states that rapidly trap charge carriers. Alkalis eliminate the trap states, helping to maintain high current density. Further to charge trapping, the interstitials accelerate charge recombination. By passivating the interstitials, alkalis make carrier lifetimes up to seven times longer than in defect‐free perovskites and up to thirty times longer than in defective perovskites. 
    more » « less
  2. Substitutionally doped transition metal dichalcogenides (TMDs) are essential for advancing TMD‐based field effect transistors, sensors, and quantum photonic devices. However, the impact of local dopant concentrations and dopant–dopant interactions on charge doping and defect formation within TMDs remains underexplored. Here, a breakthrough understanding of the influence of rhenium (Re) concentration is presented on charge doping and defect formation in MoS2monolayers grown by metal–organic chemical vapor deposition (MOCVD). It is shown that Re‐MoS2films exhibit reduced sulfur‐site defects, consistent with prior reports. However, as the Re concentration approaches ⪆2 atom%, significant clustering of Re in the MoS2is observed. Ab Initio calculations indicate that the transition from isolated Re atoms to Re clusters increases the ionization energy of Re dopants, thereby reducing Re‐doping efficacy. Using photoluminescence (PL) spectroscopy, it is shown that Re dopant clustering creates defect states that trap photogenerated excitons within the MoS2lattice, resulting in broad sub‐gap emission. These results provide critical insights into how the local concentration of metal dopants influences carrier density, defect formation, and exciton recombination in TMDs, offering a novel framework for designing future TMD‐based devices with improved electronic and photonic properties. 
    more » « less
  3. Abstract The ground state electron density — obtainable using Kohn-Sham Density Functional Theory (KS-DFT) simulations — contains a wealth of material information, making its prediction via machine learning (ML) models attractive. However, the computational expense of KS-DFT scales cubically with system size which tends to stymie training data generation, making it difficult to develop quantifiably accurate ML models that are applicable across many scales and system configurations. Here, we address this fundamental challenge by employing transfer learning to leverage the multi-scale nature of the training data, while comprehensively sampling system configurations using thermalization. Our ML models are less reliant on heuristics, and being based on Bayesian neural networks, enable uncertainty quantification. We show that our models incur significantly lower data generation costs while allowing confident — and when verifiable, accurate — predictions for a wide variety of bulk systems well beyond training, including systems with defects, different alloy compositions, and at multi-million-atom scales. Moreover, such predictions can be carried out using only modest computational resources. 
    more » « less
  4. Non-adiabatic (NA) molecular dynamics (MD) is a powerful approach for studying far-from-equilibrium quantum dynamics in photophysical and photochemical systems. Most NA-MD methods are developed and tested with few-state models, and their validity with complex systems involving many states is not well studied. By modeling intraband equilibration and interband recombination of charge carriers in MoS2, we investigate the convergence of three popular NA-MD algorithms, fewest switches surface hopping (FSSH), global flux surface hopping (GFSH), and decoherence induced surface hopping (DISH) with the number of states. Only the standard DISH algorithm converges with the number of states and produces Boltzmann equilibrium. Unitary propagation of the wave function in FSSH and GFSH violates the Boltzmann distribution, leads to internal inconsistency between time-dependent Schrödinger equation state populations and trajectory counts, and produces non-convergent results. Introducing decoherence in FSSH and GFSH by collapsing the wave function fixes these problems. The simplified version of DISH that omits projecting out the occupied state and is applicable to few-state systems also causes problems when the number of states is increased. We discuss the algorithmic application of wave function collapse and Boltzmann detailed balance and provide detailed FSSH, GFSH, and DISH flow charts. The use of convergent NA-MD methods is highly important for modeling complicated quantum processes involving multiple states. Our findings provide the basis for investigating quantum dynamics in realistic complex systems. 
    more » « less
  5. Abstract Although all‐inorganic metal halide perovskites (MHPs) have shown tremendous improvement, they are still inferior to the hybrid organic–inorganic MHPs in efficiency. Recently, a conceptually new β‐CsPbI3perovskite reached 18.4 % efficiency combined with good thermodynamic stability at ambient conditions. We use ab initio non‐adiabatic molecular dynamics to show that native point defects in β‐CsPbI3are generally benign for nonradiative charge recombination, regardless of whether they introduce shallow or deep trap states. These results indicate that MHPs do not follow the simple models used to explain defect‐mediated charge recombination in the conventional semiconductors. The strong tolerance is due to the softness of the perovskite lattice, which permits separation of electrons and holes upon defect formation, and only allows carriers to couple to the low‐frequency vibrations. Both factors decrease notably the non‐adiabatic coupling and slow down the dissipation of energy to heat. 
    more » « less