Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “potential energy surface”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Chemistry of Sugar Formation in the Gas Phase: Following the Activated Aldehyde

Sugars are produced by living organisms, and are required building blocks for life as we know it, which raises the foundational question of how sugars formed in a prebiotic environment. The abiotic formose reaction produces sugars from formaldehyde, but our understanding of its initiation step remains murky, with chemists invoking the concept of an “activated aldehyde” to seed this reaction. Singlet hydroxycarbenes, high-energy isomers of aldehydes, were recently reported to facilitate sugar formation under cold, nonaqueous conditions relevant to interstellar environments. Here, we generate singlet methylhydroxycarbene ( 1 CH 3 –C̈–OH) from the photodissociation of pyruvic acid and experimentally measure its gas-phase reaction with d 4 -acetaldehyde using multiplexed photoionization mass spectrometry. The C 4 H 4 D 4 O 2 isomer d 4 -acetoin is the sole product, which we kinetically link to the reactant CH 3 –C̈–OH, and attribute to a carbonyl-ene formation mechanism. We see no evidence of 3-hydroxybutanal, the C–H insertion product expected in carbene chemistry. Using automated exploration we calculate stationary points on the potential energy surface and report master equation rate coefficients from T = 20–600 K, providing quantitative kinetics of this fast reaction for use in chemical models. The prereactive complex in this reaction is stabilized by both hydrogen bonding and electrophilic carbene C═O interactions. These effects create a short-range dynamical bottleneck for the reaction besides the long- and midrange barrierless bottlenecks. Combined with recent reports of 1 HC̈OH production from methanol photodissociation and pyruvic acid production in cold irradiated ices, this work provides evidence that singlet hydroxycarbene + aldehyde chemistry is a feasible path to prebiotic sugar formation.

aldehydes↗
DOE OSTI2025

Acid-Catalyst is Required for Hydrogenolysis of a Late Metal Hydroxo Complex

The cationic AuIII pincer complex [(tBuPCP)AuIII -OH]OTf (AuOH, tBuPCP = 2,6-(CH2PtBu2)2C6H3) was prepared and characterized. In contrast to hydrogenolysis reactions of other late metal-hydroxides, the hydrogenolysis of AuOH to form the corresponding AuIII -H and water is shown to require an acid catalyst. The kinetic data, collected in acetone-d6, are consistent with a reaction mechanism involving initial protonation to form the AuIII aquo complex, [(tBuPCP)AuIII -OH2](OTf)2 (AuOH2), which was isolated and characterized. Following substitution of the water by coordinating solvent, the dicationic AuIII center undergoes electrophilic activation of H2 to generate the [(tBuPCP)AuIII -H]OTf product (AuH). The mechanism of hydrogenolysis is supported by experimental and computational (DFT) results. A AuIII dihydrogen complex (Au(H2)) was located along the potential energy surface for the reaction and computations indicate that deprotonation of Au(H2) by acetone solvent is barrierless. The bonding and properties of the proposed AuIII -H2 complex are consistent with predictions made for high-valent, electrophilic metal-dihydrogen complexes.

Ardon, Yotam [University of Pennsylvania] (ORCID:0↗
DOE OSTI2026

Ion-Electron Coupling-Driven Redox Behavior in Metal–Organic Frameworks

Redox-active metal–organic frameworks (MOFs) have long been proposed as electronic transport platforms, yet the microscopic origin of their conductivity remains debated. A theoretical demonstration reveals charge transport in a Zn(pyrazole–naphthalene diimide (NDI)) MOF arising not from delocalized band-like states but from redox hopping between discrete linker sites. Using ab initio molecular dynamics simulations combined with electronic structure analysis, we established a direct link among electron injection, structural reorganization, and transport. Electron accumulation proceeds sequentially and site-selectively from imide and carbonyl groups of the NDI core progressively involving pyrazole N atoms at higher reduction states, through a hierarchy of redox-active sites. In contrast, Zn nodes remain essentially redox-inactive, which confirms their structural role. Density-of-states analysis corroborates a transport regime dominated by linker-centered states with evolving p-character upon reduction, resulting in dynamically reconfigured conduction networks. Real-time trajectories reveal anisotropic linker-to-linker electron transfer modulated by counterion coordination. This cooperative ion–electron regime emerges from potential energy surface collapse into a single low-barrier transition (ΔG ‡ ≈ 45 meV), where ionic and electronic motions evolve adiabatically on the same free-energy landscape. Elucidating redox conductivity in Zn(pyrazole–NDI) MOFs provides a theoretical framework for use in neuromorphic computing and related technologies.

Charge transfer↗
DOE OSTI2021

Nonequilibrium charge-density-wave order beyond the thermal limit

The interaction of many-body systems with intense light pulses may lead to novel emergent phenomena far from equilibrium. Recent discoveries, such as the optical enhancement of the critical temperature in certain superconductors and the photo-stabilization of hidden phases, have turned this field into an important research frontier. Here, we demonstrate nonthermal charge-density-wave (CDW) order at electronic temperatures far greater than the thermodynamic transition temperature. Using time- and angle-resolved photoemission spectroscopy and time-resolved X-ray diffraction, we investigate the electronic and structural order parameters of an ultrafast photoinduced CDW-to-metal transition. Tracking the dynamical CDW recovery as a function of electronic temperature reveals a behaviour markedly different from equilibrium, which we attribute to the suppression of lattice fluctuations in the transient nonthermal phonon distribution. A complete description of the system’s coherent and incoherent order-parameter dynamics is given by a time-dependent Ginzburg-Landau framework, providing access to the transient potential energy surfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Chiral photochemistry of achiral molecules

Chirality is a molecular property governed by the topography of the potential energy surface (PES). Thermally achiral molecules interconvert rapidly when the interconversion barrier between the two enantiomers is comparable to or lower than the thermal energy, in contrast to thermally stable chiral configurations. In principle, a change in the PES topography on the excited electronic state may diminish interconversion, leading to electronically prochiral molecules that can be converted from achiral to chiral by electronic excitation. Here we report that this is the case for two prototypical examples – cis-stilbene and cis-stiff stilbene. Both systems exhibit unidirectional photoisomerization for each enantiomer as a result of their electronic prochirality. We simulate an experiment to demonstrate this effect in cis-stilbene based on its interaction with circularly polarized light. Our results highlight the drastic change in chiral behavior upon electronic excitation, opening up the possibility for asymmetric photochemistry from an effectively nonchiral starting point.

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

Step-by-step state-selective tracking of fragmentation dynamics of water dications by momentum imaging

Abstract The double photoionization of a molecule by one photon ejects two electrons and typically creates an unstable dication. Observing the subsequent fragmentation products in coincidence can reveal a surprisingly detailed picture of the dynamics. Determining the time evolution and quantum mechanical states involved leads to deeper understanding of molecular dynamics. Here in a combined experimental and theoretical study, we unambiguously separate the sequential breakup via D + + OD + intermediates, from other processes leading to the same D + + D + + O final products of double ionization of water by a single photon. Moreover, we experimentally identify, separate, and follow step by step, two pathways involving the b 1 Σ + and a 1 Δ electronic states of the intermediate OD + ion. Our classical trajectory calculations on the relevant potential energy surfaces reproduce well the measured data and, combined with the experiment, enable the determination of the internal energy and angular momentum distribution of the OD + intermediate.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Architector for high-throughput cross-periodic table 3D complex building

Abstract Rare-earth and actinide complexes are critical for a wealth of clean-energy applications. Three-dimensional (3D) structural generation and prediction for these organometallic systems remains a challenge, limiting opportunities for computational chemical discovery. Here, we introduce Architector , a high-throughput in-silico synthesis code for s-, p-, d-, and f-block mononuclear organometallic complexes capable of capturing nearly the full diversity of the known experimental chemical space. Beyond known chemical space, Architector performs in-silico design of new complexes including any chemically accessible metal-ligand combinations. Architector leverages metal-center symmetry, interatomic force fields, and tight binding methods to build many possible 3D conformers from minimal 2D inputs including metal oxidation and spin state. Over a set of more than 6,000 x-ray diffraction (XRD)-determined complexes spanning the periodic table, we demonstrate quantitative agreement between Architector-predicted and experimentally observed structures. Further, we demonstrate out-of-the box conformer generation and energetic rankings of non-minimum energy conformers produced from Architector , which are critical for exploring potential energy surfaces and training force fields. Overall, Architector represents a transformative step towards cross-periodic table computational design of metal complex chemistry.

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

Aldehyde cool-flame chemistry explains a missing source of organic acids

Combustion emission is a significant source of organic acids, impacting atmospheric chemistry and climate. Their formation mechanisms, however, remain poorly understood, leading to underestimation in kinetic models. We investigate the cool-flame oxidation of key combustion intermediates—C 1 –C 4 aldehydes and benzaldehyde. Using in-situ synchrotron vacuum ultraviolet photoionization mass spectrometry, we observe the direct conversion of aldehydes to organic acids, a process enhanced by HO 2 radicals. Quantum chemistry calculations reveal that the reaction of RC(O)O 2 with HO 2 on the singlet potential energy surface contributes to organic acids. Incorporating this pathway into a kinetic model significantly improves organic acid prediction. Despite the high-temperature nature of engine combustion, significant spatial and temporal inhomogeneities (e.g., near-wall regions and crevice volumes) lead to localized cool-flame conditions, facilitating organic acid formation and emission. Elucidating the acid formation under cool-flame conditions provides a critical mechanism for accurately modelling anthropogenic organic acid emissions and developing mitigation strategies.

SVUV-PIMS↗
DOE OSTI2022

Vibrationally resolved optical excitations of the nitrogen-vacancy center in diamond

A comprehensive description of the optical cycle of spin defects in solids requires the understanding of the electronic and atomistic structure of states with different spin multiplicity, including singlet states which are particularly challenging from a theoretical standpoint. We present a general framework, based on spin-flip time-dependent density function theory, to determine the excited state potential energy surfaces of the many-body singlet states of spin defects; we then predict the vibrationally resolved absorption spectrum between singlet shelving states of a prototypical defect, the nitrogen-vacancy center in diamond. Our results, which are in very good agreement with experiments, provide an interpretation of the measured spectra and reveal the key role of specific phonons in determining absorption processes, and the notable influence of non-adiabatic interactions. The insights gained from our calculations may be useful in defining strategies to improve infrared-absorption-based magnetometry and optical pumping schemes. The theoretical framework developed here is general and applicable to a variety of other spin defects and materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Complexity of many-body interactions in transition metals via machine-learned force fields from the TM23 data set

Abstract This work examines challenges associated with the accuracy of machine-learned force fields (MLFFs) for bulk solid and liquid phases ofd-block elements. In exhaustive detail, we contrast the performance of force, energy, and stress predictions across the transition metals for two leading MLFF models: a kernel-based atomic cluster expansion method implemented using sparse Gaussian processes (FLARE), and an equivariant message-passing neural network (NequIP). Early transition metals present higher relative errors and are more difficult to learn relative to late platinum- and coinage-group elements, and this trend persists across model architectures. Trends in complexity of interatomic interactions for different metals are revealed via comparison of the performance of representations with different many-body order and angular resolution. Using arguments based on perturbation theory on the occupied and unoccupieddstates near the Fermi level, we determine that the large, sharpddensity of states both above and below the Fermi level in early transition metals leads to a more complex, harder-to-learn potential energy surface for these metals. Increasing the fictitious electronic temperature (smearing) modifies the angular sensitivity of forces and makes the early transition metal forces easier to learn. This work illustrates challenges in capturing intricate properties of metallic bonding with current leading MLFFs and provides a reference data set for transition metals, aimed at benchmarking the accuracy and improving the development of emerging machine-learned approximations.

Chemistry↗
DOE OSTI2025

Systematic softening in universal machine learning interatomic potentials

Machine learning interatomic potentials (MLIPs) have introduced a new paradigm for atomic simulations. Recent advancements have led to universal MLIPs (uMLIPs) that are pre-trained on diverse datasets, providing opportunities for universal force fields and foundational machine learning models. However, their performance in extrapolating to out-of-distribution complex atomic environments remains unclear. In this study, we highlight a consistent potential energy surface (PES) softening effect in three uMLIPs: M3GNet, CHGNet, and MACE-MP-0, which is characterized by energy and force underprediction in atomic-modeling benchmarks including surfaces, defects, solid-solution energetics, ion migration barriers, phonon vibration modes, and general high-energy states. The PES softening behavior originates primarily from the systematically underpredicted PES curvature, which derives from the biased sampling of near-equilibrium atomic arrangements in uMLIP pre-training datasets. Our findings suggest that a considerable fraction of uMLIP errors are highly systematic, and can therefore be efficiently corrected. We argue for the importance of a comprehensive materials dataset with improved PES sampling for next-generation foundational MLIPs.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Cross-functional transferability in foundation machine learning interatomic potentials

The rapid development of foundation potentials (FPs) in machine learning interatomic potentials demonstrates the possibility for generalizable learning of the universal potential energy surface. The accuracy of FPs can be further improved by bridging the model from lower-fidelity datasets to high-fidelity ones. In this work, we analyze the challenge of this transfer learning (TL) problem within the CHGNet framework. We show that significant energy scale shifts and poor correlations between GGA and r 2 SCAN hinder cross-functional transferability. By benchmarking different TL approaches on the MP-r 2 SCAN dataset, we demonstrate the importance of elemental energy referencing in the TL of FPs. By comparing the scaling law with and without the pre-training on a low-fidelity dataset, we show that significant data efficiency can still be achieved through TL, even with a target dataset of sub-million structures. We highlight the importance of proper TL and multi-fidelity learning in creating next-generation FPs on high-fidelity data.

Huang, Xu [University of California, Berkeley, CA ↗
DOE OSTI2020

Reactants, products, and transition states of elementary chemical reactions based on quantum chemistry

Reaction times, activation energies, branching ratios, yields, and many other quantitative attributes are important for precise organic syntheses and generating detailed reaction mechanisms. Often, it would be useful to be able to classify proposed reactions as fast or slow. However, quantitative chemical reaction data, especially for atom-mapped reactions, are difficult to find in existing databases. Therefore, we used automated potential energy surface exploration to generate 12,000 organic reactions involving H, C, N, and O atoms calculated at the ω B97X-D3/def2-TZVP quantum chemistry level. We report the results of geometry optimizations and frequency calculations for reactants, products, and transition states of all reactions. Additionally, we extracted atom-mapped reaction SMILES, activation energies, and enthalpies of reaction. We believe that this data will accelerate progress in automated methods for organic synthesis and reaction mechanism generation—for example, by enabling the development of novel machine learning models for quantitative reaction prediction.

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

Dissociation dynamics of the diamondoid adamantane upon photoionization by XUV femtosecond pulses

This work presents a photodissociation study of the diamondoid adamantane using extreme ultraviolet femtosecond pulses. The fragmentation dynamics of the dication is unraveled by the use of advanced ion and electron spectroscopy giving access to the dissociation channels as well as their energetics. To get insight into the fragmentation dynamics, we use a theoretical approach combining potential energy surface determination, statistical fragmentation methods and molecular dynamics simulations. We demonstrate that the dissociation dynamics of adamantane dications takes place in a two-step process: barrierless cage opening followed by Coulomb repulsion-driven fragmentation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Ultraviolet spectroscopy of pressurized and supercritical carbon dioxide

Abstract Carbon dioxide (CO 2 ) is prevalent in planetary atmospheres and sees use in a variety of industrial applications. Despite its ubiquitous nature, its photochemistry remains poorly understood. In this work we explore the density dependence of pressurized and supercritical CO 2 electronic absorption spectra by vacuum ultraviolet spectroscopy over the wavelength range 1455-2000 Å. We show that the lowest absorption band transition energy is unaffected by a density increase up to and beyond the thermodynamic critical point (137 bar, 308 K). However, the diffuse vibrational structure inherent to the spectrum gradually decreases in magnitude. This effect cannot be explained solely by collisional broadening and/or dimerization. We suggest that at high densities close proximity of neighboring CO 2 molecules with a variety of orientations perturbs the multiple monomer electronic state potential energy surfaces, facilitating coupling between binding and dissociative states. We estimate a critical radius of ~4.1 Å necessary to cause such perturbations.

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

CHGNet as a pretrained universal neural network potential for charge-informed atomistic modelling

Large-scale simulations with complex electron interactions remain one of the greatest challenges for atomistic modelling. Although classical force fields often fail to describe the coupling between electronic states and ionic rearrangements, the more accurate ab initio molecular dynamics suffers from computational complexity that prevents long-time and large-scale simulations, which are essential to study technologically relevant phenomena. Here we present the Crystal Hamiltonian Graph Neural Network (CHGNet), a graph neural network-based machine-learning interatomic potential (MLIP) that models the universal potential energy surface. CHGNet is pretrained on the energies, forces, stresses and magnetic moments from the Materials Project Trajectory Dataset, which consists of over 10 years of density functional theory calculations of more than 1.5 million inorganic structures. The explicit inclusion of magnetic moments enables CHGNet to learn and accurately represent the orbital occupancy of electrons, enhancing its capability to describe both atomic and electronic degrees of freedom. We demonstrate several applications of CHGNet in solid-state materials, including charge-informed molecular dynamics in Li x MnO 2 , the finite temperature phase diagram for Li x FePO 4 and Li diffusion in garnet conductors. We highlight the significance of charge information for capturing appropriate chemistry and provide insights into ionic systems with additional electronic degrees of freedom that cannot be observed by previous MLIPs.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Uncertainty-driven dynamics for active learning of interatomic potentials

Machine learning (ML) models, if trained to data sets of high-fidelity quantum simulations, produce accurate and efficient interatomic potentials. Active learning (AL) is a powerful tool to iteratively generate diverse data sets. In this approach, the ML model provides an uncertainty estimate along with its prediction for each new atomic configuration. If the uncertainty estimate passes a certain threshold, then the configuration is included in the data set. Here we develop a strategy to more rapidly discover configurations that meaningfully augment the training data set. The approach, uncertainty-driven dynamics for active learning (UDD-AL), modifies the potential energy surface used in molecular dynamics simulations to favor regions of configuration space for which there is large model uncertainty. The performance of UDD-AL is demonstrated for two AL tasks: sampling the conformational space of glycine and sampling the promotion of proton transfer in acetylacetone. The method is shown to efficiently explore the chemically relevant configuration space, which may be inaccessible using regular dynamical sampling at target temperature conditions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Theoretical study of the reaction mechanism and kinetics of the phenyl + propargyl association

Potential energy surface for the phenyl + propargyl radical recombination reaction has been studied at the CCSD(T)-F12/cc-pVTZ-f12//B3LYP/6-311G** level of theory for the closed-shell singlet species and at the triplet–singlet gap CASPT2/cc-pVTZ-CCSD(T)-F12/cc-pVTZ-f12//CASSCF/cc-pVTZ level of theory for the diradical species. High-pressure limit rate constants for the barrierless channels were evaluated with variable reaction coordinate transition state theory (VRC-TST). Rice–Ramsperger–Kassel–Marcus Master Equation (RRKM-ME) calculations have been performed to assess temperature- and pressure-dependent phenomenological rate constants and product branching ratios. Here, the entrance channels of the radical association reaction produce 3-phenyl-1-propyne and phenylallene which can further dissociate/isomerize into a variety of unimolecular and bimolecular products. Theoretical evidence is presented that, at combustion relevant conditions, the phenyl + propargyl recombination provides a feasible mechanism for the addition of a second five-member ring to the first six-member aromatic ring producing the prototype two-ring species indene and indenyl. Rate expressions for all important reaction channels in a broad range of temperatures and pressures have been generated for kinetic modeling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗