Search NASA⌕ Search

SEARCH · Search NASA

Results for “Interaction energy”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Machine Learning-Based Upscaling of Finite-Size Molecular Dynamics Diffusion Simulations for Binary Fluids

Molecular diffusion coefficients calculated using molecular dynamics (MD) simulations suffer from finite-size (i.e., finite box size and finite particle number) effects. Results from finite-sized MD simulations can be upscaled to infinite simulation size by applying a correction factor. For self-diffusion of single-component fluids, this correction has been well-studied by many researchers including Yeh and Hummer (YH); for binary fluid mixtures, a modified YH correction was recently proposed for correcting MD-predicted Maxwell–Stephan (MS) diffusion rates. In this study we use both empirical and machine learning methods to identify improvements to the finite-size correction factors for both self-diffusion and MS diffusion of binary Lennard-Jones (LJ) fluid mixtures. Using artificial neural networks (ANNs), the error in the corrected LJ fluid diffusion is reduced by an order of magnitude versus existing YH corrections, and the ANN models perform well for mixtures with large dissimilarities in size and interaction energies where the YH correction proves insufficient.

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

Condensed Phase Water Molecular Multipole Moments from Deep Neural Network Models Trained on Ab Initio Simulation Data

Ionic solvation phenomena in liquids involve intense interactions in the inner solvation shell. For interactions beyond the first shell, the ion–solvent interaction energies result from the sum of many smaller-magnitude contributions that can still include polarization effects. Deep neural network (DNN) methods have recently found wide application in developing efficient molecular models that maintain near-quantum accuracy. In this work, we extend the DeePMD-kit code to produce accurate molecular multipole moments in the bulk and near interfaces. The new method is validated by comparing the DNN moments with those generated by ab initio simulations. The moments are used to compute the electrostatic potential at the center of a molecular-sized hydrophobic cavity in water. The results show that the fields produced by the DNN models are in quantitative agreement with the AIMD-derived values. These efficient methods will open the door to more accurate solvation models for large solutes such as proteins.

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

Nanoconfinement of High Hydrogen-to-Metal Ratio Lanthanum Hydrides in Functionalized Carbon Hosts

Metal hydrides with a high hydrogen content are important for materials-based hydrogen storage and high-temperature superconductivity. Nanoconfinement of metal hydrides in porous hosts is a promising strategy to tune the thermodynamic stability and control the hydrogen-to-metal ratio. However, lanthanum hydride (LaH x ) nanoconfinement in porous materials has been limited due to the challenges associated with isolating and stabilizing nanoparticles of La or La-hydrides. Here we successfully demonstrated the chemical reduction of La­(III) salts to La(0), and subsequent infiltration into pure CMK-3 and nitrogen-doped CMK-3 (NCMK-3) porous carbons. Transmission electron microscopy measurements revealed a uniform distribution of LaH x species within the carbon hosts, while X-ray absorption and photoelectron spectroscopy provided detailed information about the local chemical environment. Sieverts measurements indicate that LaH x @NCMK-3 could desorb up to 0.75 wt % hydrogen, which is higher than non-nitrogen-functionalized CMK-3 (0.43 wt % H). Density Functional Theory and ab initio molecular dynamics calculations indicate that host–guest interaction energies are favorable for porous carbon with nitrogen defects, which is supported by experimental evidence. Moreover, high-pressure synchrotron X-ray diffraction measurements were conducted using a diamond anvil cell up to 60 GPa and reveal that the nitrogen-functionalized nanoporous carbon host favors the formation of higher H:La ratios in the presence of ammonia borane compared to the pure CMK-3 host. This approach could serve as a suitable platform for developing nanoscale superconducting materials at lower pressures and temperatures compared to bulk.

Shivanna, Mohana↗
DOE OSTI2024

The Open DAC 2023 Dataset and Challenges for Sorbent Discovery in Direct Air Capture

Direct air capture (DAC) of CO 2 with porous adsorbents such as metal-organic frameworks (MOFs) has the potential to aid large-scale decarbonization. Previous screening of MOFs for DAC relied on empirical force fields and ignored adsorbed H 2 O and MOF deformation. We performed quantum chemistry calculations overcoming these restrictions for thousands of MOFs. The resulting data enable efficient descriptions using machine learning.

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

A Transferable Force Field for Predicting Adsorption and Diffusion of Water in Cationic Zeolites with Coupled Cluster Accuracy

We present a transferable force field for water in proton-exchanged, alkali (Li, Na, K, Rb, and Cs) metal-exchanged, and alkaline-earth (Mg, Ca, Sr, and Ba) metal-exchanged zeolites. The fitting methodology is based on adsorbate–adsorbent interaction energies obtained from periodic density functional theory calculations and corrected using the coupled-cluster method applied to small model clusters. To ensure an accurate prediction of both adsorption and diffusion properties of water, sets of configurations that sample both adsorption sites and intracrystalline hopping transition states were used in the fitting. The quality of the force field is assessed for a wide range of zeolites with different topologies and chemical compositions, demonstrating good agreement between theoretical predictions and experimental measurements of water adsorption and diffusion.

Adsorption↗
DOE OSTI2021

Surface-Enhanced Raman Scattering Optophysiology Nanofibers for the Detection of Heavy Metals in Single Breast Cancer Cells

Mercury(II) ions (Hg2+) and silver ions (Ag+) are two of the most hazardous pollutants causing serious damage to human health. Here, we constructed surface enhanced Raman scattering (SERS)-active nanofibers covered with 4-mercaptopyridine (4-Mpy)-modified gold nanoparticles (AuNPs) to detect Hg2+ and Ag+. Experimental evidence suggests that the observed spectral changes originate from the combined effect of (i) the coordination between the nitrogen on 4-Mpy and the metal ions and (ii) the 4-Mpy molecular orientation (from flatter to more perpendicular with respect to the metal surface). The relative intensity of a pair of characteristic Raman peaks (at ~428 and ~708 cm-1) was used to quantify the metal ion concentration, greatly increasing the reproducibility of the measurement compared to signal-on or signal-off detection based on a single SERS peak. The detection limit of this method for Hg2+ is lower than that for the Ag+ (5 nM vs 100 nM), which can be explained by the stronger interaction energy between Hg2+ and N compared to Ag+ and N, as demonstrated by density functional theory (DFT) calculations. The Hg2+ and Ag+ ions can be masked by adding ethylenediaminetetraacetate (EDTA) and Cl-, respectively, to the Hg2+ and Ag+ samples. The good sensitivity, high reproducibility and excellent selectivity of these nanosensors were also demonstrated. Furthermore, detection of Hg2+ in living breast cancer cells at the subcellular level is possible thanks to the nanometric size of the herein described SERS nanosensors, allowing high spatial resolution and minimal cell damage.

Zhao, Xingjuan↗
DOE OSTI2025

Deep Eutectic Solvents for Efficient Fractionation of Lignocellulose to Produce Uncondensed Lignin and High-Quality Cellulose

Simultaneously inhibiting lignin condensation and cellulose degradation remains a major challenge for achieving holistic valorization of lignocellulose. Here, we developed a deep eutectic solvent (DES), composed of l-cysteine (Cys) and lactic acid (LA), to fractionate both uncondensed lignin and high-quality cellulose from eucalyptus wood by leveraging the unique properties of Cys, i.e., highly nucleophilic groups (–SH) and hydrogen bond acceptor/donor groups (–NH 2 and –COOH). The nucleophilic –SH in Cys effectively quenches the benzylic carbocations (C α + ions, formed at the benzylic sites of lignin) that lead to lignin condensation. This enables the high yield of uncondensed lignin (81%) with high retention of β–O–4 bonds (up to 90%). The separated uncondensed lignin is further depolymerized to prepare monophenols in a satisfactory 43% yield, equivalent to 73% of the theoretical yield. Moreover, the –NH 2 and –COOH groups in Cys form extensive hydrogen bonds with the hydroxyl groups in cellulose, thus decreasing the interaction energy of DES on cellulose. As a result, the cellulose achieves an astonishing 99% retention and maintains a high degree of polymerization of 1160. The obtained high-quality cellulose is further conversed into cellulose nanofibers for strong and transparent films. Furthermore, this study provides new insights into the efficient separation of uncondensed lignin and high-quality cellulose from lignocellulose by a novel DES system.

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

Synthesis and Computational Analysis of Uranium(III)-Pnictogen Bonds

In this report, we describe the synthesis and characterization of two novel complexes that contain uranium(III)-pnictogen (P: 1-PMes 2 ; As: 1-AsMes 2 ) bonds to elucidate the degrees of covalency among these bonds. To our knowledge, this is the first reported uranium(III)-arsenic complex to be synthesized and characterized. Analysis of the phosphorus and arsenic bonds reveals a similar electronic environment, assessed by UV–vis NIR, that is comparable to other previously reported uranium(III) complexes. A computational analysis of these compounds and their congeners, N, Sb, and Bi, was performed to identify trends in the overall bonding character of the complexes. This analysis shows that bond covalency decreases as the pnictogen becomes heavier and that overall the interaction energy and its components decrease down the group. Here, this study provides an in-depth analysis and understanding of the nature of bonding between hard actinide and soft pnictogen centers.

Actinides↗
DOE OSTI2023

Many-Body Potential for Simulating Self-Assembly of Polymer-Grafted Nanoparticles in a Polymer Matrix

Many-body interactions between polymer-grafted nanoparticles (NPs) play a key role in promoting their assembly into low-dimensional structures within polymer melts, even when the particles are spherical and isotropically grafted. However, capturing such interactions in simulations of NP assembly is very challenging because explicit modeling of the polymer grafts and melt chains is highly computationally expensive, even using coarse-grained models. Here, we develop a many-body potential for describing the effective interactions between spherical polymer-grafted NPs in a polymer matrix through a machine-learning approach. The approach involves using permutationally invariant polynomials to fit two- and three-body interactions derived from the potential of mean force calculations. The potential developed here reduces the computational cost by several orders of magnitude, thereby, allowing us to explore assembly behavior over large length and time scales. We show that the potential not only reproduces previously known assembled phases such as 1D strings and 2D hexagonal sheets, which generally cannot be achieved using isotropic two-body potentials, but can also help discover interesting phases such as networks, clusters, and gels. We demonstrate how each of these assembly morphologies intrinsically arises from a competition between two- and three-body interactions. Our approach for deriving many-body effective potentials can be readily extended to other colloidal systems, enabling researchers to make accurate predictions of their behavior and dissect the role of individual interaction energy terms of the overall potential in the observed behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Magnetic-field-induced Wigner crystallization of charged interlayer excitons in van der Waals heterostructures

Abstract Ever since its inception, coherent excited states of semiconductors have been the focus of semiconductor materials research to evolve into a vibrant field of low-dimensional solid-state physics. The field is gaining new momentum these days due to emerging transdimensional semiconductors such as van der Waals bound layers of transition metal dichalcogenides (TMDs) of controlled thickness. Here, we develop the theory of magnetic-field-induced Wigner crystallization for charged interlayer excitons (CIE) discovered recently in TMD heterobilayers. We derive the ratio of the potential interaction energy to the kinetic energy for the many-particle CIE system in the perpendicular magnetostatic field of an arbitrary strength and predict the crystallization effect in the strong field regime. We show that magnetic-field-induced Wigner crystallization and melting of CIEs can be observed in magneto-photoluminescence experiments with TMD bilayers of systematically varied electron-hole doping concentrations. Our results advance the capabilities of this new generation of transdimensional quantum materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Ultrafast enhancement of electron-phonon coupling via dynamic quantum well states

Abstract The density of states at the Fermi surface controls many material properties by influencing the low-energy interactions of conductive electrons. Typically, external tuning knobs generate only small perturbations to the density of states in crystals, since their band structure depends strongly on the lattice potential. In contrast, quantum well states are contingent on a localized potential extrinsic to the crystal lattice and can be easily modified leading to changes in the density of states at the Fermi level. Here, we are able to control the quantum well potential on the surface of Bi 2 Se 3 with light, driving a density of states singularity below E F at ultrafast timescales, thereby triggering a reversible Lifshitz transition. We reveal a substantial ultrafast enhancement of the electron-phonon coupling with repercussions on the relaxation dynamics in the system. We argue that the relaxation dynamics are governed largely by the interplay of the dynamic density of states of the quantum wells with c -axis scattering from the $${A}_{{{{{{{{\rm{1g}}}}}}}}}^{2}$$ A 1g 2 optical phonon mode. These results demonstrate a powerful way to enhance electron-boson interaction on ultrafast timescales providing new avenues for controlling material properties and driving novel quantum phases of matter.

Ciocys, Samuel T.↗
DOE OSTI2022

Impact of secondary salts, temperature, and pH on the colloidal stability of graphene oxide in water

The stability of graphene oxide (GO) in water is extremely relevant because of its application as an adsorbent material, as well as for its fate and behavior in the environment. Zeta potential was used to study the effect of secondary salts (carbonate, sulfate, and phosphate), temperature (20 to 60 °C), and pH (5 to 9) on the stability of six different GOs produced from natural, synthetic, and amorphous graphite—with and without the use of attrition milling. Generally, GOs produced with attrition-milled graphites had lower ζ-potentials than their unmilled counterparts because of their smaller particle sizes and higher concentration of oxygen-containing functional groups. It was observed that GO produced from graphite and synthetic graphite had ζ-potential values lower than –30 mV, even at 30 °C. However, it was observed that all the GOs studied were unstable in the presence of carbonate and sulfate salts at concentrations between 170 and 1695 mg L –1 , as they reached a ζ-potential of –4.1 mV. Density-functional theory electronic structure calculations suggested that the instability of GO in the presence of carbonate and sulfate was caused by the abstraction of a proton resulting in interaction energies E int of 28.3 and 168.9 kJ mol –1 , respectively. Our results suggest that temperatures above 30 °C, as well as carbonate and sulfate salts at concentrations relevant to arid and semi-arid regions, could promote the formation of agglomerates of GO, thus limiting its use and mobility in water.

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

Feasibility study of Mg storage in a bilayer silicene anode via application of an external electric field

With the goal of developing a Si-based anode for Mg-ion batteries (MIBs) that is both efficient and compatible with the current semiconductor industry, the current research utilized classical Molecular Dynamics (MD) simulation in investigating the intercalation of a Mg 2+ ion under an external electric field (E-field) in a 2D bilayer silicene anode (BSA). First principles density functional theory calculations were used to validate the implemented EDIP potentials. Our simulation shows that there exists an optimum E-field value in the range of 0.2–0.4 V Å –1 for Mg 2+ intercalation in BSA. To study the effect of the E-field on Mg 2+ ions, an exhaustive spread of investigations was carried out under different boundary conditions, including calculations of mean square displacement (MSD), interaction energy, radial distribution function (RDF), and trajectory of ions. Our results show that the Mg 2+ ions form a stable bond with Si in BSA. The effects of E-field direction and operating temperature were also investigated. In the X–Y plane in the 0°–45° range, 15° from the X-direction was found to be the optimum direction for intercalation. The results of this work also suggest that BSA does not undergo drastic structural changes during the charging cycles with the highest operating temperature being ~300 K

25 ENERGY STORAGE↗
DOE OSTI2023

The Li + CaF → Ca + LiF chemical reaction under cold conditions

The calcium monofluoride (CaF) molecule has emerged as a promising candidate for precision measurements, quantum simulation, and ultracold chemistry experiments. Inelastic and reactive collisions of laser cooled CaF molecules in optical tweezers have recently been reported and collisions of cold Li atoms with CaF are of current experimental interest. In this paper, we report ab initio electronic structure and full-dimensional quantum dynamical calculations of the Li + CaF → LiF + Ca chemical reaction. The electronic structure calculations are performed using the internally contracted multi-reference configuration-interaction method with Davidson correction (MRCI + Q). An analytic fit of the interaction energies is obtained using a many-body expansion method. Furthermore, a coupled-channel quantum reactive scattering approach implemented in hyperspherical coordinates is adopted for the scattering calculations under cold conditions. Results show that the Li + CaF reaction populates several low-lying vibrational levels and many rotational levels of the product LiF molecule and that the reaction is inefficient in the 1–100 mK regime allowing sympathetic cooling of CaF by collisions with cold Li atoms.

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

Advancing density functional tight-binding method for large organic molecules through equivariant neural networks

Semi-empirical quantum-mechanical (QM) methods have become valuable tools for studying complex (bio)molecular systems due to their balance between computational efficiency and accuracy. A key aspect of these methods is their parameterization, which not only governs the reliability of the results but also provides an opportunity to enhance their overall performance. In our previous work [J. Phys. Chem. Lett., 2021, 11, 16], we advanced the third-order semi-empirical density functional tight-binding (DFTB3) method for computing multiple properties of small molecules by developing the machine learning (ML) potential NN rep to bridge the gap between DFTB3 electronic components and those of the hybrid DFT-PBE0 functional. To overcome the limitations of NN rep , we introduce the EquiDTB framework, which leverages physics-inspired equivariant neural networks (NN) to parameterize scalable and transferable many-body Δ TB potentials, replacing the standard pairwise DFTB repulsive potential. This advancement extends the applicability of our ML-corrected DFTB approach to larger molecules and non-covalent systems (including only C, N, O, and H atoms), going beyond the chemical space represented in the training QM datasets. The enhanced performance of EquiDTB over the standard TB methods is demonstrated by the accurate computation of the atomic forces of S66x8 molecular dimers, as well as their interaction energies. Moreover, EquiDTB can be effectively employed to explore the potential energy surfaces of large and flexible drug-like molecules—for example, to determine the minimum energy path between isomers, analyze structural transitions during dynamical simulations, compute vibrational modes, and investigate energetic rankings. The performance for single molecules slightly decreases when the DFTB electronic energy is reduced to first-order but remains superior to standard TB methods. Our work thus demonstrates that an optimal integration of an equivariant NN with QM datasets can advance the DFTB method while maintaining high efficiency, paving the way for reliable (bio)molecular simulations.

Medrano Sandonas, Leonardo [Technische Universität↗
DOE OSTI2024

Toward first principles-based simulations of dense hydrogen

Accurate knowledge of the properties of hydrogen at high compression is crucial for astrophysics (e.g., planetary and stellar interiors, brown dwarfs, atmosphere of compact stars) and laboratory experiments, including inertial confinement fusion. There exists experimental data for the equation of state, conductivity, and Thomson scattering spectra. However, the analysis of the measurements at extreme pressures and temperatures typically involves additional model assumptions, which makes it difficult to assess the accuracy of the experimental data rigorously. On the other hand, theory and modeling have produced extensive collections of data. They originate from a very large variety of models and simulations including path integral Monte Carlo (PIMC) simulations, density functional theory (DFT), chemical models, machine-learned models, and combinations thereof. At the same time, each of these methods has fundamental limitations (fermion sign problem in PIMC, approximate exchange–correlation functionals of DFT, inconsistent interaction energy contributions in chemical models, etc.), so for some parameter ranges accurate predictions are difficult. Recently, a number of breakthroughs in first principles PIMC as well as in DFT simulations were achieved which are discussed in this review. Here we use these results to benchmark different simulation methods. We present an update of the hydrogen phase diagram at high pressures, the expected phase transitions, and thermodynamic properties including the equation of state and momentum distribution. Furthermore, we discuss available dynamic results for warm dense hydrogen, including the conductivity, dynamic structure factor, plasmon dispersion, imaginary-time structure, and density response functions. We conclude by outlining strategies to combine different simulations to achieve accurate theoretical predictions that are based on first principles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Molecular dynamics based study on the effects of cation size on the local structure and diffusion in polymerized ionic liquids

We have used coarse-grained molecular dynamics simulations to understand the effect of cation size on ion diffusion in polymerized ionic liquids at temperatures well above the glass transition temperature of the polymers. We investigated dependencies of the diffusion constant on the cation radius and static dielectric constant and interpreted these results in terms of underlying structural changes, and decoupling phenomenon. We have found non-monotonic effects of the cation radius on the diffusion constant, with a maximum at an intermediate radius resulting from two different size dependent effects. Changes in the radial distribution functions characterizing spatial distribution of cations with respect to anions result in a monotonic increase in the coordination number. This yields an increase in the effective interaction energy as a function of cation radius, which can explain the decrease in the diffusion constant. This increase is counteracted at small cation radii by a corresponding increase in decoupling of cation and anion dynamics. In addition, the diffusion constant is found to increase on increasing the static dielectric constant for all radii of the cations. These results highlight the significant role of decoupling in designing polymerized ionic liquids with an enhanced diffusion constant of cations.

Gillespie, Colin [ORNL] (ORCID:0000000320287080)↗
DOE OSTI2026

Analysis of Second Target Station Target Removal Dose Rates

The Second Target Station (STS) project at Oak Ridge National Laboratory’s spallation neutron source is a crucial initiative for maintaining U.S. leadership in neutron sciences. The STS aims to create the world’s brightest pulsed cold neutron source, enabling cutting-edge research across various scientific disciplines. To ensure safe and efficient maintenance operations, understanding the effects of shutdown dose rates from activated components within the STS target systems is essential. This study establishes a computational framework for calculating decay gamma sources and subsequent shutdown dose rates utilizing advanced methods to account for all activation channels, including high-energy interactions down to thermal neutron capture. This study describes a novel integration of multiple tools and provides an effective means of analyzing activation and shutdown dose rates at spallation neutron facilities. A custom-developed script automates the decay gamma source generation process, ensuring proper sampling during the variance reduction phase, which is critical for accurate predictions of shutdown dose rates.

Transmutation↗