Search NASA⌕ Search

SEARCH · Search NASA

Results for “laser-based spectroscopic techniques”

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.

1,149 records · Page 33
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Mechanisms of Protonic Nonvolatile Memory Device

A nonvolatile memory device based on protonic transport in oxides has been proposed. The mobile H+ ions are introduced into the SiO2 layer by annealing Si/SiO2/Si structures in H2 at temperatures greater than 500 deg C. This effect has only been observed for confined oxides that have been annealed at greater than or equal to 1100 C prior to the hydrogenation anneal. This includes buried oxides such as Unibond and SIMOX as well as thermal oxides annealed with a polysilicon cap. An applied field moves the charge within the oxide and the charge stops moving when the field is removed. In a memory device, the hydrogen-annealed oxide is the gate oxide and the position of the mobile charge is sensed by the shift of the I-V curve. Much is still not understood about the motion of the charge across the buried oxide. Previous work has assumed that H+ transport and the time it takes to traverse the oxide is governed by interactions within the bulk of the oxide. Based on parameters that affect the transport time, we conclude that H+ trapping and detrapping at the Si/SiO2 interface are more important than H+ interactions within the oxide bulk. These parameters include the applied field, the H+ concentration and the oxide thickness. One consequence is that projections of device write-time based on the previous assumptions of H+ transport mechanisms may be overly optimistic.

P J Macfarlane↗
DOE OSTI2025

32 examples of LLM applications in materials science and chemistry: towards automation, assistants, agents, and accelerated scientific discovery

Abstract Large language models (LLMs) are reshaping many aspects of materials science and chemistry research, enabling advances in molecular property prediction, materials design, scientific automation, knowledge extraction, and more. Recent developments demonstrate that the latest class of models are able to integrate structured and unstructured data, assist in hypothesis generation, and streamline research workflows. To explore the frontier of LLM capabilities across the research lifecycle, we review applications of LLMs through 32 total projects developed during the second annual LLM hackathon for applications in materials science and chemistry, a global hybrid event. These projects spanned seven key research areas: (1) molecular and material property prediction, (2) molecular and material design, (3) automation and novel interfaces, (4) scientific communication and education, (5) research data management and automation, (6) hypothesis generation and evaluation, and (7) knowledge extraction and reasoning from the scientific literature. Collectively, these applications illustrate how LLMs serve as versatile predictive models, platforms for rapid prototyping of domain-specific tools, and much more. In particular, improvements in both open source and proprietary LLM performance through the addition of reasoning, additional training data, and new techniques have expanded effectiveness, particularly in low-data environments and interdisciplinary research. As LLMs continue to improve, their integration into scientific workflows presents both new opportunities and new challenges, requiring ongoing exploration, continued refinement, and further research to address reliability, interpretability, and reproducibility.

Computer Science↗
DOE OSTI2026

Achieving uniaxial magnetic anisotropy in Ce2⁢Fe17⁢N3 through Co- and Sm-substitution

Th2⁢Zn17−type structure-based permanent magnets, such as Sm2⁢Fe17⁢N3, offer strong potential as alternatives to neodymium magnets (NdFeB), but their practical use is limited by phase stability and the scarcity of Sm. Ce-based counterparts, particularly Ce2⁢Fe17⁢N3, are attractive low-cost candidates, yet their intrinsic planar magnetic anisotropy restricts permanent-magnet performance. Here, we induce uniaxial magnetic anisotropy in Ce2⁢Fe17⁢N3 through two approaches: (i) Co substitution on the Fe sublattice and (ii) partial substitution of Ce with Sm. Combined density functional theory and experimental results show that both strategies modify the 3⁢𝑑–4⁢𝑓 interactions and band filling, yielding magnetization values up to ∼1.2T and magnetocrystalline anisotropy energies exceeding 1MJ/m3 for Co-alloyed compositions, with significantly larger anisotropy achieved upon Sm substitution. In addition, the Sm-substituted Ce2⁢Fe17⁢N3 samples exhibit enhanced high-temperature stability compared to Sm2⁢Fe17⁢N3. These findings demonstrate that Ce2⁢Fe17⁢N3-based alloys can deliver magnetic performance suitable for permanent-magnet applications while reducing cost and reliance on critical rare-earth elements, and they provide practical design guidelines for rare-earth-lean magnets for energy and industrial applications.

Pokhrel, Nabaraj [ORNL] (ORCID:0000000328283076)↗
DOE OSTI2025

Evidence for large baryonic feedback at low and intermediate redshifts from kinematic Sunyaev-Zel’dovich observations with ACT and DESI photometric galaxies

Recent advances in cosmological observations have provided an unprecedented opportunity to investigate the distribution of baryons relative to the underlying matter. In this work, we show that the gas is more extended than the dark matter, and the amount of baryonic feedback at $z \lesssim 1$ disfavors low-feedback models such as that of state-of-the-art hydrodynamical simulation IllustrisTNG compared with high-feedback models such as that of the original Illustris simulation. This has important implications for bridging the gap between theory and observations and understanding galaxy formation and evolution. Furthermore, a better grasp of the baryon-dark matter link is critical to future cosmological analyses, which are currently impeded by our limited knowledge of baryonic feedback. Here, we measure the kinematic Sunyaev-Zel'dovich (kSZ) effect from the Atacama Cosmology Telescope (ACT), stacked on the luminous red galaxy (LRG) sample of the Dark Energy Spectroscopic Instrument (DESI) imaging survey. This is the first analysis to use photometric redshifts for reconstructing galaxy velocities. Due to the large number of galaxies comprising the DESI imaging survey, this is the highest signal-to-noise stacked kSZ measurement to date: we detect the signal at 13$σ$, finding strong evidence that the gas is more spread out than the dark matter, as well as a preference for larger feedback compared to some commonly used state-of-the-art hydrodynamical simulations. Here, our work opens up the possibility of recalibrating large hydrodynamical simulations using the kSZ effect. In addition, our findings highlight the importance of properly accounting for baryonic feedback with future surveys such as LSST through direct probes such as the kSZ, and shed light on long-standing enigmas in astrophysics, such as the “missing baryon” problem.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Integration of Online Cross-Section Generation Capability with Depletion and Transient Solvers in Griffin

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE)-based reactor multiphysics analysis application jointly developed by Argonne and Idaho National Laboratories under the DOENE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In FY25, an online crosssection generation capability based on the Self-Shielding Application Programming Interface (SSAPI) was demonstrated for TRISO-fueled reactor problems under steady-state conditions. This fiscal year, that capability was extended to support depletion and transient multiphysics calculations, enabling high-fidelity analyses that generate self-shielded cross sections on the fly from the actual evolving composition and temperature states rather than from pre-tabulated libraries. For depletion, a two-way coupling was established in which SSAPI computes compact-averaged self-shielded cross sections that the depletion solver then uses to advance the Bateman equations, with the updated compositions returned to SSAPI at each step; the depletion module was refactored to support both library-based and SSAPI-based cross sections, and additional logic was added to track daughter isotopes and to exclude minor isotopes for efficiency. For transient analysis, the SSAPI multigroup library was extended with the kinetics data required for time-dependent calculations, the Improved Quasi-Static (IQS) scheme was coupled with SSAPI, and several supporting capabilities were implemented, including a self-shielding treatment that lets control rods and drums move within a self-shielded model, which had previously been impossible and had ruled out rod- and drum-movement transients with on-the-fly cross sections altogether, a new mixing scheme for delayed-neutron precursor decay constants, a checkpoint-based restart workflow, and performance improvements such as pointwise cross-section interpolation and the bypassing of unnecessary Dancoff factor calculations. The implemented capabilities were verified against Serpent Monte Carlo solutions. For depletion, a prismatic pin-cell problem based on a Next Generation Nuclear Plant (NGNP) Very High Temperature Reactor benchmark showed excellent agreement, with eigenvalue differences within 200 pcm over the entire burnup range (up to 140 MWD/kgU) and fission-product and actinide inventories agreeing to within 0.8% and 2.5%, respectively; a heat-pipe microreactor assembly problem with a much higher fuel loading confirmed the same behavior and quantified the bias introduced when the multigroup equivalence effect is neglected. For transient analysis, a pin-cell problem with a step reactivity insertion and temperature feedback reproduced the analytically expected asymptotic power and showed close agreement between the direct and IQS solutions, and a two-dimensional microreactor core problem with control-drum rotation exercised the new moving-drum self-shielding treatment and demonstrated successful coupling of the online crosssection generation with both the direct and IQS transient methods. The capability was further exercised on a full-core pebble-bed problem, in which Griffin was coupled with the System Analysis Module (SAM) to simulate load-following operation of the gPBR with the Doppler feedback resolved at the TRISO fuel kernel temperature. These developments in Griffin provide a convenient, high-fidelity approach to cross-section generation for advanced thermal reactors with geometrically complex and highly heterogeneous configurations, including TRISO-fueled prismatic and pebble-bed systems, and support steady-state, depletion, and transient multiphysics calculations. They also enable self-shielded cross sections to be evaluated directly at the actual coupled state of the system, thereby establishing a foundation for high-fidelity, fully coupled multiphysics analysis of advanced reactors

Park, H.↗
NASA NTRS1984

The Telecommunications and Data Acquisition Report

Tracking and ground-based navigation; communications, spacecraft-ground; station control and system technology; capabilities for new projects; networks consolidation program; and network sustaining are described.

E C Posner↗
NASA NTRS2026

Parameterization of Nuclear Electric Propulsion Trajectories for Outer Solar System Science Missions Methodology

This manuscript addresses the methodology used to create a database of low thrust missions to outer planets. This database utilizes previous work modeling NEP systems to determine the maximum delivered mass to outer planets based on a range of mission parameters, such as time of flight, launch vehicle, and power system mass. Trajectories were selected which best utilized NEP benefits. Additionally, a discussion on the database outputs for missions to Saturn is included, such as time of flight based on trajectory type and maximum payload, given a specific launch vehicle. The purpose of this work was to create a basis for future mission design, and a tool to investigate general trends across mission options.

Outer Planets↗
DOE OSTI2024

O Corona, where art thou? eROSITA’s view of UV-optical-IR variability-selected massive black holes in low-mass galaxies

Finding massive black holes (MBHs,M BH ≈ 10 4 –10 7 M ⊙ ) in the nuclei of low-mass galaxies $\left( {{M_*}\mathop {\mathop < \limits_ }\limits_ {{10}^{10}}{M_ \odot }} \right)$ is crucial to constrain seeding and growth of black holes over cosmic time, but it is particularly challenging due to their low accretion luminosities. Variability selection via long-term photometric ultraviolet, optical, or infrared (UVOIR) light curves has proved effective and identifies lower-Eddington ratios compared to broad and narrow optical spectral lines searches. In the inefficient accretion regime, X-ray and radio searches are effective, but they have been limited to small samples. Therefore, differences between selection techniques have remained uncertain. Here, we present the first large systematic investigation of the X-ray properties of a sample of known MBH candidates in dwarf galaxies. We extracted X-ray photometry and spectra of a sample of ~200 UVOIR variability-selected MBHs and significantly detected 17 of them in the deepest available SRG/eROSITA image, of which four are newly discovered X-ray sources and two are new secure MBHs. This implies that tens to hundreds of LSST MBHs will have SRG/eROSITA counterparts, depending on the seeding model adopted. Surprisingly, the stacked X-ray images of the many non-detected MBHs are incompatible with standard disk-corona relations, typical of active galactic nuclei, inferred from both the optical and radio fluxes. They are instead compatible with the X-ray emission predicted for normal galaxies. After careful consideration of potential biases, we identified that this X-ray weakness needs a physical origin. A possibility is that a canonical X-ray corona might be lacking in the majority of this population of UVOIR-variability selected low-mass galaxies or that unusual accretion modes and spectral energy distributions are in place for MBHs in dwarf galaxies. This result reveals the potential for severe biases in occupation fractions derived from data from only one waveband combined with SEDs and scaling relations of more massive black holes and galaxies.

Astronomy & Astrophysics↗
DOE OSTI2026

Innovations Driven by Advanced Characterization to Strategize Critical Mineral Production and Beneficial Reuse from Fossil Energy Waste

Critical minerals (CM), such as rare earth elements (REE), cobalt, nickel, and lithium, have important uses in modern electronics and advanced manufacturing, yet are vulnerable to potential supply chain disruptions. Relatively abundant and readily available fossil energy (FE) wastes, such as coal combustion ash, acid mine drainage (AMD) and treatment solids (AMD solids), and Oil and Gas (O&G) drilling wastes (drill cuttings and produced waters) are under consideration as CM feedstocks. The National Energy Technology Laboratory (NETL) has studied CM resources for various FE wastes as part of the U.S. Department of Energy’s mission of bolstering the domestic CM supply, and makes the data available to the public on EDX at sites such as the NEWTS group. Advanced characterization utilizing synchrotron x-ray techniques coupled with laboratory extractions has been performed to identify CM hosting phases in these FE wastes to inform CM recoverability mechanisms. Novel methods to selectively recover CMs while co-producing other valuable byproducts have been developed. Successful examples discussed here include: (1) The identification of REE/Co/Ni/Sc binding and hosting phases in select FE waste (coal combustion ash and AMD solids), resulting in the development of a patented CM step-extraction process, (2) coupled production of functional sorbents from these extraction wastes and for CM recovery. A pilot-scale testing to evaluate the patent’s technical feasibility for extracting REE from coal ash on a barrel scale has been successfully performed. Additionally, (3) evaluation and measurements of brine geochemistry from U.S. O&G produced waters has informed a high Li recovery potential from Marcellus Shale produced water. NETL researchers have been developing tailored pre-treatment processes, an innovative and highly durable lithium sorbent, and geochemical model guided precipitation to accelerate Li production from the Marcellus Shale produced waters. These innovations driven by characterization are integral for maximizing and advancing the potential for CM recovery while offsetting the cost and environmental footprint for FE waste management.

critical mineral processing↗
NASA NTRS1973

Mariner Mars 1971 battery design, test, and flight performance

The design, integration, fabrication, test results, and flight performance of the battery system for the Mariner Mars spacecraft launched in May 1971 are presented. The battery consists of 26 20-Ah hermetically sealed nickel-cadmium cells housed in a machined magnesium chassis. The battery package weighs 29.5 kg and is unique in that the chassis also serves as part of the spacecraft structure. Active thermal control is accomplished by louvers mounted to the battery baseplate. Battery charge is accomplished by C/10 and C/30 constant current chargers. The switch from the high-rate to low-rate charge is automatic, based on terminal voltage. Additional control is possible by ground command or onboard computer. The performance data from the flight battery is compared to the data from various battery tests in the laboratory. Flight battery data was predictable based on ground test data.

Bogner, R. S.↗
DOE OSTI2024

Utilities Perspective on Protection Challenges with High IBR Penetration

Utilities have seen rapid increase in Solar, Wind, and Battery Energy Storage Resources that interconnect to their electric system through Inverters. Inverter-Based Resources (IBRs) have fault current characteristics that are unlike the fault current response of traditional rotating-machine-based generators, which is well known and repeatable. IBR’s non-traditional fault current behavior is due to the IBR control scheme, which is configured to provide a clean AC output but also protect the inverter’s sensitive power electronics devices from damage, one source of which is overcurrent. This results in low fault current magnitude, low or no negative sequence current injection, the variability of sequence component currents, the variability of voltage with respect to current angles, and the lack of inertia. The control scheme also results in a fault current response that can vary between manufactures and between models of the same manufacturer. High penetration of IBRs can adversely affect the protection schemes applied in areas with high penetration of IBRs. With the proliferation of IBRs, utilities are finding out that conventional protection schemes are not adequately equipped to protect the electric systems. This is mainly because the existing protection elements and practices have been designed based on the fault current response of conventional rotating machines. In several cases, the available literature does not provide any clear solution for the issues when the protection scheme does not operate properly near IBRs. This presentation identifies various protection challenges due to IBRs that industry is facing, from the utility perspective. Instead of facing on one issue, we are looking broadly on all the challenges that system protection has experienced with high penetration of IBRs. Based on the IBR response from various utilities during real fault events and gathering perspective from different utility SMEs via questionnaire, the presentation summarizes on gathered data and internal experiences.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRSDate not supplied

A Machine Learning Framework for Error Compensation in Radiative Transfer Calculations

Radiative heat transfer influences the amount of heat flux transferred to the surface of the hypersonic vehicle, which is essential to evaluate the performance of thermal protection systems. The radiative heat flux is found to be computationally prohibitive while accounting for the variation in spatial, angular, and spectral domains. A new methodology has been recently developed to alleviate the cost of computation in the spectral domain by constructing flow-agnostic reduced-order models (ROMs). The developed spectral ROM databases provide grouping strategies that account for non-equilibrium absorption and emission as well as interaction between disparate species due to spectral overlap in associated radiative processes. However, the developed ROMs need to be optimized for a specific combination of interacting gas species and would need to re-calibrated in case individual species are added/omitted. In this work, we use various machine learning (ML) techniques to approximate the radiative intensities determined by a ROM optimized for a specific gas mixture. The ML model relies on the ROM databases developed for a single species which ignores any spectral overlap. Thus, radiation evaluation starts with a simple summation of radiative intensities predicted using these non-calibrated ROMs for the contributing species. The ML framework then provides a correction to account for the interplay in the frequency, i.e., emission of photons by one species and absorption by another, and yields mixture-specific radiation fields. Once trained on the individual ROM databases, the ML framework offers instantaneous corrections that serves as a time/cost effective alternative to the optimization of ROMs for a specific gas mixture. The ML framework is trained on both the high fidelity and ROM evaluated line of sight (LOS) data from Orion, Stardust, and FIRE II cases to obtain a general purpose correction model for earth re-entry scenarios when radiation contributions from both atomic nitrogen and atomic oxygen are considered. A geometric length scale parameter is used in the training process to account for errors introduced in the ROM databases as a consequence of high optical thickness. The efficacy of the ML framework is underscored through extensive analysis of train and test errors with respect to all the re-entry scenarios. The applicability of such an ML framework was further corroborated by embedding it in a state-of-the-art US3D - NERO system for determining the radiative heat flux transferred to the hypersonic vehicle surface.

Radiation↗
NASA NTRS1990

Laboratory Investigations: Low Earth Orbit Environment Chemistry With Spacecraft Surfaces

Long-term space operations that require exposure of material to the low earth orbit (LEO) environment must take into account the effects of this highly oxidative atmosphere on material properties and the possible contamination of the spacecraft surroundings. Ground-based laboratory experiments at Los Alamos using a newly developed hyperthermal atomic oxygen (AO) source have shown that not only are hydrocarbon based materials effected but that inorganic materials such as MoS2 are also oxidized and that thin protective coatings such as Al2O3 can be breached, producing oxidation of the underlying substrate material. Gas-phase reaction products, such as SO2 from oxidation of MoS2 and CO and CO2 from hydrocarbon materials, have been detected and have consequences in terms of spacecraft contamination. Energy loss through gas-surface collisions causing spacecraft drag has been measured for a few select surfaces and has been found to be highly dependent on the surface reactivity.

Jon B Cross↗
DOE OSTI2025

Review on Preprocessing Strategies, Deactivation, Thermal Safety, and Future Perspectives in Lithium-Ion Battery Recycling

The rapid growth in the use of lithium-ion batteries (LIBs) in electric vehicles, consumer electronics, and renewable energy storage has made effective end-of-life management essential. Recycling LIBs is critical not only for resource recovery and environmental protection but also for ensuring safety and economic viability. This review focuses on the preprocessing technologies that precede typical recovery processes, including disassembly, sorting, discharging, electrolyte removal, dismantling, thermal treatment, separation, and flotation. These steps play a foundational role in determining the efficiency, safety, and environmental impact of LIB recycling. LIBs pose substantial fire and explosion risks due to residual charge, flammable electrolytes, and reactive materials. The conditions and successive progression of the exothermic reactions which lead to thermal runaway has been discussed. It also explores secure deactivation techniques such as external circuit discharge, saline immersion, and thermomechanical methods, alongside fire prevention strategies including the use of flame retardants, elimination of oxidants, and reduction of heat generation and accumulation. Challenges and future directions are outlined, highlighting the need for standardized designs, automation, and safer, more sustainable recycling infrastructure. Furthermore, this review is distinguished by its focused analysis of preprocessing and deactivation steps, with particular attention to the thermal safety engineering aspects of LIB recycling.

Battery deactivation↗
NASA NTRS2004

Ejecta Generation and Redistribution on 433 Eros: Modeling Ejecta Launch Conditions

The NEAR-Shoemaker mission to asteroid 433 Eros presents an unprecedented opportunity to gain fundamental new knowledge about the processes governing regolith formation and redistribution on small bodies. NEAR-Shoemaker’s high-resolution imaging of the surface of Eros makes the asteroid a valuable and heretofore unparalleled laboratory for the detailed study of impact ejecta reaccretion and regolith redistribution on low-gravity (of order 10 -3 g) objects. Regolith is produced on asteroids by impact cratering, and the existence of regolith on the smallest solar system bodies supports the view that some of the ejecta from impact events on such objects may be retained. Impact craters and retained ejecta on low-gravity objects like Eros represent valuable natural laboratories for evaluating various models of impact cratering processes, since they may present crater structures or ejecta features that either do not form or are hidden on higher-gravity bodies like the Moon. Further, quantifying the extent to which impact processes generate and redistribute regoliths on small body surfaces (excavation depths, retained fraction, turnover timescales, etc.) is pivotal to the issue of how to relate meteoritical samples to their asteroidal parent bodies when surficial processes ( i.e., “space weathering”) may disguise or cover up underlying material and confound the ability of remote sensing techniques to provide reliable mineralogical assays of the parent objects. The rich variety of data on Eros’ regolith properties and distribution returned by NEAR-Shoemaker now require detailed analysis in order to take full advantage of the clues these observations offer for elucidating details of the impact cratering process on small bodies. Complicating simple interpretations of crater and ejecta morphology are dynamical effects on ejecta emplacement resulting from Eros’ irregular shape, rapid (5.27 hr) rotation, and low gravity. Figure 1 shows the very different ejecta deposit morphology that can result if the effects of rotation alone are neglected. Considering the additional complicating factors of Eros’ irregular shape and complex gravitational field, simple calculations of the extent and thickness of ejecta blankets and the spatial distribution of ejecta blocks from basic crater scaling laws or numerical hydrocodes alone do not suffice. In order to fully interpret the suite of NEAR-Shoemaker observations of regolith features across the surface of Eros and to evaluate various impact models for specific craters on the asteroid, detailed dynamical modeling of the deposition of crater ejecta from those craters is required . Here, I describe some modifications and improvements to the dynamical model being used for these studies.

D D Durda↗
DOE OSTI2026

Extreme confinement of hydrogen gas within fullerenelike nanoporous carbon

Nanoporous carbons and carbon nanostructures can store hydrogen at cryogenic temperatures but lack the volumetric and gravimetric capacity to be industrially significant. Recent inelastic neutron scattering experiments suggest a highly dense phase of hydrogen at temperatures well above the melting point of solid hydrogen. However, it remains unclear how pore geometry and intermolecular interactions enable these dense phases to exist, with dispersion (van der Waals) or electrostatic/induction suggested to be the key effects in slit and curved pores but their relative contributions have yet to be quantified. In this paper, we perform benchmark electronic structure calculations allowing the interactions between planar and curved aromatic molecules with hydrogen to be accurately determined. Dispersion was found to dominate over electrostatic and inductive effects with some many-body charge transfer (Dobson type-A) effects needed to capture the most highly curved structures. Density functional methods that include type-A many-body effects were found to accurately describe the intermolecular interactions at a fraction of the cost of coupled-cluster simulations and these approaches were used to calculate the energies inside large carbon bowl and slit pores. The interaction energies inside the bowl pores were found to depend on the orientation of the hydrogen molecule. This rotational barrier, modeled as a quantum hindered rotor, could reproduce the peak splitting observed in inelastic neutron scattering experiments, with weak splitting arising from bowl-like fullerene pores and strong splitting from highly confining nanotubelike pores. Increasing the fraction of such curved pores in nanoporous carbons may therefore offer a pathway to enhance their hydrogen-storage capacity. Moreover, the preferential adsorption of ortho hydrogen on nanotubelike pores could enable the storage of high-density hydrogen without the need to remove heat produced during the ortho-para hydrogen conversion.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Fast Active-Set Thresholding Method for Nonnegative Least Squares

Nonnegative Least Squares (NNLS) is a fundamental constrained optimization problem encountered in many applications such as image deblurring, signal processing, nonnegative matrix factorization, magnetic microscopy, and hyperspectral imaging. Active-set based methods are a common class of algorithms for solving NNLS which identify the optimal variable set of the NNLS solution. They do so by iteratively solving a series of unconstrained least squares problems, identifying which variables violate the nonnegativity constraints, and then swapping variables in/out of consideration until the optimal set of variables is found. Several variations improving upon this method exist in the literature. In this work, we propose an active-set swap heuristic which further improves upon existing active-set based methods for NNLS. Our optimizations are based upon adding multiple variables to the passive set within a threshold of the smallest gradient value and removing variables within a similar threshold of the closest boundary constraint. We leverage these optimizations to yield a Fast Active-Set Thresholding NNLS (FAST-NNLS) algorithm which significantly outperforms the existing state-of-the-art NNLS algorithms for a wide range of problems. Rigorous convergence guarantees are proven for the proposed method. We demonstrate the effectiveness of our proposed method on multiple synthetic datasets and two realworld text analysis applications. In doing so, we present the most comprehensive NNLS solver comparison in the literature to date.

Cobb, Benjamin [Georgia Institute of Technology]↗