Search NASA⌕ Search

SEARCH · Search NASA

Results for “Nonadiabatic Molecular Dynamics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Ultrafast charge transfer coupled to quantum proton motion at molecule/metal oxide interface

Understanding how the nuclear quantum effects (NQEs) in the hydrogen bond (H-bond) network influence the photoexcited charge transfer at semiconductor/molecule interface is a challenging problem. By combining two kinds of emerging molecular dynamics methods at the ab initio level, the path integral–based molecular dynamics and time-dependent nonadiabatic molecular dynamics, and choosing CH 3 OH/TiO 2 as a prototypical system to study, we find that the quantum proton motion in the H-bond network is strongly coupled with the ultrafast photoexcited charge dynamics at the interface. The hole trapping ability of the adsorbed methanol molecule is notably enhanced by the NQEs, and thus, it behaves as a hole scavenger on titanium dioxide. The critical role of the H-bond network is confirmed by in situ scanning tunneling microscope measurements with ultraviolet light illumination. It is concluded the quantum proton motion in the H-bond network plays a critical role in influencing the energy conversion efficiency based on photoexcitation.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Nonadiabatic Excited-State Molecular Dynamics for Open-Shell Systems

Nonadiabatic Molecular Dynamics (NAMD) of excited states has been widely used in the simulation of photoinduced phenomena. However, the inability to treat bond breaking and forming processes with single-reference electronic structure methods limits their application in photochemistry for extended molecular systems. In this work, the extension of excited-state NAMD for open-shell systems is developed and implemented in the NEXMD software. We present the spin-unrestricted CIS and TDSCF formalism for the ground and excited states, analytical derivatives, and nonadiabatic derivative couplings for the respective potential energy surfaces. Next, this methodology is employed to study the photochemical reaction of three model molecules. The results demonstrate the advantage of the open-shell approach in modeling photochemical reactions, especially involving bond breaking processes. We find that the open-shell method lowers the reaction barrier at the bond-breaking limits resulting in larger calculated photochemical quantum yields compared to the respective closed-shell results. We also address problems related to spin contamination in the open-shell method, especially when molecular geometries are far from equilibrium.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Non-Adiabatic Excited State Molecular Dynamics Methodologies: comparison and convergence

Direct atomistic simulation of nonadiabatic molecular dynamics is a challenging goal that allows important insights into fundamental physical phenomena. A variety of frameworks, ranging from fully quantum treatment of nuclei to semiclassical and mixed quantum–classical approaches, were developed. These algorithms are then coupled to specific electronic structure techniques. Such diversity and lack of standardized implementation make it difficult to compare the performance of different methodologies when treating realistic systems. Here, we compare three popular methods for large chromophores: Ehrenfest, surface hopping, and multiconfigurational Ehrenfest with ab initio multiple cloning (MCE-AIMC). These approaches are implemented in the NEXMD software, which features a common computational chemistry model. The resulting comparisons reveal the method performance for population relaxation and coherent vibronic dynamics. Finally, we study the numerical convergence of MCE-AIMC algorithms by considering the number of trajectories, cloning thresholds, and Gaussian wavepacket width. Our results provide helpful reference data for selecting an optimal methodology for simulating excited-state molecular dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

An accurate, non-empirical method for incorporating decoherence into Ehrenfest dynamics

In mixed quantum-classical nonadiabatic molecular dynamics methods, the anchoring of the electronic wave function to a single nuclear geometry results in both quantitative and qualitative errors in the dynamics. In the context of both Ehrenfest and trajectory surface hopping methods, methods for incorporating decoherence are widely used to eliminate these errors. However, the accuracy of these methods often depends strongly on the parameterization of the decoherence time and/or other related quantities. Here, we present a refinement of the recently introduced collapse to a block (TAB) scheme for incorporating decoherence into Ehrenfest dynamics. The proposed approach incorporates an approximation to the history of the population dynamics and treats the coherence decay as Gaussian, rather than exponential. This method uses parameters that can be obtained from first principles, rather than empirical fitting. Application to one-dimensional models indicates excellent agreement with numerically exact simulations. We also introduce a second refinement to the TAB method: a robust linear least-squares algorithm for determining collapse probabilities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Bidentate Lewis bases are preferred for passivation of MAPbI 3 surfaces: A time-domain ab initio analysis

We report chemical passivation is a predominant approach to inhibit intrinsic defects responsible for electron–hole recombination in perovskite solar cells. Using time-domain density functional theory combined with nonadiabatic molecular dynamics, we demonstrate defect passivation by separation of electrons and holes in metal halide perovskites (MHPs), and show that bidentate ligands exhibit the best performance. Defects in traditional semiconductors create deep midgap states, and passivation eliminates these states. In contrast, common defects produce no deep midgap states in MHPs. Instead, defects localize electrons and holes around defect sites, enhancing electron–hole interaction. Defect passivation in MHPs acts to separate charges, decreasing electron–hole and charge–phonon interactions, and increasing charge lifetimes. Bidentate ligands work best, because they can passivate both unsaturated chemical bonds created due to vacancy defects. Bidentate ligands with spatially separated binding sites are preferred, since they provide better match to the sparse inorganic lattice of MHPs. Similar to the Lewis base ligands, water also acts as a ligand and extends charge lifetimes, with the oxygen atom donating its lone electron pair to the defect site. However, water accelerates chemical degradation of MHPs. The binding energy of most ligands is larger than that of water, and therefore, the ligands displace water, increase MHP stability and prolong carrier lifetimes. The established physical mechanism of defect passivation, and the specific principles guiding the choice of passivating molecules advance our understanding of the exceptional properties of MHPs and suggest routes for further improvement of MHP performance in solar energy and optoelectronic applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Real-time GW -BSE investigations on spin-valley exciton dynamics in monolayer transition metal dichalcogenide

We develop an ab initio nonadiabatic molecular dynamics (NAMD) method based on GW plus real-time Bethe-Salpeter equation (GW + rtBSE-NAMD) for the spin-resolved exciton dynamics. From investigations on MoS 2 , we provide a comprehensive picture of spin-valley exciton dynamics where the electron-phonon (e-ph) scattering, spin-orbit interaction (SOI), and electron-hole (e-h) interactions come into play collectively. In particular, we provide a direct evidence that e-h exchange interaction plays a dominant role in the fast valley depolarization within a few picoseconds, which is in excellent agreement with experiments. Moreover, there are bright-to-dark exciton transitions induced by e-ph scattering and SOI. Our study proves that e-h many-body effects are essential to understand the spin-valley exciton dynamics in transition metal dichalcogenides and the newly developed GW + rtBSE-NAMD method provides a powerful tool for exciton dynamics in extended systems with time, space, momentum, energy, and spin resolution.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dynamic signatures of electronically nonadiabatic coupling in sodium hydride: a rigorous test for the symmetric quasi-classical model applied to realistic, ab initio electronic states in the adiabatic representation

Sodium hydride (NaH) in the gas phase presents a seemingly simple electronic structure making it a potentially tractable system for the detailed investigation of nonadiabatic molecular dynamics from both computational and experimental standpoints. The single vibrational degree of freedom, as well as the strong nonadiabatic coupling that arises from the excited electronic states taking on considerable ionic character, provides a realistic chemical system to test the accuracy of quasi-classical methods to model population dynamics where the results are directly comparable against quantum mechanical benchmarks. Here, using a simulated pump–probe type experiment, this work presents computational predictions of population transfer through the avoided crossings of NaH via symmetric quasi-classical Meyer–Miller (SQC/MM), Ehrenfest, and exact quantum dynamics on realistic, ab initio potential energy surfaces. The main driving force for population transfer arises from the ground vibrational level of the D 1 Σ + adiabatic state that is embedded in the manifold of near-dissociation C 1 Σ + vibrational states. When coupled through a sharply localized first-order derivative coupling most of the population transfers between t = 15 and t = 30 fs depending on the initially excited vibronic wavepacket. While quantum mechanical effects are expected due to the reduced mass of NaH, predictions of the population dynamics from both the SQC/MM and Ehrenfest models perform remarkably well against the quantum dynamics benchmark. Additionally, an analysis of the vibronic structure in the nonadiabatically coupled regime is presented using a variational eigensolver methodology.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Mitigating Band Tailing in Kesterite Solar Absorbers: Ab Initio Quantum Dynamics

Open-circuit voltage deficits are limiting factors in kesterite solar cells. Addressing this issue by suppressing band tailing and nonradiative charge recombination is essential for enhancing the performance. We employ ab initio nonadiabatic molecular dynamics to elucidate the origin of band tailing and charge losses and propose a mitigation strategy. The simulations show that Cu–Zn disorder, associated with antisite defect clusters [Cu Zn +Zn Cu ], is a significant source of band tailing in kesterites, as evidenced by the much larger Urbach energy in disordered than ordered kesterites. Cu–Zn disorder gives rise to new sulfur-centered coordination polyhedra, increases structural inhomogeneity, changes electrostatic potential at sulfur centers, and shifts the S(3p) orbital energy. Differences in the S(3p)/Cu(3d) and S(3p)/Sn(5s) hybridization strengths and the S(3p) orbital energy shift reduce the band gap by 0.37 eV. Furthermore, Cu–Zn disorder enhances vibrational motion of sulfur anions and surrounding cations, increasing band gap fluctuations by 15 meV. The stronger electron–phonon interactions reduce charge carrier lifetimes and limit the kesterite solar cell efficiency. Partial substitution of Zn with Cd facilitates structural ordering and significantly suppresses band tailing, particularly in disordered systems. The improvement can be attributed to the larger atomic radius and mass of Cd, which weakens bonding around the anion, suppresses S-related vibrations within the covalent tetrahedra, and reduces nonadiabatic coupling, thereby increasing charge carrier lifetimes. The reported results establish the key influence of cation disorder on band tailing and reduced charge carrier lifetimes in kesterites and highlight cation disorder engineering as a strategy to achieve high-efficiency kesterite solar cells.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Elucidation of Ultrafast Decay, Vibrational Beating, and Slow Decay Processes of Excited Azulene

Azulene’s nonradiative decay dynamics and kinetics from its singlet excited-states were studied using mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) combined with trajectory surface-hopping nonadiabatic molecular dynamics (NAMD) and Fermi’s golden rule (FGR) rate theory. The NAMD dynamics reproduce experimental observations that the S 1 → S 0 decay is accelerated and exhibits a crossover from mono- to biexponential kinetics with increasing excess vibrational energy. Minimum-energy-path analyses reveal a continuous S 1 /S 0 conical-intersection seam slightly above the S 1 minimum, providing readily accessible funnels. Time-resolved normal-mode projections reveal selective energy funneling into C–C stretching modes at 1272, 1528, and 1692 cm –1 . Constructive combinations of the latter two modes appear to promote rapid access to the conical-intersection seam, whereas their beating at ∼ 160 fs imposes a natural limit on the decrease of the decay time upon a further increase of excess energies, suggesting that their interference delays the decay. Here, the FGR rate calculation data for the S 1 → S 0 transition reaffirm that the nonradiative decay proceeds mainly near or through the conical intersection, rather than via simple nonadiabatic derivative coupling around the minimum of S 1 . On the other hand, the FGR rates indicate that S 2 decays predominantly to S 0 , as the S 2 –S 1 vibronic coupling is exceptionally weak, which serves as a primary cause for azulene’s characteristic anti-Kasha emission.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Understanding and Controlling Photoexcited Molecules in Complex Environments (Final Technical Report)

Predicting the outcome of a photochemical reaction is complicated by factors outside the scope of tradition theoretical chemistry. On such small scales, fluctuations of the environment render outcomes stochastic, nonadiabatic dynamics blur the separation between nuclear and electronic motion, and ultrafast relaxation renders typical assumptions of equilibrium inappropriate. This project aimed to elucidate and control nonadiabatic molecular dynamics in condensed phases. We developed accurate, efficient simulation tools to study irreversible molecular processes computational across a range of scales. We applied these methodological advances to study electron transfer under exotic conditions where spin and topology opened new pathways for reactivity. We also brought these ideas to bare on light harvesting systems, where energy and charge transport are coupled, and where long lifetimes can be deleterious to function. These advances helped to establish general principles for chemical efficiency and guide the design of nanoscale energy materials, molecular machines, and related technologies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Charge carrier dynamics in two-dimensional hybrid perovskites: Dion–Jacobson vs. Ruddlesden–Popper phases

Two-dimensional (2D) halide perovskites are promising materials for environmentally stable next-generation optoelectronic device applications. Besides the widely investigated Ruddlesden–Popper (RP) phases, Dion–Jacobson (DJ) phases are attracting considerable attention due to their rapid emergence as efficient solar cell materials. However, there is very little atomistic understanding of the charge carrier dynamics under ambient conditions for these DJ-phases, limiting the possibilities to tune their optoelectronic performances through compositional engineering routes. In this study, by combining nonadiabatic molecular dynamics with time-domain density functional theory methods at room temperature, we compare the dominant non-radiative carrier recombination and dephasing processes in RP and DJ monolayered lead halide perovskites. Our systematic study demonstrates that performance-limiting nonradiative carrier recombination processes greatly depends on the electron–phonon interactions induced by structural fluctuations and instantaneous charge localization in these materials. The stiffer interlayer packing due to the presence of single spacer dications, which separates the lead iodide slabs, reduces the thermal fluctuations in the DJ-phase to a greater extent than that in the RP-phase 2D-perovskites. Specific electronic coupling between the closely spaced lead iodide layers enhances the delocalization of band-edge charge densities in DJ-phase systems. Compared to the RP-phase, reduced inelastic electron–phonon scattering in DJ-phase perovskites significantly limits intrinsic non-radiative recombination processes. The consequent enhancement in the photogenerated charge carrier lifetime makes DJ-phase perovskites potentially suitable for various optoelectronic devices. The computational insights gained from this study allow us to outline a set of robust design principles for DJ-phase perovskites to strategically tune their optoelectronic properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Simultaneous observation of nuclear and electronic dynamics by ultrafast electron diffraction

Simultaneous observation of nuclear and electronic motion is crucial for a complete understanding of molecular dynamics in excited electronic states. It is challenging for a single experiment to independently follow both electronic and nuclear dynamics at the same time. In this paper we show that ultrafast electron diffraction can be used to simultaneously record both electronic and nuclear dynamics in isolated pyridine molecules, naturally disentangling the two components. Electronic state changes (S 1 →S 0 internal conversion) were reflected by a strong transient signal in small-angle inelastic scattering, and nuclear structural changes (ring puckering) were monitored by large-angle elastic diffraction. Supported by ab initio nonadiabatic molecular dynamics and diffraction simulations, our experiment provides a clear view of the interplay between electronic and nuclear dynamics of the photoexcited pyridine molecule.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Excitons in Hematite Fe 2 O 3 : Short-Time Dynamics from TD-DFT and Non-Adiabatic Dynamics Theories

We present a first-principles study of the short-time dynamics of excitons in hematite Fe 2 O 3 . We used time-dependent density functional theory (TD-DFT) with an underlying DFT+U treatment of electron interactions to characterize the electronic structure of excitons and nonadiabatic molecular dynamics theory (NA-MD) to determine their recombination (electronic ground-state recovery) and relaxation dynamics. Decoherence-corrected trajectory surface hopping approaches in NA-MD simulations yielded recovery times of ~1.1 to 1.8 ns and “higher-lying” exciton relaxation times of ~60 to 70 fs, in accord with experimentally derived lifetimes. With hematite phonons in the range of ~100 to 700 cm –1 , higher-lying excitons relax within one or two oscillations of the phonons before getting trapped into an electron–hole pair Exc-3 structure on the first excited state potential energy surface. This structure resembles already a pair of polarons (electron polaron plus hole polaron) with associated lattice distortions three (3) basal planes away. On longer time scales, the electron–hole bipolaronic pair hops to structures Exc-5, then Exc-7, then Exc-9, ... with the electron polaron and hole polaron separated by 5, 7, 9, ... basal planes in a process of charge separation. The largest frequency phonon ~672 cm –1 for the Exc-3 exciton structure is associated with the electron polaron moiety of the exciton. This phonon is a good candidate for giving rise to the recently observed and reported postexcitation transient IR absorption peak.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Charge carrier nonadiabatic dynamics in non-metal doped graphitic carbon nitride

Graphitic carbon nitride (GCN) has attracted significant attention due to its excellent performance in photocatalytic applications. Non-metal doping of GCN has been widely used to improve the efficiency of the material as a photocatalyst. Using a combination of time-domain density functional theory with nonadiabatic molecular dynamics, we study the charge carrier dynamics in oxygen and boron doped GCN systems. The reported simulations provide a detailed time-domain mechanistic description of the charge separation and recombination processes that are of fundamental importance while evaluating the photovoltaic and photocatalytic performance of the material. The appearance of smaller energy gaps due to the presence of dopant states improves the visible light absorption range of the doped systems. At the same time, the nonradiative lifetimes are shortened in the doped systems as compared to the pristine GCN. In the case of boron doped at a carbon (B–C–GCN), the charge recombination time is very long as compared to the other two doped systems owing to the smaller electron–phonon coupling strength between the valence band maximum and the trap state. The results suggest B–C–GCN as the most suitable candidate among three doped systems studied in this work for applications in photocatalysis. Finally, this work sheds light into the influence of dopants on quantum dynamics processes that govern GCN performance and, thus, guides toward building high-performance devices in photocatalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Tuning charge carrier dynamics through spacer cation functionalization in layered halide perovskites: an ab initio quantum dynamics study

Dion–Jacobson (DJ) phase two-dimensional (2D) hybrid halide perovskites are promising for environmentally stable optoelectronic device applications due to their attractive photophysical properties, including long charge carrier lifetime and partially suppressed non-radiative losses. However, the atomistic details of the structure–property relationship are severely limited, which substantially restricts the strategic materials designed for these layered halide perovskites (LHPs). Here, we combine nonadiabatic molecular dynamics and time-domain density functional theory to understand the effects of spacer cation functionalization on the ground and excited-state charge carrier dynamics in DJ phase LHPs. Our in-depth study reveals that the fluorination of spacer cations considerably restricts the thermal motions of LHPs at ambient conditions. The compact structure weakens the overall electron–phonon interactions and reduces the thermal fluctuations of the band edges in real-time. These dynamic modifications partially mitigate the non-radiative recombination, prolonging the lifetime of photogenerated charge carriers. The suppressed carrier loss mechanism strongly suggests that the fluorinated spacer cation-based LHPs would exhibit enhanced performance as optoelectronic materials. These systematic simulations of excited state carrier dynamics elaborate that chemically viable functionalization of the spacer cations is a robust approach to improve the photophysical properties of 2D halide perovskites strategically. These insights also guide us to propose a few potential strategies to design highly beneficial spacer cations of LHPs that can be introduced in next-generation optoelectronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Ligand Controls Excited Charge Carrier Dynamics in Metal-Rich CdSe Quantum Dots: Computational Insights

Small metal-rich semiconducting quantum dots (QDs) are promising for solid-state lighting and single-photon emission due to their highly tunable yet narrow emission line widths. Nonetheless, the anionic ligands commonly employed to passivate these QDs exert a substantial influence on the optoelectronic characteristics, primarily owing to strong electron–phonon interactions. In this work, we combine time-domain density functional theory and nonadiabatic molecular dynamics to investigate the excited charge carrier dynamics of Cd 28 Se 17 X 22 QDs (X = HCOO – , OH – , Cl – , and SH – ) at ambient conditions. These chemically distinct but regularly used molecular groups influence the dynamic surface-ligand interfacial interactions in Cd-rich QDs, drastically modifying their vibrational characteristics. The strong electron–phonon coupling leads to substantial transient variations at the band edge states. The strength of these interactions closely depends on the physicochemical characteristics of passivating ligands. Consequently, the ligands largely control the nonradiative recombination rates and emission characteristics in these QDs. Our simulations indicate that Cd 28 Se 17 (OH) 22 has the fastest nonradiative recombination rate due to the strongest electron–phonon interactions. Conversely, QDs passivated with thiolate or chloride exhibit considerably longer carrier lifetimes and suppressed nonradiative processes. The ligand-controlled electron–phonon interactions further give rise to the broadest and narrowest intrinsic optical line widths for OH and Cl-passivated single QDs, respectively. Finally, obtained computational insights lay the groundwork for designing appropriate passivating ligands on metal-rich QDs, making them suitable for a wide range of applications, from blue LEDs to quantum emitters.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Low-frequency lattice phonons in halide perovskites explain high defect tolerance toward electron-hole recombination

Low-cost solution-based synthesis of metal halide perovskites (MHPs) invariably introduces defects in the system, which could form Shockley-Read-Hall (SRH) electron-hole recombination centers detrimental to solar conversion efficiency. Here, we investigate the nonradiative recombination processes due to native point defects in methylammonium lead halide (MAPbI 3 ) perovskites using ab initio nonadiabatic molecular dynamics within surface-hopping framework. Regardless of whether the defects introduce a shallow or deep band state, we find that charge recombination in MAPbI 3 is not enhanced, contrary to predictions from SRH theory. We demonstrate that this strong tolerance against defects, and hence the breakdown of SRH, arises because the photogenerated carriers are only coupled with low-frequency phonons and electron and hole states overlap weakly. Both factors appreciably decrease the nonadiabatic coupling. We argue that the soft nature of the inorganic lattice with small bulk modulus is key for defect tolerance, and hence, the findings are general to other MHPs.

Science & Technology - Other Topics↗
DOE OSTI2022

Multiscale Strategy for Predicting Radiation Chemistry in Polymers

A primary mode for radiation damage in polymers arises from ballistic electrons that induce electronic excitations, yet subsequent chemical mechanisms are poorly understood. We develop a multiscale strategy to predict this chemistry starting from subatomic scattering calculations. Nonadiabatic molecular dynamics simulations sample initial bond-breaking events following the most likely excitations, which feed into semiempirical simulations that approach chemical equilibrium. Application to polyethylene reveals a mechanism explaining the low propensity to cross-link in crystalline samples.

36 MATERIALS SCIENCE↗