Search NASA⌕ Search

SEARCH · Search NASA

Results for “interaction potential”

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 19 records

AP-Net: An atomic-pairwise neural network for smooth and transferable interaction potentials

Intermolecular interactions are critical to many chemical phenomena, but their accurate computation using ab initio methods is often limited by computational cost. The recent emergence of machine learning (ML) potentials may be a promising alternative. Useful ML models should not only estimate accurate interaction energies but also predict smooth and asymptotically correct potential energy surfaces. However, existing ML models are not guaranteed to obey these constraints. Indeed, systemic deficiencies are apparent in the predictions of our previous hydrogen-bond model as well as the popular ANI-1X model, which we attribute to the use of an atomic energy partition. As a solution, we propose an alternative atomic-pairwise framework specifically for intermolecular ML potentials, and we introduce AP-Net—a neural network model for interaction energies. The AP-Net model is developed using this physically motivated atomic-pairwise paradigm and also exploits the interpretability of symmetry adapted perturbation theory (SAPT). We show that in contrast to other models, AP-Net produces smooth, physically meaningful intermolecular potentials exhibiting correct asymptotic behavior. Initially trained on only a limited number of mostly hydrogen-bonded dimers, AP-Net makes accurate predictions across the chemically diverse S66x8 dataset, demonstrating significant transferability. On a test set including experimental hydrogen-bonded dimers, AP-Net predicts total interaction energies with a mean absolute error of 0.37 kcal mol−1, reducing errors by a factor of 2–5 across SAPT components from previous neural network potentials. The pairwise interaction energies of the model are physically interpretable, and an investigation of predicted electrostatic energies suggests that the model “learns” the physics of hydrogen-bonded interactions.

Glick, Zachary L. (ORCID:0000000309002849)↗

New interaction potentials for alkaline earth silicate and borate glasses

New interaction potentials were developed for molecular dynamics simulations to study the role of Mg and Ca in modifying the structure and properties of alkaline earth silicates and borates. Competition between the depolymerization of the silica network and the formation of new bonds between oxygen atoms and modifiers leads to the enhancement of the elastic moduli with increasing modifier content in alkaline earth silicate glasses. Compared with calcium silicate, the higher elastic moduli of magnesium silicate result from a higher connectivity of the overall glass network due to the incorporation of fourfold coordinated magnesium and a more rigid connection between oxygen atoms and modifiers. In contrast to the silicates, the effect of modifier on the elastic moduli of alkaline earth borates is dominated by the formation of fourfold coordinated boron (N 4 ). Calcium borate with higher N 4 shows a more rigid network structure and higher elastic moduli.

36 MATERIALS SCIENCE↗

New interaction potentials for borate glasses with mixed network formers

Here, we adapt and apply a recently developed optimization scheme used to obtain effective potentials for aluminosilicate glasses to include the network former boron into the interaction parameter set. As input data for the optimization, we used the radial distribution functions of the liquid at high temperature generated by ab initio molecular dynamics simulations, and density, coordination, and elastic modulus of glass at room temperature from experiments. The new interaction potentials are shown to reliably reproduce the structure, coordination, and mechanical properties over a wide range of compositions for binary alkali borates. Furthermore, the transferability of these new interaction parameters allows mixing to reliably reproduce the properties of various boroaluminate and borosilicate glasses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Stability of trapped solutions of a nonlinear Schrödinger equation with a nonlocal nonlinear self-interaction potential

This work focuses on the study of the stability of trapped soliton-like solutions of a (1 + 1)-dimensional nonlinear Schrödinger equation (NLSE) in a nonlocal, nonlinear, self-interaction potential of the form [|Ψ(x,t)| 2 +|Ψ(-x,t)| 2 ] κ where κ is an arbitrary nonlinearity parameter. Although the system with κ = 1 (i.e. fully integrable case) was first reported by Yang (2018 Phys. Rev. E 98 042202), here in the present work, we extend this model to the one in which κ is arbitrary. This allows us to compare the stability properties of the now trapped solutions to previously found solutions of the more usual NLSE with κ ≠ 1 which are moving soliton solutions. We show that there is a simple, one-component, nonlocal Lagrangian and corresponding action governing the dynamics of the system. Using a collective coordinate method derived from the action as well as assuming the validity of Derrick's theorem, we find that these trapped solutions are stable for 0 < κ < 2 and unstable when κ > 2. At the critical value of κ, i.e. κ = 2, the solution can either collapse or blowup linearly in time when q 0 = 0, where q 0 is the center of the initial density ρ(x, t = 0) = ψ*ψ of the solution. For q 0 ≠ 0 the displaced solution collapses. When κ > 2 initial small displacements from the origin also lead to collapse of the wave function. This phenomenon is not seen in the usual NLSE.

collective coordinates↗

Colloidal Stability of PFSA-Ionomer Dispersions. Part I. Single-Ion Electrostatic Interaction Potential Energies

Charged colloidal particles neutralized by a single counterion are increasingly important for many emerging technologies. Attention here is paid specifically to hydrogen fuel cells and water electrolyzers whose catalyst layers are manufactured from a perfluorinated sulfonic acid polymer (PFSA) suspended in aqueous/alcohol solutions. Partially dissolved PFSA aggregates, known collectively as ionomers, are stabilized by the electrostatic repulsion of overlapping diffuse double layers consisting of only protons dissociated from the suspended polymer. We denote such double layers containing no added electrolyte as "single ion". Size-distribution predictions build upon interparticle interaction potential energies from the Derjaguin-Landau-Verwey-Overbeek (DLVO) formalism. However, when only a single counterion is present in solution, classical DLVO electrostatic potential energies no longer apply. Accordingly, here a new formulation is proposed to describe how single-counterion diffuse double layers interact in colloidal suspensions. Part II (Srivastav, H.; Weber, A. Z.; Radke, C. J. Langmuir 2024 DOI: 10.1021/acs.langmuir.3c03904) of this contribution uses the new single-ion interaction energies to predict aggregated size distributions and the resulting solution pH of PFSA in mixtures of n-propanol and water. A single-counterion diffuse layer cannot reach an electrically neutral concentration far from a charged particle. Consequently, nowhere in the dispersion is the solvent neutral, and the diffuse layer emanating from one particle always experiences the presence of other particles (or walls). Thus, in addition to an intervening interparticle repulsive force, a backside osmotic force is always present. With this new construction, we establish that single-ion repulsive pair interaction energies are much larger than those of classical DLVO electrostatic potentials. The proposed single-ion electrostatic pair potential governs dramatic new dispersion behavior, including dispersions that are stable at a low volume fraction but unstable at a high volume fraction and finite volume-fraction dispersions that are unstable with fine particles but stable with coarse particles. Finally, the proposed single-counterion electrostatic pair potential provides a general expression for predicting colloidal behavior for any charged particle dispersion in ionizing solvents with no added electrolyte.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Collisional excitation of C2H− by H2: New interaction potential and scattering calculations

Interstellar anions play an important role in astrochemistry as being tracers of the physical and chemical conditions in cold molecular clouds and circumstellar gas. The local thermodynamic equilibrium is generally not fulfilled in media where anions are detected and radiative and collisional data are required to model the observed lines. The C2H− anion has not yet been detected in the interstellar medium; however, collisional data could be used for non-LTE models that would help in identifying the most intense lines. For this purpose, we have computed the first 4D potential energy surface (PES) of the C2H−–H2 complex using an explicitly correlated coupled-cluster approach. The PES is characterized by a single deep minimum with a well-depth of 924.96 cm−1. From this interaction potential, we derived excitation cross sections and rate coefficients of C2H− induced by collisions with para- and ortho-H2. The results obtained for collisions with para-H2 are compared to previous calculations performed using a 2D-PES obtained from an average over H2 rotations.

Chemistry↗

An extensive assessment of the performance of pairwise and many-body interaction potentials in reproducing ab initio benchmark binding energies for water clusters n = 2–25

We assess the performance of 7 pairwise additive (TIP3P, TIP4P, TIP4P-ice, TIP5P, OPC, SPC, SPC/E) and 8 families of many-body potentials (q-AQUA, HIPPO, AMOEBA, EFP, TTM, WHBB, MB-pol, MB-UCB) in reproducing high-level ab initio benchmark values, CCSD(T) or MP2 at the complete basis set (CBS) limit for the binding energy and the many-body expansion (MBE) of water clusters n = 2–11, 16–17, 20, 25. By including a large range of cluster sizes having dissimilar hydrogen bonding networks, we obtain an understanding of how these potentials perform for different hydrogen bonding arrangements that are mostly outside of their parameterization range. While it is appropriate to compare the results of ab initio based many-body potentials directly to the electronic binding energies (D e 's), the pairwise additive ones are compared to the enthalpies at T = 298 K, ΔH(298 K), as the latter class of force fields are parametrized to reproduce enthalpies (implicitly accounting for zero-point energy corrections) rather than binding energies. We find that all pairwise additive potentials considered overestimate the reference ΔH values for the n = 2–25 clusters by >13%. For the water dimer (n = 2) in particular, the errors are in the range 83–119% for the pairwise additive potentials studied since these are based on an effective rather than the true 2-body interaction specifically designed as a means of partially accounting for the missing many-body terms. This stronger 2-body interaction is achieved by an enhanced monomer dipole moment that mimics its increase from the gas phase monomer to the condensed phase value. Indeed, for cluster sizes n ≥ 4 the percent deviations become slightly smaller (albeit all exceeding 13%). In contrast, we find that the many-body potentials perform more accurately in reproducing the electronic binding energies (D e 's) throughout the entire cluster range (n = 2–25), all reproducing the ab initio benchmark binding energies within ±7% of the respective CBS values. We further assess the ability of a subset of the many-body potentials (MB-UCB, q-AQUA, MB-pol, and TTM2.1-F) to also reproduce the magnitude of the ab initio many-body energy terms for water cluster sizes n = 7, 10, 16 and 17. The potentials show an overall good agreement with the available benchmark values. However, we identify characteristic differences upon comparing the many-body terms at both the ab initio-optimized geometries and the respective potential-optimized geometries to the reference ab initio values. Additionally, by applying this analysis to a wide range of cluster sizes, trends in the MBE of the potentials with increasing cluster size can be identified. Finally, in an attempt to draw a parallel between the pairwise additive and many-body potentials, we report the analysis of the individual molecular dipole moments for water clusters with 1 to ~4 solvation shells with the TTM2.1-F potential. We find that the internally solvated water molecules have in general a larger molecular dipole moment ranging from 2.6–3.0 D. This justifies the use of an enhanced, with respect to the gas-phase value, molecular dipole moment for the pairwise additive potentials, which is intended to fold in the many body terms into an effective (enhanced) pairwise interaction through the choice of the charges. These results have important implications for the development of future generations of efficient, transferable, and highly accurate classical interaction potentials in both the pairwise additive and many-body categories.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Potential Skyrmion Host Fe(IO3)3: Connecting Stereoactive Lone-Pair Electron Effects to the Dzyaloshinskii-Moriya Interaction

Magnetic skyrmions, which are topologically distinct magnetic spin textures, are gaining increased attention for their unique physical properties and potential applications in spintronic devices. Here we present a design strategy for skyrmion host candidates based on combinations of magnetic spin, asymmetric building units having stereoactive lone-pair electrons, and polar lattice symmetry. To demonstrate the viability of the proposed rational design principles, we successfully synthesized a Fe(IO 3 ) 3 polycrystalline sample and single crystals by using a new simplified low-temperature pathway, which is experimentally feasible for extending materials growth of transition metal iodates. Single crystal X-ray and powder synchrotron X-ray diffraction measurements demonstrated that Fe(IO 3 ) 3 crystallizes in the polar chiral hexagonal lattice with space group P63. The combined structural features of the macroscopic electric polarization along the c-axis stemming from the coalignment of the stereoactive lone-pairs of the IO 3 – trigonal pyramid and the magnetic Fe 3+ cation residing on the 3-fold rotation axis were selected to promote asymmetric exchange coupling. We find evidence of a predicted skyrmion phase at 14 K ≤ T ≤ 16 K and 2.5 T ≤ μ 0 H ≤ 3.2 T driven by a Dzyaloshinskii–Moriya (DM) interaction, a conclusion supported by the appreciable DM exchange and the zero-field spiral antiferromagnetic ground state of Fe(IO 3 ) 3 deduced from neutron diffraction experiments. The associated magnetic modulation wavelength of the putative skyrmions is expected to be short ~18 nm, comparable to the period of the DM-driven incommensurate order. This work links stereoactive lone-pair electron effects to enhanced DM interaction, demonstrating a new approach for chemical guidelines in the search for skyrmionic states of matter.

36 MATERIALS SCIENCE↗

Crystal Nucleation Kinetics and Mechanism: Influence of Interaction Potential

Modulating liquid-to-solid transitions and the resulting crystalline structure for tailored properties is much desired. Colloidal systems are exemplary to this end, but the fundamental knowledge gaps in relating the influence of intermolecular interactions to crystallization behavior continue to hinder progress. In this study, we address this knowledge gap by studying nucleation and growth in systems with modified Lennard-Jones potential. Specifically, we study the commonly used 12-6 potential and a softer 7-6 potential. The thermodynamic state point for the study is chosen such that both systems are investigated at the same level of supercooling and pressure. Under these conditions, we find that the nucleation rate for both systems is comparable. Interestingly, the nucleation pathways and resulting crystal structures are different. In the 12-6 system, nucleation and growth occur predominantly through the FCC structure. Softening the potential alters the critical nucleus composition and introduces two distinct nucleation pathways. One pathway predominantly leads to the nucleus with a body-centered cubic (BCC) structure, while the other favors the face-centered cubic (FCC) arrangement. Our study illustrates that polymorph selection can be achieved through modifications to intermolecular interactions without impacting nucleation kinetics. The results have significant implications in designing approaches for polymorph selection and modulating self-assembly mechanisms.

36 MATERIALS SCIENCE↗

Enhancing the accuracy of density functional tight binding models through ChIMES many-body interaction potentials

Semi-empirical quantum models such as Density Functional Tight Binding (DFTB) are attractive methods for obtaining quantum simulation data at longer time and length scales than possible with standard approaches. However, application of these models can require lengthy effort due to the lack of a systematic approach for their development. In this work, we discuss the use of the Chebyshev Interaction Model for Efficient Simulation (ChIMES) to create rapidly parameterized DFTB models, which exhibit strong transferability due to the inclusion of many-body interactions that might otherwise be inaccurate. We apply our modeling approach to silicon polymorphs and review previous work on titanium hydride. We also review the creation of a general purpose DFTB/ChIMES model for organic molecules and compounds that approaches hybrid functional and coupled cluster accuracy with two orders of magnitude fewer parameters than similar neural network approaches. In all cases, DFTB/ChIMES yields similar accuracy to the underlying quantum method with orders of magnitude improvement in computational cost. In conclusion, our developments provide a way to create computationally efficient and highly accurate simulations over varying extreme thermodynamic conditions, where physical and chemical properties can be difficult to interrogate directly, and there is historically a significant reliance on theoretical approaches for interpretation and validation of experimental results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Potential Interactions of Novec 1230 and Its Thermal Degradation Products with Actinide Metals and Oxides

The environmentally friendly and mainly chemically inert fire suppressant, Novec 1230, is being increasingly used worldwide. However, no published study has considered using Novec 1230 with fires involving radioactive material including actinides; here, therefore, this research note focuses on the possible interactions between Novec 1230 and its thermal degradation products (TDPs) with some actinide compounds (e.g., actinide metal and oxides) commonly present in a radiological facility in the event of a fire. Previous studies of halogenated fire suppressants used in actinide metal fires indicate a possibility for Novec 1230 or its TDPs to chemically interact with actinide compounds at elevated temperatures; however, these reactions are highly unlikely to cause either runaway exothermic reactions or gaseous actinide release.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Collisional excitation of interstellar PN by H 2 : New interaction potential and scattering calculations

Here, rotational excitation of interstellar PN molecules induced by collisions with H 2 is investigated. We present the first ab initio four-dimensional potential energy surface (PES) for the PN–H 2 van der Waals system. The PES was obtained using an explicitly correlated coupled cluster approach with single, double, and perturbative triple excitations [CCSD(T)-F12b]. The method of interpolating moving least squares was used to construct an analytical PES from these data. The equilibrium structure of the complex was found to be linear, with H 2 aligned at the N end of the PN molecule, at an intermolecular separation of 4.2 Å. The corresponding well-depth is 224.3 cm -1 . The dissociation energies were found to be 40.19 cm -1 and 75.05 cm -1 for complexes of PN with ortho-H 2 and para-H 2 , respectively. Integral cross sections for rotational excitation in PN–H 2 collisions were calculated using the new PES and were found to be strongly dependent on the rotational level of the H 2 molecule. These new collisional data will be crucial to improve the estimation of PN abundance in the interstellar medium from observational spectra.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spliced soft-core interaction potential for filling small-scale enclosures

A molecular dynamics system employing a spliced soft-core potential (the “MD/SSCP”) facilitates studying the arrangement of particles in an enclosure and/or near an interface. In some embodiments, the MD/SSCP initializes a three-dimensional representation containing the enclosure and the particles in a first arrangement. The MD/SSCP conducts a first simulation to transition the representation to a second arrangement, during which the particles are allowed to move through a wall of the enclosure while the SSCP is unengaged. The MD/SSCP conducts a second simulation to transition the representation to a third arrangement, during which it becomes more difficult for the particles to move through the wall of the enclosure while the SSCP is gradually engaged. The MD/SSCP conducts a third simulation to transition the representation to a fourth arrangement of the particles, during which it becomes almost impossible for the particles to move through the wall of the enclosure.

Meshot, Eric↗