Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distribution function”

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

From Cation Order to Disorder: Unlocking Ion Transport Pathways in Li–Zn–Zr–Cl Halospinels

Lithium metal chloride halospinels of the general formula Li 2 MCl 4 are a promising class of earth-abundant ion conductors for all-solid-state batteries. However, poor roomtemperature ionic conductivity has historically limited their use in practical applications. Here, we substitute Zr 4+ into Li 2 ZnCl 4 along the series Li 2−2x/3 Zn 1−x Zr 2x/3 Cl 4 (x = 0, 0.1, 0.3, 0.4, 0.6, 0.9, and 1.0) to understand how cation disorder and vacancy tuning impact ion transport in “normal” halospinels. Aliovalent Zr 4+ substitution increases ionic conductivity by nearly 5 orders of magnitude, from 1.320(3) × 10 −9 S cm −1 in Li 2 ZnCl 4 to 6.74(1)× 10 −5 S cm −1 for x = 0.6. Average and local structure characterization through synchrotron X-ray diffraction (SXRD) and neutron pair distribution function (nPDF) analysis reveal that Zr 4+ redistributes the Zn 2+ and Li + sublattices into previously unoccupied interstitial sites, which form new low-energy hopping pathways that facilitate ion transport. We rationalize the dramatic rearrangement of the cation local structure by considering the coordination preferences of the cations and the potential electrostatic penalties incurred by the higher-valent Zr 4+ cations. This work delivers an atomistic understanding of substitution-induced cation disorder and ion transport properties in a new family of earth-abundant halospinels.

Cardoza, Abby M. [Colorado School of Mines, Golden↗
DOE OSTI2021

Competitive TcO4-, IO3-, and CrO42- Incorporation into Ettringite

Ettringite is a naturally occurring mineral found in cementitious matrices that is known for its ability to incorporate environmentally mobile oxyanion contaminants. To better assess this immobilization mechanism for contaminants within cementitious waste forms intended for nuclear waste storage, this work explores how mixed oxyanion contaminants compete for ettringite incorporation and influence the evolving mineralogy. Ettringite was precipitated in the presence of TcO4-, IO3-, and/or CrO42-, known contaminants of concern to nuclear waste treatment, over pre-determined precipitation periods. Solution analyses quantified contaminant removal and the collected solid was characterized using bulk and microprobe XRD coupled with PDF and microprobe XRF analyses. Results suggest that =96% IO3- is removed from solution, regardless of ettringite precipitation time or the presence of TcO4- or CrO42-. However, TcO4- removal remained <20%, was not significantly improved with longer ettringite precipitation times, and significantly decreased in the presence of IO3-. When IO3- is comingled with CrO42, gypsum is formed as a secondary mineral phase, which results in oxyanion partitioning among mineral phases, e.g., IO3- incorporation into ettringite and CrO42- incorporation into gypsum via SO42- substitution. Results from this work exemplify the importance of competitive immobilization when assessing waste form performance and environmental risk of contaminant release.

Gillispie, Elizabeth C.↗
DOE OSTI2020

Short-Range Order in VI 3

We present a detailed investigation of the crystal structure of VI 3 , a two-dimensional van der Waals material of interest for studies of low-dimensional magnetism. As opposed to the average crystal structure that features $R\bar{3}$ symmetry of the unit cell, our Raman scattering and X-ray atomic pair distribution function analysis supported by density functional theory calculations point to the coexistence of short-range ordered $P\bar{3}$1c and long-range ordered $R\bar{3}$phases. The highest-intensity peak, A 1g 3 , exhibits a moderate asymmetry that might be traced back to the spin–phonon interactions, as in the case of CrI 3 .

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

Structural Insights into the Vapochromic Behavior of Pt- and Pd-Based Compounds

Anionic complexes having vapochromic behavior are investigated: [K(H 2 O)][M(ppy)(CN) 2 ], [K(H 2 O)][M(bzq)(CN) 2 ], and [Li(H 2 O) n ][Pt(bzq)(CN) 2 ], where ppy = 2-phenylpyridinate, bzq = 7,8-benzoquinolate, and M = Pt(II) or Pd(II). These hydrated potassium/lithium salts exhibit a change in color upon being heated to 380 K, and they transform back into the original color upon absorption of water molecules from the environment. The challenging characterization of their structure in the vapochromic transition has been carried out by combining several experimental techniques, despite the availability of partially ordered and/or impure crystalline material. Room-temperature single-crystal and powder X-ray diffraction investigation revealed that [K(H 2 O)][Pt(ppy)(CN) 2 ] crystallizes in the Pbca space group and is isostructural to [K(H 2 O)][Pd(ppy)(CN) 2 ]. Variable-temperature powder X-ray diffraction allowed the color transition to be related to changes in the diffraction pattern and the decrease in sample crystallinity. Water loss, monitored by thermogravimetric analysis, occurs in two stages, well separated for potassium Pt compounds and strongly overlapped for potassium Pd compounds. Furthermore, the local structure of potassium compounds was monitored by in situ pair distribution function (PDF) measurements, which highlighted changes in the intermolecular distances due to a rearrangement of the crystal packing upon vapochromic transition.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The Role of Phonons and Oxygen Vacancies in Non-Cubic SrVO 3

Combining neutron diffraction with pair distribution function analysis, we have uncovered hidden reduced symmetry in the correlated metallic d 1 perovskite, SrVO 3 . Specifically, we show that both the local and global structures are better described using a GdFeO 3 distorted (orthorhombic) model as opposed to the ideal cubic ABO 3 perovskite type. Recent reports of imaginary phonon frequencies in the density functional theory (DFT)-calculated phonon dispersion for cubic SrVO 3 suggest a possible origin of this observed non-cubicity. Namely, the imaginary frequencies computed could indicate that the cubic crystal structure is unstable at T = 0 K. However, our DFT calculations provide compelling evidence that point defects in the form of oxygen vacancies, and not an observable symmetry breaking associated with calculated imaginary frequencies, primarily result in the observed non-cubicity of SrVO 3 . These experimental and computational results are broadly impactful because they reach into the thin-film and theoretical communities who have shown that SrVO 3 is a technologically viable transparent conducting oxide material and have used SrVO 3 to develop theoretical methods, respectively.

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

In Situ Visualization of Local Distortions in the High- T c Molecule-Intercalated Li x (C 5 H 5 N) y Fe 2– z Se 2 Superconductor

A time-resolved synchrotron X-ray total scattering study sheds light on the evolution of the different structural length scales involved during the intercalation of the layered iron–selenide host by organic molecular donors, aiming at the formation of the expanded-lattice Li x (C 5 H 5 N) y Fe 2– z Se 2 hybrid superconductor. The intercalates are found to crystallize in the tetragonal ThCr 2 Si 2 -type structure at the average level, however, with an enhanced interlayer iron–selenide spacing ( d = 16.2 Å) that accommodates the heterocyclic molecular spacers. Quantitative atomic pair distribution function (PDF) analysis at variable times suggests distorted FeSe 4 tetrahedral local environments that appear swollen with respect to those in the parent β-FeSe. Simultaneously acquired in situ synchrotron X-ray powder diffraction data disclose that secondary phases (α-Fe and Li 2 Se) grow significantly when a higher lithium concentration is used in the solvothermal reaction or when the solution is aged. These observations are in line with the strongly reducing character of the intercalation medium’s solvated electrons that mediate the defect chemistry of the expanded-lattice superconductor. In the latter, intralayer correlated local distortions indicate electron-donating aspects that reflect in somewhat enlarged Fe–Se bonds. Additionally, they also reveal a degree of relief of chemical pressure associated with a large distance between Fe and Se sheets (“taller” anion height) and a stretched Fe–Fe square planar topology. The elongation of the latter, derived from the in situ PDF study, speaks for a plausible increase in the Fe-site vacancy concentration. The evolution of the local structural parameters suggests an optimum reaction window where kinetically stabilized phases resemble the distortions of the edge-sharing Fe–Se tetrahedra, required for a high- T c in expanded-lattice iron-chalcogenides.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

KBKit: A Python Toolkit for Kirkwood–Buff Theory from Molecular Dynamics

Thermodynamic properties of liquid mixtures govern processes that range from drug delivery to energy storage, yet extracting these properties from molecular simulations remains challenging. Kirkwood–Buff (KB) theory offers a rigorous route by linking microscopic pair distribution functions to macroscopic free energies, but practical use of the theory has been hindered by two obstacles: (i) the long simulations needed to obtain well-converged Kirkwood-Buff integrals (KBIs) and (ii) the specialized corrections required to translate finite-size data to the thermodynamic limit. $\texttt{KBKit}$ is an open-source Python package that removes these barriers. It automatically computes KBIs and derived thermodynamic quantities from GROMACS input files, applies state-of-the-art finite-size corrections, and provides built-in diagnostic tools to quantify statistical uncertainty. Written with modern software-engineering practices—continuous integration, extensive unit testing, and thorough documentation—$\texttt{KBKit}$ is both reliable and easy to extend. By condensing complex KBI analysis into a few intuitive commands, $\texttt{KBKit}$ enables researchers to incorporate KB theory into routine simulation workflows and accelerate the discovery of solution-phase thermodynamics.

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

Machine Learning Potentials with the Iterative Boltzmann Inversion: Training to Experiment

Methodologies for training machine learning potentials (MLPs) with quantum-mechanical simulation data have recently seen tremendous progress. Experimental data have a very different character than simulated data, and most MLP training procedures cannot be easily adapted to incorporate both types of data into the training process. Here, we investigate a training procedure based on iterative Boltzmann inversion that produces a pair potential correction to an existing MLP using equilibrium radial distribution function data. By applying these corrections to an MLP for pure aluminum based on density functional theory, we observe that the resulting model largely addresses previous overstructuring in the melt phase. Interestingly, the corrected MLP also exhibits improved performance in predicting experimental diffusion constants, which are not included in the training procedure. The presented method does not require autodifferentiating through a molecular dynamics solver and does not make assumptions about the MLP architecture. Our results suggest a practical framework for incorporating experimental data into machine learning models to improve the accuracy of molecular dynamics simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Using Computationally-Determined Properties for Machine Learning Prediction of Self-Diffusion Coefficients in Pure Liquids

The ability to predict transport properties of liquids quickly and accurately will greatly improve our understanding of fluid properties both in bulk and complex mixtures, as well as in confined environments. Such information could then be used in the design of materials and processes for applications ranging from energy production and storage to manufacturing processes. As a first step, we consider the use of machine learning (ML) methods to predict the diffusion properties of pure liquids. Recent results have shown that Artificial Neural Networks (ANNs) can effectively predict the diffusion of pure compounds based on the use of experimental properties as the model inputs. In the current study, a similar ANN approach is applied to modeling diffusion of pure liquids using fluid properties obtained exclusively from molecular simulations. A diverse set of 102 pure liquids is considered, ranging from small polar molecules (e.g., water) to large nonpolar molecules (e.g., octane). Self-diffusion coefficients were obtained from classical molecular dynamics (MD) simulations. Since nearly all the molecules are organic compounds, a general set of force field parameters for organic molecules was used. The MD methods are validated by comparing physical and thermodynamic properties with experiment. Computational input features for the ANN include physical properties obtained from the MD simulations as well as molecular properties from quantum calculations of individual molecules. Furthermore, fluid properties describing the local liquid structure were obtained from center of mass radial distribution functions (COM-RDFs). Feature sensitivity analysis revealed that isothermal compressibility, heat of vaporization, and the thermal expansion coefficient were the most impactful properties used as input for the ANN model to predict the MD simulated self-diffusion coefficients. The MD-based ANN successfully predicts the MD self-diffusion coefficients with only a subset (2 to 3) of the available computationally determined input features required. A separate ANN model was developed using literature experimental self-diffusion coefficients as model targets. Although this second ML model was not as successful due to a limited number of data points, a good correlation is still observed between experimental and ML predicted self-diffusion coefficients.

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

Relation between the Hydrated Electron Solvation Structure and Its Partial Molar Volume

It is now generally accepted that the hydrated electron occupies a cavity in water, but the size of the cavity and the arrangements of the solvating water molecules are not fully characterized. Here, we use the Kirkwood-Buff (KB) approach to examine how the partial molar volume (V M ) provides insight into these issues. The KB method relates V M to an integral of the electron-water radial distribution function, a key measure of the hydrated electron structure. Here we have applied it to three widely-used pseudopotentials and the results show that V M is a sensitive measure of the fidelity of hydrated electron descriptions. Thus, the measured V M places constraints on the hydrated electron structure that are important in developing and evaluating model descriptions. Importantly, we find that V M does not reflect only the cavity size (and thus should not be used to infer the cavity radius), but is strongly dependent on the extended solvation structure.

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

Structural Properties of [N1888][TFSI] Ionic Liquid: A Small Angle Neutron Scattering and Polarizable Molecular Dynamics Study

In this study, we investigate the quaternary ammonium-based ionic liquid (QAIL), methyltrioctylammonium bis(trifluoromethylsulfonyl)imide, [N 1888 ][TFSI], utilizing small angle neutron scattering (SANS) measurements and polarizable molecular dynamics (MD) simulations to characterize the shortand long-range liquid structure. Scattering structure factors show signatures of three length scales in reciprocal space indicative of alternating polarity (k ~ 0.44 Å –1 ), charge (k ~ 0.75 Å –1 ), and neighboring or adjacent (k ~ 1.46 Å –1 ) domains. Excellent agreement between simulation and experimental scattering structure factors validates various simulation analyses that provide detailed atomistic characterization of the different length scale correlations. The first solvation shell structure is illustrated by obtaining radial, angular, dihedral, and combined distribution functions, where two dominant spatial motifs, N + ···N – and N + ···O – , compete for optimal packing around the polar head of the [N 1888 ] + cation. Intermediate and long-range structures are governed by the balance between local electroneutrality and octyl chain networking, respectively. By computing the charge-correlation structure factor, S ZZ , and the spatial extent of the octyl chain network using graph theory, the bulk-phase structure of [N 1888 ][TFSI] is characterized in terms of electrostatic screening and apolar domain formation length scales.

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

Atomistic Insights into Lithium–Glyme Solvate Ionic Liquids: Effects of Chain Length and Anion Coordination

For this study, mixtures of lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) in diglyme (G2), triglyme (G3) and tetraglyme (G4) at solvate ionic liquid (SIL) concentrations were investigated using classical molecular dynamics (cMD) simulations with a physically motivated force-field specifically developed for modeling these systems. The structural and dynamical properties of the mixtures were computed and analyzed. Lithium solvation shells, radial distribution functions, and X-ray structure factors were studied across the different SIL systems. Translational diffusion and rotational relaxation times were also evaluated, exhibiting similar trends with increasing glyme chain length. The results are consistent with experimental data and in good agreement with previous computational studies on G3 and G4. These findings validate the accuracy of the force field in modeling glyme systems and its use for describing the [Li(G2) 4/3 ][TFSI] mixture. Additionally, the thermal and electrochemical stability of these electrolytes were systematically examined. The thermal stability appears to be governed by cooperative interactions among glyme molecules, while the electrochemical stability is primarily influenced by Li + -anion interactions, which vary significantly with glyme chain length. Overall, the study sheds light on the crucial role of the anion in these glyme-based SILs and offers valuable insights into Li + -glyme systems at SIL concentrations, highlighting their promise as potential Li-ion battery electrolytes.

anions↗
DOE OSTI2021

Single-Molecule Kinetics of Styrene Hydrogenation on Silica-Supported Vanadium: The Role of Disorder for Single-Atom Catalysts

We report a theoretical approach for the study of supported atom catalysis is developed based on recent advances in the study of single-molecule kinetics. This view is particularly useful in exhibiting the role of disorder in single-atom and single-site catalysts on amorphous supports. The distribution of passage times (or waiting times) through a complex catalytic network originating from a set of coupled active sites is described by a probability distribution function (PDF), f(t), that reflects the local environment of the reaction center. An efficient algorithm is developed based on the linear algebra of the Markov transition matrix that produces f(t) or its moments. The kinetics of the hydrogenation reaction of styrene on an organovanadium(III) catalyst supported on amorphous silica is studied. A kinetic model consisting of three intertwined catalytic cycles emanating from three chemically distinct active sites is proposed to describe the chemistry. Density functional theory (DFT) calculations are employed to determine the free energy barriers of the reactions, which are used to construct the rate coefficient matrix. The disorder induced by the amorphous support material is divided into a low-dimensional short-range component reflecting the covalent structures near the reaction center and a weaker long-range component modeling the bulk randomness. The results are computed and analyzed for a wide range of concentration values and disorder scenarios. The unusual structure in the f(t) PDF is found to occur for certain cases that reveal the contribution of multiple catalytic pathways acting in concert.

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

Temperature Dependent Local Atomic Structure and Vibrational Dynamics of Barium Hydride and Calcium Hydride

Solid-state ionic conductors that exhibit pure ionic transport of hydride anions are rare. Here, we investigate two alkaline earth metal hydrides, barium hydride and calcium hydride, using neutron scattering techniques to understand how the local atomic environment plays a role in the diffusion of hydride ions. At high temperatures, barium hydride exhibits exceptional transport properties with ionic conductivities that are higher than those of many of the typical proton and oxide ion conductors in use today. Total neutron scattering and pair distribution function analysis reveal how a structural phase transition converts barium hydride from a modest ionic conductor into a fast ionic conductor through the introduction of disorder, deuterium site splitting, and dynamic structural fluctuations. Furthermore, neutron vibrational spectroscopy is employed to probe changes in the temperature evolution of the lattice dynamics and local energy landscape. These results improve our fundamental knowledge of the interplay between structure and dynamics governing a rare conduction process of hydride ions in solid-state materials.

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

Additive Molar Volumes in Amorphous Ca/Sr Carbonate Solid Solutions

The development of predictive models of minor element incorporation in crystalline carbonate end products requires an understanding of the fundamental controls on metastable intermediate phase composition. In this study, we used small-angle X-ray scattering (SAXS), X-ray pair distribution function (PDF), thermogravimetric analysis (TGA), inductively coupled plasma mass spectrometry (ICP-MS), and transmission electron microscopy (TEM) to determine the composition-dependent density of an amorphous calcium–strontium carbonate (ACSC) solid solution. The amorphous calcium carbonate (ACC) and strontium carbonate (ASC) measured densities were ρ ACC = 2.19 ± 0.04 g/cm 3 and ρ ASC = 2.97 ± 0.05 g/cm 3 . Throughout our work, the experiments showed a dependence of the water content of the amorphous solid solution on the Sr mole fraction. An equation that relates the molar volume to the average cation radius, the carbonate ion radius, and the water volume was parameterized for hydrated crystalline carbonates and predicted well the molar volume of the ACSC solid solution. This finding indicates that, as for hydrated crystalline carbonates, the molar volumes of amorphous carbonates are additive and that water is a structural component of ACSC. Ab initio molecular dynamics (AIMD) simulations of the ACSC solid solutions showed strong linear correlations between calculated molar volumes and Sr and H 2 O contents, thus supporting the experimental results. In conclusion, our findings highlight the need to consider the full CaCO 3 –MeCO 3 –H 2 O ternary when quantifying metal cation incorporation in ACC.

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

Molecular Origins of Temperature-Dependent X-ray Absorption in YAG:Dy

The absorption and emission of X-rays in dysprosium-doped yttrium aluminum garnet (YAG:Dy) has produced unexpected thermographic behavior, which is investigated using a combination of finite temperature ab initio molecular dynamic simulations, structural characterization, and electronic structure calculations of X-ray characteristics. Calculated average peak X-ray absorption spectra (XAS) from simulations between 300 and 600 K result in peak intensity loss due to thermalization effects, matching experimentally measured behavior of YAG:Dy. Here, investigation of atomic snapshots indicates structural factors that correlated with the X-ray behavior, with the first Y–O coordination sphere identified as the primary structural feature unique to high XAS intensity as calculated by radial and pair distribution functions.

Ab initio molecular dynamics↗
DOE OSTI2025

Ion Clusters in Multicomponent Solutions Determined from X-ray PDF and SAXS Analysis: The NaNO 2 –NaOH–H 2 O System

Here, this study explores ion cluster formation in the NaNO 2 –NaOH–H 2 O system to understand how ion cluster formation is influenced by the composition in multicomponent mixtures. X-ray pair distribution function (PDF) and small-angle X-ray scattering (SAXS) identified complex ion clusters in concentrated NaNO 2 and NaOH solutions as well as their mixtures. PDF analysis showed that the Na–O distance depended primarily on total Na + concentrations rather than anion composition, whereas the nitrite-water oxygen distance stayed the same regardless of concentration or composition. This result indicates that the nitrite ion hydration was relatively independent of composition. SAXS confirms local fluctuations and coherent clusters with notable size differences between NaOH and NaNO 2 solutions. SAXS analysis of mixed solutions shows that their clusters are an average of the individual solutions, indicating mixed clusters.

Reynolds, Jacob G. [Central Plateau Cleanup Compan↗
DOE OSTI2026

Resolving Competing Ultraviolet Photoinduced Ring-Opening and Dissociation Pathways in Iodothiophenes by Ultrafast Electron Diffraction

The ultraviolet photochemistry of halothiophenes is thought to proceed via competing pathways, prompt C–X (X = halogen) bond fission, and parent ring opening followed by fragmentation, yet direct structural evidence distinguishing these mechanisms is limited. Here, time-resolved gas-phase ultrafast electron diffraction (UED) has been used to directly follow the structural dynamics of 2- and 3-iodothiophene following photoexcitation at 247 nm. For both isomers, the time-resolved pair distribution functions reveal the appearance of new interatomic distances that cannot be explained by simple C–I bond cleavage and formation of C 4 SH 3 radical cofragments. The signatures, particularly those for 2-iodothiophene, are consistent with a rival ring-opening channel of the intact parent, which evolve on time scales distinct from those associated with the direct C–I bond cleavage. In addition, the extracted kinetics indicate that the ring-opening channel occurs faster in 2-iodothiophene than in the 3-isomer at this excitation wavelength. Furthermore, these measurements provide direct geometric evidence for competing pathways leading to different structures following photodissociation in a model heteroaromatic and further establish UED as a sensitive probe of complex excited-state chemistry.

Bond cleavage↗