Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular potentials”

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 OSTI2026

Charge accumulation and solvation in $β$-NiOOH: Surface chemistry of an OER catalyst from ML-aided simulations

Electrochemical water splitting is a key technology for a sustainable energy transition, providing a route to store surplus electricity from renewable sources. A central bottleneck is the sluggish oxygen evolution reaction (OER), which drives the search for catalysts that are active, stable, and inexpensive enough for large-scale deployment. Within this context, pure and doped NiO x H y combine high activity with low cost, making them prime candidates for alkaline OER. Yet, despite extensive study, the atomistic structure of NiOOH under operando conditions and the associated reaction mechanisms remain debated. Here, we investigate the structural complexity of pure β-NiOOH, the scaffold for its doped derivatives. We systematically investigate the oxidation of the surface adsorbates via proton-coupled electron transfer steps across relevant facets and sites, identifying the most probable sequence of deprotonation events. Our results reveal asymmetric charge accumulation on Wulff-relevant surfaces and show how applied potential can promote morphological restructuring. Explicit solvation is included through machine-learning interatomic potential molecular dynamics of the NiOOH/water interface, which allows us to resolve the hydrophobic and hydrophilic character of different surfaces and the associated interfacial water structure. Together, these insights demonstrate how surface chemistry and solvation jointly govern the stability of NiOOH and the accumulation of surface charge, with possible implications for catalytic performance.

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

RMCProfile7 : reverse Monte Carlo for multiphase systems

This work introduces a completely rewritten version of the programRMCProfile(version 7), big-box, reverse Monte Carlo modelling software for analysis of total scattering data. The major new feature ofRMCProfile7is the ability to refine multiple phases simultaneously, which is relevant for many current research areas such as energy materials, catalysis and engineering. Other new features include improved support for molecular potentials and rigid-body refinements, as well as multiple different data sets. An empirical resolution correction and calculation of the pair distribution function as a back-Fourier transform are now also available.RMCProfile7is freely available for download at https://rmcprofile.ornl.gov/.

Chemistry↗
DOE OSTI2025

BONCAT-Live for isolation and cultivation of active environmental bacteria

In diverse environments, microbes drive a myriad of processes, from geochemical and nutrient cycling to interspecies interactions, including associations with plants and animals. Their physiological state is dynamic and impacted by abiotic and biotic conditions, responding to environmental fluctuations by changes in cellular metabolism, according to their genetic potential. Molecular, cellular, and genomic approaches can identify and measure microbial responses and adaptation to environmental changes in native communities. However, isolating individual microbial cells that respond to specific changes for cultivation has been difficult. To address this, we implemented a novel bacterial isolation approach (BONCAT-Live) by integrating bio-orthogonal non-canonical amino acid tagging (BONCAT) in diverse native communities, with isolation and cultivation of cells responding to specific stimuli, at different time scales. In frozen Arctic permafrost samples, we identified and isolated dormant bacteria that become active after thawing under native or nutrient-enriched conditions. From the Populus tree rhizosphere, we isolated strains that thrive under high concentrations of root exudates that act as defense compounds and nutrients. In the human microbiome, we identified and isolated bacteria that rapidly proliferated when exposed to metabolites provided by the host or other co-occurring microbes. Further characterization of isolated bacterial strains will provide opportunities for in-depth determination of how these microbes adapt to changes in their environments, individually and as part of model communities.

Analytical Methods↗
DOE OSTI2025

Predicting Fuel Properties and Emissions for Advanced Biofuels for Diesel Engines (CRADA Final Report)

The project will investigate variations in biofuel composition and optimize performance in combustion for conventional and future compression ignition engines. It will evaluate a variety of bio-derived molecules in the diesel range that can be produced using technology in ExxonMobil’s portfolio as well as fuels that cover the range of potential molecular structures for robust model development. Changes in fuel/air premixing and stratification in advanced engines could alter the relationship between fuel properties and performance in comparison to current generation spray combustion approaches. NREL experience in fuel and combustion modeling will enable development of general rules for predicting performance of a wide range of biofuel options.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

A review of displacement cascade simulations using molecular dynamics emphasizing interatomic potentials for TPBAR components

This review explores molecular dynamics simulations for studying radiation damage in Tritium Producing Burnable Absorber Rod (TPBAR) materials, emphasizing the role of interatomic potentials in displacement cascades. Recent machine learning potentials (MLPs), trained on quantum data, enhance prediction accuracy over traditional models like EAM. We highlight temperature, PKA energy, and composition effects on damage evolution in TPBAR components, recommending suitable potentials and discussing advancements for materials in extreme radiation environments.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Hybrid Quantum Mechanical, Molecular Mechanical, and Machine Learning Potential for Computing Aqueous-Phase Adsorption Free Energies on Metal Surfaces

Performing reliable computer simulations of elementary processes occurring at metal–water interfaces is pivotal for novel catalyst design in sustainable energy applications. Computational catalyst design hinges on the ability to reliably and efficiently compute the potential energy surface (PES) of the system. Here, due to the large system sizes needed for studying processes at liquid water–metal interfaces, these systems can currently not be described using density functional theory (DFT). In this work, we used a hybrid quantum mechanical, molecular mechanical, and machine learning potential for studying the adsorption behavior of phenol, atomic hydrogen, 2-butanol, and 2-butanone on the (0001) facet of Ru under reducing conditions when Ru is not oxidized. Specifically, we describe the adsorbate and the surrounding metal atoms at the DFT level of theory. Here, we also considered the electrostatic field effect of the water molecules on adsorbate–metal interactions. Next, for the water–water and water–adsorbate interactions, we used established classical force fields. Finally, for the water–Ru surface interaction, for which no reliable force fields have been published, we used Behler–Parrinello high-dimensional neural network potentials (HDNNPs). Employing this setup, we used our explicit solvation for metal surface (eSMS) approach to compute the aqueous-phase effect on the low-coverage adsorption of selected molecules and atoms on the (0001) facet of Ru. In agreement with previous experimental and computational studies of oxygenated molecules over transition metal facets, we found that liquid water destabilizes the tested adsorbates on Ru(0001). Interestingly, our findings indicate that adsorbates on Ru are less affected by the presence of an aqueous phase than on other transition metals (e.g., Pt), highlighting the necessity of experimental investigations of Ru-based catalytic systems in liquid water.

Adsorption↗
DOE OSTI2024

Understanding Strain and Failure of a Knot in Polyethylene Using Molecular Dynamics with Machine-Learned Potentials

A neural network potential (NNP) has been developed by fitting to ab initio electronic structure data on hydrocarbons and is used to study failure of linear and knotted polyethylene (PE) chains. A linear PE chain must be highly strained before breaking as the stress is equally distributed across the chain. In contrast, the stress in a PE chain with a 31 or overhand knot, accumulates at the knot’s entrance/exit. We find the strain energy is greatest when the bond length and angle are strained simultaneously, and that the knot weakens the chain by increasing the variance of the C–C–C angle, thereby allowing rupture at lower bond strains. Here, we extend our analysis to both 51 and 52 knots and find that both break at the entrance/exit of a loop. Notably, molecular scale PE knots exhibit many of the same characteristics as knots in a macroscopic rope, with stick–slip phenomena upon tightening and similar points of failure.

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

The design space of E(3)-equivariant atom-centred interatomic potentials

Abstract Molecular dynamics simulation is an important tool in computational materials science and chemistry, and in the past decade it has been revolutionized by machine learning. This rapid progress in machine learning interatomic potentials has produced a number of new architectures in just the past few years. Particularly notable among these are the atomic cluster expansion, which unified many of the earlier ideas around atom-density-based descriptors, and Neural Equivariant Interatomic Potentials (NequIP), a message-passing neural network with equivariant features that exhibited state-of-the-art accuracy at the time. Here we construct a mathematical framework that unifies these models: atomic cluster expansion is extended and recast as one layer of a multi-layer architecture, while the linearized version of NequIP is understood as a particular sparsification of a much larger polynomial model. Our framework also provides a practical tool for systematically probing different choices in this unified design space. An ablation study of NequIP, via a set of experiments looking at in- and out-of-domain accuracy and smooth extrapolation very far from the training data, sheds some light on which design choices are critical to achieving high accuracy. A much-simplified version of NequIP, which we call BOTnet (for body-ordered tensor network), has an interpretable architecture and maintains its accuracy on benchmark datasets.

Computer Science↗
DOE OSTI2025

Computational investigation of water glasses using machine-learning potentials

The molecular origins of water’s anomalous properties have long been a subject of scientific inquiry. The liquid–liquid phase transition hypothesis, which posits the existence of distinct low-density and high-density liquid states separated by a first-order phase transition terminating at a critical point, has gained increasing experimental and computational support and offers a thermodynamically consistent framework for many of water’s anomalies. However, experimental challenges in avoiding crystallization near the postulated liquid–liquid critical point have focused attention to water’s canonical glassy states: low-density and high-density amorphous ice. Here, we use two Deep Potential machine-learning models, trained on the Strongly Constrained and Appropriately Normed density functional and the highly accurate Many-Body Polarizable potential, to conduct an investigation of water’s glassy phenomenology based on quantum mechanical calculations. Despite not being explicitly trained on amorphous ices, both models accurately capture the structure and transformation of the water glasses, including their interconversion along different thermodynamic paths. Isobaric quenching of liquid water at various pressures generates a continuum of intermediate amorphous ices and density fluctuations increase near the liquid–liquid critical pressure. The glass transition temperatures of the amorphous ices produced at different pressures exhibit two distinct branches, corresponding to low-density and high-density amorphous ice behaviors, consistent with experiment and the liquid–liquid transition hypothesis. Extrapolating transformation pressures from isothermal compressions to experimental compression rates brings our simulations into excellent agreement with data. Our findings demonstrate that machine-learning potentials trained on equilibrium phases can effectively model nonequilibrium glassy behavior and pave the way for studying long-timescale, out-of-equilibrium processes with quantum mechanical accuracy.

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

Breaking the Million-Electron and 1 EFLOP/s Barriers: Biomolecular-Scale Ab Initio Molecular Dynamics Using MP2 Potentials

The accurate simulation of complex biochemical phenomena has historically been hampered by the computational requirements of high-fidelity molecular-modeling techniques. Quantum mechanical methods, such as ab initio wave-function (WF) theory, deliver the desired accuracy, but have impractical scaling for modeling biosystems with thousands of atoms. Combining molecular fragmentation with MP2 perturbation theory, this study presents an innovative approach that enables biomolecular-scale ab initio molecular dynamics (AIMD) simulations at WF theory level. Leveraging the resolution-of-the-identity approximation for Hartree-Fock and MP2 gradients, our approach eliminates computationally intensive four-center integrals and their gradients, while achieving near-peak performance on modern GPU architectures. The introduction of asynchronous time steps minimizes time step latency, overlapping computational phases and effectively mitigating load imbalances. Utilizing up to 9,400 nodes of Frontier and achieving 59% (1006.7 PFLOP/s) of its double-precision floating-point peak, our method enables us to break the million-electron and 1EFLOP/s barriers for AIMD simulations with quantum accuracy.

Kurzak, Jakub↗
DOE OSTI2024

Shallow Rate-Redox Potential Scaling in Aqueous Molecular Oxygen Reduction Electrocatalysis Across a Family of Iron Macrocycles

Rate-overpotential scaling relationships have been employed widely to understand trends in oxygen reduction reaction (ORR) electrocatalysis by dissolved metal macrocycles in organic electrolytes. Similar scaling relationships remain unknown for surface-adsorbed ORR electrocatalysts in the acidic aqueous environments germane to proton-exchange membrane (PEM) fuel cells. Herein, we examine ORR catalysis in aqueous perchloric acid media for a structurally diverse array of iron macrocycle complexes adsorbed on Vulcan carbon black. The macrocycles encompass Fe– N 4 , Fe–N 2 N' 2 and Fe–N x C 4-x motifs bearing pyrrolic, pyridinic, and N-heterocyclic carbene (NHC) moieties in the primary ligation sphere, giving rise to a 670 mV range in Fe(III/II) redox potentials, E Fe(III/II) . Experimental Tafel data in the micropolarization regime were extrapolated to the E Fe(III/II) to furnish estimated per-site-normalized current density (j per-site ) values that span ~4.6 orders of magnitude across the family of compounds. Despite the structural diversity of this family of compounds, extrapolated j per-site values correlate with the Fe(III/II) redox potentials in a roughly log-linear fashion with a shallow scaling factor of approximately 145 mV/decade. Further, these findings highlight that negative shifts in E Fe(III/II) lead to diminishing returns in catalytic rate promotion and suggest that changes to the primary ligating environment in a macrocycle are insufficient to break fundamental rate-potential scaling relationships in aqueous ORR catalysis. Together these studies motivate the further development of higher-potential iron complexes that employ motifs beyond the equatorial ligation plane to enhance ORR catalysis.

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

The role of an intramolecular hydrogen bond in the redox properties of carboxylic acid naphthoquinones

A bioinspired naphthoquinone model of the quinones in photosynthetic reaction centers but bearing an intramolecular hydrogen-bonded carboxylic acid has been synthesized and characterized electrochemically, spectroscopically, and computationally to provide mechanistic insight into the role of proton-coupled electron transfer (PCET) of quinone reduction in photosynthesis. The reduction potential of this construct is 370 mV more positive than the unsubstituted naphthoquinone. In addition to the reversible cyclic voltammetry, infrared spectroelectrochemistry confirms that the naphthoquinone/naphthoquinone radical anion couple is fully reversible. Calculated redox potentials agree with the experimental trends arising from the intramolecular hydrogen bond. Molecular electrostatic potentials illustrate the reversible proton transfer driving forces, and analysis of the computed vibrational spectra supports the possibility of a combination of electron transfer and PCET processes. The significance of PCET, reversibility, and redox potential management relevant to the design of artificial photosynthetic assemblies involving PCET processes is discussed.

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

X‑ray Coherent Diffractive Imaging of Large Helium Nanodroplets Doped with Small Molecules

We report the first X-ray coherent diffractive imaging experiment on molecule-doped helium nanodroplets. It complements previous work, where we reported single-shot X-ray coherent diffractive imaging studies of Xe dopant clusters formed in 4He and 3He droplets. These noble gas clusters were used to visualize the impact of rotational excitation of the droplets on the spatial distribution of atomic dopants within the droplets, and to study the differences and connections between quantum and classical droplet rotational motion. Here, we expand our studies to the molecular dopants CF4, CHF3, CH3CN, and SF6, imaged with 1.5 keV photons. We find multiple Bragg spots in the diffraction patterns of molecule-doped droplets with radii of approximately 600 nm, which provide evidence that molecules form elongated clusters with preferential alignment along the angular momentum axis of the 4He droplets, in agreement with our previous results on the aggregation of Xe clusters on quantum vortices. Real-space reconstructions of molecular dopant cluster density profiles are obtained for droplets with smaller radii of approximately 300 nm. The diffuse images suggest the formation of low-density, potentially porous, molecular clusters upon aggregation at T = 0.4 K in 4He droplets. In the normal fluid 3He droplets, molecules aggregate into loose clusters on the droplets' equator, similar to previous observations for Xe atoms. Time-of-flight mass spectra reveal that the doped helium nanodroplet moieties fragment extensively into constituent atomic ions, producing only a small fraction of molecular fragment ions. The findings are discussed in the context of previously proposed schemes to use He droplets as potential tamper materials for ultrafast X-ray imaging experiments.

Feinberg, AlexandraJ↗
DOE OSTI2026

Foundational insights into the mechanical and molecular evolution of porcine skin gelatin during gelation

Gelatin is a widely used material in biomedical fields, particularly in regenerative medicine and tissue engineering, due to its biocompatibility and versatile properties. While prior research has explored methods to enhance gelatin's mechanical strength and stability, fundamental studies on gelatin, specifically its curing process, mechanical stiffness, and chemical evolution during gelation, remain limited. This study uses ultrasonic testing and Fourier Transform Infrared Spectroscopy (FTIR) to examine gelatin's stiffness and molecular changes during gelation. Samples of 175 and 300 Porcine Skin Bloom Strength Gelatin at concentrations of 2% and 6% (w/v) were analyzed. Through transmission ultrasonic testing helped identify key transition points in gelation, with higher concentrations exhibiting delayed transitions. FTIR revealed that C-N bond formation peaks early while N-H bond deformation persists. A correlation emerged between sound speed and peak absorbance, suggesting that changes in molecular mobility may contribute to the observed sound speed behavior during periods of active bond formation. However, as gelation continues, fewer bonding components may be available, potentially decreasing molecular movement and contributing to the observed increase in sound speed. These findings provide insights into gelatin's mechanical and chemical evolution, offering a framework for improved control over its gelation kinetics. Swept-Frequency Acoustic Interferometry (SFAI) was performed at the end of the curing process to measure the sound speed, enabling the calculation of the bulk moduli of the gelatin samples. The combined use of ultrasonic and FTIR testing provides a non-destructive method for characterizing gelatin and other biomaterials. This approach advances understanding of gelatin curing behavior and supports the development of safer biomaterials with tailored mechanical properties for various applications such as tissue engineering and regenerative medicine.

Biomaterials↗
DOE OSTI2024

Secondary electron emission measurements from imidazolium-based ionic liquids

The electron-induced secondary electron emission (SEE) yields of imidazolium-based ionic liquids are presented for primary electron beam energies between 30 and 1000 eV. These results are important for understanding plasma synthesis of nanoparticles in plasma discharges with an ionic liquid electrode. Due to their low vapor pressure and high conductivity, ionic liquids can produce metal nanoparticles in low-pressure plasmas through reduction of dissolved metal salts. In this work, the low vapor pressure of ionic liquids is exploited to directly measure SEE yields by bombarding the liquid with electrons and measuring the resulting currents. The ionic liquids studied are [BMIM][Ac], [EMIM][Ac], and [BMIM][BF 4 ]. The SEE yields vary significantly over the energy range, with maximum yields of around 2 at 200 eV for [BMIM][Ac] and [EMIM][Ac], and 1.8 at 250 eV for [BMIM][BF 4 ]. Molecular orbital calculations indicate that the acetate anion is the likely electron donor for [BMIM][Ac] and [EMIM][Ac], while in [BMIM][BF 4 ], the electrons likely originate from the [BMIM] + cation. The differences in SEE yields are attributed to varying ionization potentials and molecular structures of the ionic liquids. These findings are essential for accurate modeling of plasma discharges and understanding SEE mechanisms in ionic liquids.

ionic liquids↗