SEARCH · Search NASA

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Active deep kernel learning of molecular properties from structural embeddings

As vast databases of chemical identities become increasingly available, the challenge shifts to how we effectively explore and leverage these resources to study molecular properties. This paper presents an active learning approach for molecular discovery using deep kernel learning (DKL), demonstrated on the QM9 dataset. DKL links structural embeddings directly to properties, creating organized latent spaces that prioritize relevant property information. By iteratively recalculating embedding vectors in alignment with target properties, DKL uncovers concentrated maxima representing key molecular properties and reveals unexplored regions with potential for innovation. This approach underscores DKL’s potential in advancing molecular research and discovery.

Artificial neural networks↗

DOE OSTI2025

Resolving the Solvation Structure and Transport Properties of Aqueous Zinc Electrolytes from Salt-in-Water to Water-in-Salt Using Neural Network Potential

Zn Cl 2 solutions are promising electrolytes for aqueous zinc-ion batteries. Here, we report a joint computational and experimental study of the structural and dynamic properties of aqueous Zn Cl 2 electrolytes with concentrations ranging from salt-in-water to water-in-salt (WIS). By developing a neural network potential (NNP) model, we perform molecular dynamics (MD) simulations with accuracy but at much larger lengths and longer timescales. The NNP predicted structures are validated by the structure factors measured by X-ray total scattering experiments. The MD trajectories provide a comprehensive and quantitative picture of the Zn 2 + solvation shell structures. Additionally, we find that the O − H covalent bonds in water are strengthened with increasing salt concentration, thus expanding the electrochemical stability window of aqueous electrolytes. In terms of dynamic properties, the calculated and experimentally measured conductivities are in good agreement. Through the analysis of the calculated cation transference number, we propose a three-stage charge carrier transport mechanism with increasing concentration: independent ion transport, strongly correlated ion transport, and small positive charge carrier diffusion through negatively charged polymeric clusters. Our study provides fundamental atomic scale insights into the structure and transport properties of the Zn Cl 2 electrolyte that can aid the optimization and development of WIS electrolytes. Published by the American Physical Society 2025

25 ENERGY STORAGE↗

DOE OSTI2024

Bottom-up design of actinide materials from molecular clusters: Demonstration of a general-purpose simulation capability leveraging machine-learned atomic potentials

Actinide thin-film coatings such as uranium dioxide (UO 2 ) play an important role in nuclear reactors and other mission-relevant applications, but realization of their potential requires a deep fundamental understanding of the chemical vapor deposition (CVD) processes used for their growth. The slow experimental progress can be attributed, in part, to the standard safety guidelines associated with handling uranium byproducts, which are often corrosive, toxic, and radioactive. Accurate simulation techniques, when used in concert with experiment, can improve laboratory safety, material durability, and deliverable timeframes. However, state-of-the-art computational methods are either insufficiently accurate or intractably expensive. To remedy this situation, in this project we suggested a machine-learning (ML) accelerated workflow for simulating molecular clustering toward deposition. As a benchmark test case, we considered molecular clustering in steam and assessed independent components of our workflow by comparing with measured thermodynamic properties of water. After analyzing each component individually and finding no fundamental barrier to realization of the workflow, we attempted to integrate the ML component, a Sandia-developed tool called FitSNAP. As this was the first application of FitSNAP to atoms and molecules in the gas phase at Sandia, the method required more fitting data than was originally anticipated. Systematic improvements were made by including in the fit data diatomic potentials, molecular single-bond-breaking curves, and symmetry-constrained intermolecular potentials. We concluded that our strategy provides a feasible pathway toward modeling CVD and related processes, but that extensive training data must be generated before it can be of practical use.

36 MATERIALS SCIENCE↗

DOE OSTI2025

High temperature melting of dense molecular hydrogen from machine-learning interatomic potentials trained on quantum Monte Carlo

We present results and discuss methods for computing the melting temperature of dense molecular hydrogen using a machine learned model trained on quantum Monte Carlo data. In this newly trained model, we emphasize the importance of accurate total energies in the training. We integrate a two phase method for estimating the melting temperature with estimates from the Clausius–Clapeyron relation to provide a more accurate melting curve from the model. We make detailed predictions of the melting temperature, solid and liquid volumes, latent heat, and internal energy from 50 to 180 GPa for both classical hydrogen and quantum hydrogen. At pressures of roughly 173 GPa and 1635 K, we observe molecular dissociation in the liquid phase. Here, we compare with previous simulations and experimental measurements.

08 HYDROGEN↗

DOE OSTI2026

Effects of Composition and Oxidation States on the Structures of Chromium-Containing Sodium Silicate Glasses: Molecular Dynamics Simulations using Machine Learning Interatomic Potentials

Chromium represents a significant challenge for the vitrification of high-level nuclear waste into silicate and borosilicate glasses due to its low solubility and variable oxidation states, which can limit the waste loading due to promotion of crystallization or phase separation during processing. In this study, we modeled chromium containing silicate glasses using molecular dynamics simulations with three machine learning interatomic potentials (MLIPs), MACE, CHGNet, and PFP were employed, to gain insights on glass composition and oxidation states on the structures of these glasses. One of the goals is to evaluate their ability of these MLIPs to accurately represent the general structure of silicate glasses and chromium local environments as a function of chromium oxidation states. Density Functional Theory (DFT) based calculations and experimental data such as neutron structure factors were used to validate the structural models. It was found that the foundation models of all three MLIPs are able to reproduce general structural features of the sodium silicate glass structure consistent with experimental and DFT data, but only CHGNet and PFP can accurately capture the oxidation states and local environment of chromium: tetrahedral for Cr6+ and octahedral for Cr3+. Furthermore, we studied the effect of varying Cr3+/ Cr6+ (Cr3+/Crtotal) ratio and total chromium content using PFP. Our results show that Cr6+ enhances network polymerization by reducing non-bridging oxygens through Na? charge compensation required due to the formation of chromate (CrO42-) species, while Cr³? acts as a network modifier that disrupts connectivity. System size effects on the structural characteristics and chromium environments were also tested using the PFP potential. This work highlights the importance of careful validation on the precision, transferability, and potential of MLIPs for modeling glasses containing transition metal elements that can exist in multiple oxidation states. It is also encouraging to see the foundational models are all three MLFFs are able to reproduce the basic sodium silicate glass structures, while suggesting additional training or refining is needed to improve the description of more complex systems containing transition metals.

Puga, Christina L.↗

DOE OSTI2025

Modeling Equilibrium Solid–Liquid Interfaces under Effective Constant Chemical Potential Using Machine Learning Interatomic Potentials

The chemical potential (μ) of species in solution is essential for understanding various chemical processes at interfaces. Molecular dynamics (MD) simulations, constrained by fixed compositions, cannot maintain constant chemical potential with reference to a targeted concentration or chemical potential under nonequilibrium or dynamic conditions, as solute species can migrate to the interface and deplete (or enrich) the bulk due to solute-interface interactions. In this study, we introduce a simple and computationally efficient approach named iterative quasi-constant chemical potential molecular dynamics (iqCμMD) simulation, which helps simulate targeted molar concentrations of species in solution. iqCμMD overcomes the limitations of conventional MD by adjusting the number of species in the solution to reach a target bulk concentration (chemical potential), which allows simulation of the interface under the bulk conditions comparable to experiment. We demonstrate our approach using machine learning interatomic potential (MLIP)-based MD simulations of the Na 2 SO 4,aq –graphene interface, and to show the transferability of our approach, we also perform classical force field-based MD simulations of NaCl aq –air and NaCl aq –graphite interfaces, which produce comparable results to previous CμMD simulations. Our results also show that the iqCμMD approach efficiently achieves the desired bulk ion concentration within two iterations, and by utilizing MLIPs, we can achieve converged results using relatively small-scale simulations compared to previous CμMD simulations. By combining iqCμMD with MLIP-driven simulations, solid–liquid interfaces can be modeled under an effective constant chemical potential with DFT-level accuracy. Here, we show that iqCμMD offers a robust and simple computational framework for constant chemical potential simulations, as its only requirement is to be able to converge interfacial simulations with a measurable bulk region.

Chemical structure↗

DOE OSTI2025

Nuclear quantum effects of metal surface-mediated C–H activation

The nuclear quantum effects of surface-mediated C–H activation of surface CH 3 are considered for the pristine Pt(111) and Au(111) surfaces at 300 K. The kinetic barriers without nuclear quantum effects are calculated using both static density functional theory calculations and ab initio molecular dynamics. Static calculations are performed using the harmonic approximation while the free energy pathway is calculated using enhanced sampling molecular dynamics. Machine learning potentials are trained using generated datasets and validated against the ab initio molecular dynamics generated free energy pathways. The machine learning potentials are used to perform centroid molecular dynamics to consider the nuclear quantum effects of C–H activation. Nuclear quantum effects are found to have a very significant effect on the free energy pathway, with reduced importance at higher temperatures and in the CD 3 case.

Bunting, Rhys J. [Lawrence Livermore National Labo↗

DOE OSTI2026

Genetic and epigenetic signatures of Populus trichocarpa in response to abiotic stresses

Plants respond to abiotic stresses such as drought, heat, and salinity through both shared and stress-specific regulatory pathways. However, the role of epigenetic mechanisms, particularly DNA methylation, in modulating these responses is still underexplored. Here, we integrated transcriptome profiling with whole-genome bisulfite sequencing (WGBS) to investigate the dynamic relationship between gene expression and DNA methylation in Populus trichocarpa during brief early-stage exposure to abiotic stress. Each stress elicited distinct transcriptional and methylation signatures; however, a subset of conserved stress-responsive genes was commonly regulated across treatments. We found that ∼5-6% of differentially expressed genes also displayed differential methylation patterns, suggesting a coordinated role for DNA methylation in regulating gene expression. Motif enrichment analysis of differentially methylated regions revealed binding sites of key transcription factor families, including ERF, bHLH, and ABF, highlighting potential role for methylation in modulating transcription factor targeting. Furthermore, we identified stress-inducible molecular markers with potential applications in early stress detection and functional dissection of gene regulatory pathways. Together, these findings provide new insights into the coordinated genetic and epigenetic responses to abiotic stress and provide a foundation for developing biosensors and breeding strategies to enhance stress resilience woody plants.

DNA methylation↗

DOE OSTI2025

Genetic_and_epigenetic_signatures_of_Populus_Trichocarpa_in_response_to_abiotic_stress

Plants respond to abiotic stresses such as drought, heat, and salinity through both shared and stress-specific regulatory pathways. However, the role of epigenetic mechanisms, particularly DNA methylation, in modulating these responses is still underexplored. Here, we integrated transcriptome profiling with whole-genome bisulfite sequencing (WGBS) to investigate the dynamic relationship between gene expression and DNA methylation in Populus trichocarpa during brief early-stage exposure to abiotic stress. Each stress elicited distinct transcriptional and methylation signatures; however, a subset of conserved stress-responsive genes was commonly regulated across treatments. Approximately 5% of differentially expressed genes also displayed differential methylation patterns, suggesting a coordinated role for DNA methylation in regulating gene expression. Motif enrichment analysis of differentially methylated regions revealed binding sites of key transcription factor families, including ERF, bHLH, and ABF, highlighting potential role for methylation in modulating transcription factor targeting. Furthermore, we identified stress-inducible molecular markers with potential applications in early stress detection and functional dissection of gene regulatory pathways. Together, these findings provide new insights into the coordinated genetic and epigenetic responses to abiotic stress and provide a foundation for developing biosensors and breeding strategies to enhance stress resilience woody plants.

CBI↗

DOE OSTI2024

Converting Second‐Order Saddle Points to Transition States: New Principles for the Design of 4π Photoswitches

Abstract Molecular photoswitches have demonstrated potential for storing solar energy at the molecular level, with power densities comparable to commercial batteries and hydroelectric energy storage. However, development of efficient photoswitches is hindered by limitations in cyclability and optical properties of existing materials. We here demonstrate that certain limitations in photoswitches based on electrocyclizations stem from the issue of controlling competition between Woodward‐Hoffmann allowed and forbidden pathways. Our approach moves beyond the traditional view of activation barriers and reveals that second‐order saddle points are crucial in dictating the competition between disrotatory and conrotatory pathways. These insights suggest new opportunities to manipulate the competition between these pathways through geometric constraints, fundamentally altering the connectivity of the potential energy surface. Our study also emphasizes the necessity of multi‐reference methods and the need to conduct higher‐dimensional explorations for competing pathways beyond photoswitch design.

Chemistry↗

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 OSTI2026

Comparison of DeePMD, MTP, GAP, ACE and MACE Machine‐Learned Potentials for Radiation‐Damage Simulations: A User Perspective

Accurate and efficient interatomic potentials are essential for molecular dynamics (MD) simulations of radiation damage, gas diffusion, and phase stability in complex ceramics such as LiAlO 2 , especially under extreme conditions relevant to tritium production. Here, we evaluate the performance of six machine-learned interatomic potentials (MLIPs), moment tensor potential (MTP), Gaussian approximation potential, deep potential (DeePMD), atomic cluster expansion (ACE), message-passing ACE (multilayer atomic cluster expansion (MACE) pretrained) and MACE (trained from-scratch), all trained on the same density functional theory dataset with inclusion of tritium. The MLIPs are benchmarked against traditional Buckingham and ReaxFF potentials in terms of energy accuracy, density predictions, thermal equilibration behavior, threshold displacement energy (E d ), tritium diffusivity, and computational cost. Among the models, MTP shows the best overall balance between efficiency and accuracy, with low force and energy errors and realistic E d values for Li and Al. The ACE and MACE (pretrained and trained from scratch) models exhibit high E d (>200 eV) and unphysical pair interactions. DeePMD underestimates Ed due to overly repulsive behavior even at equilibrium distances. All models over-estimate tritium diffusion but the pretrained MACE model behaves well during tritium-diffusion simulations up to 500 K, maintaining diffusivities in the physically consistent 10 −11 m 2 /s range. Finally, we quantify the computational cost of each potential in large-scale atomic/molecular massively parallel simulator, finding that only MTP is more efficient than traditional empirical potentials, while others are significantly more expensive. These findings explain the trade-offs between accuracy and computational cost in MLIP development and provide essential guidance for use in high-throughput radiation damage and gas diffusion simulations in nuclear ceramics.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2026

Effect of Solvent on the Local Structure, Dynamics, and Vibrational Density of States in Sn-BEA Zeolite

Lewis acid zeolites are attractive catalysts for epoxidation and biomass valorization, as they are highly active and selective in the liquid phase and can operate at or near ambient conditions. While a rich experimental literature exists on liquid-phase Lewis acid zeolite catalysis, our understanding of the molecular organization and solvent dynamics in the vicinity of Lewis acid sites with differing metal site speciation remains limited. In this work, we investigate the molecular coordination and diffusion of two common solvents (methanol and water) around the closed and open Sn-BEA zeolite active sites using molecular dynamics simulations with a machine-learned interatomic potential trained on ab initio molecular dynamics trajectories. Molecular dynamics simulations reveal that introducing active sites significantly enhances local order in the first and second solvation shells compared to the pure silica case. For methanol, both closed and open active sites are singly coordinated, while more than two water molecules coordinate the open site. In contrast to methanol, we observed that water molecules dissociate, leading to the formation of additional Sn-OH and silanol groups away from the active site. The diffusion coefficients of water and methanol are functions of the solvent population in the pore. Here, our work provides insights into how active site speciation in Lewis acid zeolites affects solvent coordination, diffusion, and vibrational signature. This information is foundational for catalyst design and optimization of liquid-phase catalytic processes in zeolites. It also demonstrates the suitability of machine-learned interatomic potentials for modeling reactive systems, enabling sufficiently long trajectories for appropriate statistical averaging.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Development of a deep potential model for F and CF 2 etching of Si and SiO 2

An understanding of plasma-surface interactions at increasingly smaller scales is invaluable for the development of novel technologies and processing techniques. Molecular dynamics (MD) simulations can provide insights into atomic-scale interactions, though they are restricted by the availability of interatomic potentials. Machine learning methods, such as Deep Potential Molecular Dynamics (DeepMD), provide a systematic framework for the development of accurate and flexible ab initio-based models. In this work, we develop DeepMD models for the ion-enhanced etching of Si and SiO 2 by F and CF 2 radicals. We employ an active learning process to expand the data set on which the model is trained and demonstrate its effect on the model accuracy. The DeepMD results are compared to data from classical MD simulations and experiments. Physical sputtering yields of SiO 2 by Ar + ions show good agreement with previous simulation results using conventional interatomic potentials, though the predicted depth profiles are different. Etching yields are calculated as a function of ion energy and neutral to ion flux ratio for the Ar + ion-enhanced etching of SiO 2 and Si by F atoms, as well as for etching of SiO 2 by CF 2 radicals, showing reasonable agreement with experimental data. Finally, an ion-enhanced surface kinetic model is fitted to the DeepMD etch yields, and the fitted parameters are compared to quantities computed directly from DeepMD simulations. This study illustrates how molecular dynamics simulations using machine learning potentials can provide an accurate model of etching processes relevant to device manufacturing.

Kounis-Melas, Andreas [Princeton Univ., NJ (United↗

DOE OSTI2025

A cavity-enhanced spectroscopist's lens on molecular polaritons

Polariton chemistry has been hailed as a potential new route to direct molecular processes with electromagnetic fields. To make further strides, the community needs to clarify which unusual polaritonic phenomena are true hallmarks of cavity quantum electrodynamics and which can be rationalized with classical optical physics. Here, in this work, we provide a tutorial perspective on the formation, spectroscopy, and behavior of molecular polaritons using classical optics. Where possible, we draw connections to cavity-enhanced spectroscopy and recast open questions in terms that may be more familiar to the broader community of physical chemists.

Cavity quantum electrodynamics↗

DOE OSTI2025

Highly Conjugated Porphyrin Arrays Enable Optical Resolution of Ferromagnetic and Antiferromagnetic Aligned States of the Triplet Exciton and an Incorporated Stable Radical

Well-defined photogenerated molecular spin systems have potential utility in spintronics and quantum information science (QIS). Because molecular magnetic, optical, and electronic properties can be controlled by design, diverse spin systems can be prepared at modest temperatures. Photogenerated molecular spin systems often involve states prepared from the interaction of excitons and charges. Resolving the nature of electron spin alignment in photogenerated spin states described by the coupling of a triplet exciton and a stable radical commonly relies on EPR spectroscopy. Here, we describe ethyne-bridged (porphinato)metal ( PM n ) oligomers that incorporate a macrocycle-bound Cu(II) radical center. Upon photoexcitation of such PM n arrays, a singdoublet ( 2 S 1 ) state is formed; ultrafast internal conversion (IC) then produces a tripdoublet ( 2 T 1 ) state, which undergoes intersystem crossing (ISC) to produce a tripquartet ( 4 T 1 ) state, before relaxation to the ground state ( 2 S 0 ). These highly conjugated Cu(II) radical-containing PMn arrays enable direct observation of copper porphyrin 2 T 1 → 4 T 1 ISC dynamics from the biexponential decay of the near-infrared (NIR) 2,4 T 1 → 2,4 T n transient absorption manifold. Multireference n-electron valence perturbation theory (NEVPT2) computations illuminate how PM n electronic structure controls the relaxation dynamics of these long-lived (>10 ns) electronically excited multiplet states. These studies show that highly conjugated and polarizable porphyrin arrays incorporating stable spin centers provide rare π-delocalized systems where the ferromagnetic and antiferromagnetic alignment between a triplet exciton and a stable radical are both spectrally resolved and addressable using transient optical spectroscopy at wavelengths exceeding 1 μm, providing new opportunities to QIS.

14 SOLAR ENERGY↗

DOE OSTI2025

Temperature-Dependent Molecular Diffusional Properties in Deep Eutectic Solvents and Eutectogels

Eutectogels (ETGs) prepared from deep eutectic solvents (DESs), a gelator, and water have many uses in separations, catalysis, and energy storage systems. In these applications, temperature-dependent molecular diffusional properties and intermolecular interactions play a critical role in their function. Diffusional properties of Alexa Fluor 633 and ATTO 647N were measured across a range of temperatures in choline chloride:2glycerol DESs comprised of one molar equivalent of choline chloride and two molar equivalents of glycerol (also known as glyceline) as well as ETGs made from this DES, a xanthan gum gelator, and 10% w/w water or 20% w/w water. Fluorescence recovery after photobleaching (FRAP) was employed to evaluate potential changes in molecular diffusion within DESs and ETGs from 20 to 100 °C. Surprisingly, the ETGs have a larger sample viscosity but also exhibit faster molecular diffusion when compared to the dry DES. This is attributed to macroscopic properties of the ETGs (e.g., pores, three-dimensional structure). The FRAP data also show an irreversible temperature-dependent decrease in the diffusion coefficient of both fluorophores in the ETGs. This is consistent with differential scanning calorimetry data of the ETG, which shows a shift to a more negative glass transition temperature (−36 °C to −47 °C) after the first cooling/heating cycle. Raman data reveal no detectable changes in the intermolecular interactions in either the DESs or ETGs as a function of temperature or water content. Finally, the findings on temperature- and water-dependent diffusional properties of DESs and ETGs provide a foundation for optimizing the use of these materials across various applications, particularly where repeated heating/cooling cycles may be used.

alcohols↗

DOE OSTI2025

Parameterizing empirical interatomic potentials for predicting thermophysical properties via an irreducible derivative approach: the case of ThO 2 and UO 2

The accuracy of classical physical property predictions using molecular dynamics simulations is determined by the quality of the interatomic potentials. Here we introduce a training approach for empirical interatomic potentials (EIPs) which is well suited for capturing phonons and phonon-related properties. Our approach is based on direct comparisons of the second- and third-order irreducible derivatives (IDs) between an EIP and the Born–Oppenheimer potential within density functional theory (DFT) calculations. IDs fully exploit space group symmetry and allow for training without redundant information. We demonstrate the fidelity of our approach in the context of ThO 2 and UO 2 , where we optimize parameters of an embedded-atom method potential in addition to core–shell interactions. Our EIPs provide thermophysical properties in good agreement with DFT and outperform widely utilized EIPs for phonon dispersion and thermal conductivity predictions. Reasonable estimates of thermal expansion and formation energies of Frenkel pairs are also obtained.

empirical interatomic potential↗