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.
The novel synthesis of carbon spherogels in organic solvents enables energy and solvent savings and omits supercritical drying. In addition, the resulting carbon spherogels feature excellent capacitance for energy storage in a supercapacitor setup.
Recent developments in artificial intelligence (AI) and machine learning (ML), implemented through self-driving laboratories (SDLs), are rapidly creating unprecedented opportunities for the accelerated discovery and optimization of materials. This paper provides a joint analysis of SDLs from both academic and industry perspectives, highlighting the importance of integrating human intelligence in these systems. It discusses the necessity of careful planning in SDL design across physical, data, and workflow dimensions, including instrumental setup, experimental workflow, data management, and human–SDL interaction. The significance of integrating human input within SDLs, especially as the focus shifts from individual tools and tasks to the creation and management of complex workflows, is emphasized. The paper stresses on the crucial role of reward function design in developing forward-looking workflows and examines the interplay between hardware evolution, ML application across chemical processes, and the influence of reward systems in research. Ultimately, the article advocates for a future where SDLs blend human intuition in hypothesis formulation with AI's precision, speed, and data-handling capabilities.
Carbon monoxide (CO) is a versatile intermediate feedstock for many applications and can be produced through the electrochemical reduction of carbon dioxide (CO 2 ). However, current electrochemical CO production is hindered by high overall costs, primarily due to low conversion efficiencies and significant energy requirements. Herein, we report a unique way of enhancing the electrochemical reduction of CO 2 to CO using gamma (γ) photons. The γ-irradiation applied to the electrochemical cell setup induces the production of e˙ - , which results in an increased CO 2 ionization and production of excited CO 2 (CO$^{*}_{2}$) molecules via lower energy barrier. The ionized CO 2 ˙ - is quickly stabilized over a silver catalyst, providing an alternative low activation energy route for CO 2 reduction. In conclusion, a decrease in the overpotential barrier enhanced the electrochemical reduction of CO 2 to CO by 25%.
Over the past decade, there has been a sustained interest in using frequency-swept (FS) pulses for the efficient acquisition of wideline and ultra-wideline (UW) NMR powder patterns. Such experiments are typically conducted under static conditions, employing both direct- and indirect-excitation methods (i.e., WCPMG and BRAIN-CP/WCPMG, respectively). Recently, Koppe et al. demonstrated that the WCPMG pulse sequence can be used to efficiently acquire wideline and UW NMR spectra with spinning sideband (SSB) manifolds under magic-angle spinning (MAS) conditions, capitalizing on the increased signal-to-noise ratios (SNR) afforded by MAS. To date, there have been only a few instances of broadband cross-polarization (CP) experiments using FS pulses under MAS conditions and no applications to systems exhibiting wideline and/or ultra-wideline powder patterns, despite the clear advantages these experiments could offer. Herein, we demonstrate that FS pulses selectively applied to a single sideband of the S spin can be used for efficient 1 H-S polarization transfer to S = 1/2 nuclides with large anisotropic chemical shift interactions at slow to moderate MAS rates. The Hartmann–Hahn matching conditions in BRAIN-CP/WCPMG-MAS experiments bear similarity to those of standard CP sequences, yet operate over UW frequency ranges and only require low-amplitude RF pulses on the S channel. Crucial to the success of the BRAIN-CP/WCPMG-MAS experiment is careful calibration of the RF amplitude, transmitter offset, and effective frequency sweep of the FS pulse applied to the S spins at a given MAS rate. Thus, by means of numerical simulations and experimental testing, we provide recommendations for the parameterization and setup of BRAIN-CP/WCPMG-MAS experiments for their most efficient use. Results showcasing the capability of the BRAIN-CP/WCPMG-MAS pulse sequence are presented, including applications to 119 Sn, 195 Pt, and 103 Rh NMR.
Machine learning models are often used as scoring functions to predict the binding affinity of a protein–ligand complex. These models are trained with limited amounts of data with experimentally measured binding affinity values. A large number of compounds are labeled inactive through single-concentration screens without measuring binding affinities. These inactive compounds, along with the active ones, can be used to train binary classification models, while regression models are trained using compounds with binding affinities only. However, the classification and regression tasks are often handled separately, without sharing the learned feature representations. In this paper, we propose a novel model architecture that jointly performs regression and classification objectives, aiming to maximize data utilization and improve predictive performance by leveraging two complementary tasks. In our setup, the regression yields the binding affinity, whereas the classification task yields the label as active or inactive. We demonstrate our method using PDBbind, the standard 3D structure database, as well as a dataset of flavivirus protease compounds with binding affinity data. Our experiments show that the new joint training strategy improves the accuracy of the model, increasing applicability in various practical drug screening scenarios.
Lithium–sulfur (Li–S) batteries hold significant promise for electric vehicles and aviation due to their high energy density and cost-effectiveness. However, understanding the root causes of performance degradation remains a formidable challenge, as the interplay of multiple factors obscures key failure mechanisms. A major limitation has been the inability to quantify soluble sulfur species within practical detection limits accurately and to correlate electrochemical processes with associated physical inventory changes. Here, we introduce the high-performance liquid chromatography-ultraviolet spectroscopy and gas chromatography sequential characterization (HUGS) toolkit, capable of precisely quantifying seven distinct sulfur and polysulfide species at concentrations as low as 40 ppb. HUGS has been successfully applied to practical coin and pouch cells without requiring cell modification. Furthermore, our self-developed software, Dr HUGS, enhanced the data analysis speed by over 30 times, enabling multi-source data integration and delivering comprehensive analysis results within minutes. Using HUGS, we identify significant capacity losses from inactive lithium and sulfur during initial cycles and sulfide-rich solid–electrolyte interphase (SEI) formation on the anode during later cycles. Notably, our findings reveal that soluble polysulfides have minimal contributions to capacity loss, challenging long-standing assumptions. Moreover, HUGS demonstrates that constant-pressure setups in Li–S pouch cells improve compositional uniformity compared to constant-gap configurations. For sulfurized polyacrylonitrile (SPAN) cathodes, unique issues such as non-sulfide SEI formation and lithium pulverization are observed, which can be mitigated through localized high-concentration electrolytes to enhance lithium inventory retention. By enabling precise quantification of critical inventory components, HUGS provides transformative insights into failure mechanisms across various electrolytes and cathode chemistries, guiding rational design strategies for next-generation energy storage systems.
Alkaline water electrolysis (AWE) presents key advantages, including reduced material costs, enhanced operational stability, and compatibility with non-precious metal catalysts, positioning it as a scalable route for hydrogen production. In this study, we introduce a minimally invasive single-cell configuration incorporating a reference electrode via diaphragm extension to form an internal ion channel. This setup, combined with an interfaced potentiostat and auxiliary electrometer, enables real-time, independent monitoring of anode and cathode behavior, offering high-resolution electrochemical diagnostics. While it is well established that the hydrogen evolution reaction (HER) exhibits sluggish kinetics in alkaline media, our study reveals that this limitation persists even in practical AWE systems where nickel-based substrates are used as electrodes. This observation is supported by both experimental data and voltage breakdown modeling. Arrhenius-type analysis reveals that localized electric fields induced by catalysts shift the reaction kinetics from classical Butler–Volmer behavior toward a Marcus-like regime, where interfacial molecular dynamics and bimolecular charge transfer dominate. We propose a semi-empirical model and a surficial reaction mechanism to describe these dynamics. This work underscores the critical need for cathode innovation and provides a rational framework for designing advanced catalysts and electrode architectures to optimize AWE performance.
The lithium-mediated nitrogen reduction reaction (Li-NRR) is currently the most promising strategy for electrochemical ammonia synthesis. In this study, we present an operando grazing incidence wide angle X-ray scattering (GI WAXS) investigation using an improved electrochemical flow cell that enables hydrogen oxidation at the anode and thereby eliminates the need of a sacrificial proton donor. The improved cell design also increases nitrogen availability and mass transport, achieving ammonia faradaic efficiencies (FEs) up to 36%. This setup allows direct analysis of reaction intermediates and the solid electrolyte interphase (SEI) using state-of-the-art diglyme-based electrolytes. We identify lithium amide (LiNH2) as the only stable, crystalline intermediate, providing direct insight into the Li-NRR mechanism. Notably, in diglyme-based electrolytes, the SEI composition differs significantly from that in tetrahydrofuran-based systems, with reduced LiF content and the formation of previously unreported crystalline diglyme–lithium salt complexes. These species likely influence ammonia selectivity and long-term stability. Our findings highlight how the electrolyte composition and cell architecture govern Li-NRR selectivity and efficiency, offering a foundation for the rational design of next-generation SEI layers and solid electrolytes to enable scalable electrochemical ammonia synthesis.
Sequestering greenhouse gases (CO 2 and CH 4 ) in biogas into carbon nanofibers (CNF) offers a promising route to mitigate carbon emissions and create value-added solid carbon materials. Coupling non-thermal plasma with a thermocatalytic reactor in a tandem setup is a promising approach for tandem reactions of dry reforming of methane to synthesis gas and its subsequent conversion to CNF. Various parameters were studied to determine their effects on CNF growth. Decreasing the total flow rate resulted in an increase in CNF growth. Increasing the plasma power input or the plasma zone length also enhanced the production of CNF. These results illustrate that plasma-thermal tandem reactors can be used to synthesize CNF from biogas with tunable parameters that may be further optimized in future studies.
The real-time interface chemistry between the lithium cobalt oxide (LCO) working electrode and the LiClO 4 /propylene carbonate (PC) electrolyte is investigated during lithiation/delithiation using dip-and-pull ambient pressure photoelectron spectroscopy (APXPS). The APXPS results appear to exhibit the seldom discussed Co 2+ state in the LCO structure, where the operando measurements indicate electron transfer among Co 2+ , Co 3+ , and Co 4+ states. Specifically, the lithiation of LCO reduces the Co 4+ state to both Co 3+ and Co 2+ states, where, as a function of voltage, reduction to the Co 2+ state is initially more pronounced followed by Co 3+ formation. In addition, a delay in surface delithiation is observed during the reverse potential steps. This is discussed in terms of overpotential at the interface measurement position as a consequence of the dip-and-pull setup for this experiment. Finally, the shifts in the apparent binding energies of the spectral features corresponding to the electrolyte and LCO at their interface show that the electrochemical potentials at delithiation voltage steps are different from the lithiation steps at the same applied voltages. This further explains the non-responsive delithiation. The BE shift observed from the LCO surface is argued to be dominantly due to the semi-conductive nature of the sample. Overall, this article shows the importance of operando APXPS for probing non-equilibrium states in battery electrodes for understanding electron transfer in the reactions.
Efficient hydrogen evolution reaction (HER) catalysts that reduce the use of noble metals and can be synthesized on a large scale are essential for advancing anion exchange membrane water electrolyzers (AEMWEs) toward commercialization. Herein, we present a composite catalyst in which Ru nanoparticles coexist with Ru single-atom alloys (SAAs) dispersed within Ni nanoparticles (Ru-SAA/Ni), creating a highly active HER electrocatalyst. Using a one-pot and scalable synthesis method, we can tune the material composition from SAA, i.e. materials containing atomically dispersed Ru atoms (with ≤0.4 at% Ru) to composite structures in which SAAs coexist with Ru NPs. Comprehensive characterization using XPS, XAS, and TEM confirms Ru-SAA formation at a low Ru content and composite structures at higher contents. Electrochemical evaluations conducted in a three-electrode setup reveal that Ru-SAA/Ni composites achieve HER performance on par with that of Pt/C. Computational insights suggest that water dissociation is significantly faster at the Ru/Ni interface compared with that on extended surfaces. These active sites are thermodynamically as active as basal planes, preventing the excessive accumulation of reaction intermediates (H*, OH*). All these results highlight the synergistic interaction between Ru SAAs and Ru nanoparticles and their potential for large-scale applications with minimal use of precious metals. Finally, the materials are processed and tested in AEMWEs, achieving 1.85 V at 0.5 A cm −2 with a total noble metal loading of only 0.1 mg cm −2 .
As the power grid continues to evolve with advanced wide-area monitoring, protection, and control (WAMPAC) algorithms, there is an increasing need for realistic testbed environments with industry-grade software and hardware-in-the-loop (HIL) to perform verification and validation studies. Such testbed environments serve as ideal platforms to perform WAMPAC prototyping, operator training, and also to study the impacts of different types of cyberattack scenarios on the operation of the grid. In this paper, we introduce PRIME (PNNL cybeR physIcal systeMs tEstbed): the testbed that integrates real-time transmission system simulator with commercial industry grade energy management system (EMS) software and remote hardware-in-the-loop (RHIL). PRIME is an end-to-end, modular testbed that allows high-fidelity RHIL experimentation of a power system. We present two detailed case studies (fault location and clearing in transmission system, and operator training) to showcase the capabilities of our PRIME testbed. Finally, we briefly discuss some of the potential limitations of our testbed in terms of scalability and flexibility to setup larger test systems and identify directions for future work to address those limitations.
Our ability to predict the structure and evolution of stars is in part limited by complex, 3D hydrodynamic processes such as convective boundary mixing. Hydrodynamic simulations help us understand the dynamics of stellar convection and convective boundaries. However, the codes used to compute such simulations are usually tested on extremely simple problems and the reliability and reproducibility of their predictions for turbulent flows is unclear. We define a test problem involving turbulent convection in a plane-parallel box, which leads to mass entrainment from, and internal-wave generation in, a stably stratified layer. We compare the outputs from the codes FLASH, MUSIC, PPMSTAR, PROMPI, and SLH, which have been widely employed to study hydrodynamic problems in stellar interiors. The convection is dominated by the largest scales that fit into the simulation box. All time-averaged profiles of velocity components, fluctuation amplitudes, and fluxes of enthalpy and kinetic energy are within ≲3σ of the mean of all simulations on a given grid (128 3 and 256 3 grid cells), where σ describes the statistical variation due to the flow’s time dependence. They also agree well with a 512 3 reference run. The 128 3 and 256 3 simulations agree within 9% and 4%, respectively, on the total mass entrained into the convective layer. The entrainment rate appears to be set by the amount of energy that can be converted to work in our setup and details of the small-scale flows in the boundary layer seem to be largely irrelevant. Our results lend credence to hydrodynamic simulations of flows in stellar interiors. We provide in electronic form all outputs of our simulations as well as all information needed to reproduce or extend our study.
Star-forming regions may play an important role in the life cycle of Galactic cosmic rays (CRs), notably as home to specific acceleration mechanisms and transport conditions. Gamma-ray observations of Cygnus X have revealed the presence of an excess of hard-spectrum gamma-ray emission, possibly related to a cocoon of freshly accelerated particles. We seek an improved description of the gamma-ray emission from the cocoon using ~13 yr of observations with the Fermi-Large Area Telescope (LAT) and use it to further constrain the processes and objects responsible for the young CR population. We developed an emission model for a large region of interest, including a description of interstellar emission from the background population of CRs and recent models for other gamma-ray sources in the field. Thus, we performed an improved spectro-morphological characterisation of the residual emission including the cocoon. The best-fit model for the cocoon includes two main emission components: an extended component FCES G78.74+1.56, described by a 2D Gaussian of extension r 68 = 4.4° ± 0.1° -0.1° +0.1° and a smooth broken power law spectrum with spectral indices 1.67 ± 0.05 -0.01 +0.02 and 2.12 ± 0.02 -0.01 +0.00 below and above 3.0 ± 0.6 -0.2 +0.0 GeV, respectively; and a central component FCES G80.00+0.50, traced by the distribution of ionised gas within the borders of the photo-dissociation regions and with a power law spectrum of index 2.19 ± 0.03 -0.01 +0.00 that is significantly different from the spectrum of FCES G78.74+1.56. An additional extended emission component FCES G78.83+3.57, located on the edge of the central cavities in Cygnus X and with a spectrum compatible with that of FCES G80.00+0.50, is likely related to the cocoon. For the two brightest components FCES G80.00+0.50 and FCES G78.74+1.56, spectra and radial-azimuthal profiles of the emission can be accounted for in a diffusion-loss framework involving one single population of non-thermal particles with a flat injection spectrum. Particles span the full extent of FCES G78.74+1.56 as a result of diffusion from a central source, and give rise to source FCES G80.00+0.50 by interacting with ionised gas in the innermost region. For this simple diffusion-loss model, viable setups can be very different in terms of energetics, transport conditions, and timescales involved, and both hadronic and leptonic scenarios are possible. The solutions range from long-lasting particle acceleration, possibly in prominent star clusters such as Cyg OB2 and NGC 6910, to a more recent and short-lived release of particles within the last 10–100 kyr, likely from a supernova remnant. The observables extracted from our analysis can be used to perform detailed comparisons with advanced models of particle acceleration and transport in star-forming regions.
The optical depth to reionization, τ, is the least constrained parameter of the cosmological Λ cold dark matter (ΛCDM) model. To date, its most precise value is inferred from large-scale polarized cosmic microwave background (CMB) power spectra from the High Frequency Instrument (HFI) aboard the Planck satellite. These maps are known to contain significant contamination by residual non-Gaussian systematic effects, which are hard to model analytically. Therefore, robust constraints on τ are currently obtained through an empirical cross-spectrum likelihood built from simulations. In this paper, we present a likelihood-free inference of τ from polarized Planck HFI maps which, for the first time, is fully based on neural networks (NNs). NNs have the advantage of not requiring an analytical description of the data and can be trained on state-of-the-art simulations, combining the information from multiple channels. By using Gaussian sky simulations and Planck SRoll2 simulations, including CMB, noise, and residual instrumental systematic effects, we trained, tested, and validated NN models considering different setups. We inferred the value of τ directly from Stokes Q and U maps at ~4° pixel resolution, without computing angular power spectra. On Planck data, we obtained τ NN = 0.0579 ± 0.0082, which is compatible with current EE cross-spectrum results but with a ~30% larger uncertainty, which can be assigned to the inherent nonoptimality of our estimator and to the retraining procedure applied to avoid biases. While this paper does not improve on current cosmological constraints on τ, our analysis represents a first robust application of NN-based inference on real data, and highlights its potential as a promising tool for complementary analysis of near-future CMB experiments, also in view of the ongoing challenge to achieve the first detection of primordial gravitational waves.
Convection is one of the most important mixing processes in stellar interiors. Hydrodynamic mass entrainment can bring fresh fuel from neighboring stable layers into a convection zone, modifying the structure and evolution of the star. Because flows in stellar convection zones are highly turbulent, multidimensional hydrodynamic simulations are fundamental to accurately capture the physics of mixing processes. Under some conditions, strong magnetic fields can be sustained by the action of a turbulent dynamo, adding another layer of complexity and possibly altering the dynamics in the convection zone and at its boundaries. In this study, we used our fully compressible S EVEN -L EAGUE H YDRO code to run detailed and highly resolved three-dimensional magnetohydrodynamic simulations of turbulent convection, dynamo amplification, and convective boundary mixing in a simplified setup whose stratification is similar to that of an oxygen-burning shell in a star with an initial mass of 25 M ⊙ . We find that the random stretching of magnetic field lines by fluid motions in the inertial range of the turbulent spectrum (i.e., a small-scale dynamo) naturally amplifies the seed field by several orders of magnitude in a few convective turnover timescales. During the subsequent saturated regime, the magnetic-to-kinetic energy ratio inside the convective shell reaches values as high as 0.33, and the average magnetic field strength is ∼10 10 G. Such strong fields efficiently suppress shear instabilities, which feed the turbulent cascade of kinetic energy, on a wide range of spatial scales. The resulting convective flows are characterized by thread-like structures that extend over a large fraction of the convective shell. The reduced flow speeds and the presence of magnetic fields with strengths up to 60% of the equipartition value at the upper convective boundary diminish the rate of mass entrainment from the stable layer by ≈20% as compared to the purely hydrodynamic case.
Context. Third-order lensing statistics contain a wealth of cosmological information that is not captured by second-order statistics. However, the computational effort it takes to estimate such statistics in forthcoming stage IV surveys is prohibitively expensive. Aims. We derive and validate an efficient estimation procedure for the three-point correlation function (3PCF) of polar fields such as weak lensing shear. We then use our approach to measure the shear 3PCF and the third-order aperture mass statistics on the KiDS-1000 survey. Methods We constructed an efficient estimator for third-order shear statistics that builds on the multipole decomposition of the 3PCF. We then validated our estimator on mock ellipticity catalogs obtained from N -body simulations. Finally, we applied our estimator to the KiDS-1000 data and presented a measurement of the third-order aperture statistics in a tomographic setup. Results. Our estimator provides a speedup of a factor of ∼100–1000 compared to the state-of-the-art estimation procedures. It is also able to provide accurate measurements for squeezed and folded triangle configurations without additional computational effort. We report a significant detection of tomographic third-order aperture mass statistics in the KiDS-1000 data (S/N = 6.69). Conclusions. Our estimator will make it computationally feasible to measure third-order shear statistics in forthcoming stage IV surveys. Furthermore, it can be used to construct empirical covariance matrices for such statistics.
A neutron induced neutron emission experiment was conducted as the Los Alamos Neutron Science Center (LANSCE) facility at Los Alamos National Laboratory (LANL). In this experiment, a sample was placed in a well collimated neutron beam and was surrounded by an array of 28 fast neutron detectors (EJ-309). The experiment was performed with a neutron flight path of 21.5 m from the source to the sample, and 1 m from the sample to the detectors. The neutron emission from the sample was measured as a function of neutron time of flight covering an incident energy range from 0.7- 20 MeV. The samples included U-235, Pu-239, carbon (graphite), and blanks that matched the encapsulation of the sample. The measured samples were constantly cycled in and out of the neutron beam. This type of experiment measures neutron emission from all reactions occurring in the sample such as fission and elastic and inelastic scattering. Similar to the methodology previously developed at RPI, the measurements were compared with detailed simulations of the experiment using different cross section evaluations for the sample. The observed differences can be attributed to the evaluated neutron cross section and angular distributions. The carbon sample was used as a reference to validate both the experiment and simulation methodology and showed good agreement between experiments and simulations. A review of the experimental setup, analysis methods, and some of the results will be presented.