Search NASA⌕ Search

SEARCH · Search NASA

Results for “potential energy landscape”

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 OSTI2023

Dependence of simulated radiation damage on crystal structure and atomic misfit in metals

This study investigates the evolution of radiation damage in three metals in the low temperature and high radiant flux regime using molecular statics and a Frenkel pair accumulation method to simulate up to 2.0 displacements per atom. The metals considered include Fe, equiatomic CrCoNi, and a fictitious metal with similar bulk properties to the CrCoNi composed of a single atom type referred to as an A-atom. CrCoNi is found to sustain higher concentrations of dislocations than either the Fe or A-atom systems and more stacking faults than the A-atom system. The results suggest that the difference between the concentrations of vacancies and interstitials is substantially smaller for CrCoNi than the A-atom system, perhaps reflecting that the sink capture radius is smaller in CrCoNi due to the roughened potential energy landscape. A model that partitions the major contributions from defects to the stored energy is described, and serves to highlight a general need for higher fidelity approaches to point defect identification.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Shape memory effect in metallic glasses

Shape memory effect (SME), mainly present in crystalline Ti-Ni alloys, is basically missing in metallic glasses (MGs) that lack the long-range periodic order of crystals. In this work, we report experimental results of SME in annealed MGs, in which the low-energy configuration state recovery is observed by both differential scanning calorimetry and X-ray photon correlation spectroscopy. We elucidate the origin of SME in MGs under the potential energy landscape framework, i.e., after annealing, the energy of MGs enters into a deep basin and atoms are located in the low-energy configuration state. Albeit deviating from their relative stable configurations by temperature changes, atoms in the annealed MGs tend to return to the low-energy atomic configurations along certain trajectories as pre-annealing temperature is approached, These results could extend the application of MGs as functional materials by directionally manipulating their energy states via annealing and rejuvenation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Density-Matrix Based Extended Lagrangian Born–Oppenheimer Molecular Dynamics

Extended Lagrangian Born–Oppenheimer molecular dynamics [ Phys. Rev. Lett. 2008, 100, 123004] is presented for Hartree–Fock theory, where the extended electronic degrees of freedom are represented by a density matrix, including fractional occupation numbers at elevated electronic temperatures. In contrast to regular direct Born–Oppenheimer molecular dynamics simulations, no iterative self-consistent field optimization is required prior to the force evaluations. To sample regions of the potential energy landscape where the gap is small or vanishing, which leads to particular convergence problems in regular direct Born–Oppenheimer molecular dynamics simulations, an adaptive integration scheme for the extended electronic degrees of freedom is presented. The integration scheme is based on a tunable, low-rank approximation of a fourth-order kernel, K, that determines the metric tensor, T ≡ K T K, used in the extended harmonic oscillator of the Lagrangian that generates the dynamics of the electronic degrees of freedom. Here, the formulation and algorithms provide a general guide to implement extended Lagrangian Born–Oppenheimer molecular dynamics for quantum chemistry, density functional theory, and semiempirical methods using a density matrix formalism.

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

General Design Rules for Bimetallic Platinum(II) Complexes

A series of platinum(II) bimetallic complexes were studied to investigate the effects of ligands on both the geometric and electronic structure. Modulating the Pt-Pt distance through the bridging ligand architecture was found to dictate the nature of the lowest energy electronic transitions, localized in one-half of the molecule or delocalized across the entire molecule. By reducing the separation between the platinum atoms, the lowest energy electronic transitions will be dominated by the metal-metal-to-ligand charge transfer transition. Conversely, by increasing the distance between the platinum atoms, the lowest electronic transition will be largely localized metal-to-ligand charge transfer or ligand centered in nature. Additionally, the cyclometalating ligands were observed to have a noticeable stabilizing effect on the triplet excited states as the conjugation increased, arising from geometric reorientation and increased electron delocalization of the ligands. Furthermore, such stabilization of the triplet state energy has been shown to alter the excited state potential energy landscape as well as the excited state trajectory.

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

Water Network Shape-Dependence of Local Interactions with the Microhydrated –NO 2 – and –CO 2 – Anionic Head Groups by Cold Ion Vibrational Spectroscopy

We report the structural evolutions of water networks and solvatochromic response of the CH 3 NO 2 – radical anion in the OH and CH stretching regions by analysis of the vibrational spectra displayed by cryogenically cooled CH 3 NO 2 – ·(H 2 O) n=1–6 clusters. The OH stretching bands evolve with a surprisingly large discontinuity at n = 6, which features the emergence of an intense, strongly red-shifted band along with a weaker feature that appears in the region assigned to a free OH fundamental. Very similar behavior is displayed by the perdeuterated carboxylate clusters, RCO 2 – ·(H 2 O) n=5–7 (R = CD 3 CD 2 ), indicating that this behavior is a general feature in the microhydration of the triatomic anionic domain and not associated with CH oscillators. Electronic structure calculations trace this behavior to the formation of a “book” isomer of the water hexamer that adopts a configuration in which one of the water molecules resides in an acceptor–acceptor–donor (AAD) (A = acceptor, D = donor) H-bonding site. Excitation of the bound OH in the AAD site explores the local network topology best suited to stabilize an incipient –XO 2 H–OH–(H 2 O) 2 intracluster proton-transfer reaction. These systems thus provide particularly clear examples where the network shape controls the potential energy landscape that governs water network-mediated, intracluster proton transfer. The CH stretching bands of the CH 3 NO 2 – ·(H 2 O) n=1–6 clusters also exhibit strong solvatochromic shifts, but in this case, they smoothly blue-shift with increasing hydration with no discontinuity at n = 6. Furthermore, this behavior is analyzed in the context of the solute-ion polarizability response and partial charge transfer to the water networks.

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

On the Red Antenna States of Photosystem I Mutants from Cyanobacteria Synechocystis PCC 6803

To identify the molecular composition of the low-energy states in cyanobacterial Photosystem I (PSI) of Synechocystis PCC6803, we focus on high resolution (low-temperature) absorption, emission, resonant, and nonresonant hole burned spectra obtained for wild-type (WT) PSI and three PSI mutants. In the Red_a mutant, the B33 chlorophyll (Chl) is added to the B31-B32 dimer; in Red_b, histidine 95 (His95) on PsaB (which coordinates Mg in the B7 Chl within the His95-B7-A31-A32- cluster) is replaced with glutamine (Gln), while in the Red_ab mutant, both mutations are made. We show that the C706 state (B31-B32) changes to the C710 state (B31-B32-B33) in both Red_a and Red_ab mutants, while the C707 state in WT Synechocystis (localized on the His95-B7-A31-A32 cluster) is modified to C716 in both Red_b and Red_ab. Excitation energy transfer from C706 to the C714 trap in the WT PSI and Red_b mutant is hampered as reflected by a weak emission at 712 nm. Large electron-phonon coupling strength (exposed via resonant hole-burned spectra) is consistent with a strong mixing of excited states with intermolecular charge transfer states leading to significantly red-shifted emission spectra. We conclude that excitation energy transfer in PSI is controlled by fine-tuning the electronic states of a small number of highly conserved red states. Finally, we show that mutations modify the protein potential energy landscape as revealed by different shapes and shifts of the blue- and red-shifted antiholes.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Accelerated Steam Methane Reforming by Dynamically Applied Charges

Catalyst design has traditionally focused on tuning active site properties to optimally bind reaction intermediates and balance the kinetic requirements of multiple competing chemical processes, as necessitated by the Sabatier principle. It has recently been proposed that for reactions following certain potential energy landscapes, the activity limit imposed by the Sabatier principle may be overcome by using programmed oscillations of surface electron density at the timescales of surface reactions (i.e., “catalytic resonance”). Here, we use a combination of density functional theory (DFT) simulations and transient kinetic models (TKMs) to simulate the kinetics of steam methane reforming (SMR) on Ru(211) surfaces under statically and dynamically applied charges. DFT-calculated binding energies of SMR intermediates and transition states exhibit strong sensitivity to positively applied charges and follow unique scaling relationships that deviate from linear periodic trends across transition metals. Our simulations demonstrate that applying a small positive charge to Ru dramatically enhances the steady-state turnover frequency (TOF) of SMR by up to 5 orders of magnitude above the TOF observed over neutral Ru. Thus, statically charging Ru catalysts may be an effective strategy to lower the temperature requirements for SMR. Dynamic square-wave oscillations in charge resulted in SMR catalytic resonance with an onset frequency f ∼ 106 Hz and the corresponding average TOFs exceeding the statically charged Ru surface by an additional 15%. Here, based on sensitivity analyses performed for the two end points of oscillation, we propose that dynamic TOF improvement beyond the Sabatier maximum can be expected when the system oscillates between two kinetic regimes that are uniquely controlled by distinct elementary steps.

Catalysts↗
DOE OSTI2026

Ultrafast Relaxation of MLCT Excited States in Manganese Tricarbonyl Complexes: Insights from Polarization-Resolved Femtosecond X-ray Absorption Spectroscopy at the Mn and Br K-Edges

Optical absorption features that are often described as metal-to-ligand charge transfer (MLCT) bands underlie the utility of many metal coordination complexes by harnessing the energy of light to drive otherwise inaccessible chemical reactions. Excitation of these bands triggers rapid electronic dynamics that can be challenging to understand, due to complicated potential energy landscapes and highly correlated electronic and nuclear degrees of freedom in metal-containing compounds. The lowest-energy absorption bands in Mn complexes containing alpha-diimine, carbonyl, and halide ligands are particularly interesting, due to the delocalized nature of charge transfer upon excitation. Here, in this study, we report experimental ultrafast dynamics measurements for the series of compounds Mn­(CO) 3 ( R bpy)Br ( R bpy = 4,4′-disubstituted-2,2′-bipyridine; R = H, CF 3 , or NO 2 ) using polarization-resolved, femtosecond X-ray absorption spectroscopy (XAS) at both the Mn and Br K edges. The appearance of a new absorption feature in the Br pre-edge spectrum upon optical excitation reveals the instantaneous formation of an electronic hole that is partially localized on the Br atom and has a lifetime that depends on the electron withdrawing character of the R bpy ligand. For two of the complexes (R = H or CF 3 ), a large expansion of the Mn–Br distance results in a rapid redistribution of the hole and a corresponding decrease in anisotropy of the absorption feature within 50 fs, after which the absorption into the hole disappears on a time scale shorter than 300 fs. We observe a different result for the NO 2 substituted compound, for which the Mn–Br bond contracts and the absorption into the Br-centered hole persists beyond the 2 ps time scale of our measurement. The Mn pre-edge spectrum also reveals structural changes for these three complexes, but the new Mn absorption features become evident only after the nuclei respond to the initial excitation, which allows mixing of Mn 3d and 4p orbitals. The combined use of ultrafast Mn and Br K-edge spectroscopy provides unique insight into the ways in which bipyridine substitution alters the excited-state dynamics, including a very different structural response for the most strongly electron withdrawing substituent.

Otolski, Christopher J. [Univ. of Kansas, Lawrence↗
DOE OSTI2026

Engineering Structural Transitions in a Multilevel Molecular Switch via Intermolecular Coupling

Controlling molecular conformations with atomic precision is essential for advancing molecular functional electronics, as well as our understanding of molecular dynamics. While switching between bistable molecular conformers is common in nature, creating systems with multiple, addressable states remains synthetically challenging. Here, we demonstrate a bottom-up strategy in which intermolecular interactions give rise to multilevel functionality within a simple two-molecule assembly. Using low-temperature scanning tunneling microscopy, we show that a pyrrolidine dimer on Cu(100) exhibits six distinct adsorption conformations, exceeding the four expected from two independent bistable units. This unusual complexity arises from the interplay between intermolecular van der Waals attraction and steric repulsion, which reshapes the potential energy landscape and changes a single high-energy transition into a sequential two-step pathway. Each step is driven by low-energy inelastic electron excitations, achieving a switching efficiency an order of magnitude higher than that of the monomer. Here, by tuning the bias voltage and tip–molecule distance, we achieve deterministic control over multiple stable states, establishing a general design principle for on-demand engineering of collective molecular behavior and energy-efficient multilevel molecular devices.

Molecular interactions↗
DOE OSTI2025

The classical-to-quantum crossover in the strain-induced ferroelectric transition in SrTiO3 membranes

Mechanical strain presents an effective control over symmetry-breaking phase transitions. In quantum paralelectric SrTiO3, strain can induce ferroelectric order via modification of the local Ti potential energy landscape. However, brittle bulk materials can only withstand limited strain range (~0.1%). Taking advantage of nanoscopically-thin freestanding membranes, we demonstrate an in- situ strain-induced reversible ferroelectric transition in freestanding SrTiO3 membranes. We measure the ferroelectric order by detecting the local anisotropy of the Ti 3d orbital signature using x-ray linear dichroism at the Ti-K pre-edge, while the strain is determined by x-ray diffraction. With reduced thickness, the SrTiO3 membranes remain elastic with >1% tensile strain cycles. A robust displacive ferroelectricity appears beyond a temperature-dependent critical strain. Interestingly, we discover a crossover from a classical ferroelectric transition to a quantum regime at low temperatures, which enhances strain-induced ferroelectricity. Our results offer new opportunities to strain engineer functional properties in low dimensional quantum materials and provide new insights into the role of ferroelectric fluctuations in quantum paraelectric SrTiO3.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Real-space observation of the dissociation of a transition metal complex and its concurrent energy redistribution

Mechanistic insights into photodissociation dynamics of transition metal carbonyls, like Fe(CO) 5 , are fundamental for understanding active catalytic intermediates. Although extensively studied, the structural dynamics of these systems remain elusive. Using ultrafast X-ray scattering, we uncover the photochemistry of Fe(CO) 5 in real space and time, observing synchronous oscillations in atomic pair distances, followed by a prompt rotating CO release preferentially in the axial direction. This behavior aligns with simulations, reflecting the interplay between the axial Fe-C distances’ potential energy landscape and non-adiabatic transitions between metal-to-ligand charge-transfer states. Additionally, we characterize a secondary delayed CO release associated with a reduction of Fe-C steady state distances and structural dynamics of the formed Fe(CO) 4 . Our results quantify energy redistribution across vibration, rotation, and translation degrees of freedom, offering a microscopic view of complex structural dynamics, enhancing our grasp on Fe(CO) 5 photodissociation, and advancing our understanding of transition metal catalytic systems.

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

Atomistic origin of the entropy of melting from inelastic neutron scattering and machine learned molecular dynamics

The latent heat, L, is central to melting, but its atomic origin remains elusive. It is proportional to the entropy of fusion, ΔS fus = L/T m (T m is the melting temperature), which depends on changes of atom configurations, atom vibrations, and thermal electron excitations. Here, we combine inelastic neutron scattering and machine-learned molecular dynamics to separate ΔS fus into these components for Ge, Si, Bi, Sn, Pb, and Li. When the vibrational entropy of melting, ΔS vib , is zero, ΔS fus ≃ 1.2 kB per atom. This result provides a baseline for ΔSconfig and nearly coincides with “Richard’s Rule” of melting. The ΔSfus deviates from this value for most elements, however, and we show that this deviation originates with extra ΔS vib and extra ΔS config . These two components are correlated for positive and negative deviations from Richard’s rule – the extra ΔS config is consistently ~ 80% of ΔS vib . Our results, interpreted with potential energy landscape theory, imply a correlation between the change in the number of basins and the change in the inverse of their curvature for the melting of pure elements.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Fundamental data for modeling electron-induced processes in plasma remediation of perfluoroalkyl substances

Plasma treatment of per- and polyfluoroalkyl substances (PFAS) contaminated water is a potentially energy efficient remediation method. In this treatment, an atmospheric pressure plasma interacts with surface-resident PFAS molecules. Developing a reaction mechanism and modeling of plasma–PFAS interactions requires fundamental data for electron–molecule reactions. In this paper, we present results of electron scattering calculations, potential energy landscapes and their implications for plasma modelling of a dielectric barrier discharge in PFAS contaminated gases, a first step towards modelling of plasma–water–PFAS intereactions. It is found that the plasma degradation of PFAS is dominated by dissociative electron attachment with the importance of other contributing processes varying depending on the molecule. All molecules posses a large number of shape resonances – transient negative ion states – from near-threshold up to ionization threshold. These states lie in the region of the most probable electron energies in the plasma (4–5 eV) and consequently are expected to further enhance the fragmentation dynamics in both dissociative attachment and dissociative excitation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Parallel algorithms for hyperdynamics and local hyperdynamics

Hyperdynamics (HD) is a method for accelerating the timescale of standard molecular dynamics (MD). It can be used for simulations of systems with an energy potential landscape that is a collection of basins, separated by barriers, where transitions between basins are infrequent. HD enables the system to escape from a basin more quickly while enabling a statistically accurate renormalization of the simulation time, thus effectively boosting the timescale of the simulation. In [Kim, Perez, Voter, J Chem Phys, 139:144110, 2013)1, a local version of HD was formulated, which exploits the intrinsic locality characteristic typical of most systems to mitigate the poor scaling properties of standard HD as the system size is increased. In this paper, we discuss how both HD and local HD can be formulated to run efficiently in parallel. We have implemented these ideas in the LAMMPS MD code, which means HD can be used with any interatomic potential LAMMPS supports. Together, these parallel methods allow simulations of any size to achieve the time acceleration offered by HD (which can be orders of magnitude), at a cost 3-5x that of standard MD. As examples, we performed two simulations of a million-atom system to model the diffusion and clustering of Pt adatoms on a large patch of Pt(100) surface for 80 and 160 μs.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

A theory of localized excitations in supercooled liquids

Here, a new connection between the structure and dynamics in glass-forming liquids is presented. We show how the origin of spatially localized excitations, as defined by the dynamical facilitation (DF) theory, can be understood from a structure-based framework. This framework is constructed by associating excitation events in the DF theory to hopping events between energy minima in the potential energy landscape (PEL). By reducing the PEL to an equal energy well picture and applying a harmonic approximation, we develop a field theory to describe elastic fluctuations about inherent states, which are energy minimizing configurations of the PEL. We model an excitation as a shear transformation zone (STZ) inducing a localized pure shear deformation onto an inherent state. We connect STZs to T1 transition events that break the elastic bonds holding the local structure of an inherent state. A formula for the excitation energy barrier, denoted as J σ , is obtained as a function of inherent-state elastic moduli and the radial distribution function. The energy barrier from the current theory is compared to the one predicted by the DF theory where good agreement is found in various two-dimensional continuous poly-disperse atomistic models of glass formers. These results strengthen the role of structure and elasticity in driving glassy dynamics through the creation and relaxation of localized excitations.

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

Medium-density amorphous ice unveils shear rate as a new dimension in water’s phase diagram

Recent experiments revealed a new amorphous ice phase, medium-density amorphous ice (MDA), formed by ball-milling ice I h at 77 K [Rosu-Finsen et al., Science 379, 474–478 (2023)]. MDA has density between that of low-density amorphous (LDA) and high-density amorphous (HDA) ices, adding to the complexity of water’s phase diagram, known for its glass polyamorphism and two-state thermodynamics. The nature of MDA and its relation to other amorphous ices and liquid water remain unsolved. Here, we use molecular simulations under controlled pressure and shear rate at 77 K to produce and investigate MDA. Here. we find that MDA formed at constant shear rate is a steady-state nonequilibrium shear-driven amorphous ice (SDA), that can be produced by shearing ice I h , LDA, or HDA. Our results suggest that MDA could be obtained by ball-milling water glasses without crystallization interference. Increasing the shear rate at ambient pressure produces SDAs with densities ranging from LDA to HDA, revealing shear rate as a new thermodynamic variable in the nonequilibrium phase diagram of water. Indeed, shearing provides access to amorphous states inaccessible by controlling pressure and temperature alone. SDAs produced with shearing rates as high as 10 6 s −1 sample the same region of the potential energy landscape than hyperquenched glasses with identical density, pressure, and temperature. Intriguingly, SDAs obtained by shearing at ~10 8 s −1 have density, enthalpy, and structure indistinguishable from those of water “instantaneously” quenched from room temperature to 77 K over 10 ps, making them good approximants for the “true glass” of ambient liquid water.

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

Minimal cyclic behavior in sheared amorphous solids

Although jammed packings of soft spheres exist in potential energy landscapes with a vast number of minima, when subjected to cyclic shear they may revisit the same configurations repeatedly. Simple hysteretic spin models, in which particle rearrangements are represented by interacting spin flips called hysterons, capture many features of this periodic behavior. Yet it has been unclear to what extent individual rearrangements can be described by such binary objects and how such objects interact with one another. Using a particularly sensitive algorithm, we identify rearrangements in simulated jammed packings and select pairs of rearrangements that undo one another to create periodic cyclic behavior. We find that the rearrangement pairs surprisingly persist down to the smallest increments in strain, even in the smallest systems we can study. We explore the statistics of these rearrangement pairs and find that there is a relation between the amount of hysteresis and the energy drop and mean-square displacement of the particles; these results are inconsistent with the scaling found in models that treat rearrangements as localized buckling events. Finally, our analysis shows that there is no clean distinction between the particle motions that represent the identity of a single, individual rearrangement and the particle motions that lead to interactions between separated rearrangements or hysterons. These results offer insight into how complex systems such as amorphous solids can reach a limit cycle.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Formation of transfermium elements in reactions with Pb 208

Within the Langevin framework, we investigate the dynamics of the fusion process for production of transfermium elements in reactions of Ca 48 , Ti 50 , Cr 54 , and Fe 58 with Pb 208 . After the reacting nuclei have made contact, the early dynamical stage is dominated by the dissipation of the initial radial kinetic energy, while the subsequent shape evolution is diffusive. The probability for surmounting the inner barrier and forming a compound system is obtained by simulating the evolution as a Metropolis random walk in a five-dimensional potential-energy landscape. Good agreement with the available data is obtained, especially for the maximal formation probability. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗