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.
Machine learning (ML) is a powerful tool for functional material design. In this work, we combine first-principles density function theory with ML to develop perovskite database and design O<sub>2</sub> sensors for harsh environmental applications.
Here, the filterscope diagnostic uses bandpass filters and photomultiplier tubes (PMTs) to detect specific spectral emission lines. A filterscope was used to measure the W I 400.88 nm line emission as a function of ion energy on the Radio Frequency Plasma Interaction Experiment (RF PIE) for the purpose of assessing W erosion on plasma-facing components (PFCs). Different filter techniques are being explored and compared in order to effectively screen out nearby impurity lines, like the Ar II 401.39 nm line. The effectiveness of these techniques is determined by comparing the measurements to a high-resolution 1.0 m Czerny–Turner spectrometer with 0.012 nm spectral resolution. The ability to filter out nearby impurity emissions is useful when imaging PFCs in fusion devices including divertor and antenna guard limiters. Initial results with a helium plasma show little difference between the two techniques at higher bias voltages. In a helium plasma at lower dc bias voltages, a two-filter technique with filters at two separate wavelengths was shown to be more effective at screening out background signals. Data collected with an argon plasma however show the technique with overlapping filters on the line of interest is a closer match to spectrometer data at lower dc bias voltages.
High-Q superconducting resonators have been proposed and developed as detectors of light-by-light scattering mediated by the hypothesized axion or virtual electron–positron pairs in quantum electrodynamics: the Euler–Heisenberg (EH) interaction. Photon frequency and mode conversion is central to the scheme for detecting such rare events. Superconducting resonators are nonlinear devices. The Meissner screening currents that confine the electromagnetic fields to the vacuum region of a superconducting RF cavity are nonlinear functions of the EM field at the vacuum–superconducting interface, and as a result can generate source currents and frequency conversion of microwave photons in the cavity. In this report we consider photon frequency and mode conversion in superconducting resonators with high quality factors from Meissner currents in single- and dual-cavity setups proposed for axion and QED searches based on light-by-light scattering. In a single cavity with two pump modes, photon frequency conversion by the Meissner screening current dominates photon generation by the EH interaction for cavities with Q≲1012. The Meissner currents also generate background photons that limit the operation of the resonator for axion detection in three-mode, single-cavity setups. We also consider the leakage of photons from pump modes into the signal mode for both axion- and EH-mediated light-by-light scattering. Photon frequency conversion by the EH interaction can compete with Meissner and leakage radiation in ultra-high-Q cavities that are beyond the current state of the art. Meissner radiation and leakage backgrounds can be suppressed in dual-cavity setups with appropriate choices for pump and spectator modes, as well as the single-cavity setup proposed for heterodyne detection of galactic axion dark matter.
The superheavy nuclei push the periodic table of the elements and the chart of the nuclides to their limits, providing a unique laboratory for studies of the electron-nucleus interactions. The most important weak decay mode in known superheavy nuclei is electron capture (EC). In the standard calculations of EC, the lepton wave functions are usually considered in the lowest-order approximation. In this work, we investigate the sensitivity of EC rates on the choice of the electron wave functions by (i) assuming the single-particle approximation for the electron wave functions, and (ii) carrying out Dirac-Hartree-Fock (DHF) calculations. The nuclear response is generated based on the state-of-the-art quasiparticle random phase approximation employing relativistic nuclear energy density functional theory. Here, we show that using the improved lepton wave functions reduces the EC rates up to 40% in the superheavy nucleus oganesson (𝑍=118). Interestingly, because of screening effects, the difference between the EC rates obtained with the DHF and single-particle calculations is fairly small.
High-throughput density functional theory (DFT) calculations have become a vital element of computational materials science, enabling materials screening, property database generation, and training of “universal” machine learning models. While several software frameworks have emerged to support these computational efforts, new developments such as machine learned force fields have increased demands for more flexible and programmable workflow solutions. This manuscript introduces atomate2, a comprehensive evolution of our original atomate framework, designed to address existing limitations in computational materials research infrastructure. Key features include the support for multiple electronic structure packages and interoperability between them, along with generalizable workflows that can be written in an abstract form irrespective of the DFT package or machine learning force field used within them. Our hope is that atomate2's improved usability and extensibility can reduce technical barriers for high-throughput research workflows and facilitate the rapid adoption of emerging methods in computational material science.
Continuous counter‐current microfluidic liquid–liquid extraction performs separations by flowing immiscible liquids in opposing directions within a single flow channel. In principle, this flow arrangement enables a large number of theoretical separation units in a small footprint, without using interstage valving, pumping, and phase separation. Despite its potential for excellent separation performance, this microfluidic scheme rarely appears in literature due to the requirement for capillary forces to be greater than hydrodynamic forces for stable flow. We present a novel microfluidic device and flow approaches that overcome this force‐balance challenge, enabling stable, long‐duration continuous counter‐current flow. Additionally, we cover a suite of methodologies for quantifying the performance of the microfluidic device, revealing the number of theoretical equilibrium stages achieved. The enabling technologies include a woven mesh screen‐based microfluidic device architecture that is easily fabricated outside of a clean room, surface functionalization strategies to promote conjugate (organic/aqueous) wettability, flow approaches to eliminate bubbles and carryover, and computer‐aided flow automation with optical measurement of extraction performance. The reported experiments lasted for over 36 h, terminated only at experiment conclusion, where the device still exhibited good performance. Automated Raman spectroscopy was used for solute quantitation of the ternary system tert‐butanol in a toluene/water matrix, a ternary system that was specifically chosen to analyze the device's performance with a small solute partition ratio and to enable in‐line Raman measurements of solute concentrations in both phases. The microfluidic device possessed a 55 mm contact length and a 38.5 µL internal volume. During counter‐current flow, we observed approximately 37 equilibrium stages (37 ± 13) based on a best‐fit of the solute fraction remaining in the aqueous phase using a Kremser Group Method analysis.
Here, density functional theory (DFT) is routinely used to make electronic structure predictions for high-throughput screening of materials and molecules for technologically relevant areas, like the identification of better catalysts, electronic materials, and drug discovery. However, the DFT formalism is limited by (a) its poor (quadratic-to-quartic) scaling, and (b) the need to perform repeated eigenvalue computations of the electronic Hamiltonian as part of its self-consistent field (SCF) iteration procedure to obtain the converged ground state electron density, ρ (r). Approaches that directly predict ρ (r) of a structure with high accuracy can accelerate conventional SCF calculations and can also be used in linearly scaling methods such as orbital-free DFT. To this end, we present a procedure to predict the ground state electron density of molecular and periodic three-dimensional systems directly from the atomic structure with a particular emphasis on physical interpretability. In our framework, ρ (r) is modeled using many-body correlation descriptors that accurately capture the effects of local atomic arrangements in the neighborhood of a grid point. Our use of a linear regression scheme to fit to charge density data enables transparent analysis of the relative contributions of various types of local atomic correlations. By systematically including increasingly complex correlations, our model is shown to accurately predict ρ (r) for a variety of chemically and electronically diverse systems — amorphous Ge, Al(001) slab, crystalline Ga 2 O 3 , molecular benzene, and polyethylene. We then demonstrate a symbolic regression-based protocol to construct easily computable, interpretable features from lower-order correlations that significantly improves our electron density predictions with effectively no increase in the computational cost.
This study investigates the effects of ionic cores on non-isothermal plasmas using a novel ion-ion interaction potential that incorporates screening effects from both ion cores and exchange-correlation interactions. Our findings indicate that with increasing distance, the effective potential approaches a Yukawa-like screening potential, while at shorter distances, strong electron binding weakens screening. The different values of the cutoff radius ' cut and core edge steepness α significantly influence the potential behavior and radial distribution functions (RDFs). Higher coupling parameters (Γ i ) strengthen the electron-ion interactions, leading to deeper potential wells and more pronounced non-ideality corrections. Increasing Θ decreases the absolute values of non-ideality corrections, indicating fewer interactions in the system. A larger cutoff radius r cut due to weaker screening effects. As Γ i at a fixed parameter also reduces corrections increases, non-ideality corrections grow, reflecting stronger coupling. The results show the importance of taking into account the ion core effects in dense non-isothermal plasma research.
Understanding the temperature dependence of functional properties of sensing materials is vital for their applications in combustion environments. The electron-phonon coupling that derives the electronic structure change with temperatures is a key property of interest as it affects other sensing responses. Herein, we first assess the temperature dependence of band gap renormalization in sensing materials by employing Allen-Heine-Cardona (AHC) theory with density functional theory (DFT) simulations corroborated with experimental observation. As the AHC calculations are impractical for high-throughput screening of materials, we employ data-driven Gaussian process regression to predict the parameters employed in the O’Donnell empirical model from a set of physical features. To mitigate the reliability issues arising from the small size of the dataset, we apply a Bayesian technique to improve the generalizability of the data-driven models as well as to quantify the uncertainty associated with theoretical predictions. These models capture well the overall trend of the O’Donnell parameters with respect to a reduced feature set obtained by transforming the available physical features. Quantifying the associated uncertainty helps us understand the reliability of the predictions and, therefore, the variation of bandgap as a function of temperature for other novel materials. The predicted candidates from machine learning models are further validated by experiments and DFT calculations.
Abstract Removal of a catalyst after a catalytic process is a time‐consuming and challenging process. In this paper, a fullerene catalyst coated on cellulose thread is developed and applied it in a heterogeneous photoreaction screening via a thread spray mass spectrometry (MS) platform. The surface of the cellulose thread is modified with NH 2 groups, which subsequently allows immobilization of fullerene C 60 on the thread. The final fullerene functionalized thread is placed in a glass capillary and the solution of the reactant to be interrogated is added. The assembly of thread and glass capillary is placed in front of a mass spectrometer, which serves as an online photoreaction platform after a light source and direct current high voltage are applied. As an ambient ionization approach, the photocatalytic reaction screening platform enabled real‐time MS characterization of products and intermediates and using only nanomole quantities of reactants. Fullerene C 60 is identified to facilitate dehydro‐dimerization reactions from amines. Norharmane shows a unique photo‐oxygenated dimer formation. The functionalized threads immobilized with fullerene C 60 are reusable, confirming the heterogeneous nature of the platform. Thus, this system provided real‐time screening of reactions catalyzed by fullerene C 60 and a means to investigate the reaction mechanism.
Abstract High‐entropy perovskite oxides (HEPOs) have recently emerged as multifunctional catalysts. However, the HEPOs’ structural and compositional complexity hinders the easy and accurate extrapolation of activity indicators, which are essential for establishing structure‐property correlations. Here, OxiGraphX, is introduced as a novel graph neural network (GNN) model designed to capture the complex relationships among structure, composition, and atomic chemical environments for accurate prediction of oxygen vacancy formation energies (OVFEs) in HEPOs. By integrating machine learning (ML), density functional theory (DFT), and experimental validation, this work demonstrates an efficient framework for rapidly and accurately screening HEPO electrocatalysts for oxygen evolution reaction (OER). The OxiGraphX predicts OVFEs with a precision exceeding existing data, enabling the identification of compositions of higher oxygen vacancy content (OVC) and, thus, higher catalytic activity. Furthermore, the model explores latent spaces that translate effectively into experimental domains, bridging computational predictions with real‐world applications. This approach accelerates the discovery of high‐performance HEPO catalysts while providing deeper insights into their catalytic mechanisms.
A grand challenge in materials research is identifying the relationship between composition and performance. Herein, we explore this relationship for magnetic properties, specifically magnetic saturation (M s ) and magnetocrystalline anisotropy energy (K) of ferrites. Ferrites are materials derived from magnetite (chemical formula = Fe 3 O 4 ) that comprise metallic elements such as Fe, Mn, Ni, Co, Cu and Zn. Experimentally, synthesizing and characterizing ferrites is time consuming. Further, selection of compositions to achieve particular magnetic properties currently relies on intuition. To address this, in this work, density functional theory (DFT) is used to predict M s and K for 571 ferrite structures. These structures are primarily double-substituted non-stoichiometric ferrites with formulae M1 x M2 y Fe 3–x–y O 4 , where M1 and M2 can be Mn, Ni, Co, Cu and/or Zn and 0 ≤ x ≤ 1 and y = 1–x. Calculated magnetic properties for the structures in this study vary from 0.04 × 10 5 to 9.6 × 10 5 A m −1 for M s and from 0.02 × 10 5 to 14.08 × 10 5 J m −3 for K. All structures are made publicly available in a FAIR database.
Funded by the ARPA-E ULTIMATE program, this project aims to develop cost-effective refractory high entropy alloys (RHEAs) for gas turbine blades applications at 1300 degree Celsius and above. The alloy design focuses on carbides precipitation strengthening in a strong and ductile solid solution BCC matrix. The goal is to achieve balanced mechanical properties including room temperature ductility & fracture toughness and high temperature strength and creep resistance, while maintaining comparable oxidation resistance and densities to Ni-base superalloys. To accelerate alloy design, multiscale modeling including CALPHAD and first-principles density functional theory (DFT) calculations and machine learning are performed. Specifically, high throughput phase diagram calculations and screening are carried out using CALPHAD; DFT calculations are performed to predict intrinsic ductility, grain boundary strengthening potential, coefficient of thermal expansion (CTE) and temperature-dependent elastic constants. Temperature dependent yield stress, creep strain and CTE are also predicted using machine learning. Down-selected alloys are synthesized in small buttons of about 250g using arc melting for rapid evaluation on microstructure and mechanical properties before employing plasma arc melting and additive manufacturing for producing large ingots. Preliminary computational and experimental results of this project will be presented and future direction in RHEA development will be discussed.
Molecular-level insights into reactive separations are crucial for the design of new conversion pathways of carbon dioxide (CO 2 ). This work explores a postulated pathway that directs CO 2 to undergo inverse-electron-demand Diels–Alder reactions to produce heterocycles using the CO 2 chemically fixed on water-lean solvent molecules. Density functional theory calculations are applied to evaluate the lowest unoccupied molecular orbital (LUMO) energies of three types of reactants (1,3-butadiene, 1,3-cyclohexadiene, and 1,2,4,5-tetrazine) with various functional substituents. These calculations also provide a data set (5.8k data) for developing a machine learning model to efficiently predict LUMO energies. A computational screening of LUMO energies for an additional 47k diene and tetrazine candidates is performed, and a list of candidates with lowered LUMO energies by electron-withdrawing substituents is provided. These candidates are further examined by their reaction energy barriers computed from the interatomic potential or density functional theory. Two major energy barriers are identified, one for the proton transfer within the water-lean solvent and the other for the CO 2 transfer from the solvent molecule to the reactant candidate (diene or tetrazine). The functional substituents have a more significant impact on the second barrier but a very slight one on the first barrier. This exploratory work demonstrates a new possibility for guiding experimental efforts toward the chemical conversion of fixated CO 2 to value-added compounds.
The electrosynthesis of hydrogen peroxide (H 2 O 2 ) via a two-electron oxygen reduction reaction enables decentralized H 2 O 2 production. While metal-free carbon catalysts are sustainable and low-cost, their performance is hindered by poorly defined active sites and uncontrolled defect states. Here, we resolve these challenges through active site design and catalyst screening using fluorine (F) and nitrogen (N) codoped carbons as model materials. Statistical analysis combined with density functional theoretical calculations reveals that F-induced structural modification and defect passivation optimize OOH* binding, with F-doping and adjacent F atoms predominantly lowering abs ΔG(OOH*). Experimental results confirm that semi-ionic C–F bonds passivate defects in nitrogen-doped carbon, enhancing catalytic activity and durability. The resulting (N, F)-codoped carbon achieves nearly 100% H 2 O 2 selectivity at 0.5–0.65 V versus the reversible hydrogen electrode and maintains > 95% across 0.01–0.65 V versus the reversible hydrogen electrode. In an electrolyzer, (N, F)-codoped carbon exhibits an H 2 O 2 yield rate of 74.35 mol g cat. −1 h -1 and sustains 300 mA cm -2 for 105 hours with ~95% faradaic efficiency. Coupling the two-electron oxygen reduction reaction with methanol oxidation further reduces cell voltage and enhances productivity. This work provides a means to design efficient catalysts for industrial H 2 O 2 electrosynthesis.
Multivalent-ion batteries offer an alternative to Li-based technologies, with the potential for greater sustainability, improved safety, and higher energy density, primarily due to their rechargeable system featuring a passivating metal anode. Although a system based on the Ca 2+ /Ca couple is particularly attractive given the low electrochemical plating potential of Ca 2+ , the remaining challenge for a viable rechargeable Ca battery is to identify Ca cathodes with fast ion transport. In this work, a high-throughput computational pipeline is adapted to (1) discover novel Ca cathodes in a largely unexplored space of empty intercalation hosts and (2) develop material design rules for Ca-ion mobility. One candidate from the screening, W 2 O 3 (PO 4 ) 2 , is confirmed to have a low Nudged Elastic Band (NEB) barrier of 168 meV within a one-dimensional (1D) ion percolation topology. This candidate is subsequently synthesized and electrochemically tested, achieving reversible Ca cycling with a capacity of 25 mA h/g. To further accelerate the screening for promising Ca intercalation electrodes, machine learning (ML) Random Forest (RF) and Extreme Gradient Boosting (XGB) classification models are created with local environment descriptors based on a large, structurally and chemically diverse dataset of minimum energy pathways, spanning over 5,000 density functional theory (DFT) site energy calculations. Accuracies of 92% are achieved, material design metrics are quantified, ML force-fields are leveraged in an accelerated iteration of the screening, and a total of 27 novel Ca cathode materials are highlighted for further investigation.
This study investigates the facile hydride synthesis method guided by theoretical predictions to explore the K–T–Bi (T = Zn, Cd) phase spaces. Using an adaptive genetic algorithm (AGA) and density functional theory (DFT), candidate compositions are identified for experimental validation via a facile hydrides route, permitting experimental screening of K–Zn–Bi and “empty” K–Cd–Bi systems. The previously reported KZnBi and KZn 2 Bi 2 are synthesized alongside newly discovered KCdBi and KCd 2 Bi 2 . While the AGA and DFT predict the stability of these compounds, structural predictions align with the experiment only for KZnBi and KZn 2 Bi 2 . Single-crystal X-ray structure refinements confirm that KZnBi and KZn 2 Bi 2 adopt the hexagonal ZrBeSi- and tetragonal ThCr 2 Si 2 -structure types, respectively. KCdBi has tetragonal PbClF-structure type and KCd 2 Bi 2 belongs to the ThCr 2 Si 2 -structure type. A trend based on the ratio of the metal ionic radii allows to rationalize variation in the structure types within the ATBi family (A = Li–Cs), correctly identifying KCdBi as isostructural to NaZnBi. Thermal stability studied by high-temperature powder X-ray diffraction reveals that Zn-containing compounds melt at higher temperatures (821 K for KZn 2 Bi 2 ) than Cd-containing KCd 2 Bi 2 (635 K). This study highlights the efficacy of combining rapid synthesis techniques with predictive modeling, though structural predictions show some limitations in accuracy.
We present a systematic high-throughput density-functional-theory investigation of the structural and magnetic stability of 312 substitutional compounds in the magnetic kagome 𝐴𝑇 6 𝑋 6 family. Our screening confirms the stability of many previously reported structures and predicts several additional stable candidates. Within collinear spin configurations, we find that Fe-based systems predominantly adopt antiferromagnetic ground states, whereas Mn-based analogs exhibit a more balanced distribution between ferromagnetic and antiferromagnetic order. For compounds exhibiting several nearly degenerate collinear configurations, we analyze the nature of their magnetic ground states, assess the possible emergence of noncollinear order, and discuss the limitations and uncertainties inherent to standard density-functional approaches. Our electronic structure analysis further reveals that predicted ferromagnetic kagome systems display characteristic features of topological metals, with rich magnetic configurations that can be tuned by chemical substitution. Altogether, these ferromagnetic kagome compounds constitute a broad and still largely unexplored materials platform for the emergence of exciting magnetotransport phenomena.