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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2026

Investigation of the Effect of Framework Flexibility on Adsorption in SIFSIX-3-Cu using a Machine-Learned Force Field

Metal-organic frameworks (MOFs) are a promising class of adsorbents. The performance of MOF sorbents relies on high selectivity and low regeneration energy. This work focuses on the use of machine learned force fields (MLFFs) to model adsorption in a flexible MOF, SIFSIX-3-Cu. A DeePMD-based MLFF was trained to reproduce DFT (PBE+D3) energies, forces, and stresses, using an iterative sampling scheme combining sampling based on molecular dynamics, Monte Carlo, and geometry optimization to capture both attractive and repulsive regions of the potential energy surface. Flexibility of the MOF was explicitly included in this model. Hybrid Monte Carlo/molecular dynamics (MC/MD) simulations using the MLFF predicted adsorption isotherms in good agreement with experimental data for a range of pressures (40 Pa – 104 Pa) in contrast to rigid models, which overpredict CO2 adsorption at low pressures. The improvement was the result of a description of the variability of fluorine-fluorine diagonal distances at adsorption sites. This detailed description of flexibility afforded by the MLFF resulted in more accurate predictions adsorption isotherms when compared to the experimentally measured values. These results underscore the importance of including framework flexibility when modeling adsorption phenomena in MOFs, particularly for low pressure applications and provide a robust procedure for training MLFF models for MOFs.

Atomistic Simulation↗
DOE OSTI2022

Inefficient intramolecular vibrational energy redistribution for the H + HO2 reaction and negative internal energy dependence for its rate constant

Quasiclassical trajectories (QCT) and newly constructed global potential energy surfaces are used to compute thermal and nonthermal rate constants for the H + HO 2 reaction. The thermal QCTs rate constants are up to 50% smaller than transition state theory (TST) rate constants based on the same level of electronic structure theory. This reduction is demonstrated to result from inefficient intramolecular vibrational energy redistribution (IVR) in the transient H 2 O 2 well, with a significant fraction of trajectories that reach the H 2 O 2 well promptly dissociating back to reactants instead of via the heavily statistically favored 2OH channel. The nonstatistical reduction factor, κ IVR , that quantifies this effect is shown to increase in importance with temperature, with κ IVR = 0.81 at 300 K and 0.47 at 2500 K. Finally, we show that inefficient IVR causes H + HO 2 rate constants mediated by H 2 O 2 to depend inversely on the initial vibrational excitation of HO 2 .

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

Hybrid Quantum-Classical Neural Network for Calculating Ground State Energies of Molecules

We present a hybrid quantum-classical neural network that can be trained to perform electronic structure calculation and generate potential energy curves of simple molecules. The method is based on the combination of parameterized quantum circuits and measurements. With unsupervised training, the neural network can generate electronic potential energy curves based on training at certain bond lengths. To demonstrate the power of the proposed new method, we present the results of using the quantum-classical hybrid neural network to calculate ground state potential energy curves of simple molecules such as H 2 , LiH, and BeH 2 . The results are very accurate and the approach could potentially be used to generate complex molecular potential energy surfaces.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

A Theoretical Study of NH 2 Radical Reactions with Propane and Its Kinetic Implications in NH 3 -Propane Blends’ Oxidation

The reaction of NH 2 radicals with C 3 H 8 is crucial for understanding the combustion behavior of NH 3 /C 3 H 8 blends. In this study, we investigated the temperature dependence of the rate coefficients for the hydrogen abstraction reactions of C 3 H 8 by NH 2 radicals using high-level theoretical approaches. The potential energy surface was constructed at the CCSD(T)/cc-pV(T, Q)//M06-2X/aug-cc-pVTZ level of theory, and the rate coefficients were computed using conventional transition state theory, incorporating the corrections for quantum tunneling and hindered internal rotors (HIR). The computed rate coefficients showed a strong curvature in the Arrhenius behavior, capturing the experimental literature data well at low temperatures. However, at T > 1500 K, the theory severely overpredicted the experimental data. The available theoretical studies did not align with the experiment at high temperatures, and the possible reasons for this discrepancy are discussed. At 300 K, the reaction of NH 2 with C 3 H 8 predominantly occurs at the secondary C-H site, which accounts for approximately 95% of the total reaction flux. However, the hydrogen abstraction reaction at the primary C-H site becomes the dominant reaction above 1700 K. A composite kinetic model was built, which incorporated the computed rate coefficients for NH 2 + C 3 H 8 reactions. The importance of NH 2 + C 3 H 8 reactions in predicting the combustion behavior of NH 3 /C 3 H 8 blends was demonstrated by kinetic modeling.

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

Raw Data for: Ultra-Confined Environments May Restrict the Possible Configurations of Supported Metal Complexes

Raw NMR data (1D 2H spectrum and pseudo-2D 2H CODEX spectra) in Bruker topspin format. Raw DRIFTS IR data as text files. Results from the DFT potential energy surface scans (coordinates are in the supporting information of the main publication). Raw BET nitrogen physisorption isotherms and the determined more size distributions.

Perras, Frederic A [Ames National Laboratory]↗
NASA NTRS1972

Dynamics of the molecular and atomic mechanisms for the hydrogen-iodine exchange reaction.

The molecular and atomic mechanisms for the hydrogen-iodine exchange reaction are treated theoretically by means of extensive classical trajectories calculated on a reasonable potential energy surface on which the single adjustable parameter is the iodine-core effective charge. The analysis shows the molecular mechanism to be dynamically forbidden, but gives an over-all rate constant for the atomic mechanism that is in agreement with the experimental values. It is indicated that the formation of a weak H2I complex plays an important dynamical role if the atomic mechanism is limited to reactions with collision complexes involving no more than two hydrogen atoms and two iodine atoms. Excellent agreement with experiment is obtained for the rate constant for the recombination I+I+H2 yields I2+H2 and its negative temperature coefficient.

Raff, L. M.↗
NASA NTRS1972

Stabilities of nitrogen containing heterocyclic radicals and geometrical influences on non-radiative processes in organic molecules

Stabilities of nitrogen containing heterocyclic radicals were studied to detect radicals of the type R-N-R, and to theoretically rationalize their electronic structure. The computation of simple potential energy surfaces for ground and excited states is discussed along with the photophysical properties of indolizine. Methods of calculation and problems associated with the calculations are presented. Results, tables, diagrams, discussions, and references are included.

Evleth, E. M.↗
NASA NTRS1977

Rotational and vibrational transitions for Li + H2 collisions

Close coupling calculations for integral and differential cross sections have been carried out for Li + H2 collisions with an ab initio Hartree-Fock potential energy surface. Rotational, vibrational, and vib-rotational excitation cross sections are reported at 0.4336 eV, 0.7 eV, and 0.8673 eV in the center of mass system. For pure rotational excitations, which dominate the inelastic scattering, coupling with vibrational states is not very important. For vibrational transitions, the influence of large multiquantum rotational transitions is far less than that found for Li(+) + H2 collisions.

Choi, B. H.↗
NASA NTRS1978

On the possibility of negative activation energies in bimolecular reactions

The temperature dependence of the rate constants for model reacting systems was studied to understand some recent experimental measurements which imply the existence of negative activation energies. A collision theory model and classical trajectory calculations are used to demonstrate that the reaction probability can vary inversely with collision energy for bimolecular reactions occurring on attractive potential energy surfaces. However, this is not a sufficient condition to ensure that the rate constant has a negative temperature dependence. On the basis of these calculations, it seems unlikely that a true bimolecular reaction between neutral molecules will have a negative activation energy.

Jaffe, R. L.↗
NASA NTRS1978

Theoretical study of the photodissociation of HOCl

A study of the UV photodissociation of hypochlorous acid is conducted on the basis of ab initio SCF-Cl calculations. These calculations show that HOCl has only a single peak in the UV photoabsorption spectrum at 220 nm. This result implies that HOCl would have a long lifetime for photodissociation if it were to be formed in the stratosphere. The photodissociation products of HOCl have been identified as Cl+OH, based on an examination of the topographies of the excited electronic state potential energy surfaces. The results of this study indicate that HOCl could be a significant reservoir for stratospheric chlorine.

Jaffe, R. L.↗