Search NASA⌕ Search

SEARCH · Search NASA

Results for “Matrix Models”

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.

At least 523 records · Page 29

Determining the N -Representability of a Reduced Density Matrix via Unitary Evolution and Stochastic Sampling

The N-representability problem consists in determining whether, for a given p-body matrix, there exists at least one N-body density matrix from which the p-body matrix can be obtained by contraction, that is, if the given matrix is a p-body reduced density matrix (p-RDM). The knowledge of all necessary and sufficient conditions for a p-body matrix to be N-representable allows the constrained minimization of a many-body Hamiltonian expectation value with respect to the p-body density matrix and, thus, the determination of its exact ground state. However, the number of constraints that complete the N-representability conditions grows exponentially with system size, and hence, the procedure quickly becomes intractable for practical applications. This work introduces a hybrid quantum-stochastic algorithm to effectively replace the N-representability conditions. The algorithm consists of applying to an initial N-body density matrix a sequence of unitary evolution operators constructed from a stochastic process that successively approaches the reduced state of the density matrix on a p-body subsystem, represented by a p-RDM, to a target p-body matrix, potentially a p-RDM. The generators of the evolution operators follow the well-known adaptive derivative-assembled pseudo-Trotter method (ADAPT), while the stochastic component is implemented by using a simulated annealing process. The resulting algorithm is independent of any underlying Hamiltonian, and it can be used to decide whether a given p-body matrix is N-representable, establishing a criterion to determine its quality and correcting it. We apply the proposed hybrid ADAPT algorithm to alleged reduced density matrices from a quantum chemistry electronic Hamiltonian, from the reduced Bardeen–Cooper–Schrieffer model with constant pairing, and from the Heisenberg XXZ spin model. In all cases, the proposed method behaves as expected for 1-RDMs and 2-RDMs, evolving the initial matrices toward different targets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Understanding Differences in Water Adsorption Isotherms: Structural Variations, Force Fields, and Monte Carlo Simulation Approaches

Accurate prediction of water adsorption in micro- and mesoporous materials with hydrophobic pores is essential for the design and characterization of advanced adsorbent materials for separation and energy applications. Here, we assess the reproducibility and consistency of water adsorption isotherms in two microporous all-silica MFI zeolite structures (MFI-K and MFI-O) using two different zeolite force fields and three simulation approaches: grand canonical Monte Carlo (GCMC), Gibbs ensemble Monte Carlo (GEMC), and transition matrix Monte Carlo (TMMC). We demonstrate that consistent treatment of the bulk fluid phase in GCMC and TMMC simulations is critical for reconciling isotherms across methods, and we construct simulation-based equations of state for the TIP4P water model to enable rigorous fugacity-to-pressure conversions. Large shifts in the isotherms are observed for two zeolite force fields developed using different parametrization strategies, with the GCS force field representing implicitly a defect-containing all-silica zeolite, whereas the TraPPE-zeo force field accurately represents an essentially defect-free all-silica zeolite. While water in the van Koningsveld structure of MFI exhibits a first-order phase transition and condensation-like step for adsorption near room temperature, water in the Olson structure of MFI displays continuous adsorption, attributed to differences in the adsorption free energy landscapes. Structural analysis reveals that small geometric variations, particularly Si–O–Si bond angles near the strongest adsorption sites, lead to these substantial differences in adsorption behavior. Furthermore, our results highlight the sensitivity of simulated water adsorption isotherms in hydrophobic frameworks to seemingly small differences in the framework structures, force field parametrization, and simulation approaches.

36 MATERIALS SCIENCE↗

Quantum subspace expansion in the presence of hardware noise

Finding ground state energies on current quantum processing units (QPUs) using algorithms such as the variational quantum eigensolver (VQE) continues to pose challenges. Hardware noise severely affects both the expressivity and trainability of parameterized quantum circuits, limiting them to shallow depths in practice. Here, we demonstrate that both issues can be addressed by synergistically integrating VQE with a quantum subspace expansion, allowing for an optimal balance between quantum and classical computing capabilities and costs. We perform a systematic benchmark analysis of the iterative quantum-assisted eigensolver in the presence of hardware noise. We determine ground state energies of 1D and 2D mixed-field Ising spin models on noisy simulators and the IBM QPUs ibmq_quito (5 qubits) and ibmq_guadalupe (16 qubits). To maximize accuracy, we propose a suitable criterion to select the subspace basis vectors according to the trace of the noisy overlap matrix. Finally, we show how to systematically approach the exact solution by performing controlled quantum error mitigation based on probabilistic error reduction on the noisy backend fake_guadalupe.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Nonrelativistic nuclear reduction for tensor couplings in dark matter direct detection and μ → e conversion

The nonrelativistic effective field theory (NRET) is widely used in dark matter direct detection and charged-lepton flavor violation studies through μ → e conversion. However, existing literature has not fully considered tensor couplings. This study fills this gap by utilizing an innovative tensor decomposition method, extending NRET to incorporate previously overlooked tensor interactions. This development is expected to have a significant impact on ongoing experiments seeking physics beyond the Standard Model and on our understanding of the new-physics interactions. Notably, we identify additional operators in μ → e conversion that are absent in scalar and vector couplings. To support further research and experimental analyses, comprehensive tables featuring tensor matrix elements and their corresponding operators are provided. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

Developing Platinum-Group-Metal-Free Catalysts for Oxygen Reduction Reaction in Acid: Beyond the Single Metal Site

This project is to develop M (x) -N-C catalysts with dense multiple metal center (MMC) sites to meet the DOE 2025 activity target of 0.044 mA/cm 2 at 0.9 V IR-free ., as well as other goals such as durability. We have made important contributions to both catalyst development and fundamental understandings of the active sites in M-N-C catalysts in this project. We successfully made M (x) -N-C catalysts with some multiple metal center (MMC) sites by combining ionothermal carbonization with chemical vapor deposition (CVD). These catalysts, however, are not as active as the most active single-atom Fe-N-C catalysts made by the similar CVD process, owing likely to the low site density and the presence of inorganic Fe species such as iron carbides and nanoparticles. The most significant accomplishments we achieved in this project are: (1) we unraveled the formation pathway of Fe-N 4 sites during pyrolysis step-by-step and identified the trans-metalation mechanism, in collaboration with Deborah Myers and her colleagues at Argonne National Laboratory (ANL); (2) inspired by this finding, we pioneered the CVD synthesis of M-N-C catalysts (M = Mn, Fe, and Co), in which the Fe-N-C catalyst by CVD demonstrated an ORR activity of 33 mA/cm 2 at 0.9 V in H 2 -O 2 proton exchange membrane fuel cells (PEMFCs), very close to the ultimate goal of 35 mA/cm 2 of our project. This catalyst is the first Fe-N-C catalyst that contains only D1 sites without the D2 sites; whereas D1 and D2 sites have been always identified by Mossbauer in previous Fe-N-C catalysts. This finding helps to understand what the D1 and D2 sites are and their roles in catalyzing the ORR. (3) by improving the mass transport of the carbon matrix prior to the CVD process, the revised Fe-N-C catalyst made by CVD delivered a maximum power density of 0.53 W/cm 2 in H 2 -air PEMFCs. The improvement strategy was partly inspired by the computational modeling work by Adam Weber from LBNL, the Co-PI of this project, by developing, coding, and exercising a continuum level model of transport phenomena within a PGM-free catalyst layer. The model demonstrated that local resistances combined with limited site density of the PGM-free catalyst can result in limiting currents and poor polarization performance. The model also gave design guidance for impact of ECSA and overall catalyst-layer thickness. However, both Fe-N-C and Co-N-C catalysts developed by CVD showed poor durability in PEMFCS, in comparison with the traditional M-N-C catalysts synthesized via regular pyrolysis process. Consequently, we did not achieve the proposed durability targets. Despite so, we believe the FeNC-CVD catalysts with the poor durability and D1 sites only provides an excellent platform to understand the degradation mechanism of Fe-N-C in PEMFCs, the most important challenge in the development of M-N-C catalysts.

08 HYDROGEN↗

Rapid Design of Cost-Effective Refractory High Entropy Alloys Strengthened by Precipitation

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.

high temperature materials↗

Beyond leading twist: 𝜌 meson decay constants and distribution amplitudes in a self-consistent light-front quark model

In this study, we present a comprehensive analysis of decay constants and chiral-even and chiral-odd distribution amplitudes (DAs) up to twist 4 for the 𝜌 meson in the standard light-front quark model (LFQM) based on the Bakamjian-Thomas (BT) construction. For the 𝜌 meson, which possesses both longitudinal (ℎ=0) and transverse (ℎ =±1) polarizations, two types of decay constants, 𝑓$^{∥}_{𝜌}$ and 𝑓$^{⊥}_{𝜌}$, arises accordingly. We demonstrate that these decay constants can be self-consistently extracted from both local (𝑧𝜇 =0) and nonlocal (𝑧𝜇 ≠0) matrix elements ⟨0⁢|$\overline{𝑞}$(𝑧)⁢Γ⁢𝑞⁡(−𝑧)|⁢𝜌⁡(𝑃,ℎ)⟩, with Γ=(𝛾 𝜇 ,𝜎 𝜇⁢𝜈 ,𝛾 𝜇 ⁢𝛾 5 ,𝟏), in a manner independent of current components, polarizations, and reference frames. In particular, we emphasize the role of nonlocal matrix elements involving axial-vector and scalar currents, where mixing between 𝑓$^{∥}_{𝜌}$ and 𝑓$^{⊥}_{𝜌}$ occurs. We show that this mixing is consistently resolved through the BT construction, ensuring the proper extraction of these decay constants.

Electroweak interaction↗

A review of low-rank methods for time-dependent kinetic simulations

Time-dependent kinetic models are ubiquitous in computational science and engineering. The underlying integro-differential equations in these models are high-dimensional, comprised of a six–dimensional phase space, making simulations of such phenomena extremely expensive. In this article we demonstrate that in many situations, the solution to kinetics problems lives on a low dimensional manifold that can be described by a low-rank matrix or tensor approximation. We then review the recent development of so-called low-rank methods that evolve the solution on this manifold. The two classes of methods we review are the dynamical low-rank (DLR) method, which derives differential equations for the low-rank factors, and a Step-and-Truncate (SAT) approach, which projects the solution onto the low-rank representation after each time step. Thorough discussions of time integrators, tensor decompositions, and method properties such as structure preservation and computational efficiency are included. We further show examples of low-rank methods as applied to particle transport and plasma dynamics.

97 MATHEMATICS AND COMPUTING↗

Theory of strong-field sequential double ionization of polyatomic molecules

Understanding of strong-field sequential double ionization (SDI) of molecules by a highly intense infrared laser pulse could be the key to observing electron motion in molecules in the attosecond to femtosecond timescale. Based on a novel density matrix approach for SDI (DM-SDI), a recent theoretical study has shown that SDI can be used as a probe to monitor the changes in vibronic coherence in homonuclear diatomic molecules. In this article, we extend the DM-SDI model to general molecules that could possess permanent dipole moments and arbitrary symmetry. We apply the model to SDI of a water molecule and identify the formation pathway of individual dication states. We further deduce that kinetic energy release spectra from two- and three-body fragments could carry the signature of vibronic coherence between the lowest two states of the water cation. Furthermore, our results suggest that observables from the SDI probe can be interpreted intuitively with only the knowledge of electronic structures of populated ionic states, making the SDI probe to be a highly desirable probing scheme for vibronic coherence in generic molecules.

74 ATOMIC AND MOLECULAR PHYSICS↗

Partial-wave projection of relativistic three-body amplitudes

We derive the integral equations for partial-wave projected three-body scattering amplitudes, starting from the integral equations for three-body amplitudes developed for lattice QCD analyses. The results, which hold for generic three-body systems of spinless particles, build upon the recently derived partial-wave projected one-particle exchange, a primary component of the relativistic framework proven to satisfy 𝑆 matrix unitarity. We derive simplified expressions for factorizable short-distance interactions, 𝒦 3 , in two equivalent formalisms—one symmetric under particle interchange and one asymmetric. For the asymmetric case, we offer parametrizations useful for amplitude analysis. Finally, we examine toy models for 3⁢𝜋 systems at unphysically heavy pion masses with total isospins 0, 1, and 2.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Exceptional strength and wear resistance in an AA7075/TiB 2 composite fabricated via friction consolidation

The friction consolidation method successfully reinforced aluminum 7075 alloy (AA7075) with high-volume fractions (12 and 24 vol%) of titanium diboride (TiB 2 ) by high pressure and severe plastic deformation at elevated temperatures. The consolidated AMCs have a uniform dispersion of submicron- and micron-sized TiB 2 particles in the AA7075 matrix, with significant refinement of the matrix grain size and the particles. The addition of TiB 2 significantly increases hardness by up to 50 %, Young’s modulus by up to 62 %, and ultimate tensile strength by up to 28 % to 672 MPa, while reducing ductility by 80 %. Wear resistance of 7075/24 vol% TiB 2 improves seven-fold compared to baseline, making it comparable to that of carburized steels. Microstructure-based finite element modeling provided a theoretical strength limit of ~730 MPa for the composites and indicated that high triaxiality in conjunction with severe equivalent plastic strain in a narrow area between the TiB 2 particles led to early fracture initiations, limiting the ductility.

36 MATERIALS SCIENCE↗

Accelerating Floating-Point Computations with Intel AMX

Intel AMX is a built-in component of recent Intel CPU architectures, first supported by the Intel Sapphire Rapids in 2023, that enables efficient dense matrix multiplications using mixed precision with low-precision data types. The popularity of mixed-precision algorithms has grown recently, primarily due to their use on GPUs to enhance the efficiency of HPC applications, particularly for the training of large language models. The availability of mixed precision on CPUs represents a cost-effective solution for applications where high speed is not critical. This report shows how to use the Intel AMX accelerator through examples in C++ and Python. The examples will focus on mixed-precision floating-point operations obtained by the use of bfloat16 (or BF16) to accelerate code in single precision. We employ a bottom-up methodology, starting from specific register instructions (TMUL operation) to higher-level applications in libraries such as Intel MKL, PyTorch, and TensorFlow, ensuring a comprehensive understanding of the accelerator's potential. Additionally, we provide insights into the expected performance gains when leveraging the accelerator on the Kestrel HPC machine at the National Renewable Energy Laboratory.

97 MATHEMATICS AND COMPUTING↗

Rapid Bayesian High Entropy Alloy Designs Fabricated via Wire Arc Additive Manufacturing

Purpose: This project seeks to demonstrate a new high-throughput (rapid) alloy design technique applied to creating new high entropy alloys (HEAs) for extreme environments. High entropy alloys shift the design paradigm from being focused on a single principal element (e.g. nickel-based alloys) to target alloys that include high atomic fractions (X >10%) of multiple elements. These HEA materials can exhibit sluggish diffusion and enhanced corrosion resistance, ideal for potential applications in advanced ultra supercritical (A-USC) steam cycles for power generation. Scope: The addition of multiple elements in high atomic fractions creates an enormous design space that cannot easily be investigated by traditional material design strategies such as designed of experiments (DOE). This project utilizes a Bayesian machine learning algorithm that has been modified to work with calculation of phase diagrams (CALPHAD) software. This Bayesian algorithm reduces manual inputs and increase the likelihood of achieving an optimal solution. Compositional inputs to this algorithm will be assessed using existing material property models for high temperature strength and corrosion resistance. The target for alloy performance will be a 15% (~100 ⁰C) increase in allowable service temperature beyond heat-resistant stainless steels while maintaining or improving alloy cost and corrosion resistance. Haynes 230 was selected as a baseline, which is 57 wt% Ni with 22 wt% Cr 14 wt% W, and 2 wt% Mo as solid solution strengtheners. In addition to rapid design via Bayesian machine learning, the alloys were rapidly fabricated using a multi-wire arc additive manufacturing (mWAAM) technique which allows for precise control of alloy composition and assessing of alloy design “windows” to study composition effects. Build speeds for wire-arc additive processes are among the highest for additive technologies enabling rapid and reliable sample fabrication when compared to conventional methods such as arc button melting. The mWAAM samples will be rapidly characterized via instrumented indentation for room temperature modulus and strength and for elevated temperature strength via hot hardness tests. After being screened with hardness testing, potential alloys will be further evaluated with conventional microscopy techniques including scanning electron microscopy (SEM) and transmission electron microscopy (TEM) to assess agreement with modeling results. The most promising compositions will also be evaluated by printing full sized tensile specimens for mechanical behavior tests at elevated temperatures. Results: Bayesian machine learning of a single performance function was initially used to optimize five performance metrics: 1) single phase stability, 2) yield strength, 3) creep resistance (low diffusion coefficient), 4) freezing range (weldability), and 5) material cost. The single performance function was suboptimal as assumptions had to be made about the results while formulating the optimization. A goal-oriented Bayesian optimization strategy (Hanaoka, 2021) was implemented with CALPHAD for use with the five metrics above. This multi-objective Bayesian optimization (MOBO) enabled the design of NiCrCoFe alloys with V and W additions. A base composition of NiCoCr was selected as Ni provides a stable FCC matrix, Cr aids corrosion/oxidation resistance, and Co is a solid-solutions strengthener that also improves creep by increasing the activation energy. Fe helps reduce diffusion coefficients and cost. Finally, V and W were selected for their reasonable solubility and high atomic misfit to aid in solid solution strengthening. Cracking of the mWAAM specimens was an early issue, and the Easton solidification cracking model (Easton et al., 2014a) was selected for addition to the MOBO function. High performing alloys fabricated by mWAAM included Ni 28 Cr 25 Co 26 Fe 15 V 8 and Ni 62 Cr 18 Co 1 Fe 3 W 15 . It was observed that even after adapting the mWAAM process for W, the W did not fully dissolve. To fully evaluate the Ni 62 Cr 18 Co 1 Fe 3 W 15 composition, a cored wire (80-20 NiCr sheath/powder core) was manufactured and printed via WAAM, and HIP’ing was utilized to homogenize and densify the printed alloy. The V and W alloys produced met metrics 1 (solid solution), 4 (solidification cracking), and 5 (cost). However, an unmodeled mechanism of thermal stress cracking was identified in the WAAM produced materials, perhaps exacerbated by the lack of grain boundary strengthening elements (B, C). Conclusions & Recommendations: A high-throughput (rapid) alloy design technique was applied to designing and manufacturing new high entropy alloys (HEAs) for extreme environments utilizing MOBO and mWAAM. The developed process was rapid and effective in addressing the mechanisms included in the model. The lack of grain boundary strengthening element additions (e.g., B, C) was a simplification that likely produced thermal stress cracking that turned into a large part of the investigation. Additions on the order of 0.005 wt% B and 0.05 wt% C likely would have minimized thermal stress grain boundary cracking. Overall, the high throughput design strategy is promising for rapid design of metrics-driven alloys for advanced ultra supercritical (A-USC) steam cycles for power generation. The MOBO and mWAAM process could be commercialized to accelerate metrics-driven alloy design. In addition, the cored-wire process utilized for scale-up is a promising high-volume process for WAAM alloy development and scale-up.

36 MATERIALS SCIENCE↗

Time-dependent-bases with local CUR decomposition method for accelerating turbulent combustion simulations

Here, this study presents a novel reduced-order modeling framework, Time-Dependent Bases with Local CUR decomposition (TDB-L-CUR), designed to efficiently and accurately approximate the species transport equations in reacting flow simulations. The method extends the existing TDB-CUR approach for chemically reacting flows (Jung et al. Comput. Methods Appl. Mech. Engrg. 437 (2025) 117758), which leverages matrix decomposition techniques to form a global-in-space, time-dependent low-dimensional manifold. While TDB-CUR performs well in homogeneous systems, it may be less well-suited to spatially heterogeneous systems such as turbulent flames, where higher-rank approximations are typically required. The proposed TDB-L-CUR framework introduces two methodological extensions to the baseline approach. First, it applies unsupervised clustering to partition the physical domain into distinct regions, enabling spatially localized manifold construction, thereby reducing the rank required for the reduced-order representation. Second, it incorporates a computational singular perturbation (CSP)-based scheme for identifying and penalizing fast species, allowing for spatio-temporally adaptive mitigation of chemical stiffness. The proposed framework is validated on a hierarchy of test cases, including a one-dimensional premixed flame, a two-dimensional nonpremixed ignition case with vortex interaction, and a three-dimensional turbulent premixed flame. TDB-L-CUR significantly improves accuracy over TDB-CUR while further reducing computational cost. The fully on-the-fly formulation of TDB-L-CUR (i.e., requiring no offline training or prior knowledge) makes it a robust and scalable tool for reduced-order modeling of reactive flows.

Local manifold↗

Accurate numerical simulations of open quantum systems using spectral tensor trains

Decoherence between qubits is a major bottleneck in quantum computations. Decoherence results from intrinsic quantum and thermal fluctuations as well as noise in the external fields that perform the measurement and preparation processes. With prescribed colored noise spectra for intrinsic and extrinsic noise, we present a numerical method, Quantum Accelerated Stochastic Propagator Evaluation (Q-ASPEN), to solve the time-dependent noise-averaged reduced density matrix in the presence of intrinsic and extrinsic noise. Q-ASPEN is arbitrarily accurate and can be applied to provide estimates for the resources needed to error-correct quantum computations. We employ spectral tensor trains, which combine the advantages of tensor networks and pseudospectral methods, as a variational ansatz to the quantum relaxation problem and optimize the ansatz using methods typically used to train neural networks. Here, the spectral tensor trains in Q-ASPEN make accurate calculations with tens of quantum levels feasible. We present benchmarks for Q-ASPEN on the spin-boson model in the presence of intrinsic noise and on a quantum chain of up to 32 sites in the presence of extrinsic noise. In our benchmark, the memory cost of Q-ASPEN scales as a low-order polynomial in the size of the system once the number of system states surpasses the number of basis functions used in the spectral expansion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Transient fuel performance analysis for the preliminary fuel concept of general atomics fast modular reactor

This study investigates the transient fuel performance of General Atomics Fast Modular Reactor (GA-FMR) during accident scenarios, focusing on the behavior of its innovative fuel system that combines high-assay low enriched uranium dioxide (HALEUO2) fuel with SiGA® ceramic matrix composite silicon carbide cladding. The preliminary fuel design’s response was analyzed during reactivity-initiated accidents (RIA) and loss of coolant accidents (LOCA) using BISON fuel performance analysis code, which included both the diffusion enhanced and BISON-FASTGRASS coupled UO 2 models. The RIA analysis demonstrated that effective reactivity control reduced fuel temperature, though with transient fission gas release resulting in additional tensile stress state on the cladding. LOCA simulations revealed differing predictions between the two models: the BISON UO 2 model showed more transient fission gas release but minimal pellet expansion, while the BISON-FASTGRASS UO 2 model predicted less pronounced fission gas release but more fuel swelling and thermal expansion, potentially leading to pellet-cladding mechanical interaction. Here, these findings highlight critical areas for fuel design optimization and identify knowledge gaps requiring further experimental and computational investigation to advance GA-FMR fuel development.

Lee, Soon K. [Argonne National Laboratory (ANL), A↗

Explainable machine learning for incipient anomaly detection in compact molten salt heat exchanger with overlapping feature distributions

High-temperature molten salt-cooled reactors (MSCRs) are a promising next-generation nuclear technology option, offering efficient power conversion and inherent safety features. However, the reliability of these systems depends on the robust operation of heat exchangers (HXs), which are susceptible to failure due to temperature gradients and channel plugging caused by fluid freezing. Conventional monitoring methods, relying on inlet and outlet measurements, lack the spatial resolution needed to detect early-stage faults. We propose a novel design of a compact salt-to-salt matrix-type HX design consisting of interleaved arrays of parallel tubes, with integrated synthetic fiber optic distributed temperature sensing (DTS) to enable localized detection of incipient faults. To evaluate performance of this design, we generate high-fidelity synthetic data using heat transfer computational modeling to simulate channel plugging, and introduce sensor noise for realistic modeling of measurements. The dataset comprises of 97% normal operation and 3% anomaly cases, with each anomaly class representing 1% of the data. These early anomalies result in overlapping temperature profiles between normal and faulty channels, producing a non-separable dataset that challenges traditional classification techniques. We benchmark eight supervised machine learning (ML) models and demonstrate that XGBoost achieves the highest performance. To improve transparency, we develop an explainability framework combining Shapley values and partially ordered sets (POSETs) to quantify and structurally analyze feature importance. This approach identifies both dominant predictors and ambiguous feature relationships, enhancing trust and interpretability. Our results highlight the potential of combining DTS and explainable ML with intelligent feature selection to improve predictive maintenance and ensure operational resilience in advanced nuclear systems.

Prantikos, Konstantinos [Argonne National Laborato↗

The Gene Ontology knowledgebase in 2026

Abstract The Gene Ontology (GO) knowledgebase (https://geneontology.org) is a comprehensive resource describing the functions of genes. The GO knowledgebase is regularly updated and improved. We describe here the major updates that have been made in the past 3 years. The ontology and annotations have been expanded and revised, particularly in several areas of biology: cellular metabolism, multi-organism interactions (e.g. host-pathogen), extracellular matrix proteins, chromatin remodeling (e.g. the “histone code”), and noncoding RNA functions. We have released version 2 of a comprehensive set of integrated, reviewed annotations for human genes, which we call the “functionome.” We have also dramatically increased the number of GO-CAM models, with over 1500 models of metabolic and signaling pathways, primarily in human, mouse, budding and fission yeast, and fruit fly. Finally, we discuss our current recommendations and future prospects of AI in the use and development of GO.

Aleksander, Suzi A (ORCID:0000000167872901)↗