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 OSTI2023

Development of a scalar-based geometric parameterization approach for the crystal structure landscape of dithienylethene-based crystalline solids

Dithienylethenes (DTEs) are a promising class of organic photoswitches that can be used to create crystalline solids with properties controlled by light. However, the ability of DTEs to adopt multiple conformations, only one of which is photoactive, complicates the rational design of these materials. Herein, the synthesis and structural characterization of 19 crystalline solids containing a single DTE molecule are described. A novel D – D analysis of the molecular geometries obtained from rotational potential energy surface calculations and the ensemble of experimental structures were used to construct a crystal landscape for DTE. Of the 19 crystal structures, 17 contained photoinactive DTE rotamers and only 2 were photoactive. These results highlight the challenges associated with the design of these materials. Overall, the D – D analysis described herein provides rapid, effective and intuitive means of linking the molecular structure to photoactivity that could be applied more broadly to afford a general strategy for producing photoactive diarylethene-based crystalline solids.

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

Intermolecular vibrations mediate ultrafast singlet fission

Singlet fission is a spin-allowed exciton multiplication process in organic semiconductors that converts one spin-singlet exciton to two triplet excitons. It offers the potential to enhance solar energy conversion by circumventing the Shockley-Queisser limit on efficiency. We study the primary steps of singlet fission in a pentacene film by using a combination of TG and 2D electronic spectroscopy complemented by quantum chemical and nonadiabatic dynamics calculations. We show that the coherent vibrational dynamics induces the ultrafast transition from the singlet excited electronic state to the triplet-pair state via a degeneracy of potential energy surfaces, i.e., a multidimensional conical intersection. Significant vibronic coupling of the electronic wave packet to a few key intermolecular rocking modes in the low-frequency region connect the excited singlet and triplet-pair states. Along with high-frequency local vibrations acting as tuning modes, they open a new channel for the ultrafast exciton transfer through the resulting conical intersection.

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

Ab initio investigation of a hypersonic double cone experiment

This article presents a direct molecular simulation (DMS) of a reactive Mach 8.2 oxygen flow over a double cone geometry. The free stream conditions and article configuration generate a flow with thermal and chemical nonequilibrium, which are common attributes of hypersonic flight. This scenario was first studied experimentally at Calspan University of Buffalo Research Center’s test facility. DMS is a particle method that uses quantum mechanically derived interaction potentials to simulate molecular collisions within a flow field. Since these interaction potentials are the only modeling inputs used in the simulation, all flow features can solely be attributed to the ab initio potential energy surfaces. Hence, providing a comparison of a hypersonic ground test and numerical data anchored to quantum mechanics. The fundamental nature of DMS is leveraged to investigate molecular level mechanisms prevalent in the flow, and comparisons with lower fidelity simulations are presented to highlight the role of these first principles calculations as benchmark solutions.

Science & Technology - Other Topics↗
DOE OSTI2020

A Method for Dimensionally Adaptive Sparse Trigonometric Interpolation of Periodic Functions

We present a method for dimensionally adaptive sparse trigonometric interpolation of multidimensional periodic functions belonging to a smoothness class of finite order. This method targets applications where periodicity must be preserved and the precise anisotropy is not known a priori. To the authors' knowledge, this is the first instance of a dimensionally adaptive sparse interpolation algorithm that uses a trigonometric interpolation basis. The motivating application behind this work is the adaptive approximation of a multi-input model for a molecular potential energy surface (PES) where each input represents an angle of rotation. Our method is based on an anisotropic quasi-optimal estimate for the decay rate of the Fourier coefficients of the model; a least-squares fit to the coefficients of the interpolant is used to estimate the anisotropy. Thus, our adaptive approximation strategy begins with a coarse isotropic interpolant, which is gradually refined using the estimated anisotropic rates. The procedure takes several iterations where ever-more accurate interpolants are used to generate ever-improving anisotropy rates. We present several numerical examples of our algorithm where the adaptive procedure successfully recovers the theoretical “best” convergence rate, including an application to a periodic PES approximation. An open-source implementation of our algorithm resides in the Tasmanian UQ library developed at Oak Ridge National Laboratory.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Many-Body Effects in Aqueous Systems: Synergies Between Interaction Analysis Techniques and Force Field Development

Interaction analysis techniques, including the many-body expansion (MBE), symmetry-adapted perturbation theory, and energy decomposition analysis, allow for an intuitive understanding of complex molecular interactions. We review these methods by first providing a historical context for the study of many-body interactions and discussing how nonadditivities emerge from Hamiltonians containing strictly pairwise-additive interactions. We then elaborate on the synergy between these interaction analysis techniques and the development of advanced force fields aimed at accurately reproducing the Born–Oppenheimer potential energy surface. In particular, we focus on ab initio–based force fields that aim to explicitly reproduce many-body terms and are fitted to high-level electronic structure results. These force fields generally incorporate many-body effects through ( a) parameterization of distributed multipoles, ( b) explicit fitting of the MBE, ( c) inclusion of many-atom features in a neural network, and ( d) coarse-graining of many-body terms into an effective two-body term. We also discuss the emerging use of the MBE to improve the accuracy and speed of ab initio molecular dynamics.

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

Spectroscopy and Scattering Studies Using Interpolated Ab Initio Potentials

The Born–Oppenheimer potential energy surface (PES) has come a long way since its introduction in the 1920s, both conceptually and in predictive power for practical applications. Nevertheless, nearly 100 years later—despite astonishing advances in computational power—the state-of-the-art first-principles prediction of observables related to spectroscopy and scattering dynamics is surprisingly limited. For example, the water dimer, (H 2 O) 2 , with only six nuclei and 20 electrons, still presents a formidable challenge for full-dimensional variational calculations of bound states and is considered out of reach for rigorous scattering calculations. The extremely poor scaling of the most rigorous quantum methods is fundamental; however, recent progress in development of approximate methodologies has opened the door to fairly routine high-quality predictions, unthinkable 20 years ago. In this review, in relation to the workflow of spectroscopy and/or scattering studies, we summarize progress and challenges in the component areas of electronic structure calculations, PES fitting, and quantum dynamical calculations.

Chemistry↗
DOE OSTI2023

AL-ASMR: Active Learning of Atomistic Surrogate Models for Rare Events

Atomistic simulation with artificial intelligence (AI) is an emerging tool for understanding materials' properties and behaviors and predicting novel materials with optimized/targeted properties. Neural network potentials (NNPs) are outstanding in this field as they have shown a comparable accuracy to ab initio electronic structure calculations for reproducing potential energy surfaces while being several orders of magnitude faster. However, such NNPs can perform poorly outside of their training domain and typically fail catastrophically in the prediction of rare events in molecular dynamics (MD) simulations. For effective AL loops to distinguish the informative data from enhanced sampled configurations, we developed a decision engine by configurational similarity and uncertainty quantification (UQ) with data augmentation.

Jung, Gang Seob↗
DOE OSTI2024

pyKinML

SAND2024-13788O pyKinML software is used to train neural net potential energy surfaces for hydrocarbon molecules. The models study reaction kinetics but can be used for other purposes as well, such as running molecular dynamics simulation. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Najm, Habib↗
DOE OSTI2026

LBNL_DIATOMICS v1

LBNL_DIATOMICS is a software package developed by P. Zarzycki (Lawrence Berkeley National Laboratory) for computing isotopic effects in diatomic molecules using mass-dependent spectroscopic parameters obtained experimentally or from first principles. The package provides a streamlined pipeline from the potential energy surface to the equilibrium constant of an isotopic exchange reaction.

Zarzycki, Piotr [Lawrence Berkeley National Labora↗
DOE OSTI2020

Accelerating Finite-temperature Kohn-Sham Density Functional Theory\ with Deep Neural Networks

We present a numerical modeling workflow based on machine learning (ML) which reproduces the the total energies produced by Kohn-Sham density functional theory (DFT) at finite electronic temperature to within chemical accuracy at negligible computational cost. Based on deep neural networks, our workflow yields the local density of states (LDOS) for a given atomic configuration. From the LDOS, spatially-resolved, energy-resolved, and integrated quantities can be calculated, including the DFT total free energy, which serves as the Born-Oppenheimer potential energy surface for the atoms. We demonstrate the efficacy of this approach for both solid and liquid metals and compare results between independent and unified machine-learning models for solid and liquid aluminum. Our machine-learning density functional theory framework opens up the path towards multiscale materials modeling for matter under ambient and extreme conditions at a computational scale and cost that is unattainable with current algorithms.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Accelerating Finite-Temperature Kohn-Sham Density Functional Theory with Deep Neural Networks

We present a numerical modeling workflow based on machine learning (ML) which reproduces the total energies produced by Kohn-Sham density functional theory (DFT) at finite electronic temperature to within chemical accuracy at negligible computational cost. Based on deep neural networks, our workflow yields the local density of states (LDOS) for a given atomic configuration. From the LDOS, spatially-resolved, energy-resolved, and integrated quantities can be calculated, including the DFT total free energy, which serves as the Born-Oppenheimer potential energy surface for the atoms. We demonstrate the efficacy of this approach for both solid and liquid metals and compare results between independent and unified machine-learning models for solid and liquid aluminum. Our machine-learning density functional theory framework opens up the path towards multiscale materials modeling for matter under ambient and extreme conditions at a computational scale and cost that is unattainable with current algorithms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Simultaneous mitigation of density and energy errors in approximate DFT for transition metal chemistry (Final Technical Report)

There were three major goals and objectives of this project: 1) Develop tools to understand density-driven errors in transition metal complexes, 2) Evaluate and minimize energy delocalization error and static correlation error through judicious functional choice, and 3) applying this workflow to machine-learning accelerated screening of redox couples. Over the reporting period, we developed a framework for eliminating flat plane errors. We introduced fully non-empirical coefficients. We demonstrated the approach on both molecules and solids. We investigated and eliminated density driven errors and demonstrated their impact on potential energy surfaces. We trained machine learning models both in a method-dependent fashion and to predict errors in method accuracy. We built large data sets of small molecule energetics and multi-reference character.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Fission Product Evaluation (Quarterly Report)

The LLNL team has developed a theoretical framework to compute the fission fragment charge and mass distributions of odd-mass nuclei, which can be formed in neutron-induced reactions on even-even nuclei. For example, n+ 238 U leads to the compound nucleus 239 U which has 147 neutrons. Odd-mass and odd-odd nuclei are spin-polarized: the spin of the nucleus is not equal to 0 as in even-even nuclei (a consequence of pairing correlations) but is in fact set by a single particle. In practice, that single particle can occupy different quantum-mechanical configurations, each with different spin. This affects the deformation properties of the nucleus and, as a result, the distribution of fission fragments. These distributions require two inputs for the calculation: a potential energy surface (= the total energy of the nucleus as a function of the deformation of the nuclear shape, often parametrized with 2 degrees of freedom) and the collective inertia tensor associated with each point of the PES. In an odd-mass or odd-odd nucleus, each quantity depends on the configuration of the single particle. We computed of the order of 20 to 30 different such configurations to extract the PES for spins K=1/2, 3/2, 5/2, 7/2, and 9/2 of 239 U. We then computed the primary yields for each such PES, as shown in Fig. 1 below. These yields were weighed with the spin distribution associated with the reaction 238 U+n and the resulting distribution used as input to FREYA to compute the independent yields shown in Fig. 2. These yields are then used as inputs to the FIER code to extract the cumulative yields and compare with LLNL-LANL measurements at TUNL.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Advanced Electronic Structure Theories for Strongly Correlated Ground and Excited States

The aim of this project was to create widely applicable, numerically robust, and systematically improvable multireference theories based on the Driven Similarity Renormalization Group (DSRG). The DSRG is a many-body formalism recently developed in our lab that maps a complex problem involving strongly and weakly interacting electrons to a simpler one in which a few electrons interact strongly via "renormalized" interactions. During this project, we focused our efforts on developing a multireference version of the DSRG (MR-DSRG) applicable to a wide range of problems in theoretical chemistry, including the computation of global ground and excited potential energy surfaces, spin-splittings in transition metal complexes and predicting the interaction of near-degenerate electronically excited states.

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

Shadow molecular dynamics and atomic cluster expansions for flexible charge models

A shadow molecular dynamics scheme for flexible charge models is presented, where the shadow Born-Oppenheimer potential is derived from a coarse-grained approximation of range-separated density functional theory. The interatomic potential, including the atomic electronegativities and the charge-independent short-range part of the potential and force terms, are modeled by the linear atomic cluster expansion (ACE), which provides a computationally efficient alternative to many machine learning methods. The shadow molecular dynamics scheme is based on extended Lagrangian (XL) Born-Oppenheimer molecular dynamics (BOMD) [Eur. Phys. J. B 94, 164 (2021)]. XL-BOMD provides a stable dynamics, while avoiding the costly computational overhead associated with solving an all-to-all system of equations, which normally is required to determine the relaxed electronic ground state prior to each force evaluation. To demonstrate the proposed shadow molecular dynamics scheme for flexible charge models using the atomic cluster expansion, we emulate the dynamics generated from self-consistent charge density functional tight-binding (SCC-DFTB) theory using a second-order charge equilibration (QEq) model. The charge-independent potentials and electronegativities of the QEq model are trained for a supercell of uranium oxide (UO2) and a molecular system of liquid water. The combined ACE + XL-QEq dynamics are stable over a wide range of temperatures both for the oxide and the molecular systems, and provide a precise sampling of the Born-Oppenheimer potential energy surfaces. Accurate ground Coulomb energies are produced by the ACE-based electronegativity model during an NVE simulation of UO 2 , predicted to be within 1 meV of those from SCC-DFTB on average during comparable simulations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Gas-Phase Oxidation Mechanisms of Nitrogen-Containing Organic Compounds

Nitrogen-containing organic molecules play essential roles in nearly every aspect of chemistry. The importance of these species within gas-phase chemistry has increased in recent years, due to their emissions from wildfires, crude biofuels, ammonia combustion, and CO 2 capture facilities. However, there is a lack of detailed relevant mechanistic and kinetic studies of N-compounds, resulting in poor representation of these pathways within chemical models. In this report, we detail our exploration of the gas-phase oxidation mechanisms of pyrrole, imidazole, pyrrolidine, methylamine, and dimethylamine, done by employing multiplexed VUV photoionization mass spectrometry coupled with KinBot, a computational tool that automatically explores multi-well potential energy surfaces. Overall, we vastly expanded the known gas-phase oxidation mechanisms and chemical kinetics of these species. This work will improve the accuracy and completeness of atmospheric and combustion chemistry models, ultimately leading to better-informed decision-makers in the energy and environmental sectors.

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

New Methods for Predicting Non-Born-Oppenheimer Chemistry

Current methods for modeling non-adiabatic molecular dynamics face fundamental limitations when treating geometric phase effects: quantum mechanical phenomena where nuclear wavepackets acquire phase shifts when encircling conical intersections. Existing approaches either neglect these effects entirely or rely on potential energy surfaces arising from the Born-Oppenheimer approximation, which introduce artificial singularities and can overestimate geometric phase contributions. This project developed a new theoretical framework based on exact factorization (XF) methods to overcome these limitations. We derived mathematical formulations for hybrid quantum-classical XF dynamics that selectively treat critical nuclear degrees of freedom quantum mechanically while propagating others classically. This approach addresses the computational intractability that has previously limited exact methods to toy systems. Key innovations include a new approach to systematically identifying nuclear coordinates requiring quantum treatment, as well as novel implementation strategies that interface with existing quantum chemistry codes. The project also developed a proof-of-concept code for treating Jahn-Teller systems and creation of educational materials on non-adiabatic dynamics geared at the graduate level. The theoretical framework developed will enable future systematically improvable calculations of nuclear quantum effects in realistic molecular systems, filling a critical gap in non-adiabatic dynamics methods. This foundation supports future development of predictive tools for designing energy-relevant photochemical processes where quantum coherence effects may be exploited to control reaction outcomes.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

Computationally Guided and Experimentally Validated Design of Custom Chelators for Critical Mineral Recovery

Selective, high throughput separation of target critical metals from complex environments such as fly ash leachates and mining process streams presents a significant challenge for economical production. Custom chelators and sorbents are an attractive technology for selective metal extraction, however it can be difficult to predict their performance, and significant experimental efforts are often required to develop chelating technologies. Here, we present a computational strategy focused on modelling chelator-metal binding interactions and benchmark these results versus experimental data. A computational pipeline combining forcefield, semiempirical, and meta-GGA methods with a thermodynamic framework optimized for error cancellation has been developed to predict binding energies of chelator complexes towards critical mineral recovery applications. This approach, originally validated on [2.2.2] cryptates binding mono- and divalent cations, demonstrated robust predictive capabilities with an R2 of 0.850 against experimental aqueous binding energies. The workflow includes metadynamics for exploring high-dimensional potential energy surfaces and a cluster-continuum model for accurate yet computationally efficient solvation modeling. Error cancellation between solvation energies of free and chelator-coordinated ions enables faster convergence, even with finite cluster sizes. Initial studies on the cryptates revealed consistent metal-ligand coordination patterns, with systematic variations influenced by ion size and charge, highlighting key structural features linked to binding selectivity. Further studies of a proprietary chelator have resulted in identification of previously unreported selectivity towards economically significant metals, which in-house experiments have confirmed, demonstrating the feasibility of this approach. By applying this methodology to new chelators targeting critical minerals such as lithium, cobalt, nickel and other strategic metals, we aim to accelerate the discovery of next-generation chelators for efficient recovery, recycling, and separation processes. This computational framework serves as the backbone of a high-throughput design pipeline tailored for sustainable resource utilization and may be applied to a wide range of systems to meet experimental needs.

computational materials↗