SEARCH · Engineering Papers

Results for “Variational approach”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,045 records · Page 58

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Latent Neural ODE for Integrating Multi-Timescale Measurements in Smart Distribution Grids

Under a smart grid paradigm, there has been an increase in sensor installations to enhance situational awareness. The measurements from these sensors can be leveraged for real-time monitoring, control, and protection. However, these measurements are typically irregularly sampled. These measure-ments may also be intermittent due to communication bandwidth limitations. To tackle this problem, this paper proposes a novel latent neural ordinary differential equations (LODE) approach to aggregate the unevenly sampled multivariate time-series measurements. The proposed approach is flexible in performing both imputations and predictions while being computationally efficient. Simulation results on IEEE 37 bus test systems illustrate the efficiency of the proposed approach.

multi time-scale measurements↗

DOE OSTI2025

A convergence metric for counting statistics in time-resolved small angle neutron scattering

Here, this work introduces a model-independent, dimensionless metric for predicting optimal measurement duration in time-resolved small-angle neutron scattering using early-time data. Built on a Gaussian process regression framework, the method reconstructs scattering profiles with quantified uncertainty, even from sparse or noisy measurements. Demonstrated on the EQ-SANS instrument at the Spallation Neutron Source, the approach generalizes to general SANS instruments with a two-dimensional detector. A key result is the discovery of a dimensionless convergence metric revealing a universal power-law scaling in profile evolution across soft matter systems. When time is normalized by a system-specific characteristic time t*, the variation in inferred profiles collapses onto a single curve with an exponent between −2 and −1. This trend emerges within the first ten time steps, enabling early prediction of measurement sufficiency. The method supports real-time experimental optimization and is especially valuable for maximizing efficiency in low-flux environments such as compact accelerator-based neutron sources.

Tung, Chi-Huan [Oak Ridge National Laboratory (ORN↗

DOE OSTI2025

Detecting axion dynamics on the surface of magnetic topological insulators

Axions, initially proposed to solve the strong CP problem, have recently gained attention in condensed matter physics, particularly in topological insulators. However, detecting axion dynamics has proven challenging, with no experimental confirmations to date. In this paper, we identify the surface of magnetic topological insulators as an ideal platform for observing axion dynamics. The vanishing bulk gap at the surface allows for order O(1) variations in the axion field, making the detection of axionlike phenomena more feasible. In contrast, these phenomena are strongly suppressed in the bulk due to the small magnetic exchange gap. We investigate two-photon decay as a signature of axion dynamics and calculate the branching ratio using a perturbative approach. Our findings reveal that the photon flux emitted from the surface is in plane and orders of magnitude larger than that from the bulk, making it detectable with modern microwave technology. We also discuss potential material platforms for detecting axion two-photon decay and strategies to enhance the signal-to-noise ratio.

Gao, Zhi-Qiang↗

DOE OSTI2020

Learning to Optimize Variational Quantum Circuits to Solve Combinatorial Problems

Quantum computing is a computational paradigm with the potential to outperform classical methods for a variety of problems. Proposed recently, the Quantum Approximate Optimization Algorithm (QAOA) is considered as one of the leading candidates for demonstrating quantum advantage in the near term. QAOA is a variational hybrid quantum-classical algorithm for approximately solving combinatorial optimization problems. The quality of the solution obtained by QAOA for a given problem instance depends on the performance of the classical optimizer used to optimize the variational parameters. In this paper, we formulate the problem of finding optimal QAOA parameters as a learning task in which the knowledge gained from solving training instances can be leveraged to find high-quality solutions for unseen test instances. To this end, we develop two machine-learning-based approaches. Our first approach adopts a reinforcement learning (RL) framework to learn a policy network to optimize QAOA circuits. Our second approach adopts a kernel density estimation (KDE) technique to learn a generative model of optimal QAOA parameters. In both approaches, the training procedure is performed on small-sized problem instances that can be simulated on a classical computer; yet the learned RL policy and the generative model can be used to efficiently solve larger problems. Furthermore, extensive simulations using the IBM Qiskit Aer quantum circuit simulator demonstrate that our proposed RL- and KDE-based approaches reduce the optimality gap by factors up to 30.15 when compared with other commonly used off-the-shelf optimizers.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Status of Additive Manufacturing Capabilities for Processing Refractory Alloys Under the Mo-99 Program

The enriched 100 Mo target, denoted as aMo, that Northstar Medical Radioisotopes is developing to produce 99 Mo from an accelerator-based method is a disk-type geometry with a 0.5–0.715 mm thickness and a 29 mm outer diameter. A press and sinter method was adopted to fabricate these samples from aMo powder feedstock. Inherent porosity from the press and sinter method makes the disks more easily dissolvable during 99 Mo postprocessing. However, the disk-type geometry will be subject to thermal stresses due to temperature variations once the geometry is in line with the accelerated electron beam, causing warping and possible failure. Additive manufacturing (AM) processes such as laser powder bed fusion (LPBF) enable the manufacturing of custom, complex geometries that may otherwise be difficult to produce with more conventional methods (e.g., press and sinter approach). Therefore, the National Nuclear Security Administration’s Molybdenum-99 Program invested in developing an AM facility dedicated to processing and printing refractory powders. This report details the LPBF process, how irregular refractory powder can be prepared for LPBF, and the equipment that is available to support the program.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Performance Evaluation of Intelligent Solar Control Software Through Hardware-in-the-Loop (CRADA Final Report)

Recent research has highlighted the potential for solar to act as a zero-marginal-cost and zero-emission flexibility resource on the bulk power system when operated with advanced control systems. To increase the performance of these systems, leading technologies, including machine learning (ML) and hierarchical inverter set point allocation, have been developed by Latimer Controls, Inc. to estimate the headroom of large PV plants for grid operation and control; however, these technologies lack comprehensive validation under real-world application scenarios. Latimer Controls, Inc. received two voucher awards for research at a national laboratory from the Department of Energy American Made Solar Prize Round 6. The National Renewable Energy Laboratory (NREL) was selected to collaborate with Latimer staff to conduct a performance evaluation of Latimer PV control software. The NREL team will develop a hardware-in-the-loop (HIL) testbed to perform testing and validation of the Latimer PV control technology in a de-risked yet realistic testbed environment. Latimer and NREL worked together to analyze the test data, draw conclusions from the results, and disseminate the resulting scientific findings. In this CRADA work, we propose to test and validate the real-world application of the Latimer Control solution in an HIL environment. We evaluate the performance of different flexible solar technologies in responding to automatic generation control signals in a closed-loop fashion. In particular, a data-driven potential high limit (PHL) estimation is developed for large solar plants to accurately estimate their headroom so that they have fast and short-time regulation and control capability to participate in grid services and respond to grid signals in real time (e.g., AGC). This PHL estimation algorithm is embedded in a hardware power plant controller (PPC) and tested with an IEEE-39 bus system model developed in RTDS. To account for the varying cloud conditions and diverse inverter dispatches, we developed a 135-MW PV plant with detailed modeling of 27 individual PV modules and inverters using RTDS. The real-world communications used in such big plants, such as ModBus TCP/IP for inverter level and DNP3 for plant level, were developed to emulate the real-world applications in big PV plants. The ML-based PHL estimation method is tested under nine separate weather scenarios against the ‘reference-control’ solution, hereafter referred to as the baseline solution. The baseline method reserves a subset of inverters (reference group) to operate at their PHL at all times and dispatches only the remaining inverters (control group) at curtailed levels to fulfill the flexibility need. Despite being successfully piloted by NREL in California in 2017 and Chile in 2020, there exist two gaps in the state of the art to fully unlock the flexibility of PV plants: a. There is a trade-off between the PHL estimation accuracy and the flexibility range. b. There lacks granularity in the PHL estimation to capture the variation across inverters. The Latimer solution seeks to address these gaps by applying machine learning methods to improve PHL estimation accuracy while accounting for variability at every inverter. Performance metrics were taken from the 2023 Georgia Power CARES utility-scale RFP. The results demonstrate that the ML-based approach outperforms the traditional baseline method in PHL estimation accuracy for 7 of 9 scenarios. The average PHL error across the nine scenarios was 7.40% for the ML-based method, 2.06% less than the 9.46% PHL error average across scenarios that was exhibited by the baseline method. Additionally, the PHL error was below 5% for at least 95% of the testing interval for 3 of 9 tested intervals with the ML approach, whereas it did not achieve this metric for any of the baseline tests. Overall, simulation results indicate the superior performance of an ML-based approach compared to the conventional baseline reference-control approach, showcasing its potential to support grid stability and operational efficiency. This laboratory HIL testing using real PPC, representative power system simulation models in real-time with detailed PV plant and inverter models, and real-world communication protocols gives us confidence that this machine learning based PHL estimation algorithm works well in the hardware PPC and therefore de-risks future field commissioning. The end goal of this project is to advance grid technology to address the grid operation challenges brought by solar plant’s variability and uncertainties in power generation.

14 SOLAR ENERGY↗

DOE OSTI2022

Evaluation of hydrograph separation techniques with uncertain end-member composition

Hydrograph separation is one of many approaches used to analyse shifts in source water contributions to stream flow resulting from climate change in remote watersheds. Understanding these shifts is vital, as shifts in source water contributions to a stream can shape water management decisions. Because remote watersheds are often inaccessible and have poorly characterized contributing water sources, or end-members, it is critical to understand the implications of using different hydrograph separation techniques in these data-limited environments. To explore the uncertainty associated with different techniques, results from two hydrograph separation techniques, mass balance and principle component analysis, were compared using 3 years of aqueous geochemical data from the East River watershed located in the Elk Mountains of Central Colorado. Solute concentrations of the end-members were characterized by both a limited set of direct chemical measurements of different sources and detailed seasonal instream chemistry to examine the influences of uncertain end-member compositions in a data-limited environment. Annual volumetric end-member contributions to stream flow had relatively good agreement across separation techniques. Large variations in time were observed in the hydrograph separations, depending on the end-member type, and estimated flow contributions varied between the selected solutes. End-member concentrations characterized by stream chemistry showed several limitations including a reduced number of distinguishable end-members and differences in timing of flow contributions. Here the results highlight the benefits of using multiple hydrograph separation techniques by providing a ‘weight-of-evidence’ approach to environments with limited end-member concentration data.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Marcus–Hush–Chidsey kinetics at electrode–electrolyte interfaces

Electrochemical kinetics at electrode–electrolyte interfaces limit the performance of devices including fuel cells and batteries. While the importance of moving beyond Butler–Volmer kinetics and incorporating the effect of electronic density of states of the electrode has been recognized, a unified framework that incorporates these aspects directly into electrochemical performance models is still lacking. In this work, we explicitly account for the density functional theory-calculated density of states numerically in calculating electrochemical reaction rates for a variety of electrode–electrolyte interfaces. We first show the utility of this for two cases related to Li metal electrodeposition and stripping on a Li surface and a Cu surface (anode-free configuration). The deviation in reaction rates is minor for cases with flat densities of states such as Li, but is significant for Cu due to nondispersive d-bands creating large variation. Finally, we consider a semiconducting case of a solid-electrolyte interphase consisting of LiF and Li2CO3 and note the importance of the Fermi level at the interface pinned by the redox reaction occurring there. We identify the asymmetry in reaction rates as a function of discharge/charge naturally within this approach.

Kurchin, Rachel (ORCID:0000000221474809)↗

DOE OSTI2024

Dynamic Adaptive Relaying for Distribution Grids with High Inverter-Based Generation

The penetration of distributed energy resources (DERs) in distribution systems is rapidly increasing. This can lead to large variations in the operating load and the fault current seen by a relay, which can impact the reliability of overcurrent protection. Adaptive overcurrent relaying (AOCR) can address this concern by enabling the relay to adjust its pickup settings based on the system operating conditions. This paper presents a novel dynamic AOCR approach that allows the relays to locally estimate the pickup settings and dynamically modify them. The approach is based on a two-stage estimation algorithm that uses local measurements and system information to estimate the worst-case fault current and the relay pickup settings. The relay pickup settings are then dynamically modified based on the estimated fault current. The performance of the proposed approach is evaluated on multiple configurations of the EPRI J1 feeder and the IEEE 13-bus system. The results show that the proposed approach can effectively improve the reliability of overcurrent protection in distribution systems with DERs.

adaptive overcurrent protection↗

DOE OSTI2023

Spatiotemporal Variability of Electric System Reliability Metrics

As urbanization continues to grow, urban density also rises, placing significant stress on critical infrastructure. Cities often face challenges in investing in new infrastructure or upgrading existing ones. One area of infrastructure that has come under strain is the electric grid, leading to prolonged power outages. According to data from the U.S. Energy Information Administration, the average duration of power outages has steadily doubled from 2013 to 2021, primarily due to extreme weather events. The IEEE has developed a series of reliability metrics to assess the reliability of the electric system. However, these metrics are only applied at the utility level, making it challenging to comprehend the variations in these indices at finer spatial resolutions and, consequently, evaluate reliability at the non-utility level. In this study, we leveraged electric outage data collected at 15-minute intervals at the county level to estimate electric system reliability over a six-year period using a novel data analysis approach. Our findings reveal a gradual decline in reliability in the U.S. Southwest compared to national averages. Furthermore, the duration of outages per customer per event has significantly decreased in large portions of California and Florida. These discoveries can be valuable for utility companies as they seek to pinpoint regional anomalies and plan for prioritized remediation efforts.

Singh, Nagendra↗

DOE OSTI2021

Integrated reference cavity with dual-mode optical thermometry for frequency correction

Photonic integrated resonators have advantages over traditional benchtop cavities in terms of size, weight, and cost with the potential to enable applications that require spectrally pure light. However, integrated resonators suffer from temperature-dependent frequency variations and are sensitive to external environmental perturbations, which hinders their usage in precision frequency applications. One solution is to use interrogation of the cavity temperature through dual-mode optical thermometry (DMOT) by measuring the shift of the resonance frequency difference between two polarization or optical frequency modes. Yet this approach has only been demonstrated in bulk-optic whispering gallery mode and fiber resonators. In this paper, we implement dual-mode optical thermometry in an ultra-high Q integrated silicon nitride resonator. A dual-mode resonance frequency difference temperature sensitivity of 188 ± <#comment/> 15 M H z / K is measured. We demonstrate feedforward DMOT frequency correction that, under an applied external temperature ramp, is able to reduce the optical frequency change to 0.31 kHz/s as compared to an uncorrected 10.03 kHz/s, a factor of 30 × <#comment/> reduction. These results show promise for on-chip frequency correction solutions for quantum, metrology, atomic, and coherent optical communications applications.

Zhao, Qiancheng (ORCID:0000000206968667)↗

DOE OSTI2024

Contact-dependent growth inhibition (CDI) systems deploy a large family of polymorphic ionophoric toxins for inter-bacterial competition

Contact-dependent growth inhibition (CDI) is a widespread form of inter-bacterial competition mediated by CdiA effector proteins. CdiA is presented on the inhibitor cell surface and delivers its toxic C-terminal region (CdiA-CT) into neighboring bacteria upon contact. Inhibitor cells also produce CdiI immunity proteins, which neutralize CdiA-CT toxins to prevent auto-inhibition. Here, we describe a diverse group of CDI ionophore toxins that dissipate the transmembrane potential in target bacteria. These CdiA-CT toxins are composed of two distinct domains based on AlphaFold2 modeling. The C-terminal ionophore domains are all predicted to form five-helix bundles capable of spanning the cell membrane. The N-terminal "entry" domains are variable in structure and appear to hijack different integral membrane proteins to promote toxin assembly into the lipid bilayer. The CDI ionophores deployed by E. coli isolates partition into six major groups based on their entry domain structures. Comparative sequence analyses led to the identification of receptor proteins for ionophore toxins from groups 1 & 3 (AcrB), group 2 (SecY) and groups 4 (YciB). Using forward genetic approaches, we identify novel receptors for the group 5 and 6 ionophores. Group 5 exploits homologous putrescine import proteins encoded by puuP and plaP, and group 6 toxins recognize di/tripeptide transporters encoded by paralogous dtpA and dtpB genes. Finally, we find that the ionophore domains exhibit significant intra-group sequence variation, particularly at positions that are predicted to interact with CdiI. Accordingly, the corresponding immunity proteins are also highly polymorphic, typically sharing only ~30% sequence identity with members of the same group. Competition experiments confirm that the immunity proteins are specific for their cognate ionophores and provide no protection against other toxins from the same group. The specificity of this protein interaction network provides a mechanism for self/nonself discrimination between E. coli isolates.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Perspectives and Challenges in Bolide Infrasound Processing and Interpretation: A Focused Review with Case Studies

Infrasound sensing plays a critical role in the detection and analysis of bolides, offering passive, cost-effective global monitoring capabilities. Key objectives include determining the timing, location, and yield of these events. Achieving these goals requires a robust approach to detect, analyze, and interpret rapidly moving elevated sources such as bolides (also re-entry). In light of advancements in infrasonic methodologies, there is a need for a comprehensive overview of the characteristics that distinguish bolides from other infrasound sources and methodologies for bolide infrasound analysis. This paper provides a focused review of key considerations and presents a unified framework to enhance infrasound processing approaches specifically tailored for bolides. Three representative case studies are presented to demonstrate the practical application of infrasound processing methodologies and deriving source parameters while exploring challenges associated with bolide-generated infrasound. These case studies underscore the effectiveness of infrasound in determining source parameters and highlight interpretative challenges, such as variations in signal period measurements across different studies. Future research should place emphasis on improving geolocation and yield accuracy. This can be achieved through rigorous and systematic analyses of large, statistically significant samples of such events, aiming to resolve interpretative inconsistencies and explore the causes for variability in signal periods and back azimuths. The topic described here is also relevant to space exploration involving planetary bodies with atmospheres, such as Venus, Mars, and Titan.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

The Current State of Light-Duty Electric Vehicle Supply Equipment Costs: An Assessment of Contemporary Understanding

This study uses a hybrid meta-analysis and literature-review approach to understand the current state of knowledge regarding the costs of electric vehicle supply equipment (EVSE). We present a novel way to consider, categorize, and label measures of cost and show cost measure estimates from a sample of 13 recent studies. We find that in general, there is too much variation and too few commonly represented EVSE cost measures to reasonably provide aggregate figures for these measures. We propose a convention for presenting EVSE cost measures that includes the application (commercial or residential), the power level (Level 1, Level 2, DCFC [with further distinction based on rated power capacity]), and the type of cost measure (hardware, installation, operation, and total cost). We contend that providing researchers with standard cost measures will help to advance our knowledge of EVSE costs by ensuring that future work will use common metrics. Establishing common metrics will enable conventional meta-analyses that will make assessments of EVSE costs even more accessible. Additionally, common metrics will make tracking costs more reliable as the technology continues to evolve and become more ubiquitous.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2024

Microstructure and mechanical properties of friction stir weld performed on neutron-irradiated 304L steel with helium

As nuclear power plants (NPPs) approach or exceed 40–60 years of service, it may become necessary to repair damaged neutron-irradiated components to prolong longevity. However, fusion welding repair of irradiated steels and metallic materials is challenging because of helium-induced cracking. This study used friction stir welding (FSW) to address helium-related issues, such as helium bubble formation and grain boundary cracking. The microstructure, helium-induced degradation, and mechanical properties of a friction stir weld produced on neutron-irradiated 304L stainless steel with approximately 5.2 appm of helium were characterized. The analysis focused on variations in grain size, texture, and the morphology of helium-induced damage in the friction stir weld. Mechanical properties were characterized on the irradiated base metal (BM) and metallurgical zones of the friction stir weld: the stir zone (SZ), the thermo-mechanically affected zone (TMAZ), and the heat affected zone (HAZ). Only minor scattered porosity in the SZ and TMAZ and a few short microcracks (below 20 μm in length) in the TMAZ were observed, indicating limited helium-induced degradation. Tensile tests revealed good mechanical properties and fractography analysis demonstrated predominantly ductile fracture. In conclusion, the results highlight the immediate and substantial benefits of the FSW approach for repairing or joining helium-containing irradiated materials in NPPs.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Path-BigBird: An AI-Driven Transformer Approach to Classification of Cancer Pathology Reports

PURPOSE Surgical pathology reports are critical for cancer diagnosis and management. To accurately extract information about tumor characteristics from pathology reports in near real time, we explore the impact of using domain-specific transformer models that understand cancer pathology reports. METHODS We built a pathology transformer model, Path-BigBird, by using 2.7 million pathology reports from six SEER cancer registries. We then compare different variations of Path-BigBird with two less computationally intensive methods: Hierarchical Self-Attention Network (HiSAN) classification model and an offthe-shelf clinical transformer model (Clinical BigBird). We use five pathology information extraction tasks for evaluation: site, subsite, laterality, histology, and behavior. Model performance is evaluated by using macro and micro F 1 scores. RESULTS We found that Path-BigBird and Clinical BigBird outperformed the HiSAN in all tasks. Clinical BigBird performed better on the site and laterality tasks. Versions of the Path-BigBird model performed best on the two most difficult tasks: subsite (micro F 1 score of 72.53, macro F 1 score of 35.76) and histology (micro F 1 score of 80.96, macro F 1 score of 37.94). The largest performance gains over the HiSAN model were for histology, for which a Path-BigBird model increased the micro F 1 score by 1.44 points and the macro F 1 score by 3.55 points. Overall, the results suggest that a Path-BigBird model with a vocabulary derived from wellcurated and deidentified data is the best-performing model. CONCLUSION The Path-BigBird pathology transformer model improves automated information extraction from pathology reports. Although Path-BigBird outperforms Clinical BigBird and HiSAN, these less computationally expensive models still have utility when resources are constrained.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2023

S4PST: Sustainability for Programming Systems and Tools: May Workshop Report

The US Department of Energy (DOE) Exascale Computing Project (ECP) has fostered and strengthened the use of modern software engineering practices for developing applications and libraries, and this effort has resulted in the coordinated and interoperable E4S1 and xSDK2 ecosystems. Although this approach is cost-effective, it relies on robust programming systems and tools (PST) as the underlying foundation for our HPC software. At present, our primary PST stack consists of traditional high-performance computing (HPC) languages, namely Fortran, C, C++, and the popular Python language for data analysis and AI workflows. These languages support various programming frameworks and run-time abstractions that enable parallelism and concurrency across multiple node architectures and thousands of nodes through a variety of interconnect systems. However, to accommodate users’ diverse needs, certain aspects of the HPC ecosystem are delegated to vendor-specific or third-party implementations that extend beyond a particular scientific domain. This broader scope results in a multitude of specifications and variations, which leads to a complex orchestration of many-ecosystems. Unfortunately, this complexity in the ecosystem imposes additional overhead costs on consumers during the latter stages of the development cycle. In addition to the software ecosystem challenge, the upcoming conclusion of the ECP by December 2023 has raised significant concerns within the HPC programming systems community, from both the economic and social perspectives. The ECP has implemented a management structure for software development and funding decisions across all ECP participants by following a conventional hierarchical and centralized approach. However, this structure has prompted certain considerations within the community, particularly in anticipation of the Software Sustainability initiative by the DOE’s Advanced Scientific Computing Research Program (ASCR). For the success of this new initiative, it is of utmost importance to secure consistent funding and foster close engagement with researchers and core developers of existing programming-system products. This collaboration is vital to maintaining the critical capabilities of the current software during the transition phase while proactively adapting to future technology and workforce trends. The community recognizes the significance of adapting to emerging trends and is aware of the inherent fragility of the HPC software ecosystem, particularly in relation to programming systems that cater to all users. The ability to adapt and evolve is essential to staying relevant and effectively addressing these technical, economic, and social challenges. The S4PST team, which represents one of the six ASCR Software Sustainability seedling projects, is dedicated to tackling these challenges through community-based approaches that go beyond the scope of the DOE. This involves collaboration between national laboratories with academia, non-DOE institutions, hardware and system vendors, and international partners. By fostering these partnerships, we aim to create a robust and sustainable HPC software ecosystem that can effectively meet the needs of the community. This new community effort, driven by the eight DOE labs, will take on the responsibility of guiding funding decisions for programming-systems development and maintenance with transparency and consistency across all decisions. Additionally, the team will offer common technical services to the programming systems community, irrespective of their funding situations, and facilitate community-wide incubation to proactively nurture the software ecosystem. By actively engaging with stakeholders and employing a collaborative approach, we can collectively shape the future of programming systems and ensure a robust and thriving HPC software landscape. On May 11–12, 2023, the S4PST team conducted its inaugural kick-off workshop at the Innovative Computing Laboratory (ICL) in the University of Tennessee, Knoxville, hosted by Hartwig Anzt. The workshop encompassed various sessions dedicated to presentations and discussions, with the aim of comprehending the team members’ perspectives on the vision of software sustainability. Additionally, the workshop aimed to identify the technical, economic, and social requirements for sustaining the programming-systems community in the field of HPC. This report provides a summary of the S4PST effort by highlighting five major thrust areas discussed during the workshop: (i) community, (ii) technical support, (iii) training and diversity, (iv) verification, validation and correctness, and (v) emerging technologies. It also encompasses an overview of the presentations and discussions held throughout the event, our views and potential synergies with other seedling efforts, along with the outcomes and key takeaways from our initial discussions.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

CAFQA: A Classical Simulation Bootstrap for Variational Quantum Algorithms

Classical computing plays a critical role in the advancement of quantum frontiers in the NISQ era. In this spirit, this work uses classical simulation to bootstrap Variational Quantum Algorithms (VQAs). VQAs rely upon the iterative optimization of a parameterized unitary circuit (ansatz) with respect to an objective function. Since quantum machines are noisy and expensive resources, it is imperative to classically choose the VQA ansatz initial parameters to be as close to optimal as possible to improve VQA accuracy and accelerate their convergence on today’s devices. This work tackles the problem of finding a good ansatz initialization, by proposing CAFQA, a Clifford Ansatz For Quantum Accuracy. The CAFQA ansatz is a hardware-efficient circuit built with only Clifford gates. In this ansatz, the parameters for the tunable gates are chosen by searching efficiently through the Clifford parameter space via classical simulation. The resulting initial states always equal or outperform traditional classical initialization (e.g., Hartree-Fock), and enable high-accuracy VQA estimations. CAFQA is well-suited to classical computation because: a) Clifford-only quantum circuits can be exactly simulated classically in polynomial time, and b) the discrete Clifford space is searched efficiently via Bayesian Optimization. For the Variational Quantum Eigensolver (VQE) task of molecular ground state energy estimation (up to 18 qubits), CAFQA’s Clifford Ansatz achieves a mean accuracy of nearly 99% and recovers as much as 99.99% of the molecular correlation energy that is lost in Hartree-Fock initialization. CAFQA achieves mean accuracy improvements of 6.4x and 56.8x, over the state-of-the-art, on different metrics. Here, the scalability of the approach allows for preliminary ground state energy estimation of the challenging chromium dimer (Cr2) molecule. With CAFQA’s high-accuracy initialization, the convergence of VQAs is shown to accelerate by 2.5x, even for small molecules. Furthermore, preliminary exploration of allowing a limited number of non-Clifford (T) gates in the CAFQA framework, shows that as much as 99.9% of the correlation energy can be recovered at bond lengths for which Clifford-only CAFQA accuracy is relatively limited, while remaining classically simulable.

bayesian optimization↗