SEARCH · Engineering Papers
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.
Neutron diffraction, thermogravimetry, and high temperature solution calorimetry measurements reveal insights into proton and oxygen-ion kinetics in BaCo 0.4 Fe 0.4 Zr 0.2− X Y X O 3− δ triple ionic–electronic conducting ceramics.
Nickel sulfides (NiS x ) show promise for a range of energy and other applications, but their (in)stability under processing and operating conditions is scarcely studied. Herein, we have developed a new NiS x atomic layer deposition process using an easily synthesized NiCl 2 (TMPDA) precursor (TMPDA = N,N,N′,N′-tetramethyl-1,3-propanediamine) with H 2 S. Thin films deposited at 165–225 °C consist mostly of the β-NiS phase and display low resistivity (∼40–120 μΩ cm), high purity (<3 at% impurities), and a rough morphology. The thermal stability of the NiS x thin films is studied using high-temperature X-ray diffraction, revealing that structural and compositional changes occur in reducing, inert, and oxidizing atmospheres at approximately 300–400 °C. Under electrochemical water splitting conditions, the films are unstable in acid due to dissolution, especially at oxidizing potentials. In an alkaline electrolyte, we do not observe Ni dissolution, but β-NiS transforms to Ni 3 S 2 under HER conditions, possibly supplemented with Ni and/or Ni(OH) 2 species. Under alkaline OER, all sulfur is lost and NiOOH is formed. In addition to offering an attractive, scalable route to the synthesis of NiS x thin films, our work highlights the importance of thermal and electrochemical (in)stability of sulfides as a crucial step for understanding and engineering materials for energy and other applications.
BaTaO 2 N stands out among oxynitride photocatalysts because of its ability to capture visible light and to drive the photoelectrochemical water oxidation reaction. However, its solar energy conversion performance is limited by electron–hole recombination at Ta( IV ) defects in the material. These defects are formed by overreduction of the Ta(v) oxide precursor by excess ammonia under the high temperature conditions during ammonolysis. Here we show for the first time that Ta( IV ) defect concentrations can be lowered by conducting the ammonolysis reaction in mixed NH 3 /N 2 gas. The obtained BaTaO 2 N samples crystallize in the cubic CaTiO 3 structure type and form 200–300 nm faceted nanocrystals, based on X-ray diffraction, scanning electron microscopy, and HRTEM. Electron paramagnetic resonance spectra observe the Ta( IV ) defects at g = 1.999 and confirm an 11-fold reduction for the product synthesized in mixed (0.13 : 1.0 vol) NH 3 /N 2 gas, equivalent to 1.14 × 10 16 cm −3 Ta( IV ) ions. This optimized BaTaO 2 N has nearly twice the photocatalytic oxygen evolution activity (AQE of 6.78% at 400 nm) of a reference material made with 1.0 atm ammonia and 78% higher photoelectrochemical water oxidation photocurrent (0.9 mA cm −2 at 1.23 V vs. RHE) under simulated sunlight. According to X-ray photoelectron spectroscopy, remaining Ta( IV ) defects are concentrated in the surface region of the BaTaO 2 N particles, where >50% of all Ta ions are found in the +4 oxidation state. This surface Ta( IV ) population can be directly observed in Vibrating Kelvin Probe Surface Photovoltage Spectra (VK-SPV) via its 1.2–1.4 eV photovoltage onset. Here, it suggests that the surface Ta( IV ) ions contribute empty d-states 0.5–0.7 eV below the BaTaO 2 N conduction band edge. These findings highlight how the energetics and concentrations of Ta( IV ) defects influence the photoelectrochemical water oxidation ability of BaTaO 2 N. Additionally, the work establishes ammonolysis with diluted NH 3 as a new tool to minimize defects in BaTaO 2 N and to raise its solar energy conversion efficiency toward its theoretical limit. Because of its simplicity, the reduced ammonia pressure strategy will likely be applicable to other oxynitrides, which generally suffer from overreduction problems during ammonolysis.
Energy production is a fundamental pillar of modern society, but it remains heavily reliant on fossil fuel combustion, which leads to environmental degradation and resource depletion. Electrocatalytic reduction reactions (ERRs) offer a cleaner alternative by functioning under mild conditions, eliminating the need for high temperatures or pressures. However, their widespread adoption is hindered by challenges such as low efficiency, inadequate selectivity, and sluggish kinetics. Here, the introduction of chiral-induced spin selectivity (CISS) offers a promising strategy to address these limitations by reducing energy barriers in reactions involving spin-sensitive intermediates or transition states. By leveraging CISS, ERR pathways can be manipulated to favor desired reactions, enhance product yield, and suppress unwanted byproducts, improving selectivity. Beyond reducing the impact of competing side reactions, CISS also enhances the kinetics of ERRs by accelerating reaction rates and lowering the energy barriers associated with multi-electron transfer steps. Thanks to its straightforward and versatile application, CISS shows significant potential in optimizing well-studied ERRs, such as the hydrogen evolution reaction (HER) and carbon dioxide reduction reaction (CO2RR), facilitating the future transition of these technologies from research laboratories to commercial energy plants worldwide.
OpenStar Technologies is pursuing the levitated dipole (LD) as a highly modular, loosely-coupled system that leverages their expertise in high temperature superconductor (HTS) technology. The next generation experiment at OpenStar, Tahi (Ma¯ori for “first”), will demonstrate the generation and confinement of fast ions in a levitated dipole for the first time. Ion cyclotron range of frequency (ICRF) heating is a leading candidate for energetic ion formation in Tahi. A frequency in the 10 MHz range will be used for H minority heating or D majority heating with waves launched from an antenna located above the floating coil. Unlike a tokamak, where the targeted cyclotron resonance is typically a vertical path through the center of the plasma, in a dipole the resonance location follows a C-shaped path from the separatrix to the center of the plasma. The value of B also varies significantly within the confined plasma resulting in a large number of cyclotron harmonics present in the low field region. Furthermore, levitated dipoles contain a “first closed flux surface” surrounding the floating coil, in addition to the traditional separatrix /last closed flux surface. Simulation e ff orts using full-wave ICRF codes show that ICRF heating of a levitated dipole reactor is feasible using a pair of toroidal current straps phased to launch the appropriate parallel ( i.e. poloidal) refractive index.
Here, solutions and minerals containing sulfate (SO 4 2- ), and Ca 2+ and Mg 2+ cations, are ubiquitous throughout the lithosphere and are significant components of seawater, thus presenting a prototypical system for the study of strong electrolytes and crystal nucleation mechanisms. However, despite their relative abundance, key questions remain unanswered about the most fundamental atomic-level steps of their mineralization pathways and aqueous dynamics. Here, we carry out enhanced sampling multi-level molecular dynamics (MD) embedded correlated wavefunction theory simulations to elucidate ion-pairing mechanisms for Mg–SO 4 and Ca–SO 4 in concentrated aqueous solution, accurately capturing effects arising from both structural dynamics and electron exchange–correlation. We predict contact-ion-pair formation to be barrierless and highly exoergic for Ca–SO 4 , in agreement with its minimal solubility, whereas for Mg–SO 4 , solvent-shared and contact ion pairs have similar free energies, qualitatively consistent with its higher solubility. Finally, we demonstrate that brief high-temperature pre-equilibration may be utilized to accelerate convergence of free energies in blue-moon-ensemble enhanced-sampling MD.
Knowing the thermophysical properties of high-temperature melts can aid the design of melt processes and validate atomic structural models, such as those used in studying glass formation. Property measurements on such melts are challenging, however, due to container-related contamination and heterogeneous nucleation. Containerless processing techniques that employ levitation can be used to avoid these obstacles. In that context, we demonstrate here the application of silhouette imaging to measure the density of molten oxides in an aero-acoustic levitation instrument (AAL). The AAL combines gas jet levitation with actively controlled acoustic positioning to enable full optical access to samples ca. 2–4 mm in diameter, which are laser beam heated and melted. The cross sections of molten drops are imaged using a monochromatic light source and narrowband-filtered camera. Melt volume is calculated from the cross sections and used to find density at several temperatures ranging 1530–1920 K, including up to 350 K of supercooling. We report densities for CaAl 2 O 4 , Ca 12 Al 14 O 33 , CaSiO 3 , their Fe 2 O 3 -bearing analogs, and 83TiO 2 -17RE 2 O 3 (RE = La or Nd). These provide important benchmarks of the capabilities, measurement uncertainties, and future outlook for this technique.
We investigate electron transport in the uniform electron gas using ring-polymer molecular dynamics (RPMD). Working in the weakly coupled, non-degenerate regime, we use RPMD to probe how the onset of quantum diffraction effects at high temperature reshapes electron–electron collisions and leads to a classical-to-quantum crossover in macroscopic transport properties. Static thermodynamics obtained with RPMD are consistent with the weak-coupling equation of state, confirming correct quantum Boltzmann sampling. Real-time transport extracted from mean square displacements exhibits the expected ballistic-to-diffusive transition and a systematic reduction of the electronic self-diffusivity as quantum effects strengthen, due to quantum diffraction modifying electron–electron collisions. Direct ring-polymer scattering simulations reveal diffractive “softening” of binary deflections, providing a micro-to-macro link between collision physics and diffusion. The present study establishes RPMD as a quantitative, trajectory-based tool for electron transport across the classical–quantum crossover and furnishes benchmarks for improving Coulomb-log interpolation models. We outline extensions to multi-component plasmas and a path to incorporate Fermi–Dirac statistics within path-integral dynamics.
In an inertially confined fusion experiment, a capsule of fusion fuel is compressed to high temperatures and densities using lasers, creating a hot plasma in which thermonuclear fusion reactions take place. To achieve efficient conversion of laser energy to fuel heating, the implosion must be as spherically symmetric as possible. One common source of asymmetry that degrades performance is the stalk and glue spot that hold the capsule in place before the experiment. The exact nature and magnitude of this perturbation had not been previously quantified. In this work, inertially confined fusion experiments were executed with an additional stalk on the capsule to measure and quantify the effect of the stalk and glue spot on the symmetry and performance. It is found that the stalk drives low-mode asymmetry by shadowing the laser drive on the part of the shell near the glue spot, enhancing the hot-spot velocity toward the stalk by 40–80 km/s and elongating the hot-spot along the stalk axis. It is also found that localized mix due to the additional stalk has a negligible effect on implosion performance. Two-dimensional radiation-hydrodynamic simulations using the xRAGE code reproduce the direction of the experimentally observed change to hot-spot velocity and implosion performance metrics, but not the observed elongation, suggesting that this feature results from an interaction between stalk perturbations and other preexisting asymmetry seeds, such as laser drive asymmetry.
Sm2Fe17C3 phase (2:17) is metastable and exhibits excellent intrinsic hard magnetic properties. Doping elements such as Ga facilitate the formation of a single-phase 2:17 structure in arc-melted Sm2Fe17Cx alloys, which opens a promising route for fabricating fully dense bulk Sm2Fe17Cx magnets via high-temperature techniques such as melting and sintering. First-principles electronic structure calculation indicates that Ga prefers to partially replace Fe at the 9d and 18h crystallographic sites in Sm2Fe17C3 and Sm2Fe17, respectively. This difference in site preference is attributed to the distinct chemical environments surrounding the Fe atoms in the two compounds. Ga substitution favors the Sm–Ga bonding formation while avoiding Ga–C interactions. Doped Ga atoms result in more negative formation energy in Sm2(Fe, Ga)17C3, indicating improved structural stability. Crystal Orbital Hamilton Population analysis reveals that carbon insertion weakens the bonding of Sm-Fe (18h) and Sm-Fe (18f) in Sm2Fe17C3. Ga doping facilitates electron redistribution across chemical bonds, thereby reinforcing Fe(18h)–Sm and Fe(18f)–Sm interactions and stabilizing the carbon-centered octahedral local structure. This synergistic effect contributes significantly to the observed enhancement in phase stability of Sm2(Fe, Ga)17Cx. These findings suggest that chemical bond engineering through the selective doping of Ga can enhance phase stability and facilitate the synthesis of Sm2Fe17C3, providing a viable strategy for developing advanced magnets.
The pseudogap phenomena have been a long-standing mystery of the cuprate high-temperature superconductors. The pseudogap in the electron-doped cuprates has been attributed to band folding due to antiferromagnetic (AFM) long-range order or short-range correlation. We performed an angle-resolved photoemission spectroscopy study of the electron-doped cuprates Pr 1.3−x La 0.7 Ce x CuO 4 showing spin-glass, disordered AFM behaviors, and superconductivity at low temperatures and, by measurements with fine momentum cuts, found that the gap opens on the unfolded Fermi surface rather than the AFM Brillouin zone boundary. The gap did not show a node, following the full symmetry of the Brillouin zone, and its magnitude decreased from the zone-diagonal to (π,0) directions, opposite to the hole-doped case. These observations were reproduced by cluster dynamical-mean-field-theory calculation, which took into account electron correlation precisely within a (CuO 2 ) 4 cluster. In conclusion, the present experimental and theoretical results are consistent with the mechanism that electron or hole doping into a Mott insulator creates an in-gap band that is separated from the upper or lower Hubbard band by the pseudogap.
The quantum acoustic framework has recently emerged as a nonperturbative, coherent approach to electron–lattice interactions, uncovering rich physics often obscured by perturbative methods with incoherent scattering events. Here, we model the strongly coupled dynamics of electrons and acoustic lattice vibrations within this framework, representing lattice vibrations as coherent states and electrons as quantum wave packets, in a manner distinctively different from tight-binding or discrete hopping-based approaches. We derive and numerically implement electron backaction on the lattice, providing both visual and quantitative insights into electron wave packet evolution and the formation of acoustic polarons. We investigate polaron binding energies across varying material parameters and compute key observables—including mean square displacement, kinetic energy, potential energy, and vibrational energy—over time. Our findings reveal the conditions that favor polaron formation, which is enhanced by low temperatures, high deformation potential constants, slow sound velocities, and high effective masses. Additionally, we explore the impact of external electric and magnetic fields, showing that while polaron formation remains robust under moderate fields, it is weakly suppressed at higher field strengths. These results deepen our understanding of polaron dynamics and pave the way for future studies into nontrivial transport behavior in quantum materials.
Large Language Models (LLMs) show great promise as a powerful tool for scientific literature exploration. However, their effectiveness in providing scientifically accurate and comprehensive answers to complex questions within specialized domains remains an active area of research. Using the field of high-temperature cuprates as an exemplar, we evaluate the ability of LLM systems to understand the literature at the level of an expert. We construct an expert-curated database of 1,726 scientific papers that covers the history of the field, and a set of 67 expert-formulated questions that probe deep understanding of the literature. We then evaluate six different LLM-based systems for answering these questions, including both commercially available closed models and a custom retrieval-augmented generation (RAG) system capable of retrieving images alongside text. Experts then evaluate the answers of these systems against a rubric that assesses balanced perspectives, factual comprehensiveness, succinctness, and evidentiary support. Among the six systems, two using RAG on curated literature outperformed existing closed models across key metrics, particularly in providing comprehensive and well-supported answers. We discuss promising aspects of LLM performances as well as critical short-comings of all the models. The set of expert-formulated questions and the rubric will be valuable for assessing expert level performance of LLM based reasoning systems.
The accelerated fuel qualification (AFQ) framework has been used for the initial development of multiscale modeling of silicon carbide (SiC) fiber reinforced composite (SiC-SiC). The AFQ framework provides a methodology to leverage physics-informed multiscale modeling along with a reduced set of empirical test data to reduce the time and cost of licensing and qualification of new nuclear fuel systems while maintaining the overall nuclear power plant safety case. SiC-SiC is being proposed for in-core applications, most notably fuel cladding, for current and next-generation nuclear reactors because of its high temperature stability, irradiation tolerance, and ability to withstand many accident conditions. As these composites exhibit multiscale architectures and complex microstructure-based fracture mechanics, it is an appealing use case for the AFQ methodology. While the end goal of this work is a single multiscale model that can be used for predictive in-core performance, current focus is on the individual various length scale models. Four individual models have been initially developed from microscale to engineering system level to capture key physics-based effects across different length scales. These models include a microscale homogenized tow model, a mesoscale fast Fourier transform–based weave model that integrates the homogenized tow model, a mesoscale finite element–based weave model, and a system-level BISON fuel performance model. Results of these models have undergone an initial comparison with separate-effects test data showing a good match to experimental results. By using the AFQ framework during model development, several near-term benefits have been secured including a reduction in development time for the SiC-SiC cladding, more targeted irradiation testing, and a better understanding of uncertainty.
Advanced high-temperature fluid reactors, such as sodium-cooled fast reactors (SFRs) and molten salt–cooled reactors (MSCRs), require coolant purification systems to prevent fluid contamination and local freezing that can lead to plugging. Liquid sodium purification can be achieved with a cold trap, where the sodium temperature is reduced to a near-freezing point to precipitate out impurities. Automation of monitoring of the cold trap performance with machine learning algorithms can aid in early detection of incipient anomalies. An efficient approach to loss-of-coolant–type anomaly detection in a cold trap monitored with more than two dozen thermal-hydraulic sensors consists of a long short-term memory (LSTM) autoencoder. This work develops the uncertainty quantification of the LSTM autoencoder performance for cold trap anomaly detection using the Monte Carlo (MC) dropout method. The MC dropout methodology creates a distribution of sister distributions that all slightly differ from each other because of random neurons being turned off for testing. The variances of the sister network distributions are used to make an uncertainty interval. Our analysis shows that the uncertainty in the autoencoder performance is largest near the peak of the anomaly signal. Using the MC dropout method, we investigate the uncertainty in the anomaly detection with missing sensor inputs. This capability allows the reactor operator to evaluate resilience of the anomaly detection system and to make informed decisions about continuity of operation in the event of sensor failure.
Manifold-based models offer a computationally efficient alternative to directly transporting the thermochemical state in computational simulations of turbulent reacting flows, projecting the high-dimensional thermochemical state-space onto a low-dimensional manifold. Recent efforts have yielded a manifold-based model applicable to multi-modal combustion, enabling reconstruction of the thermochemical state from solutions to two-dimensional manifold equations in mixture fraction and generalized progress variable that are parameterised by three scalar dissipation rates. In coarse-grained simulations such as Large Eddy Simulation (LES), closure of the multi-modal manifold equations and subfilter variances/covariance requires closure of three filtered scalar dissipation rates. Here, the present work adopts a data-based approach, providing closure for the three filtered scalar dissipation rates via deep neural networks (DNNs). High-fidelity datasets corresponding to an autoigniting n-dodecane jet flame and a bluff body swirl-stabilized confined lifted spray flame of two aviation fuels (Jet-A and C1) with different ignition propensities are leveraged to generate training data that spans a diverse range of thermodynamic conditions and combustion modes, including low- and high-temperature ignition regimes in addition to premixed and nonpremixed behaviour. A final DNN model is trained to enforce inherent physical constraints by learning nonlinear functional transformations of the three filtered scalar dissipation rates. The generalizability of this constrained DNN model is demonstrated a priori via conditional statistics evaluated on the lifted spray flame with C1–a configuration that had not been included in the training data. Excellent DNN agreement with conditional DNS statistics is observed, and integrated gradients are computed to identify the most sensitive input variables. The similarity of the marginal PDFs of the most informative input variables and outputs across configurations are quantified via the Wasserstein metric, demonstrating that data-based models may successfully generalize to unseen parametric conditions so long as the most informative input variables share similar distributions across training and testing datasets.
The report summarizes DOE INFUSE-supported work between HB11 Energy and the University of Rochester’s TriForce Institute to improve computational modeling of advanced fusion fuels, especially proton–boron-11 (p- 11 B). The project extended the TriForce particle-in-cell/Monte Carlo collision code to include physics needed for dense, high-temperature p- 11 B burn studies, including p- 11 B fusion reactions, three-alpha-particle reaction products, relativistic Coulomb collisions, large-angle nuclear scattering, bremsstrahlung radiation, inverse bremsstrahlung absorption, and photon transport. The upgraded models were verified against focused physics tests and against known deuterium–tritium burn behavior. The study then used one-dimensional spherical simulations to estimate the conditions required for thermonuclear burn propagation in isochoric p- 11 B fuel. The calculations found that burn propagation is possible in the model, but only under very extreme hot-spot conditions, such as about 7000 g/cm 3 at 500 keV or 9000 g/cm 3 at 300 keV for a 20-micron hot spot. These conditions are much more demanding than current demonstrated inertial confinement fusion hot spots. The report concludes that the INFUSE collaboration successfully advanced kinetic and radiation modeling capabilities for p- 11 B fusion and provided useful estimates of ignition requirements. However, the simulated fuel gains remain below what would be needed for practical inertial fusion energy, and further work is needed to reconcile differences among kinetic, radiation-hydrodynamic, and analytic models and to identify more achievable target designs.
High energy x-ray diffraction (HEXRD) has emerged as a pivotal technique for examining the structural properties of molten salts, which play vital roles in various high-temperature applications, such as energy storage, nuclear reactors, and metallurgy. This review provides a comprehensive analysis of recent progress in HEXRD studies of various molten salt systems, including molten halide salts, their mixtures, and salts with molecular anions, focusing on the insights into their atomic-scale arrangements. Additionally, we discuss methodological advancements that have improved the resolution and accuracy of HEXRD measurements, as well as recent developments in computational simulations, such as machine learning interatomic potentials, for interpreting results. This review aims to serve as a valuable resource for researchers in this field, offering a detailed overview of the progress of HEXRD studies of molten salts and suggesting pathways for further exploration.