Search NASA⌕ Search

SEARCH · Search NASA

Results for “surface hopping”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Stable excited dication: trapping on the S 1 state of formaldehyde dication after strong field ionization

Combined theoretical and experimental work examines the dynamics of dication formaldehyde produced by strong field ionization. Trajectory surface hopping dynamics on the first several singlet electronic states of the formaldehyde dication are used to examine the relaxation pathways and dissociation channels, while kinetic energy distributions after strong field ionization of formaldehyde and deuterated formaldehyde are used to confirm the theoretical predictions. Here, we find that the first excited state of the formaldehyde dication is stable, neither decays to the ground state nor dissociates, even though the ground state and higher lying states are directly dissociative. The stability of the first excited state is explained by its symmetry which does not allow for radiative or nonradiative transitions to the ground state and by large barriers to dissociate on the excited state surface.

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

A Theory of Ultrafast Charge Transfer Relaxation with Non-Innocent Solvent Molecules

In this work we revisit the photodynamics of tetracyanoethylene-hexamethylbenzene (TCNE- HMB), the molecular complex studied by Hochstrasser et al. [J. Chem. Phys. 100, 4797–4810, 1994] that has long challenged the applicability of Marcus’s theory of elec- tron transfer for predicting photochemical reactions. Using a novel black-box electronic structure algorithm (time-dependent density functional theory with one double, TD- DFT-1D) to efficiently run molecular dynamics that can treat charge recombination, we run ab initio surface hopping molecular dynamics and confirm that, for a polar solvent, charge recombination rates can be incredibly fast (indeed faster than the sol- vent relaxation time); for non-polar solvents, the rate is much slower. We demonstrate that, although Marcus theory cannot be directly applied, these nonequilibrium (and sometimes incredibly fast) photoexcited dynamics can be effectively explained within a two-state model without any evidence of a transition through a conical intersection. Most importantly, for this paradigmatic model system, we are able to identify two nuclear coordinates of interest (rather than the single coordinate predicted by Marcus or a full set of internal quantum modes studied by Bixon and Jortner): the solvent relaxation in the first shell (that strongly modulates the energies of the charge trans- fer state and differentiates time scales for relaxation) and a nuclear displacement in the TCNE-HMB complex arising from a handful of vibrations that induces non-Born Oppenheimer motion and eventually facilitates an abrupt electronic transition to the ground state. Altogether, these findings suggest a tractable generalization of Marcus theory for future simulations of photochemistry with non-innocent solvent environ- ments in the spirit of a Hamiltonian suggested by Stuchebrukhov (J. Chem. Phys. 107, 3821, 1997).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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

Nonadiabatic Molecular Dynamics Study of the Relaxation Pathways of Photoexcited Cyclooctatetraene

In the current study, we present nonadiabatic (NAMD) and adiabatic molecular dynamics simulations of the transition-state dynamics of photoexcited cyclooctatetraene (COT). The equilibrium-state structure and absorption spectra are analyzed using the semiempirical Austin Model 1 potential. The NAMD simulations are obtained by a surface-hopping algorithm. We analyzed in detail an active excited to ground state relaxation pathway accompanied by an S 2 /S 3 (D 2d ) → S 1 (D 8h ) → S o (D 4h ) → S o (D 2d ) double-bond shifting mechanism. The simulated excitation lifetime is in good agreement with experiment. The first excited singlet state S1 plays a crucial role in the photochemistry. Here, the obtained critical molecular conformations, energy barrier, and transition-state lifetime results will provide a basis for further investigations of the bond-order inversion and photoswitching process of COT.

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

Real-Space Constrained Density Functional Theory Investigation of Site-Specific, Interfacial Charge Recombination Dynamics Across the Au Nanoparticle/TiO 2 Heterojunction

Au nanoparticle (NP)/TiO 2 heterojunction is a representative system to study interfacial charge transfer in photocatalysis and photovoltaics, where suppressing recombination from TiO 2 to Au can enhance hot carrier extraction. We apply real-space constrained density functional theory (CDFT) with Marcus theory to quantify charge recombination time scales across Au/TiO 2 . This approach enables direct control and visualization of charge-separated states, aligning with site-specific probes like time-resolved X-ray photoelectron spectroscopy (trXPS). We find that the charge-separated state features a bipolaron, with recombination dominated by TiO 2 LUMO to Au HOMO transitions, primarily at interfacial Au sites. Marcus rate predictions are benchmarked with surface hopping methods, quantifying differences in time scales and computational efficiency. Lastly, we examine how the Au cluster size affects the free energy change (ΔG) and reorganization energy (λ), explaining trends in closed-shell systems and highlighting challenges for open-shell extrapolations. Overall, CDFT + Marcus theory provides efficient, mechanistically transparent interfacial charge transfer modeling, and we clearly defined its applicability and limitation.

Glenna, Drew M. [Univ. of Idaho, Idaho Falls, ID (↗
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

Ultrafast Proton-Coupled Electron Transfer Reactions at Unit Quantum Yield

Proton-coupled electron transfer (PCET) reactions mediate essential biological energy-conversion processes, yet dynamical experiments that distinguish concerted electron/proton motion from sequential electron-transfer/proton-transfer pathways are unprecedented on ultrafast time scales. Here we report a tightly coupled donor–spacer–acceptor assembly, NDI-Sp-PhOH/BI, that enables direct interrogation of light-triggered and ultrafast PCET dynamics. Photoexcitation of the naphthalene diimide chromophore drives charge separation forming NDI radical anion and benzimidazolium phenoxy radical (PhO●/BIH+) products with a 440 fs time constant, followed by thermal charge recombination on a 5 ps time scale. Transient visible/near-infrared and mid-infrared spectroscopies independently track electron and proton motion and reveal coincident formation and decay of the charge-separated PCET products, excluding detectable phenol radical cation or phenoxide intermediates expected for sequential electron transfer-proton transfer (ET-PT) pathways. Surface-hopping quasiclassical trajectory simulations demonstrate rapid proton relocation coupled to evolution from locally excited to charge-transfer states, congruent with a concerted mechanism. These combined spectroscopic and computational results establish NDI-Sp-PhOH/BI as an archetype for ultrafast, concerted PCET occurring with unit quantum yield; further, because this platform provides distinct electronic and vibrational spectroscopic handles in its ground, electronically excited, and charge-separated states, it provides new opportunities to characterize mechanistic crossover in biomimetic PCET systems.

10 SYNTHETIC FUELS↗
DOE OSTI2024

Capturing the generation and structural transformations of molecular ions

Molecular ions are ubiquitous and play pivotal roles in many reactions, particularly in the context of atmospheric and interstellar chemistry. However, their structures and conformational transitions, particularly in the gas phase, are less explored than those of neutral molecules owing to experimental difficulties. A case in point is the halonium ions, whose highly reactive nature and ring strain make them short-lived intermediates that are readily attacked even by weak nucleophiles and thus challenging to isolate or capture before they undergo further reaction. Here we show that mega-electronvolt ultrafast electron diffraction (MeV-UED), used in conjunction with resonance-enhanced multiphoton ionization, can monitor the formation of 1,3-dibromopropane (DBP) cations and their subsequent structural dynamics forming a halonium ion. We find that the DBP + cation remains for a substantial duration of 3.6 ps in aptly named ‘dark states’ that are structurally indistinguishable from the DBP electronic ground state. The structural data, supported by surface-hopping simulations and ab initio calculations, reveal that the cation subsequently decays to iso-DBP + , an unusual intermediate with a four-membered ring containing a loosely bound bromine atom, and eventually loses the bromine atom and forms a bromonium ion with a three-membered-ring structure. We anticipate that the approach used here can also be applied to examine the structural dynamics of other molecular ions and thereby deepen our understanding of ion chemistry.

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

Electronic relaxation dynamics in [Au 25 (SR) 18 ] -1 (R = CH 3 , C 2 H 5 , C 3 H 7 , MPA, PET) thiolate-protected nanoclusters

Here we investigate the excited electron dynamics in [Au 25 (SR) 18 ] -1 (R = CH3, C 2 H 5 , C 3 H 7 , MPA, PET) [MPA = mercaptopropanoic acid, PET = phenylethylthiol] nanoparticles to understand how different ligands affect the excited state dynamics in this system. The population dynamics of the core and higher excited states lying in the energy range 0.00–2.20 eV are studied using a surface hopping method with decoherence correction in a real-time DFT approach. All of the ligated clusters follow a similar trend in decay for the core states (S 1 –S 6 ). The observed time constants are on the picosecond time scale (2–19 ps), which agrees with the experimental time scale, and this study confirms that the time constants observed experimentally could originate from core-to-core transitions and not from core-to-semiring transitions. In the presence of higher excited states, R = H, CH 3 , C 2 H 5 , C 3 H 7 , and PET demonstrate similar relaxations trends whereas R = MPA shows slightly different relaxation of the core states due to a smaller gap between the LUMO+1 and LUMO+2 gap in its electronic structure. The S 1 (HOMO → LUMO) state gives the slowest decay in all ligated clusters, while S 7 has a relatively long decay. Furthermore, separate electron and hole relaxations were performed on the [Au 25 (SCH 3 ) 18 ] -1 nanocluster to understand how independent electron and hole relaxations contribute to the overall relaxation dynamics.

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

Pathways to fluorescence via restriction of intramolecular motion in substituted tetraphenylethylenes

The design of materials with enhanced luminescence properties is a fast-developing field due to the potential applicability of these materials as light-emitting diodes or for bioimaging. A transparent way to enhance the emission properties of interesting molecular candidates is blocking competing and unproductive non-radiative relaxation pathways by the restriction of intramolecular motions. Rationalized functionalization is an important possibility to achieve such restrictions. Using time-dependent density functional theory (TD-DFT) based on the ωB97XD functional and the semiempirical tight-binding method including long-range corrections (TD-LC-DFTB), this work investigates the effect of functionalization of the paradigmatic tetraphenylethylene (TPE) on achieving restricted access to conical intersections (RACI). Photodynamical surface hopping simulations have been performed on a larger set of compounds including TPE and ten functionalized TPE compounds. Functionalization has been achieved by means of electron-withdrawing groups, bulky groups which block the relaxation channels via steric hindrance and groups capable of forming strong hydrogen bonds, which restrict the motion via the formation of hydrogen bond channels. Most of the investigated functionalized TPE candidates show ultrafast deactivation to the ground state due to their still existing structural flexibility, but two examples, one containing –CN and –CF 3 groups and a second characterized by a network of hydrogen bonds, have been identified as interesting candidates for creating efficient luminescence properties in solution.

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

Dissociative electron attachment to 5-bromo-uracil: non-adiabatic dynamics on complex-valued potential energy surfaces

Electron induced dissociation reactions are relevant to many fields, ranging from prebiotic chemistry to cancer treatments. However, the simulation of dissociation electron attachment (DEA) dynamics is very challenging because the auto-ionization widths of the transient negative ions must be accounted for. We propose an adaptation of the ab initio multiple spawning (AIMS) method for complex-valued potential energy surfaces, along the lines of recent developments based on surface hopping dynamics. Our approach combines models for the energy dependence of the auto-ionization widths, obtained from scattering calculations, with survival probabilities computed for the trajectory basis functions employed in the AIMS dynamics. The method is applied to simulate the DEA dynamics of 5-bromouracil in full dimensionality, i.e., taking all the vibrational modes into consideration. The propagation starts on the $π_{2}^{*}$ resonance state and describes the formation of Br - anions mediated by non-adiabatic couplings. The potential energies, gradients and non-adiabatic couplings were computed with the fractional-occupancy molecular orbital complete-active-space configuration-interaction method, and the calculated DEA cross section are consistent with the observed DEA intensities.

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

UV photodissociation dynamics of the acetone oxide Criegee intermediate: experiment and theory

Here, the photodissociation dynamics of the dimethyl-substituted acetone oxide Criegee intermediate [(CH 3 ) 2 COO] is characterized following electronic excitation to the bright 1 ππ* state, which leads to O ( 1 D) + acetone [(CH 3 ) 2 CO, S 0 ] products. The UV action spectrum of (CH 3 ) 2 COO recorded with O ( 1 D) detection under jet-cooled conditions is broad, unstructured, and essentially unchanged from the corresponding electronic absorption spectrum obtained using a UV-induced depletion method. This indicates that UV excitation of (CH 3 ) 2 COO leads predominantly to the O ( 1 D) product channel. A higher energy O ( 3 P) + (CH 3 ) 2 CO (T 1 ) product channel is not observed, although it is energetically accessible. In addition, complementary MS-CASPT2 trajectory surface-hopping (TSH) simulations indicate minimal population leading to the O ( 3 P) channel and non-unity overall probability for dissociation (within 100 fs). Velocity map imaging of the O ( 1 D) products is utilized to reveal the total kinetic energy release (TKER) distribution upon photodissociation of (CH 3 ) 2 COO at various UV excitation energies. Simulation of the TKER distributions is performed using a hybrid model that combines an impulsive model with a statistical component, the latter reflecting the longer-lived (>100 fs) trajectories identified in the TSH calculations. The impulsive model accounts for vibrational activation of (CH 3 ) 2 CO arising from geometrical changes between the Criegee intermediate and the carbonyl product, indicating the importance of CO stretch, CCO bend, and CC stretch along with activation of hindered rotation and rock of the methyl groups in the (CH 3 ) 2 CO product. Detailed comparison is also made with the TKER distribution arising from photodissociation dynamics of CH 2 OO upon UV excitation.

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

Different flavors of exact-factorization-based mixed quantum-classical methods for multistate dynamics

Here, the exact factorization approach has led to the development of new mixed quantum-classical methods for simulating coupled electron–ion dynamics. We compare their performance for dynamics when more than two electronic states are occupied at a given time, and analyze: (1) the use of coupled versus auxiliary trajectories in evaluating the electron–nuclear correlation terms, (2) the approximation of using these terms within surface-hopping and Ehrenfest frameworks, and (3) the relevance of the exact conditions of zero population transfer away from nonadiabatic coupling regions and total energy conservation. Dynamics through the three-state conical intersection in the uracil radical cation as well as polaritonic models in one dimension are studied.

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

Ultrafast structural dynamics of UV photoexcited cis , cis -1,3-cyclooctadiene observed with time-resolved electron diffraction

Conjugated diene molecules are highly reactive upon photoexcitation and can relax through multiple reaction channels that depend on the position of the double bonds and the degree of molecular rigidity. Understanding the photoinduced dynamics of these molecules is crucial for establishing general rules governing the relaxation and product formation. Here, in this study, we investigate the femtosecond time-resolved photoinduced excited-state structural dynamics of cis,cis-1,3-cyclooctadiene, a large-flexible cyclic conjugated diene molecule, upon excitation with 200 nm using mega-electron-volt ultrafast electron diffraction and trajectory surface hopping dynamics simulations. We tracked the photoinduced structural changes from the Franck–Condon region through the conical intersection seam to the ground state. Our findings revealed a novel primary reaction coordinate involving ring distortion, where the ring stretches along one axis and compresses along the perpendicular axis. The nuclear wavepacket remains compact along this reaction coordinate until it reaches the conical intersection seam, and it rapidly spreads as it approaches the ground state, where multiple products are formed.

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

Influence of substitution pattern on the dynamics of internal conversion and intersystem crossing in thiopyridone isomers

We report a combined experimental and theoretical investigation of the ultrafast internal conversion (IC) and intersystem crossing (ISC) dynamics of two thiopyridone (TP) isomers in solution. Our study used ultrafast transient X-ray absorption spectroscopy (XAS) at the sulfur K-edge, in conjunction with electronic excited state surface hopping molecular dynamics and simulations of the excited state XAS, to investigate the impact of the functional group substitution pattern and solvent on the dynamics of IC and ISC. The combination of the localized X-ray probe and the simulation results enables, in part, the differentiation between ππ* and nπ* character excited states, as well as singlet and triplet states. Access to nπ* character excitations has particular value since they often prove challenging to assess with optical spectroscopy. For 2-TP, the photoexcited S 2 (ππ*) state rapidly undergoes IC to the S 1 (nπ*) state below the instrument response time, followed by ISC to the T 1 (ππ*) state on a timescale of 600 fs in acetonitrile. For 4-TP, the timescale of S 2 to S 1 IC increases to 330 fs and the timescale of ISC increases to more than 10 ps. The differences between isomers are rationalized by considering the key role of the, nπ* intermediates in mediating the intersystem crossing of these systems. Varying the substitution pattern of the molecule can stabilize or destabilize these intermediates leading to the increase in ISC rate in the ortho isomer as compared to the para isomer, while changing the solvent from acetonitrile to water had minimal effect on the electronic excited state relaxation mechanism.

Garratt, Douglas [SLAC National Accelerator Labora↗
DOE OSTI2022

Nonadiabatic Molecular Dynamics Simulations Based on Time-Dependent Density Functional Tight-Binding Method

Nonadiabatic excited-state molecular dynamics underpin many photophysical and photochemical phenomena, such as exciton dynamics, charge separation and transport. In this work, we present an efficient nonadiabatic molecular dynamic (NAMD) simulation method based on time-dependent density functional tight-binding (TDDFTB) theory. Specifically, the adiabatic electronic structure, an essential NAMD input, is described at the TDDFTB level. Here, the nonadiabatic effects originating from the coupled motions of electrons and nuclei are treated by the trajectory surface hopping algorithm. To improve the computational efficiency, nonadiabatic couplings between excited states within the TDDFTB method are derived and implemented using an analytical approach. Further, the time-dependent nonadiabatic coupling scalars are calculated based on the overlap between molecular orbitals rather than the Slater determinants to speed up the simulations. In addition, the electronic decoherence scheme and a state reassigned unavoided crossings algorithm, which has been implemented in the NEXMD software, are used to improve the accuracy of the simulated dynamics and handle trivial unavoided crossings. Finally, the photoinduced nonadiabatic dynamics of a benzene molecule are simulated to demonstrate our implementation. The results for excited state NAMD simulations of benzene molecule based on TDDFTB method compare well that obtained with numerically expensive time-dependent density functional theory. The proposed methodology provides an attractive theoretical simulation tool for predicting the photophysical and photochemical properties of complex materials.

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

Ultrafast photochemistry and electron-diffraction spectra in n → (3s) Rydberg excited cyclobutanone resolved at the multireference perturbative level

We study the ultrafast time evolution of cyclobutanone excited to the singlet n → Rydberg state through non-adiabatic surface-hopping simulationsperformed at extended multi-state complete active space second-order perturbation (XMS-CASPT2) level of theory. These dynamics predict relaxation to the ground-state with a timescale of 822 ± 45 fs with minimal involvement of the triplets. The major relaxation path to the ground-state involves a three-state degeneracy region and leads to a variety of fragmented photoproducts. We simulate the resulting time-resolved electron-diffraction spectra, which track the relaxation of the excited state and the formation of various photoproducts in the ground state.

Chemistry↗
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↗