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 55 records · Page 3

Energetics of Base–Acid Pairs for the Design of High-Temperature Fuel Cell Polymer Electrolytes

The interaction energy of base–acid plays a key role in acid retention of phosphoric acid (PA)-doped polymer electrolytes under fuel cell operating conditions. Here, we investigate the energetics of proton-accepting and hydroxide-donating organic bases using density functional theory calculations. Because of their weak basicity, proton-accepting organic bases such as benzimidazole have relatively low interaction energy with the acid in the absence of water (15.3–28.0 kcal mol –1 ). Energetics of the proton-accepting base–PA complex increases by adding water, indicating that the interactions in the base–acid complex strengthen in the presence of water. On the other hand, hydroxide-donating organic bases, such as tetramethylammonium hydroxide, have high interaction energy with PA (~110 kcal mol –1 ), which remains high in the presence of water. The chemical shifts of 31P NMR support the energetics of the base–acid complexes. Overall, this study further discusses the benefit of incorporating hydroxide-donating organic bases into the polymeric structure over proton-accepting bases as a way to increase acid retention.

25 ENERGY STORAGE↗

Molecular interactions from the density functional theory for chemical reactivity: Interaction chemical potential, hardness, and reactivity principles

In the first paper of this series, the authors derived an expression for the interaction energy between two reagents in terms of the chemical reactivity indicators that can be derived from density functional perturbation theory. While negative interaction energies can explain reactivity, reactivity is often more simply explained using the “|dμ| big is good” rule or the maximum hardness principle. Expressions for the change in chemical potential (μ) and hardness when two reagents interact are derived. A partial justification for the maximum hardness principle is that the terms that appear in the interaction energy expression often reappear in the expression for the interaction hardness, but with opposite sign.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theoretical tools for neutrino scattering: interplay between lattice QCD, EFTs, nuclear physics, phenomenology, and neutrino event generators

Maximizing the discovery potential of increasingly precise neutrino experiments will require an improved theoretical understanding of neutrino-nucleus cross sections over a wide range of energies. Low-energy interactions are needed to reconstruct the energies of astrophysical neutrinos from supernovae bursts and search for new physics using increasingly precise measurement of coherent elastic neutrino scattering. Higher-energy interactions involve a variety of reaction mechanisms including quasi-elastic scattering, resonance production, and deep inelastic scattering that must all be included to reliably predict cross sections for energies relevant to DUNE and other accelerator neutrino experiments. This white paper discusses the theoretical status, challenges, required resources, and path forward for achieving precise predictions of neutrino-nucleus scattering and emphasizes the need for a coordinated theoretical effort involved lattice QCD, nuclear effective theories, phenomenological models of the transition region, and event generators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Systematic Evaluation of Counterpoise Correction in Density Functional Theory

A widespread belief persists that the Boys–Bernardi function counterpoise (CP) procedure “overcorrects” supramolecular interaction energies for the effects of basis-set superposition error. To the extent that this is true for correlated wave function methods, it is usually an artifact of low-quality basis sets. The question has not been considered systematically in the context of density functional theory, however, where basis-set convergence is generally less problematic. We present a systematic assessment of the CP procedure for a representative set of functionals and basis sets, considering both benchmark data sets of small dimers and larger supramolecular complexes. The latter include layered composite polymers with ~150 atoms and ligand–protein models with ~300 atoms. Provided that CP correction is used, we find that intermolecular interaction energies of nearly complete-basis quality can be obtained using only double-ζ basis sets. Furthermore, this is less expensive as compared to triple-ζ basis sets without CP correction. CP-corrected interaction energies are less sensitive to the presence of diffuse basis functions as compared to uncorrected energies, which is important because diffuse functions are expensive and often numerically problematic for large systems. Our results upend the conventional wisdom that CP “overcorrects” for basis-set incompleteness. In small basis sets, CP correction is mandatory in order to demonstrate that the results do not rest on error cancellation.

74 ATOMIC AND MOLECULAR PHYSICS↗

Evaluating the active site-substrate interplay between x-ray crystal structure and molecular dynamics in chorismate mutase

Designing realistic quantum mechanical (QM) models of enzymes is dependent on reliably discerning and modeling residues, solvents, and cofactors important in crafting the active site microenvironment. Interatomic van der Waals contacts have previously demonstrated usefulness toward designing QM-models, but their measured values (and subsequent residue importance rankings) are expected to be influenceable by subtle changes in protein structure. Using chorismate mutase as a case study, this work examines the differences in ligand-residue interatomic contacts between an x-ray crystal structure and structures from a molecular dynamics simulation. Select structures are further analyzed using symmetry adapted perturbation theory to compute ab initio ligand-residue interaction energies. The findings of this study show that ligand-residue interatomic contacts measured for an x-ray crystal structure are not predictive of active site contacts from a sampling of molecular dynamics frames. In addition, the variability in interatomic contacts among structures is not correlated with variability in interaction energies. However, the results spotlight using interaction energies to characterize and rank residue importance in future computational enzymology workflows.

Chemistry↗

Studying light–matter interactions and energy transfer at the nanoscale with a trapped–ion quantum computer

A grand challenge in materials sciences is to control energy transduction, transfer and dissipation at the nanoscale. Meeting this challenge requires progress in nanoscale fabrication and synthesis, but equally importantly, also progress in our understanding of how nanoscale structure influences these phenomena. Although detailed quantum mechanical models exist that capture the structure and physics of materials and molecular species, solving these models to calculate the relevant emergent phenomena, which often occur over vastly different time and energy scales, has so far proven difficult. In this project, we studied how to model energy transfer and light-matter interactions using trapped-ion quantum computers. Our studies focused on emulation of critical processes in natural systems, such as the well-known photosynthetic pigment-protein complexes responsible for capturing and transducing light to electrons. In addition, we implemented proof-of-principle experiments demonstrating that such calculations are possible on trapped-ion quantum computing hardware. These studies and experiments pave the way to investigate energy transport and transduction in a detail not accessible with classical computing technologies. We expect that our studies will allow to extract important guiding principles on how to improve the efficiency of nano-scale light-harvesting and transduction devices.

36 MATERIALS SCIENCE↗

Density Functionals for Hydrogen Storage: Defining the H2Bind275 Test Set with Ab Initio Benchmarks and Assessment of 55 Functionals

Efficient and high-capacity storage materials are indispensable for a hydrogen-based economy. In silico tools can accelerate the process of discovery of new adsorbent materials with optimal hydrogen adsorption enthalpies. Density functional theory is well-poised to become a very useful tool for enabling high-throughput screening of potential materials. In this work, we have identified density functional approximations that provide good performance for hydrogen binding applications following a two-pronged approach. First, we have compiled a data set (H2Bind275) that comprehensively represents the hydrogen binding problem capturing the chemical and mechanistic diversity in the binding sites encountered in hydrogen storage materials. We have also computed reference interaction energies for this data set using coupled-cluster theory. Second, we have assessed the performance of 55 density functional approximations for predicting H 2 interaction energies and have identified two hybrid density functionals (ωB97X-V and ωB97M-V), two double hybrid density functionals (DSD-PBEPBE-D3(BJ) and PBE0-DH), and one semilocal density functional (B97M-V) as the best performing ones. We have recommended the addition of empirical dispersion corrections to systematically underbinding density functionals such as revPBE, BLYP, and B3LYP for improvements in performance at negligible additional cost. We have also recommended the usage of the def2-TZVPP basis set as it represents a good compromise between accuracy and cost, limiting the finite basis set errors to less than 1 kJ/mol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Toward an Understanding of Linear Scaling Relations through Energy Decomposition Analysis

The discovery of linear scaling relations has fundamentally changed the field of heterogeneous catalysis. The scaling relations have been rationalized based on the d-band theory, specifically a separation of sp and d electron contributions to adsorption energies. Within the framework of energy decomposition analysis, a full understanding of such a separation would require one to further break down the adsorption energy into distinct energy components such as electrostatics, polarization, charge transfer, and van der Waals interactions, and to examine the sp and d contributions to each of them. As a step in this direction, we analyzed the interaction energy between CH x (x = 1–4) adsorbates and fcc(100) transition metal surfaces (M = Cu, Ag, Au, Rh, and Pt), with the surfaces represented both as slabs in plane-wave density functional theory (pw-DFT) calculations and as atomic clusters in atomic-orbital basis density functional theory (ao-DFT) calculations. Through an absolutely localized molecular orbital (ALMO) based energy decomposition analysis of the ao-DFT adsorption energy, each of the interaction energy components (electrostatics, polarization, van der Waals, and charge transfer) was found to follow its own scaling relations, with an intricate interplay among these energy components yielding the overall scaling relations for the total adsorption energies. Using the recently introduced ALMO-based polarization and charge-transfer analysis schemes, we further dissected polarization into metal surface and adsorbate contributions, and charge transfer into metal → adsorbate and adsorbate → metal contributions. The contributions from the sp and d electrons of the metal to these terms were further quantified, and the dominant role of the metal d electrons was reaffirmed. These results shed light on how CHx adsorbates interact with metal surfaces and further reveal the physical origin of the scaling relations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spectroscopic Probe Molecule Selection Using Quantum Theory, First-Principles Calculations, and Machine Learning

Probe molecule vibrational spectra have a long history of being used to characterize materials including metals, oxides, metal-organic frameworks, and even human proteins. Furthermore, recent advances in machine learning have enabled computationally generated spectra to aid in detailed characterization of complex surfaces with probe molecules. Despite widespread use of probe molecules, the science of probe molecule selection is underdeveloped. Here, we develop physical concepts, including orbital interaction energy and the energy overlap integral, to explain and predict the ability of probe molecules to discriminate structural descriptors. We resolve the crystal orbital overlap population (COOP) to specific molecular orbitals and quantify their bonding character, which directly influences vibrational frequencies. Using only a single adsorbate calculation from density function theory (DFT), we compute the interaction energy of individual adsorbate molecular orbitals with adsorption site atomic orbitals across many different sites. Combining the molecular orbital resolved COOP and changes in orbital interaction energy enables probe molecule selection for improved discrimination of various sites. We demonstrate these concepts by comparing the predicted effectiveness of carbon monoxide (CO), nitric oxide (NO), and ethylene (C 2 H 4 ) to probe Pt adsorption sites. Finally, using a previously developed machine learning framework, we show that models trained on hundreds of thousand C 2 H 4 spectra, computed from DFT, which regress surface binding-type and generalized coordination number (GCN), outperform those trained using CO and NO spectra. Lastly, a python package, pDOS_overlap, for implementing the electron density based analysis on any combination of adsorbates and materials, is also made available.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Equilibrium Thermodynamic Properties of Aqueous Solutions of Ionic Liquid 1-Ethyl-3-Methylimidazolium Methanesulfonate [EMIM][MeSO3]

Abstract The ionic liquid 1-ethyl-3-methylimidazolium methanesulfonate ([EMIM][MeSO 3 ]) has been considered as a promising alternative desiccant to triethylene glycol and lithium bromide commonly used in the industry. In this paper, the water activity coefficient of this binary system was measured from 303 K to 363 K with water concentration from 18% to 92%. The interaction energies between the ionic liquid molecules ( $${g}_{22}$$ g 22 ) and between the ionic liquid and water molecules ( $${g}_{12}$$ g 12 ) for the [EMIM][MeSO 3 ]/water binary system were determined from the water activity coefficient data using the Non-Random Two-Liquid (NRTL) model. The magnitude of the interaction energy between the [EMIM][MeSO 3 ] and water molecules ( $${g}_{12}$$ g 12 ) was found to be in the range of 45~49 kJ/mol, which was about 20% larger than that between the water molecules ( $${g}_{11}$$ g 11 ) in the [EMIM][MeSO 3 ]/water system. The large ( $${g}_{12}$$ g 12 ) can explain many observed macroscopic thermodynamic properties such as strong hygroscopicity in the ionic liquid [EMIM][MeSO 3 ]. These interaction energies were used to determine the heat of desorption of the [EMIM][MeSO 3 ]/water system, and the obtained heat of desorption was in good agreement with that calculated from the conventional Clausius-Clapeyron Equation.

42 ENGINEERING↗

Phonon gas model for thermal conductivity of dense, strongly interacting liquids

Developing predictive thermal property models for liquids based on microscopic principles has been elusive. The difficulty is that liquids have gas-like and solid-like attributes that are at odds when considering the frameworks of microscopic models: Models for gases are simple due to randomness and low density, whereas models for crystalline solids rely on symmetry and long-range order for easier calculation. The short-range order in liquids does, however, provide structure to neighboring molecules similar to amorphous solids, and there have been recent advances indicating that collective vibrational modes store heat in liquids. Models combining Debye approximations from solid-state physics and Frenkel’s theory of liquids can accurately predict the heat capacity of liquids. Phonon-like dispersions in liquids have also been widely observed in neutron scattering experiments. These developments motivate us to propose a model where high-frequency vibrational modes, which travel at the speed of sound and have a mean free path on the order of the average intermolecular distance, conduct heat in liquids. We use this liquid phonon gas model to calculate the thermal conductivity of liquids with varying intermolecular interaction energies from strongest to weakest—Coulomb, hydrogen-bonding, Keesom, and London dispersion energy. Generally, the model is more accurate as the intermolecular interaction energy and density of liquids increase. The calculated thermal conductivity of Coulombic-bound molten sodium nitrate and hydrogen-bonded water is within 1.46% and 2.98% of the experimentally measured values, respectively, across their entire temperature ranges. Further modal analysis of the velocity and the mean free path of collective vibrations could establish the liquid phonon gas model as an accurate model for weakly interacting liquids as well.

36 MATERIALS SCIENCE↗

Predicting and Understanding Non-Covalent Interactions Using Novel Forms of Symmetry-Adapted Perturbation Theory

Although sometimes derided as “weak” interactions, non-covalent forces play a critical role in ligand binding and crystal packing and in determining the conformational landscape of flexible molecules. Symmetry-adapted perturbation theory (SAPT) provides a framework for accurate ab initio calculation of intermolecular interactions and furnishes a natural decomposition of the interaction energy into physically meaningful components: semiclassical electrostatics (rigorously obtained from monomer charge densities), Pauli or steric repulsion, induction (including both polarization and charge transfer), and dispersion. This decomposition helps to foster deeper understanding of non-covalent interactions and can be used to construct transferable, physics-based force fields. Separability of the SAPT interaction energy also provides the flexibility to construct composite methods, a feature that we exploit to improve the description of dispersion interactions. Furthermore, these are challenging to describe accurately because they arise from nonlocal electron correlation effects that appear for the first time at second order in perturbation theory but are not quantitatively described at that level.

74 ATOMIC AND MOLECULAR PHYSICS↗

Beyond Price Taker: Optimizing Integrated Energy Systems Considering Market/Grid Interactions

Integrated Energy Systems (IES) combine two or more processes to increase the efficiency, flexibility of operation, and the overall reliability. However, analyzing IESs in volatile electricity markets is challenging, since the volatility in electricity prices makes the conventional levelized cost-type analysis less realistic. This work presents two approaches to address the challenge: price-taker and a surrogates-based approach for incorporating market interactions. The price-taker approach formulates a multiperiod optimization problem that takes the time-varying electricity prices into account, and solves the optimization problem to determine the optimal operational schedule that maximizes the chosen economic metric. This approach is successfully applied to investigate the performance of flexible power and hydrogen co-production systems. The market surrogates approach trains a machine learning model to predict the market behavior as a function of the characteristics of the IES. The trained surrogate model is used to optimize the design and operation of the given IES in an electricity market. This approach is demonstrated on a case study involving a nuclear power plant retrofitted with a low-temperature electrolysis unit to co-produce power and hydrogen.

beyond price taker↗

Nuclear Data for High Energy Ion Interactions and Secondary Particle Production

Understanding the harmful effects of galactic cosmic rays (GCRs) on space exploration requires a substantial amount of nuclear data. Specifically, the interaction of energetic GCR charged particles with spacecraft materials generates secondary radiations that, through energy deposition, can harm astronauts and electronic systems. By identifying the gaps in our knowledge of the relevant nuclear data – interaction cross sections – and identifying ways to fill those gaps – with measurements, compilations, evaluations, disseminations, reaction modeling, sensitivity studies, and uncertainty quantification – the safety and viability of space exploration can be improved. This work surveys the state of the art in this interdisciplinary field and identifies promising collaborative research topics that have significant potential to advance our understanding of the effects of the space radiation environment on space exploration.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Low-energy electronic interactions in ferrimagnetic Sr 2 CrReO 6 thin films

We reveal in this study the fundamental low-energy landscape in the ferrimagnetic Sr 2 CrReO 6 double perovskite and describe the underlying mechanisms responsible for the three low-energy excitations below 1.4 eV. Here, based on resonant inelastic x-ray scattering and magnetic dynamics calculations, and experiments collected from both Sr 2 CrReO 6 powders and epitaxially strained thin films, we reveal a strong competition between spin-orbit coupling, Hund’s coupling, and the strain-induced tetragonal crystal field. Here, we also demonstrate that a spin-flip process is at the origin of the lowest excitation at 200 meV, and we bring insights into the predicted presence of orbital ordering in this material. We study the nature of the magnons through a combination of ab initio and spin-wave theory calculations, and show that two nondegenerate magnon bands exist and are dominated either by rhenium or chromium spins. The rhenium band is found to be flat at about 200 meV (±25 meV) through X-L-W-U high-symmetry points and is dispersive toward Γ.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Density Functional Theory Solute-Dislocation Interactions

DFT-optimized structures of dislocations with a solute at different sites. The dataset include the VASP input files (except POTCAR) and output files OSZICAR and OUTCAR. There are also text files containing the DFT-computed interaction energies. The README file contains expressions that can be used to compute the interaction energies and chemical misfit energies using the volumetric size misfit for any solute species in Mg.

36 MATERIALS SCIENCE↗

Density Functional Theory Solute-Dislocation Interactions

DFT-optimized structures of dislocations with a solute at different sites. The dataset include the VASP input files (except POTCAR) and output files OSZICAR and OUTCAR. There are also text files containing the DFT-computed interaction energies. The README file contains expressions that can be used to compute the interaction energies and chemical misfit energies using the volumetric size misfit for any solute species in Mg.

Trinkle, Dallas↗

Density Functional Theory Solute-Dislocation Interactions

DFT-optimized structures of dislocations with a solute at different sites. The dataset include the VASP input files (except POTCAR) and output files OSZICAR and OUTCAR. There are also text files containing the DFT-computed interaction energies. The README file contains expressions that can be used to compute the interaction energies and chemical misfit energies using the volumetric size misfit for any solute species in Mg.

Trinkle, Dallas↗