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

Conformer-Specific Dissociation Dynamics in Dimethyl Methylphosphonate Radical Cation

The dynamics of the dimethyl methylphosphonate (DMMP) radical cation after production by strong field adiabatic ionization have been investigated. Pump-probe experiments using strong field 1300 nm pulses to adiabatically ionize DMMP and a 800 nm non-ionizing probe induce coherent oscillations of the parent ion yield with a period of about 45 fs. The yields of two fragments, PO 2 C 2 H 7 + and PO 2 CH 4 + , oscillate approximately out of phase with the parent ion, but with a slight phase shift relative to each other. We use electronic structure theory and nonadiabatic surface hopping dynamics to understand the underlying dynamics. The results show that while the cation oscillates on the ground state along the P=O bond stretch coordinate, the probe excites population to higher electronic states that can lead to fragments PO 2 C 2 H 7 + and PO 2 CH 4 + . The computational results combined with the experimental observations indicate that the two conformers of DMMP that are populated under experimental conditions exhibit different dynamics after being excited to the higher electronic states of the cation leading to different dissociation products. These results highlight the potential usefulness of these pump-probe measurements as a tool to study conformer-specific dynamics in molecules of biological interest.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Exact Factorization Adventures: A Promising Approach for Non-Bound States

Modeling the dynamics of non-bound states in molecules requires an accurate description of how electronic motion affects nuclear motion and vice-versa. The exact factorization (XF) approach offers a unique perspective, in that it provides potentials that act on the nuclear subsystem or electronic subsystem, which contain the effects of the coupling to the other subsystem in an exact way. We briefly review the various applications of the XF idea in different realms, and how features of these potentials aid in the interpretation of two different laser-driven dissociation mechanisms. We present a detailed study of the different ways the coupling terms in recently-developed XF-based mixed quantum-classical approximations are evaluated, where either truly coupled trajectories, or auxiliary trajectories that mimic the coupling are used, and discuss their effect in both a surface-hopping framework as well as the rigorously-derived coupled-trajectory mixed quantum-classical approach.

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

Unusual magnetotransport in twisted bilayer graphene from strain-induced open Fermi surfaces

Anisotropic hopping in a toy Hofstadter model was recently invoked to explain a rich and surprising Landau spectrum measured in twisted bilayer graphene away from the magic angle. Suspecting that such anisotropy could arise from unintended uniaxial strain, we extend the Bistritzer–MacDonald model to include uniaxial heterostrain and present a detailed analysis of its impact on band structure and magnetotransport. We find that such strain strongly influences band structure, shifting the three otherwise-degenerate van Hove points to different energies. Coupled to a Boltzmann magnetotransport calculation, this reproduces previously unexplained nonsaturating B 2 magnetoresistance over broad ranges of density near filling ν = ± 2 and predicts subtler features that had not been noticed in the experimental data. In contrast to these distinctive signatures in longitudinal resistivity, the Hall coefficient is barely influenced by strain, to the extent that it still shows a single sign change on each side of the charge neutrality point—surprisingly, this sign change no longer occurs at a van Hove point. The theory also predicts a marked rotation of the electrical transport principal axes as a function of filling even for fixed strain and for rigid bands. More careful examination of interaction-induced nematic order versus strain effects in twisted bilayer graphene could thus be in order.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Versatile stochastic model for predictive KMC simulation of fcc metal nanostructure evolution with realistic kinetics

Stochastic lattice-gas models provide the natural framework for analysis of the surface diffusion-mediated evolution of crystalline metal nanostructures on the appropriate time scale (often 10 1 –10 4 s) and length scale. Model behavior can be precisely assessed by kinetic Monte Carlo simulation, typically incorporating a rejection-free algorithm to efficiently handle the broad range of Arrhenius rates for hopping of surface atoms. The model should realistically prescribe these rates, or the associated barriers, for a diversity of local surface environments. However, commonly used generic choices for barriers fail, even qualitatively, to simultaneously describe diffusion for different low-index facets, for terrace vs step edge diffusion, etc. We introduce an alternative Unconventional Interaction–Conventional Interaction formalism to prescribe these barriers, which, even with few parameters, can realistically capture most aspects of behavior. Here, the model is illustrated for single-component fcc metal systems, mainly for the case of Ag. It is quite versatile and can be applied to describe both the post-deposition evolution of 2D nanostructures in homoepitaxial thin films (e.g., reshaping and coalescence of 2D islands) and the post-synthesis evolution of 3D nanocrystals (e.g., reshaping of nanocrystals synthesized with various faceted non-equilibrium shapes back to 3D equilibrium Wulff shapes).

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

Influence of tungsten doping on nonradiative electron–hole recombination in monolayer MoSe 2 with Se vacancies

Two-dimensional transition metal dichalcogenides (TMDs) are receiving significant attention due to their excellent electronic and optoelectronic properties. The material quality is greatly affected by defects that are inevitably generated during material synthesis. Focusing on chalcogenide vacancies, which constitute the most common defect, we use the state-of-the-art simulation methodology developed in our group to demonstrate that W doping of MoSe 2 with Se vacancies reduces charge carrier losses by two mechanisms. First, W doping makes the formation of double Se vacancies unfavorable, while it is favorable in undoped MoSe 2 . Second, if a Se vacancy is present, the charge carrier lifetimes are extended in the W-doped MoSe 2 . Combining ab initio real-time time-dependent density functional theory with nonadiabatic molecular dynamics, the simulations show that the nonradiative carrier losses in the presence of Se vacancies proceed by sub-10 ps electron trapping and relaxation down the manifold of trap states, followed by a 100 ps recombination of trapped electrons with free holes. The electron–vibrational energy exchange is driven by both in-plane and out-of-plane vibrational motions of the MoSe 2 layer. Additionally, the atomistic studies advance our understanding of the influence of defects on charge carrier properties in TMDs and guide improvements of material quality and development of TMD applications.

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

Case studies of the time-dependent potential energy surface for dynamics in cavities

Here, the exact time-dependent potential energy surface driving the nuclear dynamics was recently shown to be a useful tool to understand and interpret the coupling of nuclei, electrons, and photons in cavity settings. Here, we provide a detailed analysis of its structure for exactly solvable systems that model two phenomena: cavity-induced suppression of proton-coupled electron-transfer and its dependence on the initial state, and cavity-induced electronic excitation. We demonstrate the inadequacy of simply using a weighted average of polaritonic surfaces to determine the dynamics. Such a weighted average misses a crucial term that redistributes energy between the nuclear and the polaritonic systems, and this term can in fact become a predominant term in determining the nuclear dynamics when several polaritonic surfaces are involved. Evolving an ensemble of classical trajectories on the exact potential energy surface reproduces the nuclear wavepacket quite accurately, while evolving on the weighted polaritonic surface fails after a short period of time. The implications and prospects for application of mixed quantum-classical methods based on this surface are discussed.

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

Non-Kasha fluorescence of pyrene emerges from a dynamic equilibrium between excited states

Pyrene fluorescence after a high-energy electronic excitation exhibits a prominent band shoulder not present after excitation at low energies. The standard assignment of this shoulder as a non-Kasha emission from the second-excited state (S 2 ) has been recently questioned. To elucidate this issue, we simulated the fluorescence of pyrene using two different theoretical approaches based on vertical convolution and nonadiabatic dynamics with nuclear ensembles. To conduct the necessary nonadiabatic dynamics simulations with high-lying electronic states and deal with fluorescence timescales of about 100 ns of this large molecule, we developed new computational protocols. The results from both approaches confirm that the band shoulder is, in fact, due to S 2 emission. We show that the non-Kasha behavior is a dynamic-equilibrium effect not caused by a metastable S 2 minimum. However, it requires considerable vibrational energy, which can only be achieved in collisionless regimes after transitions into highly excited states. This strict condition explains why the S 2 emission was not observed in some experiments.

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

A windowed mean trajectory approximation for condensed phase dynamics

We propose a trajectory-based quasi-classical method for approximating dynamics in condensed phase systems. Building upon the previously developed optimized mean trajectory approximation that has been used to compute linear and nonlinear spectra, we borrow some ideas from filtering trajectory methods to obtain a novel semiclassical method for the dynamical propagation of density matrices. This new approximation is tested rigorously against standard multistate electronic models, spin-boson models, and models of the Fenna–Matthews–Olson complex. For dissipative systems, the current method is significantly better or as good as many other semiclassical methods available, especially at low temperatures and for off-diagonal density matrix elements, whereas for scattering models, the current method bears similar limitations as mean-field propagation schemes. All results are tested against the numerically exact hierarchical equations of motion method. In conclusion, the new method shows excellent agreement across various parameter regimes with numerically exact results, highlighting the robustness and accuracy of our approach.

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

Shape Stability of Truncated Octahedral fcc Metal Nanocrystals

Metallic nanocrystals (NCs) can be synthesized with tailored nonequilibrium shapes to enhance desired properties, e.g., octahedral fcc metal NCs optimize catalytic activity associated with {111} facets. However, maintenance of optimized properties requires stability against thermal reshaping. Thus, we analyze the reshaping of truncated fcc metal octahedra mediated by surface diffusion using a stochastic atomistic-level model with energetic input parameters for Pd. Furthermore, the model describes NC thermodynamics by an effective nearest-neighbor interaction and includes a realistic treatment of diffusive hopping for undercoordinated surface atoms. Kinetic Monte Carlo simulation reveals that the effective barrier, E eff , for the initial stage of reshaping is strongly tied to the degree of truncation of the vertices in the synthesized initial octahedral shapes. This feature is elucidated via exact analytic determination of the energy variation along the optimal reshaping pathway at low-temperature (T), which involves transfer of atoms from truncated {100} vertex facets to form new layers on {111} side facets. Deviations from predictions of the low-T analysis due to entropic effects are more prominent for higher T and larger NC sizes.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Reshaping of Truncated Pd Nanocubes: Energetic and Kinetic Analysis Integrating Transmission Electron Microscopy with Atomistic-Level and Coarse-Grained Modeling

Stability against reshaping of metallic fcc nanocrystals synthesized with tailored far-from-equilibrium shapes is key to maintaining optimal properties for applications such as catalysis. Yet Arrhenius analysis of experimental reshaping kinetics, and appropriate theory and simulation, is lacking. As such, we use TEM to monitor the reshaping of Pd nanocubes of ~25 nm side length between 410 °C (over ~4.5 h) and 440 °C (over ~0.25 h), extracting a high effective energy barrier of E eff ≈ 4.6 eV. We also provide an analytic determination of the energy variation along the optimal pathway for reshaping that involves transfer of atoms across the nanocube surface from edges or corners to form new layers on side {100} facets. The effective barrier from this analysis is shown to increase strongly with the degree of truncation of edges and corners in the synthesized nanocube. Theory matches experiment for the appropriate degree of truncation. In addition, we perform simulations of a stochastic atomistic-level model incorporating a realistic description of diffusive hopping for undercoordinated surface atoms, thereby providing a visualization of the initial reshaping process.

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

Enhanced Selectivity for C 2 H 4 Production from C 2 H 6 on Partially Chlorinated IrO 2 (110) Surfaces

Modifying metal oxide surfaces to limit their oxidizing activity can provide a means of improving catalytic selectivity toward the partial oxidation of light alkanes. Here, in this study, we investigated the oxidation of C 2 H 6 on Cl-modified IrO 2 (110) surfaces using temperature-programmed reaction spectroscopy (TPRS) and first-principles microkinetic modeling. We find that substituting Cl for O in the IrO 2 (110) surface enhances the selectivity for C 2 H 6 conversion to C 2 H 4 during TPRS by suppressing extensive oxidation to CO x products, while also either enhancing C 2 H 4 production or altering it to a lesser extent, depending on the initial C 2 H 6 coverage. The C 2 H 4 selectivity increased with increasing C 2 H 6 and Cl coverage, but reached a limiting value below 50%. The Cl coverage changed negligibly during C 2 H 6 oxidation, and the surface reactivity decreased only marginally for Cl coverages up to 0.5 ML (monolayer). TPRS simulations using a microkinetic model predict C 2 H 4 and CO x product yields as a function of the Cl coverage that agree closely with the experimental results. According to the simulations, C 2 H 6 conversion to C 2 H 4 occurs on Cl-IrO 2 (110) by the hydrogenation of C 2 H 3 * species adsorbed in blocked states, in which neighboring sites are occupied only by unreactive HO and Cl species. The microkinetic modeling shows that H-hopping away from surface HO groups provides a relatively efficient route for C 2 H 3 * to escape blocked configurations and dehydrogenate, and that this process can limit the C 2 H 4 selectivity on Cl-IrO 2 (110) under the conditions studied. Overall, our results demonstrate that Cl-substitution into IrO 2 (110) enhances the selectivity for C 2 H 4 production from C 2 H 6 and provides insights into the reaction mechanism that can guide strategies to further improve the C 2 H 4 selectivity.

IrO2↗
DOE OSTI2022

Estimates of Quantum Tunneling Effects for Hydrogen Diffusion in PuO 2

We detail the estimation of activation energies and quantum nuclear vibrational tunneling effects for hydrogen diffusion in PuO 2 based on Density Functional Theory calculations and a quantum double well approximation. We find that results are relatively insensitive to choice of exchange correlation functional. In addition, the representation of spin in the system and use of an extended Hubbard U correction has only a small effect on hydrogen point defect formation energies when the PuO 2 lattice is held fixed at the experimental density. We then compute approximate activation energies for transitions between hydrogen interstitial sites seeded by a semi-empirical quantum model and determine the quantum tunneling enhancement relative to classical kinetic rates. Our model indicates that diffusion rates in H/PuO 2 systems could be enhanced by more than one order of magnitude at ambient conditions and that these effects persist at high temperature. The method we propose here can be used as a fast screening tool for assessing possible quantum nuclear vibrational effects in any number of condensed phase materials and surfaces, where hydrogen hopping tends to follow well defined minimum energy pathways.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Deep reinforcement learning for predicting kinetic pathways to surface reconstruction in a ternary alloy

Abstract The majority of computational catalyst design focuses on the screening of material components and alloy composition to optimize selectivity and activity for a given reaction. However, predicting the metastability of the alloy catalyst surface at realistic operating conditions requires an extensive sampling of possible surface reconstructions and their associated kinetic pathways. We present CatGym, a deep reinforcement learning (DRL) environment for predicting the thermal surface reconstruction pathways and their associated kinetic barriers in crystalline solids under reaction conditions. The DRL agent iteratively changes the positions of atoms in the near-surface region to generate kinetic pathways to accessible local minima involving changes in the surface compositions. We showcase our agent by predicting the surface reconstruction pathways of a ternary Ni 3 Pd 3 Au 2 (111) alloy catalyst. Our results show that the DRL agent can not only explore more diverse surface compositions than the conventional minima hopping method, but also generate the kinetic surface reconstruction pathways. We further demonstrate that the kinetic pathway to a global minimum energy surface composition and its associated transition state predicted by our agent is in good agreement with the minimum energy path predicted by nudged elastic band calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Vacancy surface migration mechanisms in dilute nickel-chromium alloys

We investigate the unitary mechanisms related to the surface migration of vacancies in dilute Ni-Cr alloys via first-principle calculations. We survey a complete set of surface and sub-surface migration paths for vacancies near the (100) free surface and calculate the corresponding migration barriers. Our results show that a vacancy migrating towards the free surface will face lower energy barriers to migrate via atomic exchange with a neighboring Cr atom rather than with a Ni atom. Once a vacancy reaches the free surface, it will be trapped there. Our results also reveal that, when a Cr atom sits in the atomic plane just below the free surface, any in-plane vacancy hopping jump that would result in the vacancy sitting on top a subsurface Cr atom is energetically unfavorable. Taken together, these fundamental unitary surface migration mechanisms offer insights into the complex interactions between surface segregation and vacancy migration phenomena in Ni-Cr-based alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Surface coverage dynamics for reversible dissociative adsorption on finite linear lattices

Dissociative adsorption onto a surface introduces dynamic correlations between neighboring sites not found in non-dissociative absorption. We study surface coverage dynamics where reversible dissociative adsorption of dimers occurs on a finite linear lattice. We derive analytic expressions for the equilibrium surface coverage as a function of the number of reactive sites, N, and the ratio of the adsorption and desorption rates. Using these results, we characterize the finite size effect on the equilibrium surface coverage. For comparable N’s, the finite size effect is significantly larger when N is even than when N is odd. Moreover, as N increases, the size effect decays more slowly in the even case than in the odd case. The finite-size effect becomes significant when adsorption and desorption rates are considerably different. These finite-size effects are related to the number of accessible configurations in a finite system where the odd-even dependence arises from the limited number of accessible configurations in the even case. We confirm our analytical results with kinetic Monte Carlo simulations. We also analyze the surface-diffusion case where adsorbed atoms can hop into neighboring sites. As expected, the odd-even dependence disappears because more configurations are accessible in the even case due to surface diffusion.

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

Investigation of the electrical conduction mechanisms in P-type amorphous germanium electrical contacts for germanium detectors in searching for rare-event physics

For the first time, electrical conduction mechanisms in the disordered material system is experimentally studied for p-type amorphous germanium (a-Ge) used for high-purity Ge detector contacts. The localization length and the hopping parameters in a-Ge are determined using the surface leakage current measured from three high-purity planar Ge detectors. The temperature dependent hopping distance and hopping energy are obtained for a-Ge fabricated as the electrical contact materials for high-purity Ge planar detectors. As a result, we find that the hopping energy in a-Ge increases as temperature increases while the hopping distance in a-Ge decreases as temperature increases. The localization length of a-Ge is on the order of 2.13$^{-0.05}_{+0.07}$A° to 5.07$^{0.83}_{+2.58}$A°,depending on the density of states near the Fermi energy level within bandgap. Using these parameters, we predict that the surface leakage current from a Ge detector with a-Ge contacts can be much smaller than one yocto amp (yA) at helium temperature, suitable for rare-event physics searches.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Hot-lines topology and the fate of the spin resonance mode in three-dimensional unconventional superconductors

In the quasi-two-dimensional (quasi-2D) copper- and iron-based superconductors, the onset of superconductivity is accompanied by a prominent peak in the magnetic spectrum at momenta close to the wave-vector of the nearby antiferromagnetic state. Such a peak is well described in terms of a spin resonance mode, i.e., a spin-1 exciton theoretically predicted for quasi-2D superconductors with a sign-changing gap. The same theories, however, indicate that such a resonance mode should be absent in a three-dimensional (3D) system with a spherical Fermi surface. This raises the question of the fate of the spin resonance mode in layered unconventional superconductors that are not strongly anisotropic, such as certain heavy-fermion compounds and potentially the newly discovered nickelate superconductor NdNiO 2 . Here, we use the random phase approximation to calculate the dynamical spin susceptibility of 3D superconductors with a d x 2 –y 2 -wave gap symmetry and corrugated cylindrical-like Fermi surfaces. By varying the out-of-plane hopping anisotropy t z /t, we demonstrate that the appearance of a spin resonance mode is determined by the topology of the hot lines, i.e., lines on the Fermi surface that are connected by the magnetic wave vector. For an in-plane antiferromagnetic wave vector, the hot lines undergo a topological transition from open lines to closed loops at a critical t z /t value. The closed hot lines cross the nodal superconducting lines, making the spin resonance mode overdamped and incoherent. In contrast, for an out-of-plane antiferromagnetic wave vector, the hot lines remain open and the spin resonance mode remains sharp. We discuss the experimental implications of our results for the out-of-plane dispersion of the spin resonance mode and, more generally, for inelastic neutron scattering experiments on unconventional superconductors.

36 MATERIALS SCIENCE↗