Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Interatomic potentials”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Phase Boundaries, Isotope Effect and Superconductivity of Lithium Under Hydrostatic Conditions

We present theoretical and experimental studies of superconductivity and low temperature structural phase boundaries in lithium. We mapped the structural phase diagram of 6Li and 7Li under hydrostatic conditions between 5 – 55 GPa and within the temperature range of 15 – 75 K, observing the FCC → hR1 → cI16 phase transitions. 6Li and 7Li show some differences at the structural boundaries, with a potential shift of the phase boundaries of 6Li to lower pressures. Density functional theory calculations and topological analysis of the electron density elucidates the superconducting properties and interatomic interactions within these phases of lithium.

Deemyad, Shanti↗
NASA NTRS2002

Equivalent-Continuum Modeling With Application to Carbon Nanotubes

A method has been proposed for developing structure-property relationships of nano-structured materials. This method serves as a link between computational chemistry and solid mechanics by substituting discrete molecular structures with equivalent-continuum models. It has been shown that this substitution may be accomplished by equating the vibrational potential energy of a nano-structured material with the strain energy of representative truss and continuum models. As important examples with direct application to the development and characterization of single-walled carbon nanotubes and the design of nanotube-based devices, the modeling technique has been applied to determine the effective-continuum geometry and bending rigidity of a graphene sheet. A representative volume element of the chemical structure of graphene has been substituted with equivalent-truss and equivalent continuum models. As a result, an effective thickness of the continuum model has been determined. This effective thickness has been shown to be significantly larger than the interatomic spacing of graphite. The effective thickness has been shown to be significantly larger than the inter-planar spacing of graphite. The effective bending rigidity of the equivalent-continuum model of a graphene sheet was determined by equating the vibrational potential energy of the molecular model of a graphene sheet subjected to cylindrical bending with the strain energy of an equivalent continuum plate subjected to cylindrical bending.

Odegard, Gregory M.↗
NASA NTRS1982

Determination of accurate dissociation limits and interatomic interactions at large internuclear distances

This paper illustrates (using the molecules LiH, Li2 and Na2) how laser-induced fluorescence can be used to greatly expand the range of observed vibrational levels in ground electronic states. This expanded vibrational range leads to the determination of virtually the full well of the potential energy curve. This also leads to improved determination of the dissociation limit and serves as a severe test for highly accurate ab initio calculations now available for many small molecules.

Stwalley, W. C.↗
NASA NTRSDate not supplied

Machine-Learned Committor Functions for Reactive Molecular Dynamics

Reactive molecular dynamics (MD) is a powerful tool for atomistic-scale modeling of a diverse range of chemical processes. However, scaling these simulations to large systems and long times scales remains a challenge because of the complexity of the potential energy function required. The authors previously developed a heuristic approach, called REACTER, that incorporates reactivity in MD simulations in a less general but much more computationally efficient manner. REACTER uses standard, fixed valence force fields as the underlying potentialenergy surface for describing all interatomic interactions but adds a procedure for enforcing user-defined reactions that occur when certain geometric constraints on relative atomic positions are satisfied. Further, these bonding changes can be accepted or rejected with a probability related tothe local thermal energy. This work seeks to generalize this approach by replacing the set of user defined geometric constraints and energetic criteria with a committor function that specifies the probability of a reaction occurring on the basis of the local atomic configuration. The committor function is a useful mathematical tool for modeling rare events but, unfortunately, is very difficult to compute for realistic systems in a general way. This work describes a method for approximating the committor function using a machine learning approach, specifically a deep neural network trained with data from reactive MD and DFT-based dynamics simulations. This network is coupled to the existing REACTER protocol, as implemented in the LAMMPS MD package, and used to make on-the-fly predictions of reaction probabilities without the more extensive user input previously required. The new method is demonstrated using the polymerization of polystyrene as a case study. Although very dependent on the quality and quantity of training data, machine-learned committor functions show promise as a method for incorporating reaction probability from higher level calculations into highly scalable MD simulations.

polymer simulations↗
DOE OSTI2025

Experimental confirmation of first-principles thermal conductivity in Zirconium-doped ThO 2

The degradation of thermal conductivity in advanced nuclear fuels due to the accumulation of fission products and irradiation-induced defects is inevitable, and must be considered as part of safety and efficiency analyses of nuclear reactors. Here, this study examines the thermal conductivity of a zirconium-doped ThO 2 crystal, synthesized via the hydrothermal method using a spatial domain thermoreflectance technique. Zirconium is one of the soluble fission products in oxide fuels that can effectively scatter heat-carrying phonons in the crystalline lattice of fuel. Thus, thermal property measurements of zirconium-doped ThO 2 single crystals provide insights into the effects of substitutional zirconium doping, isolated from extrinsic factors such as grain boundary scattering. The experimental results are compared with first-principles calculations of the lattice thermal conductivity of ThO 2 , employing an iterative solution of the Peierls-Boltzmann transport equation. Additionally, the non-perturbative Green's function methodology is utilized to compute phonon-point defect scattering rates, accounting for local distortions around point defects, including mass difference changes, interatomic force constants, and structural relaxation. The congruence between the predicted results from first-principles calculations and the measured temperature-dependent thermal conductivity validates the computational methodology. Furthermore, the methodologies employed in this study enable systematic investigations of thermal conductivity reduction by fission products, potentially leading to the development of more accurate fuel performance codes.

36 - MATERIALS SCIENCE↗
DOE OSTI2022

Effect of Redox-Inactive Metal Ion–Nickel(III) Interactions on the Redox Properties and Proton-Coupled Electron Transfer Reactivity

We report Mononuclear nickel(II) and nickel(III) complexes of a bisamidate-bisalkoxide ligand, (NMe 4 ) 2 [Ni- II (HMPAB)] (1) and (NMe 4 )[Ni III (HMPAB)] (2), respectively, have been synthesized and characterized by various spectroscopic techniques including X-ray crystallography. The reaction of redox-inactive metal ions (M n+ = Ca 2+ , Mg 2+ , Zn 2+ , Y 3+ , and Sc 3+ ) with 2 resulted in 2-M n+ adducts, which was assessed by an array of spectroscopic techniques including X-ray absorption spectroscopy (XAS), electron paramagnetic resonance (EPR), and reactivity studies. The X-ray structure of Ca 2+ coordinated to Ni(III) complexes, 2Ca 2+ T, was determined and exhibited an average Ni-Ca distance of 3.1253 angstrom, close to the metal ions' covalent radius. XAS analysis of 2-Ca 2+ and 2-Y 3+ in solution further revealed an additional coordination to Ca and Y in the 2-M n+ adducts with shortened Ni-M distances of 2.15 and 2.11 Å, respectively, implying direct bonding interactions between Ni and Lewis acids (LAs). Such a short interatomic distance between Ni(III) and M is unprecedented and was not observed before. EPR analysis of 2 and 2-M n+ species, moreover, displayed rhombic signals with g av > 2.12 for all complexes, supporting the +III oxidation state of Ni. The Ni III /Ni II redox potential of 2 and 2-M n+ species was determined, and a plot of E 1/ 2 of 2-M n+ versus pK a of [M(H 2 O) n ] m+ exhibited a linear relationship, implying that the Ni III /Ni II potential of 2 can be tuned with different redox-inactive metal ions. Reactivity studies of 2 and 2-M n+ with different 4-X-2,6-ditert-butylphenol (4-XDTBP) and other phenol derivatives were performed, and based on kinetic studies, we propose the involvement of a proton-coupled electron transfer (PCET) pathway. Analysis of the reaction products after the reaction of 2 with 4-OMe-DTBP showed the formation of a Ni(II) complex (1a) where one of the alkoxide arms of the ligand is protonated. A pK a value of 24.2 was estimated for 1a. The reaction of 2-M n+ species was examined with 4-OMe-DTBP, and it was observed that the k 2 values of 2-M n+ species increase by increasing the Lewis acidity of redox-inactive metal ions. However, the obtained k 2 values for 2-M n+ species are much lower compared to the k 2 value for 2. Such a variation of PCET reactivity between 2 and 2-M n+ species may be attributed to the interactions between Ni(III) and LAs. Our findings show the significance of the secondary coordination sphere effect on the PCET reactivity of Ni(III) complexes and furnish important insights into the reaction mechanism involving high-valent nickel species, which are frequently invoked as key intermediates in Ni-mediated enzymatic reactions, solar-fuel catalysis, and biomimetic/synthetic transformation reactions.

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

Periodic GFN1-xTB Tight Binding: A Generalized Ewald Partitioning Scheme for the Klopman–Ohno Function

A novel formulation is presented for the treatment of electrostatics in the periodic GFN1-xTB tight-binding model. Periodic GFN1-xTB is hindered by the functional form of the second-order electrostatics, which only recovers Coulombic behavior at large interatomic distances and lacks a closed-form solution for its Fourier transform. We address this by introducing a binomial expansion of the Klopman–Ohno function to partition short- and long-range interactions, enabling the use of a generalized Ewald summation for the solution of the electrostatic energy. This approach is general and is applicable to any damped potential of the form |R n + c| –m . Benchmarks on the X23 molecular crystal dataset and a range of prototypical bulk semiconductors demonstrate that this systematic treatment of the electrostatics eliminates unphysical behavior in the equation of state curves. In the bulk systems studied, we observe a mean absolute error in total energy of 35 meV/atom, comparable to the machine-learned universal force field, M3GNet, and sufficiently precise for structure relaxation. These results highlight the promising potential of GFN1-xTB as a universal tight-binding parametrization.

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

Local Atomic Configuration in Pristine and A-Site Doped Silver Niobate Perovskite Antiferroelectrics

Antiferroelectrics have attracted increasing research interests in recent years due to both their great potential in energy storage applications and intriguing structural characteristics. However, the links between the electrical properties and structural characteristics of distorted perovskite antiferroelectrics are yet to be fully deciphered. Here, we adopt local-structure methods to elucidate the nanoscale atomic structure of AgNbO3-based antiferroelectrics and their structural evolution upon La doping. The local structural features including interatomic distance distributions and atomic displacements have been analyzed using neutron small-box pair distribution function (PDF) refinement in conjunction with large-box Reverse Monte Carlo modelling. Our results highlight the correlation of cation displacements in AgNbO 3 and its disruption by the incorporation of La, apparently in corroboration with the observed anomalous dielectric properties. Spatial ordering of cation vacancies is observed in La-doped AgNbO 3 samples, which coordinates with oxygen octahedral tilting to relieve lattice strain. These results provide renewed insights into the atomic structure and antiferroelectric phase instabilities of AgNbO 3 and relevant perovskite materials, further lending versatile opportunities for enhancing their functionalities.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Complexity of many-body interactions in transition metals via machine-learned force fields from the TM23 data set

Abstract This work examines challenges associated with the accuracy of machine-learned force fields (MLFFs) for bulk solid and liquid phases ofd-block elements. In exhaustive detail, we contrast the performance of force, energy, and stress predictions across the transition metals for two leading MLFF models: a kernel-based atomic cluster expansion method implemented using sparse Gaussian processes (FLARE), and an equivariant message-passing neural network (NequIP). Early transition metals present higher relative errors and are more difficult to learn relative to late platinum- and coinage-group elements, and this trend persists across model architectures. Trends in complexity of interatomic interactions for different metals are revealed via comparison of the performance of representations with different many-body order and angular resolution. Using arguments based on perturbation theory on the occupied and unoccupieddstates near the Fermi level, we determine that the large, sharpddensity of states both above and below the Fermi level in early transition metals leads to a more complex, harder-to-learn potential energy surface for these metals. Increasing the fictitious electronic temperature (smearing) modifies the angular sensitivity of forces and makes the early transition metal forces easier to learn. This work illustrates challenges in capturing intricate properties of metallic bonding with current leading MLFFs and provides a reference data set for transition metals, aimed at benchmarking the accuracy and improving the development of emerging machine-learned approximations.

Chemistry↗
DOE OSTI2023

Reconsidering short-range order in complex concentrated alloys

Abstract The seemingly contradictory state of research on short-range order in many-component alloys is addressed through a critical review of the characterization of face-centered-cubic 3 d systems. Despite the paucity of direct observations, the ordering of many widely studied alloys is argued to be primarily interesting for its potential ubiquity. To clarify this situation, future research directions are proposed with reference to historical results, including a review of the fundamental principles of ordering and clustering. Graphical abstract

36 MATERIALS SCIENCE↗
NASA NTRS1990

The relation between ferroelasticity and superconductivity

The high-temperature superconductivity is explained widely by the layered crystal structure. The one- and two-dimensional subsystems and their interaction are investigated here. It is assumed that the high-T(sub c) superconductivity takes place in the two-dimensional subsystem and the increase of the phase transition temperature from 60 K up to 90 K is the consequence of turning on the influence of one-dimensional chains. The interaction between the two subsystems is transferred along the c axis by the phonons of breathing mode, which causes the hybridization of the electronic bonds between these subsystems. The experimental works indicate that the existence of both the chains Cu(1)-O and their interaction with the superconducting plane of Cu(2)-O modify the temperature of the transition to the superconducting state. It is seen from the neutron scattering data that the rates of the interatomic distance dependencies on temperature are changed around 140 K and 90 K. The 'zig-zag' order in Cu(1)-O chains has been postulated but, on the other hand, the vibrations with a large amplitude only were reported. The bi-stabilized situation of the oxygen ions can be caused by the change of distance between these ions and the Ba ions. It leads to the appearance of a two-well potential. Its parameters depend on temperature and the dynamics of the oxygen ions' movement. They can induce the antipolar order, which can be, however, more or less chaotic. The investigation of the ferroelastic properties of Y-Ba-Cu-O samples lead to the conclusion that they are related to jumps of ions inside the given chain and not to a diffusion between different sites in the ab plane. Researchers deduce thus that the fluctuating oxygen ions from these chains create dipoles in the ab plane. They can be described with the pseudo-spin formalism/ - Pauli matrices/. The system can be described with the Ising model. The pseudo-spins interact with phonons and influence the superconductivity in the second subsystem.

Molak, A.↗
NASA NTRS1991

The relation between ferroelasticity and superconductivity

The high-temperature superconductivity is explained widely by the layered crystal structure. The one- and two-dimensional subsystems and their interaction are investigated here. It is assumed that the high-T(sub c) superconductivity takes place in the two-dimensional subsystem and the increase of the phase transition temperature from 60 K up to 90 K is the consequence of turning on the influence of one-dimensional chains. The interaction between the two subsystems is transferred along the c axis by the phonons of breathing mode, which causes the hybridization of the electronic bonds between these subsystems. The experimental works indicate that the existence of both the chains Cu(1)-O and their interaction with the superconducting plane of Cu(2)-O modify the temperature of the transition to the superconducting state. It is seen from the neutron scattering data that the rates of the interatomic distance dependencies on temperature are changed around 240 K and 90 K. The 'zig-zag' order in Cu(1)-O chains has been postulated but, on the other hand, the vibrations with a large amplitude only were reported. The bi-stabilized situation of the oxygen ions can be caused by the change of distance between these ions and the Ba ions. It leads to the appearance of a two-well potential. Its parameters depend on temperature and the dynamics of the oxygen ions' movement. They can induce the antipolar order, which can be, however, more or less chaotic. The investigation of the ferroelastic properties of Y-Ba-Cu-O samples lead to the conclusion that they are related to jumps of ions inside the given chain and not to a diffusion between different sites in the ab plane. Researchers deduce, thus, that the fluctuating oxygen ions from these chains create dipoles in the ab plane. They can be described with the pseudo-spin formalism (- Pauli matrices). The system can be described with the Ising model. The pseudo-spins interact with phonons and influence the superconductivity in the second subsystem.

Molak, A.↗
DOE OSTI2020

Ambient and High Pressure CuNiSb 2 : Metal-Ordered and Metal-Disordered NiAs-Type Derivative Pnictides

The mineral Zlatogorite, CuNiSb 2 , was synthesized in the laboratory for the first time by annealing elements at ambient pressure (CuNiSb 2 -AP). Rietveld refinement of synchrotron powder X-ray diffraction data indicates that CuNiSb 2 -AP crystallizes in the NiAs-derived structure ( P 3 m 1, #164) with Cu and Ni ordering. The structure consists of alternate NiSb 6 and CuSb 6 octahedral layers via face-sharing. The formation of such structure instead of metal disordered NiAs-type structure ( P 6 3 / mm c, #194) is validated by the lower energy of the ordered phase by first-principle calculations. Interatomic crystal orbital Hamilton population, electron localization function, and charge density analysis reveal strong Ni-Sb, Cu-Sb, and Cu-Ni bonding and long weak Sb-Sb interactions in CuNiSb 2 -AP. The magnetic measurement indicates that CuNiSb 2 -AP is Pauli paramagnetic. First-principle calculations and experimental electrical resistivity measurements reveal that CuNiSb 2 -AP is a metal. The low Seebeck coefficient and large thermal conductivity suggest that CuNiSb 2 is not a potential thermoelectric material. Single crystals were grown by chemical vapor transport. The high pressure sample (CuNiSb 2 -8 GPa) was prepared by pressing CuNiSb 2 -AP at 700 °C and 8 GPa. However, the structures of single crystal and CuNiSb 2 -8 GPa are best fit with a disordered metal structure in the P 3 m 1 space group, corroborated by transmission electron microscopy.

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

Minimization of Disorder as a Key Design Principle for Natural Sizes of Light Harvesting 2 Complexes

The light harvesting 2 (LH2) complex of purple bacteria has excellent energy conversion efficiency. Clarifying the design principle behind such efficiency at the atomistic level is crucial for understanding its structure–function relationship and can be utilized for the design of artificial light harvesting systems. To this end, we conducted comprehensive computational investigation of the dynamical and statistical nature of electronic excited states of pigment molecules in a natural LH2 complex with 9-fold symmetry and its two non-natural in silico analogues with 6- and 12-fold symmetries. To ensure reliable and efficient all-atomistic molecular dynamics simulations, we combined a well established interpolation approach for the construction of the potential energy surface with a neural network machine learning approach. Outcomes of these calculations clarify that non-natural forms of LH2-type complexes have significantly larger quasistatic disorder than those for the natural one. In addition, non-natural systems have more disruptions of the hydrogen bonding, underscoring its crucial role for reducing the disorder. On the other hand, local environmental dynamics are relatively insensitive to the structural changes although there is moderate enhancement in the anharmonic or interatomic components for the synthetic ones. These findings based on all-atomistic simulations provide direct computational evidence that the structure and sizes of natural LH2 complexes are designed to minimize the energetic disorder. We analyze quantitative implications of these for the energy transferring capability of the LH2 complex.

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

Multitarget Rydberg gates via spatial blockade engineering

Multi-target gates offer the potential to reduce gate depth in syndrome extraction for quantum error correction. Although neutral-atom quantum computers have demonstrated native multi-qubit gates, existing approaches that avoid additional control or multiple atomic species have been limited to single-target gates. We propose single-control-multi-target CZ^n gates on a single-species neutral-atom platform that require no extra control and have gate durations comparable to standard CZ gates. Our approach leverages tailored interatomic distances to create an asymmetric blockade between the control and target atoms. Using a GPU-accelerated pulse synthesis protocol, we design smooth control pulses for CZZ and CZZZ gates, achieving fidelities of up to 99.55% and $99.24\%$, respectively, even in the presence of simulated atom placement errors and Rydberg-state decay. Our approach is most effective for N=2 (CZZ) and N=3 targets (CZZZ); for larger N, increasing spatial crowding of the targets introduces significant challenges for maintaining the required blockade asymmetry. This work presents a practical path to implementing low-overhead multi-target gates in single-species neutral-atom systems, significantly reducing the resource overhead for syndrome extraction. To motivate the impact of these gates, we apply a greedy scheduling algorithm and we demonstrate that our proposed gates can reduce the number of atom reconfiguration costs by up to 50% for color code syndrome extraction of code distances greater than 5.

Stein, Samuel A.↗
DOE OSTI2023

Architector for high-throughput cross-periodic table 3D complex building

Abstract Rare-earth and actinide complexes are critical for a wealth of clean-energy applications. Three-dimensional (3D) structural generation and prediction for these organometallic systems remains a challenge, limiting opportunities for computational chemical discovery. Here, we introduce Architector , a high-throughput in-silico synthesis code for s-, p-, d-, and f-block mononuclear organometallic complexes capable of capturing nearly the full diversity of the known experimental chemical space. Beyond known chemical space, Architector performs in-silico design of new complexes including any chemically accessible metal-ligand combinations. Architector leverages metal-center symmetry, interatomic force fields, and tight binding methods to build many possible 3D conformers from minimal 2D inputs including metal oxidation and spin state. Over a set of more than 6,000 x-ray diffraction (XRD)-determined complexes spanning the periodic table, we demonstrate quantitative agreement between Architector-predicted and experimentally observed structures. Further, we demonstrate out-of-the box conformer generation and energetic rankings of non-minimum energy conformers produced from Architector , which are critical for exploring potential energy surfaces and training force fields. Overall, Architector represents a transformative step towards cross-periodic table computational design of metal complex chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2001

Equivalent-Continuum Modeling of Nano-Structured Materials

A method has been developed for modeling structure-property relationships of nano-structured materials. This method serves as a link between computational chemistry and solid mechanics by substituting discrete molecular structures with an equivalent-continuum model. It has been shown that this substitution may be accomplished by equating the vibrational potential energy of a nano-structured material with the strain energy of representative truss and continuum models. As an important example with direct application to the development and characterization of single-walled carbon nanotubes, the model has been applied to determine the effective continuum geometry of a graphene sheet. A representative volume element of the equivalent-continuum model has been developed with an effective thickness. This effective thickness has been shown to be similar to, but slightly smaller than, the interatomic spacing of graphite.

Odegard, Gregory M.↗
NASA NTRS1989

Correlation of molecular valence- and K-shell photoionization resonances with bond lengths

The relationship between the interatomic distance and the positions of valence-shell and K-shell sigma(asterisk) photoionization resonances is investigated theoretically for the molecules C2, F2, N2, O2, CO, NO, C2H2, C2H4, C2H6, HCN, H2CO, N20, CO2, and C2N2. The results of molecular-orbital computations are presented in three-dimensional diagrams, which are shown to be similar to the wave functions of a particle in a cylindrical well, confirming the validity of free-electron molecular-orbital (FEMO) approximations for modeling the potential along the symmetry axis. FEMO orbital energies and resonance positions are found to be in good agreement with previous theoretical and experimental results. Also included is a Feshbach-Fano analysis of the relevance of virtual-valence orbitals to the appearance of single-channel resonances in molecular photoionization cross sections.

Sheehy, J. A.↗