Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time shift”

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 613 records · Page 34

Real-time sensor measurements to evaluate impacts of load cycling on high-temperature fire-side corrosion in a PC boiler

The main objective of this research program is to design, manufacture, and demonstrate a miniaturized, multi-process, monitoring system (mMPMS) for boiler condition management and easy system deployment to obtain a higher spatial resolution. This system will facilitate a Condition-Based Maintenance (CBM) philosophy that actively monitors the health of assets to predict and prevent failures and maximize availability and generating capacity at a reduced cost. CBM systems can provide boiler data that the advanced process control (APC) system can utilize for plant performance optimization, which is increasingly relevant as coal power plants shift from predominantly base-load operation to predominantly transient operation involving large load swings. The mMPMS is based on an electrochemical sensor that provides a real-time indication of the risk of damage to key locations in the radiant or convective section of a coal-fired boiler such as metal loss rates, heat flux, metal surface temperature, and deposit thickness. These indications can be utilized to optimize boiler performance as well as improve boiler availability in conjunction with corresponding operating conditions. This monitoring system was developed and tested in the high-temperature regions of coal-fired utility boilers in this project but can be applied to many other industries and applications as well. This project leveraged the existing electrochemical noise-based monitoring system and the new sensor design is small enough to be installed through the webbing of the waterwalls without the need for long shut-downs to bend tubes and to make it feasible to obtain high spatial resolution in the boiler. Data is transferred to the plant distributed control system (DCS) and any other control system. The sensor body that houses the sensor assembly was designed to ensure good conductive contact with boiler tubes to ensure the sensor is held at an identical temperature to the tube surface temperature. The data acquisition and signal conditioning modules were redesigned into a small footprint with optimized cooling of the module. System software was developed specifically for the new signal conditioning module and is compatible with plant PLCs. After the preliminary testing at a pilot-scale facility, the three mMPMS were installed at a full-scale pulverized coal-fired plant, PacifiCorp’s Hunter 3. The systems were demonstrated over 20,000 hours at the plant during regular operation. Also, the sensor data was fed to the plant’s advanced process control system to evaluate the corrosion control by the operation changes and utilized to understand the impacts of load cycling with different ramping up and down speeds. At the end of the project, the systems were converted to the permanent installation at the power plant to be used with the advanced process control system installed at the plant.

Advanced Sensors, Corrosion, Ash Deposition, Optim↗

Fluid structural interaction simulation of the MP-1 plate performance irradiated in the advanced test reactor

Due to the high neutron fluxes they generate and increased heat transfer performance, plate type fuels are used in the U.S. high-performance research reactors. During irradiation, a significant amount of fission energy (i.e., ~200 MeV per atom) is released by the U-235 chain reaction that is carried away by the coolant. Through thermal hydraulic analysis, the coolant’s heat transfer capability is investigated to ensure that the desired fuel temperature can be maintained. At high temperatures, the plate undergoes elastic/plastic deformation, creep, and swelling as a result of both the temperature gradients and the fission gas production within the fuel. These effects are studied via fuel performance analysis. When the plate deformation is small enough that the coolant flow remains relatively unchanged, conducting these two types of analyses independently will suffice. But at high fuel burnups, the swollen plates may encroach into the coolant channels that separate the fuel plates from each other and cause these channels to narrow. Large reductions in channel gap size imperil cooling performance, causing fuel temperatures to rise. Further, if the plate deformation is asymmetric, the fuel centerline will shift toward one side of the channel, causing an uneven reduction in coolability. In addition, a boehmite (oxide) layer will, over time, grow on the plate surface, further obstructing the heat transfer process. To precisely predict fuel temperatures/deformation, a complete coupled analysis that considers coolant flow, heat transfer, oxide growth, elastic/plastic deformation, creep, and swelling is needed; however, this type of analysis method is not available in the literature. To fill the gap, this research developed a fluid structure interaction (FSI) approach to the fuel plate analysis, then successfully applied it to the Mini-Plate (MP)-1 experiment, which was irradiated in the Advanced Test Reactor (ATR) for both one and two cycles. The complete analysis coupled STAR-CCM+, a computational fluid dynamic (CFD) software for calculating flow, with Abaqus, a finite element analysis code for calculating plate deformation. Improvements in the results were found when comparing the fully coupled analysis to the independently conducted analyses but they were not significant due to the miniature size of the plates and the relatively short irradiation time. In the future, the fully coupled approach presented herein will be applied to full-size fuel plates with longer irradiation cycles once additional experiments become available.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Electrical injection and transport in Teflon-diluted hole transport materials

Diluting small molecule organic semiconductor hole transport materials with Teflon AF (TAF) dramatically increases their thermal/morphological stability while counterintuitively also increasing the current density in hole-only diodes. Here, we explore the origin of this current enhancement in co-evaporated blends of TAF and N,N'-Di(1-naphthyl)-N,N'-diphenyl-(1,1'-biphenyl)-4,4'-diamine (NPD) and find that it results from improved hole injection and not from a mobility enhancement due to semiconductor dilution effects. Whereas electroabsorption measurements indicate that an interface dipole shift of up to 0.5 eV (depending on TAF concentration) aides hole injection at the indium-tin-oxide anode, time-of-flight and temperature-dependent current-voltage measurements show that the mobility of 25 vol% TAF:NPD blends decreases by roughly an order of magnitude relative to neat NPD despite having slightly lower energetic disorder. For comparison, diluting with 25 vol% high density polyethylene achieves a similar reduction in energetic disorder but only decreases the mobility by a factor of two. Furthermore, these results highlight the tradeoff between reducing energetic disorder and frustrating percolation in diluted organic semiconductors, and represent a step toward thermally-evaporated polymer-small molecule composites that maintain high mobility while offering improved morphological stability and lower refractive index for organic light emitting diodes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Harnessing Quantum Computing for Energy Materials: Opportunities and Challenges

Developing high-performance materials is critical for diverse energy applications to increase efficiency, improve sustainability and reduce costs. Classical computational methods have enabled important breakthroughs in energy materials development, but they face scaling and time-complexity limitations, particularly for high-dimensional or strongly correlated material systems. Quantum computing (QC) promises to offer a paradigm shift by exploiting quantum bits with their superposition and entanglement to address challenging problems intractable for classical approaches. This Perspective discusses the opportunities in leveraging QC to advance energy materials research and the challenges QC faces in solving complex and high-dimensional problems. We present cases on how QC, when combined with classical computing methods, can be used for the design and simulation of practical energy materials. We also outline the outlook for error-corrected, fault-tolerant QC capable of achieving predictive accuracy and quantum advantage for complex material systems.

Algorithms↗

Wisconsin In Situ Penning (WISP) gauge: A versatile neutral pressure gauge to measure partial pressures in strong magnetic fields

A new type of in-vessel Penning gauge, the Wisconsin In Situ Penning (WISP) gauge, has been developed and successfully operated in the Wendelstein 7-X (W7-X) island divertor baffle and vacuum vessel. The capacity of the quantitative measurements of the neutral reservoir for light impurities, in particular, helium, is important for tokamaks as well as stellarator divertors in order to avoid fuel dilution and radiative energy loss. Penning gauges assisted by spectroscopy are a powerful tool to obtain the total neutral pressure as well as fractional neutral pressures of specific impurities. The WISP gauge is a miniaturized Penning gauge arrangement, which exploits the ambient magnetic field of magnetic confinement fusion experiments to establish the Penning discharge. Then, in situ spectroscopy is conducted to separate the fractional neutral pressures of hydrogen, helium, and possibly also other impurities. The WISP probe head was qualified using the magnetic field of the Magnetized Dusty Plasma Experiment at Auburn University between 0.25 T and 3.5 T [E. Thomas et al ., J. Plasma Phys. 81, 345810206 (2015)]. The in-depth quantitative evaluation for hydrogen and helium will be shown as well as an exploration of nitrogen, argon, and neon. A power law scaling between current I and pressure p, I = f (Gas,V) p n⁽Gas, B⁾ , was shown. The factor f is gas and anode potential dependent, while n is gas and magnetic field strength dependent. Pressure measurements from 0.1 mbar and down to 1 × 10 -5 mbar were achieved, demonstrating a reliable operating range for relevant pressure levels in the divertor and main vessel regions in current and future fusion devices, with a time resolution of up to 1 kHz. The lowest achievable pressure measurement increases with an increase in B and can be shifted with the anode potential V . At W7-X, the WISP probe head was mounted on an immersion tube setup that passes through the cryostat and places the probe head close to the plasma. Two probe heads were positioned in different divertor pump gaps, top and bottom, and one close to the plasma on the midplane in one module. The gauges were in situ calibrated together with the ASDEX pressure gauges [G. Haas and H.-S. Bosch, Vacuum 51, 39 (1998)]. Data were taken during the entire operation phase 1.2b, and measurements were coherent with other neutral gas pressure gauges. For the spectroscopic partial pressure measurements, channels of a spectroscopic detection system based on photo-multipliers, a so-called filterscope [R. J. Colchin et al ., Rev. Sci. Instrum. 74, 2068 (2003)], provided by the Oak Ridge National Lab were used.

47 OTHER INSTRUMENTATION↗

Distributed Fiber Sensors With High Spatial Resolution in Extreme Radiation Environments in Nuclear Reactor Cores

This paper is a comprehensive experimental report on the neutron radiation effects of distributed optical fiber sensors with enhanced Rayleigh scattering profiles in an in-pile environment. Femtosecond laser direct writing was used to inscribe Type-II modifications in standard telecom fibers and radiation-hardened fibers with fluorine-doped cores. Rayleigh backscattering signals were enhanced for continuous 1.5 m. In-pile lead-out sensors tests were carried out at the MIT Research Reactor for two months, which was operated at a nominal power of 5.7 MW with fast neutron (>0.1 MeV) flux of 1.29 × 1014 n/cm2/s and an in-core temperature of up to 560°C. Using the Optical Frequency Domain Reflectometry technique, the backscattering profiles of fiber sensors were interrogated with a 3-cm spatial resolution to monitor the temperature profile of the reactor. Results show that laser inscribed Type-II modifications in the form of nanogratings are highly stable against extreme temperature and ionizing radiation. Both standard telecom fibers and radiation-hardened fibers with laser-enhanced Rayleigh profiles can continuously perform distributed temperature measurements over the entire duration of the in-pile testing. Temperature coefficients of sensors and spectral shift quality were studied as functions of total radiation fluence. To the best of our knowledge, we present for the first time, the temperature profile of an operating nuclear reactor core with 3-cm spatial resolution, enabled by distributed fiber sensors with laser-enhanced Rayleigh scattering profiles. The high spatial resolution measurements can provide valuable data for the design and validation of digital twin and virtual reality of nuclear energy systems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Biochemical properties of glycerol kinase from the hypersaline-adapted archaeon Haloferax volcanii

ABSTRACT Extremophilic microorganisms are promising candidates for industrial and analytical biocatalysis.Haloferax volcanii, a halophilic archaeon that prefers glycerol over glucose, channels this substrate into central metabolism through glycerol kinase (GK). Here, we report the biochemical properties ofH. volcaniiGK and its potential for biotechnological applications. An N-terminal His-tagged GK was functionalin vivoand yielded 3 mg/L culture—4.5 times more enzyme than a C-terminal StrepII-tagged version. Size exclusion chromatography revealed a glycerol-induced oligomeric shift from homodimer to a dimer-dominant state with detectable tetramer. The purified enzyme showed robust activity across broad pH and salinity ranges, with optimal activity at 100 mM NaCl and 50°C–60°C. It retained catalytic activity in 5%–10% dimethyl sulfoxide (DMSO) and crude glycerol containing methanol. His-GK was freeze-thaw stable and thermotolerant in 2 M NaCl buffers. In the absence of ligands, the enzyme’s melting temperature (T m ) was 80°C. Glycerol increased the T m to 85°C, and combinations with MgCl₂ (84°C) or ATP (88°C) provided further stabilization. The highest T m (89°C) occurred with all three ligands, suggesting a cumulative stabilizing effect. Kinetic analyses revealed positive cooperativity for glycerol, ATP, and Mg² + ; Mn² + and Co² + also supported the activity.H. volcaniiGK is the first known GK to exhibit positive cooperativity with glycerol and ATP. Its high stability and substrate flexibility support its use in biodiesel waste valorization,in vitrobiocatalysis, and biosensor development—applications demanding robust, specific, and stable enzymes. IMPORTANCE This study reveals thatH. volcaniiGK exhibits positive cooperativity for glycerol, ATP, and Mg² + , a kinetic feature not previously reported for glycerol kinases. This behavior enables steep, switch-like responses to small substrate changes, offering unique advantages for biosensor design. Importantly,H. volcaniiGK also maintains high activity under extreme salinity, temperature, broad pH, and solvent conditions that typically limit enzyme use in industrial and environmental applications. These traits make this GK an ideal candidate for enzyme-based biosensors, which often suffer from poor tolerance to pH, solvent, and thermal stress. Its robustness supports its use in cross-linked enzyme crystals, an immobilization method that enhances enzyme stability and reusability under harsh conditions. Moreover, GKs are already employed in Mg² + detection kits; however,H. volcaniiGK’s ability to tolerate and respond to diverse divalent cations (e.g., Co² + , Mn² + ) broadens their potential for pollutant detection and environmental monitoring. These features collectively positionH. volcaniiGK as a valuable biocatalyst for biosensing,in vitrodiagnostics, and biotechnological applications requiring both precision and durability.

Biotechnology & Applied Microbiology↗

Fast and Scalable FFT-Based GPU-Accelerated Algorithms for Block-Triangular Toeplitz Matrices with Application to Linear Inverse Problems Governed by Autonomous Dynamical Systems

In this work, we present an efficient and scalable algorithm for performing matrix-vector multiplications (matvecs) for block Toeplitz matrices. Such matrices, which are shift-invariant with respect to their blocks, arise in the context of solving inverse problems governed by autonomous systems, and time-invariant systems in particular. In this article, we consider inverse problems that infer unknown parameters from observational data of a linear time-invariant dynamical system given in the form of partial differential equations (PDEs). Matrix-free Newton-conjugate-gradient methods are often the gold standard for solving these inverse problems, but they require numerous actions of the Hessian on a vector. Matrix-free adjoint-based Hessian matvecs require solution of a pair of linearized forward/adjoint PDE solves per Hessian action, which may be prohibitive for large-scale inverse problems. Time invariance of the forward PDE problem leads to a block Toeplitz structure of the discretized parameter-to-observable (p2o) map defining the mapping from inputs (parameters) to outputs (observables) of the PDEs. This block Toeplitz structure enables us to exploit two key properties: (1) compact storage of the p2o map and its adjoint, and (2) efficient fast Fourier transform–based Hessian matvecs. The proposed algorithm is mapped onto large multi-GPU clusters and achieves more than 80% of peak bandwidth on NVIDIA A100 GPUs. Excellent weak scaling is shown for up to 48 A100 GPUs. For the targeted problems, the implementation executes Hessian matvecs within fractions of a second, which is orders of magnitude faster than can be achieved by conventional matrix-free Hessian matvecs via forward/adjoint PDE solves.

97 MATHEMATICS AND COMPUTING↗

Single-shot observation of nonlinear pulse splitting in a Kerr medium

We report single-shot, time-resolved observation of self-steepening and temporal splitting of near-infrared, 50 fs, micro-joule pulses propagating nonlinearly in flint (SF11) glass. A coherent, smooth-profiled, 60-nm-bandwidth probe pulse that propagated obliquely to the main pulse through the Kerr medium recorded a time sequence of longitudinal projections of the main pulse’s induced refractive index profile in the form of a phase-shift “streak,” in which frequency–domain interferometry recovered with ∼10 fs temporal resolution. A three-dimensional simulation based on a unidirectional pulse propagation equation reproduced observed pulse profiles.

Chang, Yen-Yu↗

American-Made Challenges Round 2 Voucher: Orison Enables Solar

This report documents the technical assistance provided to Orison as part of the American-Made Challenges Round 2 Voucher: Orison Enables Solar project. A testbed was developed and used to demonstrate, using modeling and simulation, the capability of controlled behind-the-meter energy storage to (1) reduce net load variability caused by appliances and distributed generation and (2) enable customers to respond to the time-of-use pricing. We demonstrate that load leveling can shave peak loads or limit solar photovoltaic export and that load shifting enables customers to reduce load during “peak” pricing intervals. In the scenarios that we simulated, load leveling was effective up to the charging and discharging limits of the storage systems and the effectiveness of load shifting was a function of the energy capacity of the storage systems. This report describes the co-simulation testbed and the scenarios simulated, including controller setup and simulation results. Considerations for a multi-objective controller are also discussed.

14 SOLAR ENERGY↗

RapidClave® Technology Demonstrations – II Hat Stiffener

This project sought to evaluate the cost and performance of both polymeric and metallic tooling for use with the RapidClave® composite process. Teijin’s rapid cure carbon/epoxy prepreg, Q183, was selected to fabricate demonstration components, of aerospace design, which currently are fabricated via compression molding using matched steel tooling. The ability to fabricate components using this prepreg in the RapidClave®, and using low cost single sided tooling, could significantly reduce the cost of part manufacture and shorten production lead times for tooling fabrication. The component demonstration article was a “hat stiffener” geometry typically used to reinforce fuselage or wing skins of an aircraft. A rapid preforming process was demonstrated and rapid process cycle times similar to compression molding were achieved by placing the preform on hot tooling, thus avoiding the time and energy associated with typical thermal cycling. A 90% reduction in tooling cost was achieved through the use of single sided aluminum or polymeric additive tooling when compared to compression molding matched steel tooling. In addition to rapid cycle times, RapidClave® offers the ability to change out tooling in minutes such that multiple part geometries can be fabricated in one shift.

36 MATERIALS SCIENCE↗

Optimizing Control for Efficient Load Shifting with Thermal Energy Storage in Existing HVAC Systems

This project developed the integration of a direct-contact heat exchanger (DCHX) based thermal energy storage (TES) system with a chiller–air handling unit (AHU) plant to evaluate its potential for reducing building energy costs. Detailed physical models of the TES unit, building envelope, and HVAC components were developed alongside simplified control-oriented models to support both high-fidelity simulation and real-time optimization. Two control strategies were implemented and compared: a rule-based control (RBC) aligned with utility time-of-use (ToU) rates, and a model predictive control (MPC) framework leveraging forecasts of building load, weather, and internal gains.Simulation results show that the RBC strategy reduced daily electricity costs by around 30% by shifting cooling production from on-peak to off-peak hours. In contrast, the MPC strategy achieved significantly greater performance, reducing daily operating costs by up to 44% and peak-hour costs by more than 60%. Both strategies maintained indoor thermal comfort within acceptable limits, with MPC further improving load distribution and reducing equipment cycling.The outcomes confirm that TES integration, particularly when coordinated with advanced predictive control, can provide substantial cost savings and on-peak demand reduction. . These findings directly support the U.S. Department of Energy’s goals for grid-interactive efficient buildings and demonstrate the potential of TES-enabled HVAC systems for scalable deployment across the commercial building.

99 GENERAL AND MISCELLANEOUS↗

Accurate in vivo tumor detection using plasmonic-enhanced shifted-excitation Raman difference spectroscopy (SERDS)

For the majority of cancer patients, surgery is the primary method of treatment. In these cases, accurately removing the entire tumor without harming surrounding tissue is critical; however, due to the lack of intraoperative imaging techniques, surgeons rely on visual and physical inspection to identify tumors. Surface-enhanced Raman scattering (SERS) is emerging as a non-invasive optical alternative for intraoperative tumor identification, with high accuracy and stability. However, Raman detection requires dark rooms to work, which is not consistent with surgical settings. Methods: Herein, we used SERS nanoprobes combined with shifted-excitation Raman difference spectroscopy (SERDS) detection, to accurately detect tumors in xenograft murine model. Results: We demonstrate for the first time the use of SERDS for in vivo tumor detection in a murine model under ambient light conditions. We compare traditional Raman detection with SERDS, showing that our method can improve sensitivity and accuracy for this task. Conclusion: Our results show that this method can be used to improve the accuracy and robustness of in vivo Raman/SERS biomedical application, aiding the process of clinical translation of these technologies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Modeling the Effect of the Interplanetary Magnetic Field on Cosmic Ray Muon Shadows

Cosmic ray (CR) muon shadows were located using muon data from the MINOS Far Detector in Northern Minnesota. The shadows were located independently across 3 time periods; near solar minimum, solar maximum, and over the entire 13 year span of the data. A muon distribution for each shadow was then sampled to simulate CR motions near the Sun using the Parker Spiral model of the Sun’s IMF and a dipole model of the Earth’s GMF. The resulting particle distributions were then compared to their position with respect to the Sun. An analytical relationship between the primary CR cutoff energy E0 and the associated shadow deflections was observed. The minimum cutoff energy of the primary CR energy spectrum consistent with the Parker Spiral model was determined to be $E_0 = 22.8 \pm 1.2$~TeV. This value is 2-3 times higher than the experimentally determined value. Complexities that can be added to the Parker Spiral Model and affect the overall shift of the the Sun’s cosmic ray shadow are discuss ed to explain discrepancies. Results show that the Parker Spiral model is most consistent with the the solar minimum shadow and least consistent with the solar maximum shadow, as expected.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Submodule-Level Decentralized Control for Cascaded Quadruple Active Bridge-Based Medium-Voltage Converters

High-power, low-voltage DC (LVDC) loads, such as AI data centers, electric vehicle fast-charging hubs, and hydrogen electrolyzers, are accelerating the deployment of solid-state transformers (SSTs) that directly interface medium-voltage AC (MVAC) to LVDC, streamlining power delivery. This paper proposes a submodule-level decentralized control strategy for an SST comprising a cascaded H-bridge (CHB) active rectifier on the MVAC side and per-module quadruple-active-bridge (QAB) converters on the LVDC side with local controllers. The proposed strategy eliminates the need for high-fidelity, real-time cross-isolation communication by synchronizing all active bridges to a common heartbeat signal and enabling the primary- and secondary-side controllers to generate independent local control signals (phase shifts). The resultant phase shifts naturally govern average and double-line frequency power transfer between active bridges, enabling each side to meet its control objectives using only local measurements. Detailed controller design considerations are presented based on a stability analysis. Comprehensive electromagnetic-transient simulations demonstrate tight DC-link voltage regulation, high-quality grid currents, and robust operation under load transients and grid disturbances, without real-time communication among the submodule controllers. These results indicate that the proposed scheme can significantly enhance the reliability, scalability, and flexibility of cascaded bridge-based medium-voltage converters.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Submodule-Level Decentralized Control for Cascaded Quadruple Active Bridge-Based Medium-Voltage Converters: Preprint

High-power, low-voltage DC (LVDC) loads, such as AI data centers, electric vehicle fast-charging hubs, and hydrogen electrolyzers, are accelerating the deployment of solid-state transformers (SSTs) that directly interface medium-voltage AC (MVAC) to LVDC, streamlining power delivery. This paper proposes a submodule-level decentralized control strategy for an SST comprising a cascaded H-bridge (CHB) active rectifier on the MVAC side and per-module quadruple-active-bridge (QAB) converters on the LVDC side with local controllers. The proposed strategy eliminates the need for high-fidelity, real-time cross-isolation communication by synchronizing all active bridges to a common heartbeat signal and enabling the primary- and secondary-side controllers to generate independent local control signals (phase shifts). The resultant phase shifts naturally govern average and double-line frequency power transfer between active bridges, enabling each side to meet its control objectives using only local measurements. Detailed controller design considerations are presented based on a stability analysis. Comprehensive electromagnetic-transient simulations demonstrate tight DC-link voltage regulation, high-quality grid currents, and robust operation under load transients and grid disturbances, without real-time communication among the submodule controllers. These results indicate that the proposed scheme can significantly enhance the reliability, scalability, and flexibility of cascaded bridge-based medium-voltage converters.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Direct estimation of the density of states for fermionic systems

Simulating time evolution is one of the most natural applications of quantum computers and is thus one of the most promising prospects for achieving practical quantum advantage. Here, we develop quantum algorithms to extract thermodynamic properties by estimating the density of states (DOS), which is a central object in quantum statistical mechanics. We introduce several key innovations that significantly improve the practicality and extend the generality of previous techniques. First, our approach allows one to estimate the DOS only for a specific subspace of the full Hilbert space. This is crucial for fermionic systems, since both canonical and grand canonical ensemble thermal equilibrium properties depend on subspaces of fixed number. Second, in our approach, by time evolving very simple, random initial states, such as randomly chosen computational basis states, we can exactly recover the DOS on average. Third, due to circuit-depth limitations, we only reconstruct the DOS up to a convolution with a Gaussian window—thus all imperfections that shift the energy levels by less than the width of the convolution window will not significantly affect the estimated DOS. For these reasons, we find the approach is a promising candidate for early quantum advantage as even short-time, noisy dynamics can yield a semiquantitative reconstruction of the DOS (convolution with a broad Gaussian window), while early fault-tolerant devices will likely enable higher-resolution DOS reconstruction through longer time evolutions. We demonstrate the practicality of our approach in representative Fermi-Hubbard and spin models and indeed find that our approach is highly robust against algorithmic errors in the time evolution and against gate noise. We further demonstrate that our approach is compatible with noisy intermediate-scale quantum (NISQ) computing NISQ-friendly variational techniques, introducing and leveraging a technique for variational time evolution.

97 MATHEMATICS AND COMPUTING↗

Automated Hybrid Variance Reduction on Advanced Architectures in the Shift Monte Carlo Code

Monte Carlo transport methods are the most accurate schemes for solving problems with complex energy and spatial features, but they come with a high computational cost. Although hybrid methods have enabled the use of Monte Carlo transport for a large class of problems, they still require significant computing resources. Modern multicore CPUs with large numbers of compute cores and graphical processing units (GPUs) provide opportunities to optimize the memory and run-time costs of hybrid Monte Carlo methods. This paper documents the development and analysis of three Monte Carlo transport algorithms that support hybrid transport using the consistent adjoint-driven importance sampling (CADIS) and forward-weighted CADIS methods in the Shift Monte Carlo code: history-based transport using static and dynamic threading on multicore CPUs and event-based transport enabling weight window tracking on GPUs. The results are shown for two challenging hybrid problems on the Frontier supercomputer at the Oak Ridge Leadership Computing Facility. The results show that all three methods yield good performance and enable solutions of difficult fixed-source transport problems in less than 2 min on 20 nodes of Frontier. Dynamic threading was observed to give up to 20% better scaling behavior than static threading. Moreover, the AMD Instinct 250X GPU was found to give 9 to 11 times greater throughput per graphics compute die than the best CPU performance. In conclusion, additional opportunities for optimization of hybrid transport on GPUs are discussed.

Denovo↗