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.
A new family of Isothermal, redox-activated CO 2 sorbents were successfully developed using a high-throughput combinatorial approach to facilitate the generation of green hydrogen from biogenic carbonaceous feedstocks.
Advancing thermal/environmental barrier coating (TEBC) materials with integrated thermal-mechanical functions is paramount for safeguarding SiC-based ceramic matrix composites (CMCs) in high-efficiency gas turbines. Herein, we employ a synergistic approach, combining density functional theory (DFT) methods and combinatorial chemistry techniques, to design high-performance and low-cost RE 2 Si 2 O 7 (RE = rare earth elements) TEBC materials tailored for enhanced compatibility with SiC-based CMCs. Expanding on phase stability of alloying pure RE 2 Si 2 O 7 , the investigation extends to the mechanical and thermal properties of solid solution systems, including Er 1/2 Y 3/4 Yb 3/4 Si 2 O 7 , Gd 1/4 Er 1/4 Y 3/4 Yb 3/4 Si 2 O 7 , and Eu 1/4 Er 1/4 Y 3/4 Yb 3/4 Si 2 O 7 . The solid solution systems exhibit a major reduction in lattice thermal conductivity relative to their pure counterparts, achieving ultralow values of 0.25 to 0.39 W m −1 K −1 at 1500 K. Furthermore, the coefficients of thermal expansion (CTE) of these solid solutions are precisely tuned within the desired range for SiC (4.4 to 5.5 × 10 −6 K −1 ), while maintaining good mechanical properties. Here, in particular, the addition of Eu 2 Si 2 O 7 demonstrates to be an important variable to the tuning of CTE and lattice thermal conductivity by leveraging its strong anharmonicity, presenting a pioneering avenue for fine-tuning material properties. In summary, this research not only identifies promising TEBC materials with superior thermal properties, but also introduces a valuable computational material design methodology for the rapid discovery of complex materials for harsh environments.
Ionic liquids (ILs), which are a class of materials with versatile nature and growing popularity, are facing impediments toward widespread usage as electrolytes due to various factors such as low ionic conductivity, high viscosity, high market price etc. One of the ways these limitations can be addressed is by mixing ILs with a molecular solvent. In a combinatorial sense, there exists an immense number of specific IL–solvent combinations. An exhaustive experimental or even simulation-based investigation of the chemical space spanned by such combinations can be extremely time-consuming, expensive, and nearly impossible. An alternative approach is to employ machine learning-based models developed from available databases. Although there exists prior literature that integrates machine learning to investigate mixtures of specific solvents with ILs, these models lack generalization necessitating development of a large number of ML models to handle various solvents. To remedy this shortcoming, as a part of designing green electrolytes with high ionic conductivity that can have potential applications in next-generation batteries and solar cells, this work aims to develop a unified machine learning model to predict ionic conductivity of any IL–solvent mixture system. In this regard, three models, namely, Random Forest, extreme gradient boosting (XGBoost), and artificial neural network (ANN) were formulated using the NIST ILThermo database. The dataset contained 549 unique ionic liquids from 16 cation families and 81 unique solvents, representing a total of 23 712 datapoints. SHAPLEY additive explanation (SHAP) method was used to assess the impact of various features on model prediction and their significance was compared with literature to gain physical insight about the model behavior. Finally, using the developed models, approximately 2.5 million IL–solvent mixtures at five different compositions were screened at room temperature. The high-throughput screening yielded nearly 19 000 IL–solvent mixtures for which ionic conductivity was found to exceed the ionic conductivity of conventional Li-ion battery electrolyte.
Covering: 2020 to 2025Natural products remain indispensable sources of therapeutic and bioactive compounds, yet traditional discovery strategies are constrained by compound rediscovery. Modular biosynthetic enzymes, such as type I polyketide synthases (PKSs) and type A non-ribosomal peptide synthetases (NRPSs), offer promising platforms for combinatorial biosynthesis owing to their programmable architectures. However, practical implementation is frequently limited by inter-modular incompatibility and domain-specific interactions. This review highlights recent advances in modular enzyme assembly enabled by synthetic interfaces-including cognate docking domains, synthetic coiled-coils, SpyTag/SpyCatcher, and split inteins-which function as orthogonal, standardized connectors to facilitate post-translational complex formation. These interfaces support rational investigations into substrate specificity, module compatibility, and pathway derivatization as well as general enzyme clustering applications beyond PKS and NRPS systems. Synthetic interfaces can be integrated with computational tools to support a more systematic and scalable framework for modular enzyme engineering by providing predictive insights into domain compatibility and interface design. These approaches within iterative design-build-test-learn workflows can accelerate the programmable assembly of biosynthetic systems and expand the accessible chemical space for natural products.
Photoelectrochemical (PEC) carbon dioxide reduction reaction (CO 2 RR) has been considered as a promising route to convert and store solar energy into chemical fuels. It is crucial to find suitable photoelectrode materials that are photo-catalytically active and exhibit excellent photochemical stability. One of the promising contenders is ZnTe with the ∼2.26 eV band gap and prolonged stability under CO 2 RR PEC conditions. Herein, a new telluride based thin-film ZnGa 2 Te 4 photocathode with lower band gap and stronger visible light absorption compared to ZnTe is synthesized and characterized using a combinatorial sputtering technique. A two-step annealing method with excess Te supply is implemented to synthesize nearly stoichiometric ZnGa 2 Te 4 absorber material with a zincblende-derived tetragonal crystal structure confirmed by synchrotron X-ray and electron diffraction. Theoretical calculations show that ZnGa 2 Te 4 has suitable direct bandgap (∼1.86 eV) and high absorption coefficient ∼10 5 cm −1 , in agreement with experimentally prepared films. Transient absorption spectroscopy reveals the biexponential decay dynamics, with time constants, τ 1 ∼ 0.04, and τ 2 ∼ 0.65 μs in microsecond time scales and provides the optical transition pathways for this semiconductor thin film. PEC measurements show that the ZnGa 2 Te 4 photocurrent densities are comparable to the widely investigated ZnTe photocathodes or even surpass it under simulated sunlight condition. ZnGa 2 Te 4 samples demonstrate promising photoelectrochemical stability, maintaining consistent performance under illumination. The inclusion of diaryliodonium additive substantially increases its CO 2 RR selectivity to ∼60%. These findings open a new avenue for the synthesis of telluride-based thin-film photocathodes for further exploration and will motivate future research to integrate this potential photocathode material into PEC devices.
Accurately predicting solid electrolyte interphase (SEI) formation requires explicitly resolving the electric double layer (EDL) structure, which deviates significantly from that of the bulk electrolyte. Although an established molecular dynamics (MD) and Density Functional Theory (DFT) framework can model SEI formation by evaluating reduction reactions of local clusters in the EDL, it suffers from a combinatorial computational bottleneck. To overcome this limitation, we introduce a machine-learning-accelerated simulation workflow (MD–DFT–ML), integrating a gradient-boosted regression model trained on EDL composition data to efficiently predict reduction potentials. We apply this framework to seven fluorinated electrolytes comprising fluorinated anions, a fluorinated ester solvent, two types of diluent (ion-solvating ester vs. non-solvating ether), and an FEC additive. The analysis shows that the EDL selectively accumulates cation-binding species; consequently, the non–cation-binding ether diluent rarely enters the EDL and makes minimal contributions to SEI formation. DFT calculations on statistically representative EDL clusters provide reduction potentials and fluorine-release pathways, while the ML model, which substantially reduces the DFT workload, predicts cluster reduction energies with a mean absolute error of 0.1 eV. The combined MD–DFT–ML approach also quantifies contributions from different sources to LiF formation in the SEI. This methodology establishes a generalizable route for multiscale modeling electrolyte and interphase design for next-generation electrochemical energy-storage systems.
Thin-film material libraries in the ternary and quaternary metal oxide systems Fe–V–O, Cu–V–O, and Cu–Fe–V–O were synthesized using combinatorial reactive co-sputtering with subsequent annealing in air. Their compositional, structural, and functional properties were assessed using high-throughput characterization methods. Prior to the investigation of the quaternary system Cu–Fe–V–O, the compositions (Fe 61 V 39 )O x and (Cu 52 V 48 )O x with promising photoactivity were identified from their ternary subsystems Fe–V–O and Cu–V–O, respectively. Two Cu–Fe–V–O material libraries with (Cu 29-72 Fe 4-27 V 22-57 )O x and (Cu 11-55 Fe 27-73 V 12-34 )O x composition spread were investigated. Seven mixed ternary and quaternary phase regions were identified: I (α-Cu 3 FeV 6 O 26 /FeVO 4 ), II (Cu 5 V 2 O 10 /FeVO 4 /α-Cu 3 Fe 4 V 6 O 26 ), III (Cu 5 V 2 O 10 ), IV (Cu 5 V 2 O 10 /FeVO 4 , V (FeVO 4 /γ-Cu 2 V 2 O 7 /α-Cu 3 Fe 4 V 6 O 26 ), VI (β-Cu 2 V 2 O 7 /α-Cu 3 Fe 4 V 6 O 26 /FeVO 4 ), and VII (β-Cu 3 Fe 4 V 6 O 26 /FeVO 4 ). Furthermore, in the investigated composition range, two photoactive regions, (Cu 53 Fe 7 V 40 )O x and (Cu 45 Fe 21 V 34 )O x , were identified, exhibiting 103 μ A/cm 2 and 108 μ A/cm 2 photocurrent density for the oxygen evolution reaction at 1.63 V vs reversible hydrogen electrode, respectively. The highest photoactive region (Cu 45 Fe 21 V 34 )O x comprises the dominant α-Cu 3 Fe 4 V 6 O 24 phase and minor FeVO 4 phase. This photoactive region corresponds to having an indirect bandgap of 1.87 eV and a direct bandgap of 2.58 eV with an incident photon-to-current efficiency of 30% at a wavelength of 310 nm.
Scanning Thermal Microscopy (SThM) provides efficient thermal property measurement with micro- or nanoscale spatial resolution. However, the sensitivity and accuracy of state-of-the-art thermal probes have been limited by excessive thermal contact resistance between the probe and sample. Introduced herein is a robust thermal microprobe that can increase the probe-sample contact force by more than two orders of magnitude, thereby reducing the probe-sample thermal contact resistance by as much as 96% and increasing measurement sensitivity by more than 240% compared to a commercial thermal probe with the same dimensions and measurement principle. The relationship between the probe-sample thermal contact resistance, thermal exchange radius, and sample thermal conductivity is determined experimentally. Simultaneous measurement of thermal conductivity and Seebeck coefficient with unprecedented sensitivity is demonstrated using the enhanced scanning thermal microprobe on samples of an extended range of thermal conductivity up to 18 W/m K, increasing the range of samples applicable to SThM when compared to the conventional commercial probe with diminished measurement sensitivity above ~10 W/m K. The probe is further demonstrated by simultaneously mapping thermal conductivity and Seebeck coefficient as a function of depth from the irradiated surface of an ion-irradiated bulk nanostructured thermoelectric material. In addition to enabling microscale thermal conductivity and Seebeck coefficient measurement of materials previously not applicable to SThM, the probe can also facilitate high-throughput characterization of combinatorial materials to aid the rapid discovery of compositions and processing conditions that yield highly desired thermal and thermoelectric properties.
In recent years, there is growing interest in using quantum computers for solving combinatorial optimization problems. In this work, we developed a generic, machine learning-based framework for mapping continuous-space inverse design problems into surrogate quadratic unconstrained binary optimization (QUBO) problems by employing a binary variational autoencoder and a factorization machine. The factorization machine is trained as a low-dimensional, binary surrogate model for the continuous design space and sampled using various QUBO samplers. Using the D-Wave Advantage hybrid sampler and simulated annealing, we demonstrate that by repeated resampling and retraining of the factorization machine, our framework finds designs that exhibit figures of merit exceeding those of its training set. We showcase the framework’s performance on two inverse design problems by optimizing (i) thermal emitter topologies for thermophotovoltaic applications and (ii) diffractive meta-gratings for highly efficient beam steering. This technique can be further scaled to leverage future developments in quantum optimization to solve advanced inverse design problems for science and engineering applications.
Surface termination on a graphitic surface and the type of electrolytes in lithium-ion batteries (LIBs) play an important part in determining the structure, composition, and thus, the quality of the emergent solid electrolyte interphase. In this paper, we analyze the structure and dynamics of electrolyte molecules in multi-component electrolyte with varying species compositions combinatorially paired with four different graphitic surfaces terminated with hydrogen, hydroxyl, carbonyl, and carboxyl to explore the interplay between surface chemistry and electrolyte dynamics at electrode/electrolyte interfaces. Addition of dimethyl carbonate and fluoroethylene carbonate brought substantial changes in the ethylene carbonate (EC) and LiPF6 surface population density for hydroxyl and carbonyl surfaces. Strong density oscillation and drastic slowing of the dynamics of the electrolyte molecules at the interface are reported for all the systems. While these observations are universal, carboxyl surfaces have the strongest local and long-range effects. Characterization of the average dipole direction at the interface shows strong orientational preferences of ethylene carbonate molecules. EC molecules are preferred to be oriented either almost parallel or perpendicular to the hydroxyl surface, are tilted between parallel and perpendicular with a higher angle of incidence of the dipole vs surface normal on the carbonyl surface than on the hydroxyl surface, and are oriented perpendicularly against the carboxyl surface. These differences highlight the significant effect of graphite surface termination on the dynamics of the electrolytes and provide insight into the complex interplays between electrolyte species and graphite anode in LIBs.
Here a novel dual-energy fast neutron imaging technique is presented using short-pulse laser-driven neutron sources to leverage their inherent adaptive spectral control to enable 3D volume segmentation and reconstruction. Laser-accelerated ion beams incident onto secondary targets create directional, broadband, MeV-class neutrons. Synthetic radiographs are produced of multi-material objects using ion and neutron spectra derived from analytic and numerical models. It is demonstrated that neutron images generated from small changes to the neutron spectra, controlled by altering the initial laser conditions, are sufficient to isolate materials with differing attenuation coefficients. This is first demonstrated using a simplistic combinatorial isolation method and then by employing more advanced reconstruction algorithms to reduce artifacts and generate a segmentation volume of the constituent materials.
Autonomous experimentation is an emerging area of research, primarily related to autonomous vehicles, scientific combinatorial discovery approaches in materials science and drug discovery, and iterative research loops of planning, experimentation, and analysis. However, autonomous approaches developed in these contexts are difficult to apply to high-dimensional mapping technologies, such as scanning hyperspectral imaging of biological systems, due to sample complexity and heterogeneity. We briefly cover the history of adaptive sampling algorithms and surrogate modeling in order to define autonomous adaptive data acquisition as an objective-based, flexible building block for future biological imaging experimentation driven by intelligent infrastructure. We subsequently summarize the recent implementations of autonomous adaptive data acquisition (AADA) for scanning hyperspectral imaging, assess how these address the difficulties of autonomous approaches in hyperspectral imaging, and highlight the AADA design variation from a goal-oriented perspective. Finally, we present a modular AADA architecture that embeds AADA-driven flexible building blocks to address the challenge of time resolution for high-dimensional scanning hyperspectral imaging of nonequilibrium dynamical systems. In our example research-driven experimental design case, we propose an AADA infrastructure for time-resolved, noninvasive, and label-free scanning hyperspectral imaging of living biological systems. This AADA infrastructure can accurately target the correct state of the system for experimental workflows that utilize subsequent expensive, high-information-content analytical techniques.
Oxide ionotronics is an interdisciplinary field in which systems and devices rely on the migration of ions/ionic defects to alter or drive functionality. In this perspective, we focus on epitaxial oxide heterostructures and the contributing roles of oxygen vacancies and interfaces in ionotronics. We begin with a description of oxygen vacancy behavior, with a focus on vacancy ordering and the effects of interfaces and electric fields on particular epitaxial oxide systems. We then emphasize the use of synchrotron x-ray techniques for investigating system structure and dynamics in situ at interfaces and surfaces. Finally, an outlook on the future of epitaxial oxide ionotronics is provided, and several key areas for research are identified, such as freestanding heterostructures, combinatorial synthesis and machine learning, and next-generation synchrotron x-ray characterization.
Nitride perovskites ABN3 are an emerging and highly underexplored class of materials that are of interest due to their intriguing calculated ferroelectric, optoelectronic, and other functional properties. Incorporating novel A-site cations is one strategy to tune and expand such properties; for example, Gd 3+ is compelling due to its large magnetic moment, potentially leading to multiferroic behavior. However, the theoretically predicted ground state of GdWN 3 was a non-perovskite monoclinic structure. Here, we experimentally show that GdWN 3-y crystallizes in a perovskite structure. High-throughput combinatorial sputtering with activated nitrogen is employed to synthesize thin films of Gd 2-x W x N 3-y O y with oxygen content y < 0.05. Ex situ annealing crystallizes a polycrystalline perovskite phase in a narrow composition window near x = 1. LeBail fits of synchrotron grazing incidence wide angle x-ray scattering data are consistent with a perovskite ground-state structure. Refined density functional theory calculations that included antiferromagnetic configurations confirm that the ground-state structure of GdWN 3 is a distorted Pnma perovskite with antiferromagnetic ordering, in contrast to prior predictions. Initial property measurements find that GdWN 3-y is paramagnetic down to T = 2 K with antiferromagnetic correlations and that the absorption onset depends on cation stoichiometry. This work provides an important path toward both the rapid expansion of the emerging family of nitride perovskites and understanding their potential multiferroic properties.
Physical imaging is a foundational characterization method in areas from condensed matter physics and chemistry to astronomy and spans length scales from atomic to universe. Images encapsulate crucial data regarding atomic bonding, materials microstructures, and dynamic phenomena such as microstructural evolution and turbulence, among other phenomena. The challenge lies in effectively extracting and interpreting this information. Variational Autoencoders (VAEs) have emerged as powerful tools for identifying the underlying factors of variation in image data, providing a systematic approach to distilling meaningful patterns from complex data sets. However, a significant hurdle in their application is the definition and selection of appropriate descriptors reflecting local structures. Here, we introduce the scale-invariant VAE approach (SI-VAE) based on the progressive training of the VAE with the descriptors sampled at different length scales. The SI-VAE allows the discovery of the length scale-dependent factors of variation in the system. Here, we illustrate this approach using the ferroelectric domain images and generalize it to the movies of the electron-beam induced phenomena in graphene and topography evolution across combinatorial libraries. This approach can further be used to initialize the decision making in automated experiments including structure–property discovery and can be applied across a broad range of imaging methods. This approach is universal and can be applied to any spatially resolved data including both experimental imaging studies and simulations, and can be particularly useful for exploration of phenomena such as turbulence and scale-invariant transformation fronts.
Highly alloyed (Al,Gd)N is of potential interest in a variety of applications, including neutron detection and in devices such as non-volatile memory. Gd has been shown to have very low equilibrium solubility in AlN at room temperature; however, non-equilibrium deposition techniques such as sputtering are able to deposit thin films, which incorporate large amounts of Gd. Here, we characterize a highly-alloyed (Al,Gd)N combinatorial thin film grown by RF sputtering on a GaN substrate, looking for any evidence of chemical or phase segregation or structural disorder in the films. Compositions with between 13% and 32% Gd (on a cation basis) were studied. No evidence was found for chemical or phase segregation in any studied composition. Higher degrees of Gd incorporation led to greater structural disorder in the film and a tendency toward amorphization; however, electron diffraction shows that the film does not become fully amorphous at any of the studied compositions, instead retaining textured local order even at 32% Gd. Electron energy loss spectra suggest that the material retains a locally wurtzite-like tetrahedral bonding environment at all studied compositions.
The evolution of the atomic structures of the combinatorial library of Sm-substituted thin film BiFeO 3 along the phase transition boundary from the ferroelectric rhombohedral phase to the non-ferroelectric orthorhombic phase is explored using scanning transmission electron microscopy. Localized properties, including polarization, lattice parameter, and chemical composition, are parameterized from atomic-scale imaging, and their causal relationships are reconstructed using a linear non-Gaussian acyclic model. This approach is further extended to explore the spatial variability of the causal coupling using the sliding window transform method, which revealed that new causal relationships emerged at both the expected locations, such as domain walls and interfaces, and at additional regions forming clusters in the vicinity of the walls or spatially distributed features. While the exact physical origins of these relationships are unclear, they likely represent nanophase-separated regions in the morphotropic phase boundaries. Overall, we posit that an in-depth understanding of complex disordered materials away from thermodynamic equilibrium necessitates understanding not only the generative processes that can lead to observed microscopic states but also the causal links between multiple interacting subsystems.
Quantum computers and simulators may offer significant advantages over their classical counterparts, providing insights into quantum many-body systems and possibly improving performance for solving exponentially hard problems, such as optimization and satisfiability. Here, we report the implementation of a low-depth Quantum Approximate Optimization Algorithm (QAOA) using an analog quantum simulator. We estimate the ground-state energy of the Transverse Field Ising Model with long-range interactions with tunable range, and we optimize the corresponding combinatorial classical problem by sampling the QAOA output with high-fidelity, single-shot, individual qubit measurements. We execute the algorithm with both an exhaustive search and closed-loop optimization of the variational parameters, approximating the ground-state energy with up to 40 trapped-ion qubits. We benchmark the experiment with bootstrapping heuristic methods scaling polynomially with the system size. We observe, in agreement with numerics, that the QAOA performance does not degrade significantly as we scale up the system size and that the runtime is approximately independent from the number of qubits. We finally give a comprehensive analysis of the errors occurring in our system, a crucial step in the path forward toward the application of the QAOA to more general problem instances.