Search NASA⌕ Search

SEARCH · Search NASA

Results for “Inverse Optimization”

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 127 records · Page 7

Understanding the Origin of Negative Temperature Dependence and Activity of N-Coordinated Cobalt Sites During Ethylene Dimerization

The on-demand production of short-chain linear alpha olefins (LAOs; C4-C8) via C2H4 dimerization and oligomerization is industrially attractive, prompting extensive research on designing active, selective, and stable catalysts for industrial use. Cobalt supported on ammoniated carbon (Co(NH3)x/C) catalysts have shown remarkable activity and selectivity in this process. However, critical aspects such as the active phase, active site structure, the role of the catalyst support, cobalt loading effects, and the inverse correlation of the reaction rate with temperature remain inadequately understood. This study systematically explores these factors using a combination of steady-state differential catalytic tests, in situ molecular characterization including diffuse reflectance UV-Vis (DR-UV-Vis), Infrared, and Raman spectroscopies, and ex situ X-ray diffraction (XRD) and high annular aberration-corrected dark field transmission electron microscopy (HAADF-STEM). Various supports (SiO2, Al2O3, NH4-ZSM-5, g-C3N4, and C) and cobalt loadings (1.0-3.0 Co nm-2) were studied to determine the optimal catalyst composition and identify the active phase and sites. Carbon-supported catalysts uniquely produce C4-8 LAOs during C2H4 dimerization, with site-time-yield remaining constant (~10-3 s-1) for 1.0-4.0 Co nm-2 at prolonged reaction times (24-48?h time-on-stream). At higher loadings of 6.0 Co nm-2, the formation of crystalline CoO and Co3O4 phases reduces catalytic activity and LAO selectivity. Our findings show that active catalysts lack crystalline cobalt oxides and instead feature dispersed Co2+ sites, tetra-coordinated to a mix of N/NH3 and O/H2O ligands, which catalyze C2H4 dimerization via the Cossee-Arlman mechanism, exhibiting 1st order dependence on C2H4 concentration. The observed inverse rate-temperature correlation is attributed to compensation effects (i.e., presence of Cremer-Constable relationship) linked to changes in adsorption enthalpic and entropic factors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Optimal sensing on an asymmetric exceptional surface

We study the connection between exceptional points (EPs) and optimal parameter estimation, in a simple system consisting of two counterpropagating traveling wave modes in a microring resonator. The unknown parameter to be estimated is the strength of a perturbing cross-coupling between the two modes. Partially reflecting the output of one mode into the other creates a non-Hermitian Hamiltonian that exhibits a family of EPs, creating an exceptional surface (ES). We use a fully quantum treatment of field inputs and noise sources to obtain a quantitative bound on the estimation error by calculating the quantum Fisher information (QFI) in the output fields, whose inverse gives the Cramér-Rao lower bound on the mean-squared error of any unbiased estimator. We determine the bounds for two input states, namely, a semiclassical coherent state and a highly nonclassical NOON state. We find that the QFI is enhanced in the presence of an EP for both of these input states and that both states can saturate the Cramér-Rao bound. We then identify idealized yet experimentally feasible measurements that achieve the minimum bound for these two input states. We also investigate how the QFI changes for parameter values that do not lie on the ES, finding that these can have a larger QFI, suggesting alternative routes to optimize the parameter estimation for this problem.

Exceptional points↗

Inducing a tunable skyrmion-antiskyrmion system through ion beam modification of FeGe films

Abstract Skyrmions and antiskyrmions are nanoscale swirling textures of magnetic moments formed by chiral interactions between atomic spins in magnetic noncentrosymmetric materials and multilayer films with broken inversion symmetry. These quasiparticles are of interest for use as information carriers in next-generation, low-energy spintronic applications. To develop skyrmion-based memory and logic, we must understand skyrmion-defect interactions with two main goals—determining how skyrmions navigate intrinsic material defects and determining how to engineer disorder for optimal device operation. Here, we introduce a tunable means of creating a skyrmion-antiskyrmion system by engineering the disorder landscape in FeGe using ion irradiation. Specifically, we irradiate epitaxial B20-phase FeGe films with 2.8 MeV Au 4+ ions at varying fluences, inducing amorphous regions within the crystalline matrix. Using low-temperature electrical transport and magnetization measurements, we observe a strong topological Hall effect with a double-peak feature that serves as a signature of skyrmions and antiskyrmions. These results are a step towards the development of information storage devices that use skyrmions and antiskyrmions as storage bits, and our system may serve as a testbed for theoretically predicted phenomena in skyrmion-antiskyrmion crystals.

74 ATOMIC AND MOLECULAR PHYSICS↗

MTUQ: a framework for estimating moment tensors, point forces, and their uncertainties

SUMMARY We introduce MTUQ, an open-source Python package for seismic source estimation and uncertainty quantification, emphasizing flexibility and operational scalability. MTUQ provides MPI-parallelized grid search and global optimization capabilities, compatibility with 1-D and 3-D Green’s function database formats, customizable data processing, C-accelerated waveform and first-motion polarity misfit functions, and utilities for plotting seismic waveforms and visualizing misfit and likelihood surfaces. Applicability to a range of full- and constrained-moment tensor, point force, and centroid inversion problems is possible via a documented application programming interface, accompanied by example scripts and integration tests. We demonstrate the software using three different types of seismic events: (1) a 2009 intraslab earthquake near Anchorage, Alaska; (2) an episode of the 2021 Barry Arm landslide in Alaska; and (3) the 2017 Democratic People’s Republic of Korea underground nuclear test. With these events, we illustrate the well-known complementary character of body waves, surface waves, and polarities for constraining source parameters. We also convey the distinct misfit patterns that arise from each individual data type, the importance of uncertainty quantification for detecting multimodal or otherwise poorly constrained solutions, and the software’s flexible, modular design.

58 GEOSCIENCES↗

Photoredox Catalysis with Spin Magnetic Field Effects: A Scheme for Chiral Resolution

Quantities that break both mirror symmetry and time-reversal symmetry, such as the orbital angular momentum, are known to connect molecular chirality with an applied magnetic field. This concept has led to observations such as magneto-chiral dichroism and chirality-induced spin selectivity (CISS). However, being small, these effects often require additional amplification procedures such as flow chemistry to achieve bulk enantioseparation. In this work, we demonstrate how the magnetic field effect on photogenerated radical pairs, which also breaks time-reversal symmetry, can be harnessed for enantiopurification. Fundamental to this process is the collective decay of the singlet and triplet radical-pair states made possible by an applied magnetic field. Because opposite enantiomers exhibit spin-orbit coupling matrix elements of opposite signs, the singlet and triplet decay channels interfere constructively in one enantiomer. Meanwhile, molecules of the other enantiomer are funnelled into the first enantiomer through excited-state chirality inversion, achieving (dynamic kinetic) chiral resolution. Using an axially chiral binaphthyl derivative and a borane photosensitizer as prototype, we predict an appreciable enantiomeric excess (e.e.) of 90% to be possible at steady state, attained within hundreds of milliseconds when irradiated by a laser. Importantly, our analytical results showcase regimes of perfect enantioselectivity (100% e.e.), accessible by further chemical optimization of the photosensitizer for which general strategies are discussed. Altogether, this work illustrates a so-far untapped but powerful control knob for photoredox catalysis based on spin chemistry principles.

Magnetic properties↗

Conjugation effect of amine molecules in non-aqueous Mg redox flow batteries

Nonaqueous magnesium redox flow batteries (Mg RFBs) are attractive for low-cost, high-energy-density and long-cycle-life stationary energy storage applications. However, state-of-the-art cathode redox-active molecules suffer from low solubility and low redox potential. Herein, we screened a range of cathode redox-active molecules and identified amine molecules as optimal to couple with the Mg anode. The properties of amine derivatives and their performances were collected to establish the correlation between molecular structures and electrochemical performances. The redox potential and solubility of these amine molecules were influenced by the π-conjugated and non-conjugated structures of amine derivatives. Density functional theory (DFT) simulations and the inverse aromatic fluctuation index (FLU −1 ) verified that conjugation had an important role in stabilizing the molecule and increasing its redox potential. Notably, tris[4-(diethylamino)phenyl]amine (TDPA) achieved the highest theoretical energy density (∼120 Wh L −1 ) due to its high solubility (∼0.9 M) and voltage (∼2.5 V vs. Mg/Mg 2+ ). We also demonstrated that ether solvents were crucial for stable, high-solubility catholytes, while bulk anions did not affect the redox potential of these p-type molecules. In a Mg-amine RFB configuration, the battery delivered 2.50 V, a specific discharge capacity of 106.5 mAh g −1 , an initial coulombic efficiency of 90.74%, and a capacity retention of 93.88% after 150 cycles.

Qin, Yunan [University of Utah, Salt Lake City, UT↗

High-throughput spin-bath characterization of spin defects in semiconductors

Detailed knowledge of the local environments of spin defects in semiconductors, such as nitrogenvacancy (NV) centers in diamond or divacancies in silicon carbide, is crucial for optimizing control and entanglement protocols in quantum sensing and information applications. However, at present a direct experimental characterization of individual defect environments is not scalable, as conventional spin-bath measurements are time consuming and difficult to automate. Achieving high-throughput characterization requires short experiments to probe the spin bath. However, with fewer and noisier measurements, the inverse problem of recovering spin-bath properties from measured data becomes ill posed, with multiple spin baths having a high likelihood of yielding the same data. In this work, we present a set of computational tools to resolve the ill-posed inverse problem of recovering the atomic positions and hyperfine couplings of random nuclei surrounding spin defects from sparse, noisy experimental coherence data, which can be obtained in hours. Here, we use a trans-dimensional Bayesian approach that incorporates ab initio data to yield full posterior distributions over nuclear spin environments, enabling robust recovery from limited data. We also provide practical tools and guidelines to determine the limits of detectability for hyperfine couplings under specific dynamical decoupling sequences and sampling conditions. In addition, we demonstrate how the tools developed here, in combination with ab initio simulations of spin baths, can guide the design of efficient experimental protocols for application-specific high-throughput screening. To showcase the utility of our approach, we apply it to design fast dynamical decoupling experiments to characterize the spin baths often individual NV centers in diamond. While the primary focus is on accelerating spin-bath characterization of spin defects, this Bayesian approach also lays the foundation for digital-twin studies of spin defects, where a virtual model of the spin-defect system evolves in real time with ongoing experimental measurements. Together, the set of tools we designed and applied paves the way for scalable deployment of spin defects in semiconductors for quantum sensing and information applications.

Bayesian methods↗

Nondestructive Evaluation of Concrete: Elastic Property Imaging Through Full-Waveform Inversion

Concrete is a vital material in construction—especially in the nuclear industry, where it is used in critical structures such as containment vessels. Over time, concrete can degrade due to harsh operational and environmental conditions, necessitating that its elastic properties be accurately evaluated to ensure structural integrity and safety. Traditional nondestructive evaluation methods such as ultrasound-based techniques often rely on simplifying assumptions that may not hold true for concrete. This paper presents an advanced ultrasound-based method that uses elastic full-waveform inversion (EFWI) to create detailed images of concrete’s mechanical properties. By accurately modeling wave behaviors such as scattering and reflection, we aim to overcome the limitations of conventional ultrasonic-based methods. In this work, the imaging problem involved reconstructing the various elastic properties of a heterogenous concrete block with three steel rebars embedded in it. The ultrasonic measurements were synthetically generated from multiple sources and receivers, and the reconstruction process was performed using a gradient-based optimization algorithm. Our approach leveraged EFWI to reconstruct high-resolution images of the pressure wave speed, shear wave speed, and density. Multiple misfit functions—including L2-norm, cross-correlation (CC), and L1-norm—combined with total variation (TV) regularization and parameter constraints using a Sigmoid function—were explored for the reconstruction. The results demonstrated that using the L1-norm misfit function in conjunction with TV regularization and Sigmoid constraints significantly improved the reconstruction quality in comparison to traditional methods. This approach provided clearer images with fewer artifacts and better captured background heterogeneity. Our findings highlight that, when properly designed, EFWI carries great potential for providing comprehensive, more accurate, and more reliable assessments of concrete conditions, as is crucial for the maintenance and safety of nuclear power plant structures.

97 - MATHEMATICS AND COMPUTING↗

Determining the Solubility Behavior of Kogarkoite in Simulated Nuclear Waste

Kogarkoite (Na 3 FSO 4 ) is a sparingly soluble fluoride–sulfate double salt that has been identified in high level nuclear waste sludge at the Hanford Site and, more recently, in sludge batch compilation samples at the Savannah River Site (SRS). Due to its complex dissolution behavior, which exhibits an inverse dependence on sodium ion activity, the presence of this mineral poses significant challenges to waste retrieval and processing. Incomplete dissolution during sludge washing can lead to the retention of fluoride and sulfate in the high-level waste feed, potentially causing the formation of corrosive, immiscible molten salt layers, known as "glass gall,” in vitrification melters. Current efforts to optimize flowsheet parameters and wash-water volumes are hindered by the absence of a commercially available, certified reference material, which prevents the accurate calibration of analytical methods and the verification of dissolution kinetics. To address this critical gap, this research focuses on the laboratory synthesis of pure Kogarkoite to serve as a standard for comprehensive solubility and washing performance testing. A coupled synthesis and simulant campaign was executed using an evaporative crystallization protocol designed to replicate the dynamic concentration effects observed in tank farm operations. Thirteen simulant matrices were prepared by dissolving systematically varied ratios of sodium fluoride (NaF) and sodium sulfate (Na 2 SO 4 ) in deionized water under three distinct caustic regimes: 0.0 g (control), 4.0 g (~1 M), and 12.0 g (~3 M) sodium hydroxide (NaOH). While thermodynamic equilibrium models suggest that high-caustic environments should favor the stability of the double salt7, results from this evaporative study at 25 0 C revealed a distinct kinetic divergence. Simulants with high hydroxide loading predominantly yielded large, blocky crystals of sodium sulfate decahydrate (Na 2 SO 4 .10H 2 O). Successful synthesis of pure Kogarkoite was achieved exclusively in specific NaOH-free compositional windows, where the precipitate manifested as fine, opaque granular aggregates. Ion chromatography (IC) analysis confirmed phase purity through the simultaneous stoichiometric depletion of both fluoride and sulfate from the supernatant. This successful synthesis establishes a reproducible route to generate bulk Kogarkoite, enabling the subsequent phase of quantitative dissolution testing using inhibited water to optimize sludge-batch assembly.

Sarker, Md Sharif [Florida International Univ. (FI↗

Position-Enhanced Gradient Attack (PEGA) on Medical Language Models

Federated Learning (FL) enables collaborative training of language models on sensitive clinical notes without sharing the data. However, this paradigm is vulnerable to gradient inversion attacks that can reconstruct private data from shared gradients. We find that state-of-the-art attacks are less effective in the medical domain, failing to overcome the unique challenges posed by its specialized vocabulary and unstructured format. To address this, we introduce the Position-Enhanced Gradient Attack (PEGA), a novel attack that makes gradients position-aware by optimizing token and position embeddings simultaneously. PEGA employs two key innovations: a periodic sorting of positional embeddings to resolve token order ambiguity and a late-stage embedding replacement strategy to correct hard-to-recover critical tokens. To evaluate the leakage of sensitive data more directly, we also propose the Unified PHI-Recall (UPHI), a new metric measuring the recovery of Protected Health Information. Experiments on the MIMIC-III dataset show that PEGA significantly outperforms leading attacks like TAG and LAMP, particularly in its ability to reconstruct identifiable patient information, exposing a more severe and nuanced privacy risk in federated medical NLP.

Xu, Nuo [University of Minnesota]↗

On the engulfment of antifreeze proteins by ice

Antifreeze proteins (AFPs) are remarkable biomolecules that suppress ice formation at trace concentrations. To inhibit ice growth, AFPs must not only bind to ice crystals, but also resist engulfment by ice. The highest supercooling, ΔT*, for which AFPs are able to resist engulfment is widely believed to scale as the inverse of the separation, L, between bound AFPs, whereas its dependence on the molecular characteristics of the AFP remains poorly understood. By using specialized molecular simulations and interfacial thermodynamics, here, we show that in contrast with conventional wisdom, ΔT* scales as L –2 and not as L –1 . We further show that ΔT* is proportional to AFP size and that diverse naturally occurring AFPs are optimal at resisting engulfment by ice. By facilitating the development of AFP structure–function relationships, we hope that our findings will pave the way for the rational design of AFPs.

Science & Technology - Other Topics↗

Modeling hydrogen markets: Energy system model development status and decarbonization scenario results

Hydrogen can be used as an energy carrier and chemical feedstock to reduce greenhouse gas emissions, especially in difficult-to-decarbonize markets such as medium- and heavy-duty vehicles, aviation and maritime, iron and steel, and the production of fuels and chemicals. Significant literature has been accumulated on engineering-based assessments of various hydrogen technologies, and real-world projects are validating technology performance at larger scales and for low-carbon supply chains. While energy system models continue to be updated to track this progress, many are currently limited in their representation of hydrogen, and as a group they tend to generate highly variable results under decarbonization constraints. Here, the present work provides insights into the development status and decarbonization scenario results of 15 energy system models participating in study 37 of the Stanford Energy Modeling Forum (EMF37), focusing on the U.S. energy system. The models and scenario results vary widely in multiple respects: hydrogen technology representation, scope and type of hydrogen end-use markets, relative optimism of hydrogen technology input assumptions, and market uptake results reported for 2050 under various decarbonization assumptions. Most models report hydrogen market uptake increasing with decarbonization constraints, though some models report high carbon prices being required to achieve these increases and some find hydrogen does not compete well when assuming optimistic assumptions for all advanced decarbonization technologies. Across various scenarios, hydrogen market success tends to have an inverse relationship to success with direct air capture (DAC) and carbon capture and storage (CCS) technologies. While most model-scenario combinations predict modest hydrogen uptake by 2050 – <10 million metric tons (MMT) – aggregating the top 10 % of market uptake results across sectors suggests an upper range demand potential of 42–223 MMT. The high degree of variability across both modeling methods and market uptake results suggests that increased harmonization of both input assumptions and subsector competition scope would lead to more consistent results across energy system models. The wide variability in results indicates strongly divergent conclusions on the role of hydrogen in a decarbonized energy future.

08 HYDROGEN↗

Adjustable collimation for dark-field proton radiography contrast enhancement

A dark-field imaging condition established through pre- and post-object beam collimation has been tested for a variety of different settings at the proton radiography facility (pRad) at Los Alamos National Laboratory. This technique facilitates imaging of thin samples and materials with small areal density changes with 800 MeV protons that would otherwise appear nearly transparent for high energy protons. In the configuration that was examined, the pre-object (upstream) collimator is utilized to remove protons with large scattering angles prior to the target object. The post-object collimation entails a combination of a conventional collimator, which eliminates protons with large scattering angles after target interaction, and an inverse collimator. The inverse collimator removes unscattered protons traveling on the beam axis and those with smaller scattering angles up to a desired angle. This establishes a dark-field condition by permitting solely protons that exhibit scattering angles from target interaction within a designated range to contribute to the image. This study explores the potential of combining different sized collimators for pre- and post-object collimation. The objective is to optimize contrast for a range of samples with minimal areal density variations. A series of dark-field conditions are explored employing a remote-controlled collimator wheel and the contrast-to-noise ratio, spatial resolution, and proton transmission are compared for those combinations. It is found that the contrast-to-noise ratio can be substantially enhanced for thin samples and low areal density objects, thereby facilitating capturing high-contrast images of objects that were previously invisible to high energy protons.

43 PARTICLE ACCELERATORS↗

Evaluating the feasibility of using downwind methods to quantify point source oil and gas emissions using continuously monitoring fence-line sensors

The dependable reporting of methane (CH 4 ) emissions from point sources, such as fugitive leaks from oil and gas infrastructure, is important for profit maximization (retaining more hydrocarbons), evaluating climate impacts, assessing CH 4 fees for regulatory programs, and validating CH 4 intensity in differentiated gas programs. Currently, there are disagreements between emissions reported by different quantification techniques for the same sources. It has been suggested that downwind CH 4 quantification methods using CH 4 measurements on the fence line of production facilities could be used to generate emission estimates from oil and gas operations at the site level, but it is currently unclear how accurate the quantified emissions are. To investigate the accuracy of downwind methods, this study uses fence-line simulated data collected during controlled-release experiments as input for a non-standard closed-path eddy covariance (EC), the Gaussian plume inverse model (GPIM), and the backward Lagrangian stochastic (bLs) model in a range of atmospheric conditions. This study's EC attempt was unsuccessful due to data collection and instrumentation issues, resulting in invalid results characterized by underestimated emissions, large negative fluxes, and cospectra/ogives that deviated from their ideal shapes. Consequently, the EC results could not be compared with the GPIM and bLS model. The bLs model demonstrated the highest accuracy for single-release single-point emissions, though it exhibited greater uncertainty than GPIM under multi-release conditions. Across the GPIM and bLs model, the most reliable quantification was achieved with 15 min averaging and a narrow 5° wind sector range. Although EC was limited in this context, future studies should consider employing a standard EC system and further optimizing GPIM and bLs approaches – particularly for complex multi-source scenarios – to enhance quantification accuracy and reduce uncertainty.

03 NATURAL GAS↗

Seismic Event Characterization Using Full Moment Tensors on the Hypersphere

Moment tensor solutions provide insights into the deformation that has occurred in the source region of a seismic event and are therefore of great value in identifying different types of seismic sources, such as when monitoring for underground nuclear tests. Despite this utility, inversion of waveforms recorded by seismometers for their full seismic moment tensor is not yet routine, and development of robust methods to classify events based on this information is in its infancy. Here, we assemble an inventory of 1405 full moment tensor solutions that include explosive, earthquake, and collapse events, and investigate the use of anisotropic probability distribution functions on the 5D hypersphere to discriminate between these sources. Using a Bayesian classifier, we obtain optimal success rates of 98.4% across all events and demonstrate that modification of the prior probabilities provides a natural way to alter the balance between not missing desirable events (such as explosions) versus misclassifying large numbers of undesired events (such as earthquakes). The approach is specifically designed to progress from traditional, bipolar event screening metrics to more generalized event identification across multiple types of seismic sources. Despite current databases containing insufficient numbers of events to definitively demonstrate at present, we also find intriguing evidence of subgroupings within individual source populations on the hypersphere, for example, between chemical and nuclear explosions, raising the potential possibility of discriminating between these event types in the future.

Geosciences↗

Inverse design of a pyrochlore lattice of DNA origami through model-driven experiments

Sophisticated statistical mechanics approaches and human intuition have demonstrated the possibility of self-assembling complex lattices or finite-size constructs. However, attempts so far have mostly only been successful in silico and often fail in experiment because of unpredicted traps associated with kinetic slowing down (gelation, glass transition) and competing ordered structures. Theoretical predictions also face the difficulty of encoding the desired interparticle interaction potential with the experimentally available nano- and micrometer-sized particles. To overcome these issues, we combine SAT assembly (a patchy-particle interaction design algorithm based on constrained optimization) with coarse-grained simulations of DNA nanotechnology to experimentally realize trap-free self-assembly pathways. In this paper, we use this approach to assemble a pyrochlore three-dimensional lattice, coveted for its promise in the construction of optical metamaterials, and characterize it with small-angle x-ray scattering and scanning electron microscopy visualization.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Electron Cyclotron Emission–Based Separatrix Identification in ITER with OMFIT Synthetic Modeling

Accurate determination of the separatrix location is essential for understanding edge plasma behavior and optimizing confinement in tokamaks, especially in next-generation devices such as ITER. In this study, a synthetic microwave diagnostics module was developed and implemented in the OMFIT framework to assess the feasibility of an electron cyclotron emission–based separatrix detection method in ITER plasmas. Simulations were carried out using ITER H-mode equilibrium scenarios with different plasma density profiles and different pedestal widths. Here, the results show that the electron emission temperature profiles consistently exhibit an inversion pattern near the edge, with a well-defined minimum point that could serve as a proxy for the separatrix location. However, unlike in the DIII-D, the minimum point in ITER is systematically offset by approximately 2 cm into the scrape-off layer, independent of density or pedestal width, which is within the radial resolution range (2 to 5 cm) determined by the 500-MHz channel spacing. While the method does not provide the exact separatrix location, it offers a reliable indicator of the boundary region and has potential applications for real-time boundary monitoring in ITER and other future fusion devices.

Electron cyclotron emission↗

Sparse measurement medical CT reconstruction using multi-fused block matching denoising priors

A major challenge for medical X-ray CT imaging is reducing the number of X-ray projections to lower radiation dosage and reduce scan times without compromising image quality. However these under-determined inverse imaging problems rely on the formulation of an expressive prior model to constrain the solution space while remaining computationally tractable. Traditional analytical reconstruction methods like Filtered Back Projection (FBP) often fail with sparse measurements, producing artifacts due to their reliance on the Shannon-Nyquist Sampling Theorem. Consensus Equilibrium, which is a generalization of Plug and Play, is a recent advancement in Model-Based Iterative Reconstruction (MBIR), has facilitated the use of multiple denoisers are prior models in an optimization free framework to capture complex, non-linear prior information. However, 3D prior modelling in a Plug and Play approach for volumetric image reconstruction requires long processing time due to high computing requirement. Instead of directly using a 3D prior, this work proposes a BM3D Multi Slice Fusion (BM3D-MSF) prior that uses multiple 2D image denoisers fused to act as a fully 3D prior model in Plug and Play reconstruction approach. Our approach does not require training and are thus able to circumvent ethical issues related with patient training data and are readily deployable in varying noise and measurement sparsity levels. In addition, reconstruction with the BM3D-MSF prior achieves similar reconstruction image quality as fully 3D image priors, but with significantly reduced computational complexity. We test our method on clinical CT data and demonstrate that our approach improves reconstructed image quality.

Hossain, Maliha [ORNL]↗