SEARCH · Search NASA
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.
This report includes 1) Advanced Fuels Campaign, 2) Material Recovery and Waste Forms Development, 3) MPACT Campaign, 4) Systems Analysis and Integration (SA&I) Campaign), 5) Joint Fuel Cycle Study Activities, and 6) AFCI_HQ Program Support.
Explore the source record for details and available documents.
Motivation: Predicting microbial gene fitness across environmental conditions remains a central challenge for predictive phenomics and autonomous experimentation. Fitness assays generate large volumes of genotype–phenotype measurements difficult to integrate with experimental metadata and biological function in a form that supports mechanistic reasoning. Knowledge graphs offer a semantic framework for unifying modalities and enabling context-aware inference. Results: We build GIMME (Graph Inference for Microbial Metabolism Exploration), a semantically grounded knowledge graph that unifies gene fitness measurements spanning 10 Pseudomonas species with experimental metadata and biological context. Media are decomposed into chemical components and experiments carry structured links to natural-language descriptions. The resulting graph supports two inference modes: (1) symbolic graph traversal to surface candidate gene–environment and gene–chemical associations, and (2) learned inference using heterogeneous graph neural networks that propagate information across neighborhoods. We formulate link regression over (gene, media, experiment) triplets, combining learned gene embeddings with pretrained LLM sourced text embeddings of node descriptions to predict gene fitness. We then augment a baseline MLP with an auxiliary message-passing encoder (GraphSAGE/GAT) that propagates information over gene–protein–function and media–chemical subgraphs, and fuse the two pathways with a gated residual connection. This approach produces strong agreement with held-out fitness measurements (GraphSAGE Pearson r 0.74) while also highlighting inference challenges in extreme-fitness regimes. We aggregate GAT edge-attention weights by relation type and layer to estimate which biological and environmental relations most influence fitness predictions. Conclusion: This work explores using knowledge graphs as “context graphs” for microbial phenotype prediction. They provide a rich substrate which enables explainable retrieval of supporting evidence, and provides a natural bridge to autonomous workflows that prioritize the next experiment.
Microwave connectors with a planar geometry have been demonstrated using additive manufacturing techniques. This planar form factor is expected to enable connectors to be integrated directly with other components onto a printed circuit board. The connectors were developed with multiple signal lines that will increase the density of signal connections to PCBs, eliminating the physical constraints of single-signal COTS microwave connectors. This work developed electromagnetic models (using Ansys HFSS) of these planar connectors to predict performance and optimize connector design. Materials and printing processes were developed to fabricate the connectors based on high-temperature PEEK thermoplastics. The multiple signal connectors were designed to operate over the 1-6GHz band and are built upon the single signal printed planar geometry connector introduced in a previous paper [1].
As the integration of renewable energy resources such as wind and solar generation continues to grow, the use of hybrid power plants with both synchronous machines and power electronic inverters operating in parallel becomes more common. This paper studies the impact on small-signal stability for these pairs of devices, first as an isolated two-device system and then as part of a 9-bus test system. Each operational scenario was studied with both a grid-following inverter and a grid-forming inverter. In each scenario, the state variables associated with the synchronous generator's excitation system were found to be responsible for the inception of instability. These results underscore the importance of accurate modeling of excitation systems when modeling power systems with hybrid power plants.
Characterizing structural ensembles of intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs) of proteins is essential for studying structure–function relationships. Due to the different neutron scattering lengths of hydrogen and deuterium, selective labeling and contrast matching in small-angle neutron scattering (SANS) becomes an effective tool to study dynamic structures of disordered systems. However, experimental timescales typically capture measurements averaged over multiple conformations, leaving complex SANS data for disentanglement. We hereby demonstrate an integrated method to elucidate the structural ensemble of a complex formed by two IDRs. We use data from both full contrast and contrast matching with residue-specific deuterium labeling SANS experiments, microsecond all-atom molecular dynamics (MD) simulations with four molecular mechanics force fields, and an autoencoder-based deep learning (DL) algorithm. From our combined approach, we show that selective deuteration provides additional information that helps characterize structural ensembles. We find that among the four force fields, a99SB-disp and CHARMM36m show the strongest agreement with SANS and NMR experiments. In addition, our DL algorithm not only complements conventional structural analysis methods but also successfully differentiates NMR and MD structures which are indistinguishable on the free energy surface. Finally, we present an ensemble that describes experimental SANS and NMR data better than MD ensembles generated by one single force field and reveal three clusters of distinct conformations. Our results demonstrate a new integrated approach for characterizing structural ensembles of IDPs.
This TCF project developed a thin-oxide (SiO 2 )/polycrystalline silicon (poly-Si) passivated contact solar cell on CubicPV's (formally 1366 Technologies, Inc.) Direct Wafer® Product (DWP) kerfless wafers. The project used two NREL-developed technologies described in U.S. Patent No. 9,911,873, Hydrogenation of Passivated Contacts and U.S. Patent Application Serial No. 15/890,172, Doped Passivated Contacts . The project was motivated by a potential higher efficiency cell (compared to a PERC cell) using passivated contacts on the ultra-low cost kerfless wafers grown using the Direct Wafer process. The hope was to accelerate market adoption of the cell and wafer by delivering the lowest LCOE in the PV industry. The project tested both n-type and p-type SiO 2 /poly-Si passivated contacts grown by thermal oxidation and plasma enhanced chemical vapor deposition (PECVD) of the poly-Si layer on DWP with varying wafer resistivities. Both deposition techniques are industry standards and thus economically viable methods for commercializing the contacts. The results indicated that both n-type and p-type poly-Si passivated contacts can be formed on polycrystalline DWP wafers, but implied open-circuit voltages (i Voc ) were limited to below 0.65 mV (compared with ~ 730 mV on n-Cz wafers). Diffusion of H to the Si/SiO 2 /poly-Si interface was key to obtaining high i Voc values. In this study, H was diffused from a high-temperature SiN x layer deposited over the poly-Si layer during a high-temperature firing step, similar to one used for screen printed metals. The study concluded that poly-Si passivated contacts on DWP wafers passivated the surface of the wafers as well as PERC passivated surfaces, which use a less expensive dielectric layer stack. The project showed that Direct Wafer Product wafers grown by CubicPV could produce high i Voc values (~0.647 mV), which could produce a cell over 20% efficient with proper processing and metallization. These cells, though not economically viable in 2024 as a stand-alone cell, could be integrated with a wide-bandgap top solar cell to form a two-junction tandem cell that could be viable under certain circumstances. This is because the bottom cell of a 30%, two-terminal tandem only needs to be a 20% cell under one-sun conditions. Thus, the DWP could be an ideal low-cost wafer for tandems. The project also revealed that a TOPCon type cell could be formed on a p-type DWP wafer using a P-diffused emitter and a p-type poly-Si contact. In fact, the p-type version of the poly-Si contact out-performed the n-type version for a variety of wafer resistivities, from highly doped to lowly doped. This curiosity requires more work to understand because on Cz wafers, the n-type poly-Si contact is of much higher quality than the p-type version.
A guided cold-atom inertial sensor system comprises an atom trap integrated platform, a laser system, a magnetic field system, a control system, and a computing system. The laser system and magnetic field system are adapted to form a magneto-optical trap (MOT) about a suspended waveguide of the atom trap integrated platform made of membrane integrated photonics. After loading cold atoms from a MOT, the photonic atom trap integrated platform generates one-dimensional guided atoms with an evanescent field optical dipole trap (EF-ODT) along the optical waveguide to create guided atomic accelerometers/gyroscopes. Motion of atomic wavepackets in a superposition state is created along the guided atom geometry by way of state-dependent momentum kicks. The light-pulse sequence of guided atom interferometry splits, redirects, and recombines atomic wavepackets, which allows measurement of atom interference fringes sensitive to inertial forces via a probe laser.
Numerical diffuse scattering cross-section calculations are used to establish a rigorous basis for determining the concentration and size distribution of dislocation loops in irradiated single crystals from integral X-ray diffuse scattering (XRDS) measurements. Differential XRDS intensities for prismatic {111} type dislocation loops are numerically calculated as a function of loop radius R and wavevectors q relative to Bragg reflections in tungsten. The results show the well known 1/q 2 Huang scattering form at small q that transitions to a ∼1/q 4 dependence associated with the Stokes–Wilson approximation for q ≳ 1/R. More importantly, they show further that the 1/q 4 falloff is not the asymptotic large-q form of the diffuse scattering for small loops (R < 200 Å) as has often been assumed. Rather, for loop sizes as small as R ≃ 5 Å with strong curvature, the calculations show definitively that the scattering transitions to a robust 1/q 5 falloff at larger q that arises due to the local strains near the dislocation core defining the circumference of the dislocation loops. The presence of this 1/q 5 asymptotic form for both small and large loops is experimentally confirmed using an integral XRDS measurement around the 110 reflection on self-ion-irradiated tungsten combined with numerically calculated integral XRDS cross-sections. Accordingly, the historical two-region theoretical treatment of the cross-sections for integral XRDS is extended to a three-region model that has direct sensitivity to the (first-moment) dislocation line lengths of dislocation loops. In conclusion, these developments enable the use of both numerical and analytically modeled cross-sections to make accurate integral XRDS determinations of dislocation loop sizes and concentrations using modest-intensity laboratory X-ray sources.
Engineered heterostructures derive distinct properties from materials integration and interface formation. Two-dimensional crystals have been combined to form vertical stacks and lateral heterostuctures with covalent line interfaces. While thicker vertical stacks have been realized, lateral heterostructures from multilayer van der Waals crystals, which could bring the benefits of high-quality interfaces to bulk-like layered materials, have remained much less explored. In this work, we demonstrate the integration of anisotropic layered Sn and Ge monosulfides into complex heterostructures with seamless lateral interfaces and tunable vertical design using a two-step growth process. The anisotropic lattice mismatch at the lateral interfaces between GeS and SnS is relaxed via dislocations and interfacial alloying. Nanoscale optoelectronic measurements by cathodoluminescence spectroscopy show the characteristic light emission of joined high-quality van der Waals crystals. Spectroscopy across the lateral interface indicates valley-selective luminescence in the bulk SnS component that arises due to anisotropic electron transfer across the interface. The results demonstrate the ability to realize high-quality lateral heterostructures of multilayer van der Waals crystals for diverse applications, e.g., in optoelectronics or valleytronics.
We consider the problem of minimizing multi-modal loss functions with a large number of local optima. Since the local gradient points to the direction of the steepest slope in an infinitesimal neighborhood, an optimizer guided by the local gradient is often trapped in a local minimum. To address this issue, we develop a novel nonlocal gradient to skip small local minima by capturing major structures of the loss’s landscape in black-box optimization. The nonlocal gradient is defined by a directional Gaussian smoothing (DGS) approach. The key idea of DGS is to conducts 1D long-range exploration with a large smoothing radius along d orthogonal directions in Rd, each of which defines a nonlocal directional derivative as a 1D integral. Such long-range exploration enables the nonlocal gradient to skip small local minima. The d directional derivatives are then assembled to form the nonlocal gradient. We use the Gauss-Hermite quadrature rule to approximate the d 1D integrals to obtain an accurate estimator. The superior performance of our method is demonstrated in three sets of examples, including benchmark functions for global optimization, and two real-world scientific problems.
The gravitational path integral is usually implemented with a covariant action by analogy with other gauge field theories, but the gravitational case is different in important ways. A key difference is that the integrand has an essential singularity, which occurs at zero lapse where the spacetime metric degenerates. The lapse integration contour required to impose the local time reparametrization constraints must run from − ∞ to + ∞ , yet must not pass through zero. This raises the question: for an application—such as a partition function—where the constraints should be imposed, what is the correct integration contour, and why? We study that question by starting with the reduced phase space path integral, which involves no essential singularity. We observe that if the momenta are to be integrated before the lapse, to obtain a configuration space path integral, the lapse contour should pass below the origin in the complex lapse plane. This contour is also consistent with the requirement that quantum field fluctuation amplitudes have the usual short distance vacuum form, and with obtaining the Bekenstein-Hawking horizon entropy from a Lorentzian path integral. Published by the American Physical Society 2025
The project Site Characterization of the Highest-Priority Geologic Formations for CO2 Storage in Wyoming is one of 9 site characterization projects that were implemented as part of ARRA (American Recovery and Reinvestment Act). Data from this project was used to improve resolution of data in NATCARB in the area of study. Data related to this study has already been incorporated in NATCARB Atlas. The Wyoming Carbon Underground Storage Project (WY-CUSP) consisted of CO2 storage site characterization and evaluation, focusing on Wyoming’s most promising CO2 storage reservoirs (the Pennsylvanian Weber/Tensleep Sandstone and Mississippian Madison Limestone) and premier CO2 storage site (Rock Springs Uplift). Results from the WY-CUSP project suggest the two reservoirs could store up to 17,000 million tons of CO2. The WY-CUSP team drilled a stratigraphic test well and acquired a 3-D seismic survey covering 25 square miles of the Rock Springs Uplift site. The team retrieved 916 feet of core from the 12,810-foot-deep well, along with a complete log suite, borehole images, fluid samples, and other data. Project partners (1) provided continuous visual documentation of the core, including grain size, mineralogy, facies distribution, and porosity; (2) performed continuous permeability and velocity scans of selected reservoir intervals; and (3) chemically analyzed the fluid samples. WY-CUSP scientists integrated seismic attributes with observations from log suites, a VSP survey, core, fluid samples, and laboratory analyses, including continuous permeability scans. From these integrations, researchers constructed 3-D spatial distribution volumes of reservoir and seal properties that represent geological heterogeneity at the targeted CO2 storage site. The WY-CUSP team used this data to perform new CO2 plume migration simulations. Baker Hughes, Inc., completed a series of small-scale, in-situ water injectivity measurements. A database was formed when observations, analyses, and experiments from the stratigraphic test well were integrated. Correlation of these data allowed petrophysical parameters to be extrapolated from the test well out into the storage domain (5x5 mile 3-D seismic survey volume). This resulted in an improved, realistic understanding of performance assessments for potential CO2 storage scenarios. The WY-CUSP team worked on (1) improving CO2 storage resource estimates, (2) establishing long-term integrity and permanence of confining layers, (3) designing a profitable strategy for pressure management, and (4) evaluating the utilization of stored CO2 at the Rock Spring Uplift. Finally, Baker Hughes developed a microseismic baseline for the test site using in-bore geophones to complete field operations.
We report that recently, there has been a significant interest in detonation-based combustion systems, such as rotating detonation engines (RDE), due to potential performance advantages in propulsion and energy applications. This includes improvements in thermodynamic efficiency, inherently high gas velocities, and the possibility of achieving an increase in total pressure (i.e. pressure gain) through the combustion process. Due to approximating a constant volume combustion process, these devices also tend to exhibit extremely high local gas temperatures relative to a comparable constant pressure combustor. These advantages overlap well with desirable performance characteristics of direct power extraction (DPE) technologies, such as a magnetohydrodynamic (MHD) generator. Typically, in DPE systems hot combustion products are seeded with an easily ionizable material such as potassium carbonate (K 2 CO 3 ) in order to boost the electrical conductivity. However, due to the short gas residence times within an RDE, it was unclear whether forming an electrically conductive combustion plasma would be feasible for integration with a downstream MHD generator. A model is presented which describes the heating, decomposition, and ionization of solid K 2 CO 3 particles and aqueous solutions of K 2 CO 3 in water, for a given initial particle diameter. This model was combined with available computational fluid dynamics (CFD) data for an oxygen-methane RDE in a one-way coupled Eulerian-Lagrangian framework to predict particle trajectories and the corresponding heating, decomposition, and ionization histories. Electrical conductivities were computed using a previously developed model, and a method was proposed to determine an equivalent average electrical conductivity. Results show that particle sizes below ~30 µm are able to fully decompose before reaching the exit of the RDE. While the one-way coupled nature of the simulations precluded rigorous evaluation of the effects of seed material on detonability, a substantial temperature reduction is expected at the detonation wave front due to heating and decomposition. A preliminary comparison is presented between the RDE and an equivalent constant pressure adiabatic combustor, showing a potential performance advantage for the RDE.
Membrane technology lies at the heart of many industrial gas separation processes and applications. Molecular sieving membranes that break the Robeson limit are desirable for energy-efficient gas separation. In this paper, we report a facile strategy of directly integrating ionic liquids (ILs) into porous membranes. Particularly, the ILs form an ultra-thin layer on the carbon molecular sieve (CMS) membranes rather than penetrating into the pores, acting as a smart gate for gas entry to boost the selectivity. The hybrid membrane exhibits CO 2 permeability >600 barrer and enhanced CO 2 /N 2 selectivity >50, which surpasses the Robeson limit and shows potential in CO 2 /N 2 separation process. Molecular dynamics simulations confirm the gating effect of the IL layer of molecular thickness. In conclusion, this work demonstrates a universal strategy to improve CMS membrane performance by creating an IL-membrane interface and tuning the ion-pore interaction.
Here, we present a novel discontinuous Galerkin algorithm for the solution of a class of Fokker–Planck collision operators. These operators arise in many fields of physics, and our particular application is for kinetic plasma simulations. In particular, we focus on an operator often known as the ‘Lenard–Bernstein’ or ‘Dougherty’ operator. Several novel algorithmic innovations, based on the concept of weak equality, are reported. These weak equalities are used to define weak operators that compute primitive moments, and are also used to determine a reconstruction procedure that allows an efficient and accurate discretization of the diffusion term. We show that when two integrations by parts are used to construct the discrete weak form, and finite velocity-space extents are accounted for, a scheme that conserves density, momentum and energy exactly is obtained. One novel feature is that the requirements of momentum and energy conservation lead to unique formulas to compute primitive moments. Careful definition of discretized moments also ensure that energy is conserved in the piecewise linear case, even though the kinetic-energy term, $v^2$ is not included in the basis set used in the discretization. A series of benchmark problems is presented and shows that the scheme conserves momentum and energy to machine precision. Empirical evidence also indicates that entropy is a non-decreasing function. The collision terms are combined with the Vlasov equation to study collisional Landau damping and plasma heating via magnetic pumping.
Bilayered vanadium oxides are attractive for energy storage due to their high initial specific capacities, which could be stabilized by integrating the bilayers with conductive nanoflakes often produced in a form of aqueous dispersions. Therefore, exfoliation of the bilayered vanadium oxides in water with high yield is desirable. Here, this work introduces the first aqueous exfoliation of chemically prelithiated bilayered vanadium oxide (i.e., δ-Li x V 2 O 5 ·nH 2 O or LVO) followed by vacuum filtration to produce a free-standing film, exhibiting a lamellar stacking of the bilayered vanadium oxide nanoflakes as evidenced by scanning electron microscopy. Due to the hydrated nature of bilayered vanadium oxides, the relationship between interlayer water content and the vacuum drying temperature (105 °C vs 200 °C) was studied using X-ray diffraction, thermogravimetric analysis, and Raman spectroscopy. It was found that vacuum drying the LVO nanoflakes at 200 °C enabled more efficient removal of crystallographic water than drying at 105 °C, and did not induce a phase transformation. Scanning transmission electron microscopy confirmed the layered structure of the samples, which was more well-ordered in the 200 °C case and had no clear boundaries between flakes at the atomic scale. Furthermore, electrochemical testing in nonaqueous Li-ion cells revealed that vacuum drying at 200 °C led to improvements in ion storage capacity and electrochemical stability. Improvements in electrochemical charge storage properties of the electrodes obtained via LVO exfoliation and free-standing film assembly in water dried at 200 °C reveal that conventional battery electrode drying protocols need to be revised as new electrochemically active materials are synthesized, such as hydrated layered oxides with expanded interlayer regions. The remaining capacity fading can be attributed to the structural LVO degradation, dissolution of vanadium oxide in electrolyte, and parasitic effects of the remaining interlayer water molecules. Our results establish an environmentally friendly and safe approach to obtain two-dimensional (2D) bilayered vanadium oxide nanoflakes and create a pathway to constructing novel 2D heterostructures for improved performance in energy storage applications.
Cooperativity is used by living systems to circumvent energetic and entropic barriers to yield highly efficient molecular processes. Cooperative structural transitions involve the concerted displacement of molecules in a crystalline material, as opposed to typical molecule-by-molecule nucleation and growth mechanisms which often break single crystallinity. Cooperative transitions have acquired much attention for low transition barriers, ultrafast kinetics, and structural reversibility. However, cooperative transitions are rare in molecular crystals and their origin is poorly understood. Crystals of 2-dimensional quinoidal terthiophene (2DQTT-o-B), a high-performance n-type organic semiconductor, demonstrate two distinct thermally activated phase transitions following these mechanisms. Here we show reorientation of the alkyl side chains triggers cooperative behavior, tilting the molecules like dominos. Whereas, nucleation and growth transition is coincident with increasing alkyl chain disorder and driven by forming a biradical state. We establish alkyl chain engineering as integral to rationally controlling these polymorphic behaviors for novel electronic applications.