Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reduced order systems”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

2023 Critical Materials Strategy

The global effort to curb carbon emissions is accelerating demand for clean energy technologies and the materials they rely on. Demand for these materials will only continue to grow, especially as some nations aim to achieve net zero emissions by 2050. While some major materials like steel, copper, and aluminum are already powering the fossil fuel economy, others are more minor materials with potential supply risks. These risks could jeopardize the ability to reduce greenhouse gas emissions within the desirable timeframe to avoid significant climate change. In some cases, it may be necessary to take action to improve the resilience of material supply chains and mitigate supply risks. Understanding the importance of individual materials to clean energy and the supply risks associated with them is necessary to identify which materials may serve as potential roadblocks to a clean energy future. The U.S. Department of Energy (DOE) issued a series of 13 supply chain deep dive assessment reports on various energy technologies in 2022 in response to President Biden’s Executive Order on America’s Supply Chains (E.O. 14017). These reports emphasized that supply chain bottlenecks can occur at any stage of the value chain from mining and refining to component and even sub-system manufacturing. The bottlenecks are a combination of factors such as material availability, equipment availability, work force availability and quality, logistics, regulatory framework, and market conditions. These bottlenecks were worsened during the global Covid-19 pandemic. Its lingering impacts have hindered capacity expansion for material supply chains and prevented product lead-time recovery. One approach to reduce supply chain risks for the United States is to have a strong domestic manufacturing sector with a diverse set of producers. Boosting responsible domestic production would require leveraging the latest science not only in material extraction but also in developing substitutes, recycling, reuse, and remanufacturing. This report is an updated analysis of previous Critical Materials Strategy (CMS) reports published by the DOE in 2010, 2011, and 2019 based on national and global priorities, technology advancement, and technology adoption trends. Like the CMS reports, this analysis presents the results of a formal material criticality assessment to identify which materials are critical to the continued deployment of clean energy technologies globally. The analysis in this report leveraged the DOE supply chain deep dive assessments to develop the initial list of materials to evaluate. This DOE Critical Materials Assessment (CMA) is conducted independently of criticality assessments performed by other U.S. government agencies, such as that conducted by the U.S. Geological Survey (USGS). This analysis complements the USGS critical minerals determination in three aspects. First, the DOE assessment is performed from a global perspective, while the USGS analysis focusses on the importance of minerals to the U.S. economy. Second, this report focuses on the importance of materials to clean energy technologies, rather than to the economy in general. Lastly, this study is forward looking to 2035 based on clean energy deployment scenarios, whereas the USGS assessment is retrospective. Materials evaluated in this report that do not appear in the USGS Critical Minerals List include copper, uranium, electrical steel, and SiC. A draft version of this report received ~80 public comments related to supporting data and methodological improvement. Those comments have been incorporated as much as possible where appropriate. Highlights of findings from this 2023 CMA include: Rare earth materials (neodymium, praseodymium, dysprosium, and terbium) used in magnets in electric vehicle (EV) motors and wind turbine generators continue to be critical. While dysprosium (Dy) and terbium (Tb) are both heavy rare earth elements that serve the same function in magnets, the criticality of Tb is slightly lower than that for Dy in the short term due to the widespread use of Dy in high-grade magnets and Tb’s present role as a substitute. Similarly, praseodymium (Pr) is critical in the medium term but only near critical in the short term because it is more substitutable in magnets than neodymium (Nd); Materials used in batteries for EVs and stationary storage are now considered to be critical. While cobalt (Co) was found to be critical in this and previous reports, lithium (Li) becomes critical in the medium term due to its broader use in various battery chemistries and the rampant growth of the EV industry. Natural graphite is a new addition in this assessment and is also found to be critical; Platinum group metals used in hydrogen electrolyzers, such as platinum (Pr) and iridium (Ir), are critical due to an increased focus on hydrogen technologies to achieve net zero carbon emissions, while those used in catalytic converters, such as rhodium (Rh) and palladium (Pd), were screened out due to the decreased importance of catalytic converters in the medium term; Gallium (Ga) continues to be critical due to its use in light-emitting diodes (LEDs). In addition, the use of Ga has increased in magnet manufacturing and in semiconductor in forms such as gallium arsenide (GaAs) or gallium nitride (GaN); Major materials like Aluminum (Al), copper (Cu), nickel (Ni), and silicon (Si) move from noncritical in the short term to near critical in the medium term due to their importance in electrification; Electrical steel is near critical due to its use in transformers for the grid and electric motors in EVs.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

End-Use Savings Shapes Measure Documentation: Window Film

This documentation focuses on a single end-use savings shape measure - window film. The window film studied in this analysis, called solar control film, is a passive retrofit solution for windows that does not involve window replacement. This type of film is composed of transparent, tinted, or metalized laminated polyester layers and can be attached to an existing window surface (either on the exterior or interior side of the window). The properties of the window film are designed to shift the thermal and optical performances of the overall glazing system in order to serve various needs the customer would have (e.g., heat, glare). While the practical goal of purchasing and installing a window film varies widely in the real market, this study only focuses on the goal of energy savings, and highlighting the corresponding emissions. Other important aspects that customers typically consider include visual comfort, privacy, aesthetics, ultraviolet protection, etc. Thus, in practice, customers often choose a window film product not only to save energy (or cost) but also to mitigate issues around glare, excessive light, daytime privacy, or inconsistent appearance of the building. Window film products that were modeled in this analysis significantly reduced the solar heat gain coefficient of the overall glazing system, resulting in better energy savings for buildings in hot climate regions. However, significantly reducing the solar heat gain coefficient that can block unfavorable heat during the summer can actually harm blocking favorable heat during the winter. By applying window films on a stock of buildings covering various load and weather conditions, this analysis highlights when (e.g., time of day) and where (e.g., geospacial location) we can save energy with window films.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

A parallel p ‐adaptive discontinuous Galerkin method for the Euler equations with dynamic load‐balancing on tetrahedral grids

Abstract A novel p ‐adaptive discontinuous Galerkin (DG) method has been developed to solve the Euler equations on three‐dimensional tetrahedral grids. Hierarchical orthogonal basis functions are adopted for the DG spatial discretization while a third order TVD Runge‐Kutta method is used for the time integration. A vertex‐based limiter is applied to the numerical solution in order to eliminate oscillations in the high order method. An error indicator constructed from the solution of order and is used to adapt degrees of freedom in each computational element, which remarkably reduces the computational cost while still maintaining an accurate solution. The developed method is implemented with under the Charm++ parallel computing framework. Charm++ is a parallel computing framework that includes various load‐balancing strategies. Implementing the numerical solver under Charm++ system provides us with access to a suite of dynamic load balancing strategies. This can be efficiently used to alleviate the load imbalances created by p ‐adaptation. A number of numerical experiments are performed to demonstrate both the numerical accuracy and parallel performance of the developed p ‐adaptive DG method. It is observed that the unbalanced load distribution caused by the parallel p ‐adaptive DG method can be alleviated by the dynamic load balancing from Charm++ system. Due to this, high performance gain can be achieved. For the testcases studied in the current work, the parallel performance gain ranged from 1.5× to 3.7×. Therefore, the developed p ‐adaptive DG method can significantly reduce the total simulation time in comparison to the standard DG method without p ‐adaptation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Harnessing Cation Disorder for Enhancing Ionic Conductivity in Lithium Inverse Spinel Halides

Halides are promising solid-state electrolytes for all-solid-state lithium batteries due to their exceptional oxidation stability, high Li-ion conductivity, and mechanical deformability. However, their practicality is limited by the reliance on rare and expensive metals. This study investigates the Li 2 MgCl 4 inverse spinel system as a cost-effective alternative. Molecular dynamics simulations reveal that lithium disordering at elevated temperatures significantly reduces the activation energy in Li 2 MgCl 4 . To stabilize this disorder at lower temperatures, we experimentally explored the Li x Zr 1–x/2 Mg x/2 Cl 4 system and found that Zr doping induces both Zr and Li disorder at the 16c site at room temperature (RT). This leads to a 2 order-of-magnitude increase in ionic conductivity for the Li 1.25 Zr 0.375 Mg 0.625 Cl 4 composition, achieving 1.4 × 10 –5 S cm –1 at RT, compared to pristine Li 2 MgCl 4 . By deconvoluting the role of lithium vacancies and dopants, we reveal that cation disordering to the 16c site predominantly enhances ionic conductivity, whereas lithium vacancy concentration has a very limited effect.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Upgrade and characterization of the fast-channels of Imaging Neutral Particle Analyzers in the DIII-D tokamak

This publication presents the characterization of the recently upgraded fast (MHz) channels of the DIII-D Imaging Neutral Particle Analyzers. The new system allows for a phase space resolution of < 10 keV and < 10 cm. The main source of noise is due to the background radiation (neutrons and gammas) present in the diagnostic lab during experiments, which can hit the photomultipliers and generate large spikes in signal. Good linear correlations (r 2 > 0.9) are found between the average rate of noise spikes in the fast channels and the average neutron rate. Additional noise is caused by the scintillation of the optical fibers used to guide the signal to the acquisition system, which is induced by neutron and gamma impacts. A correlation factor r 2 > 0.6 between the baseline noise caused by this scintillation and the signals in the neutron diagnostic was found. Possible mitigation strategies, such as the use of shielding and filters, are discussed. These allow to reduce the noise spike rate by an order of magnitude and the noise baseline by a factor 4. Despite the noise, it is shown that the upgraded fast channels provide enough signal-to-noise ratio to observe coherent fluctuations of the energetic particle population.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Methods for Estimating Dose from Galactic Cosmic Rays and Jovian Radiation Belts on Reduced-Order Monte Carlo Geometries

This document outlines the methods and tools developed to estimate dose rates from galactic cosmic rays and the trapped radiation belts surrounding Jupiter using Monte Carlo particle transport software. Furthermore, the techniques described here are specific for translating high-resolution MAVRIC models to analogous but lower resolution Monte Carlo N-Particle Transport (MCNP) code models to estimate charged particle effects on models originally developed to study neutral particles only. This capability was initially developed to predict doses to various components of dynamic radioisotope power systems for spaceflight applications, but the tools and methods can be easily applied to other problem geometries and expanded to other space radiation environments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Facility-Level Industry Representation for Decarbonization Modeling [Slides]

The largest facilities of energy-intensive materials processing industries are disproportionate contributors to U.S. greenhouse gas (GHG) emissions. In general, industrial energy system modeling for the United States currently represents industrial demands at a much lower resolution than other end-use sector modeling. Therefore, characterizing even a subset of energy-intensive materials processing facilities will capture a significant portion of industrial GHG emissions. In order to further the development of publicly-available data to support modeling of industrial decarbonization, we summarize a set of approaches and results for characterizing the location, energy intensity and mix, process emissions intensity, and general production technology of existing clinker, ammonia, and iron and steel facilities in the United States. We also characterize facilities that represent options for reducing GHG emissions from each industry.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Thermal design analysis for SuperTruck II lithium-titanate battery pack

Here, this paper presents a systematic thermal management analysis for a new lithium-titanate-oxide battery pack to be installed in a SuperTruck II, Class 8 hybrid truck. The authors investigate the feasibility of mounting the battery pack inside the vehicle and air-cooling it with fans supplying conditioned air from the cabin. Moreover, the cells within each module are to be immersed in a heat-transfer fluid to improve temperature homogeneity. A multi-stage thermal analysis is performed to ensure adequate thermal regulation of the battery under the proposed design, considering installation and operational constraints. The authors perform thermal and electrical characterization tests on a single cell. Results inform the development of computational fluid dynamics models of the cells, modules, and pack. Module-level analysis shows that the cell temperatures can be maintained below their upper operational limit of 55 °C with module wall-to-air heat transfer coefficients between 10 and 20 Wm -2 K -1 . Pack-level analysis of distinct configurations determines effective airflow paths for adequate heat transfer and delivers a final battery pack design that achieves sufficient cooling and temperature uniformity. A reduced-order electrothermal model is developed to rapidly predict the transient battery performance and to develop a temperature control strategy.

25 ENERGY STORAGE↗
DOE OSTI2023

Automated High‐Throughput Fatigue Testing of Freestanding Thin Films

Abstract Mechanical testing at small length scales has traditionally been resource‐intensive due to difficulties with meticulous sample preparation, exacting load alignments, and precision measurements. Microscale fatigue testing can be particularly challenging due to the time‐intensive, tedious repetition of single fatigue experiments. To mitigate these challenges, this work presents a new methodology for the high‐throughput fatigue testing of thin films at the microscale. This methodology features a microelectromechanical systems‐based Si carrier that can support the simultaneous and independent fatigue testing of an array of samples. To demonstrate this new technique, the microscale fatigue behavior of nanocrystalline Al is efficiently characterized via this Si carrier and automated fatigue testing with in situ scanning electron microscopy. This methodology reduces the total testing time by an order of magnitude, and the high‐throughput fatigue results highlight the stochastic nature of the microscale fatigue response. This manuscript also discusses how this initial capability can be adapted to accommodate more samples, different materials, new geometries, and other loading modes.

Barrios, Alejandro↗
DOE OSTI2024

Metabolic flux optimization of iterative pathways through orthogonal gene expression control: Application to the β-oxidation reversal

Balancing relative expression of pathway genes to minimize flux bottlenecks and metabolic burden is one of the key challenges in metabolic engineering. This is especially relevant for iterative pathways, such as reverse β-oxidation (rBOX) pathway, which require control of flux partition at multiple nodes to achieve efficient synthesis of target products. Here, we develop a plasmid-based inducible system for orthogonal control of gene expression (referred to as the TriO system) and demonstrate its utility in the rBOX pathway. Leveraging effortless construction of TriO vectors in a plug-and-play manner, we simultaneously explored the solution space for enzyme choice and relative expression levels. Remarkably, varying individual expression levels led to substantial change in product specificity ranging from no production to optimal performance of about 90% of the theoretical yield of the desired products. We obtained titers of 6.3 g/L butyrate, 2.2 g/L butanol and 4.0 g/L hexanoate from glycerol in E. coli, which exceed the best titers previously reported using equivalent enzyme combinations. Since a similar system behavior was observed with alternative termination routes and higher-order iterations, we envision our approach to be broadly applicable to other iterative pathways besides the rBOX. Here, considering that high throughput, automated strain construction using combinatorial promoter and RBS libraries remain out of reach for many researchers, especially in academia, tools like the TriO system could democratize the testing and evaluation of pathway designs by reducing cost, time and infrastructure requirements.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Laser-Induced Enhancement of Acoustic Mode Quality Factor Revealed by Correlated Single-Particle Optical and Electron Microscopy

Acoustic modes in plasmonic nanostructures provide fundamental insights into their optomechanical behavior at the nanoscale, enabling emerging applications in plasmon-enhanced optomechanics, ultrasensitive sensing, and nanoscale energy transduction. Here we explore the modulation of acoustic phonon dynamics in lithographically fabricated gold nanodisks via laser-induced photothermal annealing. Using a correlated approach that utilizes both single-particle transient extinction spectroscopy and advanced electron microscopy, we directly link nanoscale structural transformations to changes in mechanical properties as probed through the coherence of the excited acoustic modes. Specifically, ultrafast pump–probe microscopy reveals an enhancement in the acoustic mode quality factor of annealed gold nanodisks, indicative of reduced damping and improved vibrational coherence. Structural characterization via scanning electron microscopy and electron backscatter diffraction confirms that photoinduced annealing results in smoother surface morphology and overall enhanced crystallinity. The improved crystalline order reduces defect and crystal boundary scattering, which we suggest as the reason underlying the lower quality factor before annealing. Furthermore, these findings demonstrate that targeted structural engineering at the nanoscale offers a powerful strategy for optimizing the optomechanical performance of plasmonic nanostructures, with broad implications for the design of next-generation nanophotonic and optomechanical systems.

Annealing (metallurgy)↗
DOE OSTI2024

Guided combinatorial synthesis and automated characterization expedites the discovery of hard, electrically conductive Pt x Au 1-x films

Sputter-deposited Pt-Au thin films have been reported to develop a hard, stable, nanocrystalline structure, yet little is known about how these characteristics vary with Pt x Au 1-x composition and process conditions. Toward this end, this document describes an extensive, combinatorial Pt-Au thin film library including characterized film compositions, structure, and properties. Complemented by kinematic Monte Carlo simulations of codeposition, a broad range of Pt x Au 1-x compositions (from x ~ 0.02 to 0.93) was first established by sputtering with varied magnetron powers and gun tilt angles. Further, the produced films were subsequently interrogated using automated nanoindentation, x-ray reflectivity, x-ray diffraction, atomic force microscopy, surface profilometry, four-point probe sheet resistance techniques, and wavelength dispersive spectroscopy in order to determine how hardness, modulus, density, surface roughness, structure, and resistivity vary with film stoichiometry and process parameters. Combinatorial films displayed an assortment of properties with the hardness of some films exceeding values reported previously for this material system. High hardness, high modulus, and low resistivity were generally attained when using increased deposition energy and reduced angle-of-incidence processes. Overall, the research identified promising, new Pt x Au 1-x compositions for future study and pinpointed strategies for improved deposition.

36 MATERIALS SCIENCE↗
DOE OSTI2025

LLRF System Analysis for the Fermilab PIP-II LINAC

Developing long-lived quantum processing units (QPUs) capable of supporting high-fidelity quantum operations is a crucial challenge on the path toward fault-tolerant quantum computing. TESLA-shaped superconducting RF (SRF) cavities, known for photon relaxation times on the order of seconds, provide an excellent foundation for 3D QPUs and quantum memory. This talk presents a novel design that leverages TESLA cavity modes coupled to ancillary transmon qubits, optimized to preserve coherence and control. By carefully engineering the package geometry, optimizing Hamiltonian parameters, and minimizing lossy participation ratios, we achieve photon relaxation times of over 16 ms and 20 ms for the two cavity modes, representing a significant improvement over previous multimode quantum memories. Despite the reduced coupling between the qubit and cavity modes, which is necessary to preserve long lifetimes, the platform supports robust and universal control schemes that are not limited by low coupling strength. We will also discuss how this architecture can lead to scalable, modular quantum computing systems.

Varghese, P. [Fermilab]↗
DOE OSTI2023

A Review of Distributed Secondary Control Architectures in Islanded-Inverter-Based Microgrids

The increasing energy demand, the shortage of energy resources, and the environmental challenges faced by conventional power-generation systems are some of the ongoing challenges faced by modern power systems. Therefore, many efforts have been made by the scientific community to develop comprehensive solutions to overcome these issues. For instance, current technological advances have allowed the integration of distributed generators into the power systems, promoting the use of microgrids to overcome these issues. However, the use of renewable distributed generators have introduced new challenges to the traditional control system schemes. To overcome these challenges, a hierarchical control approach has been proposed for distributed renewable sources. In other words, the control scheme have been divided into three hierarchical levels, primary, secondary, and tertiary, to overcome the new challenges present in modern power systems. Due to extensiveness of this topic, this overview is focused on secondary control systems, mainly for AC isolated microgrids. To improve the power quality of modern systems, several secondary control schemes have been proposed to overcome the well-known problem of frequency and voltage deviation. Some of these schemes have also introduced adequate active/reactive power sharing techniques to optimize the utilization of resources. Additionally, other secondary control schemes have also focused on reducing the communication load, to lower the network cost and adding robustness against communication problems. This article presents an insight of the different control techniques used to overcome power quality and communication problems. A comprehensive overview of distributed secondary control techniques for islanded microgrids is presented. In addition, the implementation of these techniques is explained in an orderly and sequential manner.

14 SOLAR ENERGY↗
DOE OSTI2024

Cobalt-60, Barium-133, Cesium-137, and Europium-152 migration from cementitious sources through sediment under field conditions

Safe and effective storage of radioactive waste is essential to protect human and environmental health. Due to the potential for accidental releases and the severity of the associated risks, it is imperative to further understand radionuclide transport should an accident occur. This study was the second set of measurements conducted in 2022 of an ongoing experiment that has analyzed the vadose zone migration of radionuclides from cementitious wasteforms at the Savannah River Site over the last ten years. The radionuclides introduced within the sources are prominent constituents of radioactive waste or analogs for other groups or series of radionuclides. Lysimeters were first analyzed in 2016 using a collimated high-purity germanium gamma-ray spectrometer to non-destructively measure the concentration of each radionuclide in the sediment column as a function of depth. Following these measurements, the lysimeters were redeployed for another 4 years. All radionuclides in all lysimeters were observed to transport further during the redeployment period; however, the extent of migration varied with the material used for introduction. Except for 137Cs, migration through the sediment control system increased with decreasing ionic potential (ionic charge/radius); migration order: 152Eu<137Cs<60Co<133Ba. Overall, the cementitious wasteforms were observed to decrease radionuclide migration extent relative to natural vadose zone conditions. In both cementitious wasteforms, the migration extent increased in the order 152Eu<133Ba<60Co<137Cs. However, less migration was measured when the radionuclides were incorporated into a reducing grout wasteform. The novelty of this paper is the demonstration of a technique capable of creating non-destructive measurements over decade time scales. Ultimately, this work provides insight into the long-term migration of alkali, alkali earth, divalent transition metal, and trivalent (e.g., lanthanide and actinide element) isotopes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Evaluation of Minimum NOx Emission From Ammonia Combustion

Abstract Ammonia (NH3) is being explored as a hydrogen carrier with no carbon emissions. However, if burned directly as NH3, rather than being completely decomposed back to N2/H2, the fuel-bound nitrogen comes with a potentially significant NOx emissions penalty. Indeed, several existing studies are showing ammonia combustion NOx emissions that exceed current natural gas fueled, DLN technologies by one to two orders of magnitude. Therefore, it is important to establish the theoretical minimum NOx emissions for an ammonia combustor, to determine how much NOx levels can be reduced via further technology development. In other words, the purpose of this work is not to analyze the performance of a specific combustor but, rather, the fundamental limits of what is achievable. This study quantifies this minimum NOx level for a two-stage combustor system for a given combustor exit temperature and residence time, with a constraint on unburned fuel levels. As expected, the optimum configuration is a rich front end combustor to burn and crack ammonia with significant H2 production, followed by an NO relaxation reactor, followed by a lean stage that consumes the remaining H2. The optimum residence time and stoichiometry of each zone are determined in the fast mixing limit, which essentially balances between NOx production in the primary and secondary zones. These results show minimum NOx levels are in 200–400 ppm range at 1 bar, but drop to levels of ∼25 ppm at 20 bar. These NOx emissions are dominated by NOx production in the primary stage which relaxes to equilibrium levels quite slowly. As processes controlling NOx relaxation to equilibrium in the primary stage dominate overall NO emission levels, combustor NOx sensitivities are essentially opposite that of natural gas fired, DLN systems. Specifically, NOx values drop with increased combustor residence time, increased pressure, and increased combustor exit temperature. These results also suggest that the most important strategy for NOx minimization is to provide sufficient relaxation time after the primary zone for NOx to approach equilibrium—this can be done via kinetic means to accelerate this relaxation rate, such as enhancing pressure or temperature, or increasing residence times. Indeed, this work shows that low pressure combustors specifically optimized for ammonia will have residence times that are one to two orders of magnitude larger than current natural gas systems. By doing so, NOx levels below 10 ppm may be achievable. Finally, we discuss the sensitivity of these values to uncertainties in ammonia kinetics.

Engineering↗
DOE OSTI2024

SPARC-X: Quantum simulations at extreme scale - reactive dynamics from first principles

We have developed the massively parallel electronic structure code SPARC-X: a computational framework for performing Kohn-Sham Density Functional Theory (DFT) calculations that can scale linearly with the number of atoms in the system, while being able to leverage petascale and emerging exascale parallel computers to study chemical phenomena at unprecedented length and time scales. SPARC-X exploits a recent breakthrough in electronic structure methodologies: systematically improvable, strictly local, orthonormal, discontinuous real-space bases that efficiently and systematically capture the local chemistry of the system. With further adaptation using new machine-learning techniques and the use of the massively parallel Spectral Quadrature (SQ) electronic structure method, the algorithmic complexity and prefactor associated with DFT calculations involving semilocal as well as hybrid functionals are dramatically reduced. Using petascale computational resources, SPARC-X enables quantum mechanical simulations at length and time scales previously accessible only by empirical approaches, e.g., 1,000,000 atoms for a few picoseconds using semilocal functionals or 1,000 atoms for a few picoseconds using hybrid functionals. Using exascale resources, the sizes and times targeted are two orders of magnitude larger. Such a capability has applications in a wide variety of chemical sciences, including reactive interfaces where large length- and/or long time-scales are needed and traditional force fields fail. This is particularly important in dynamic catalysis, where bond breaking and formation must be understood in detail. We developed, tested, and employed the SPARC-X framework to understand the photocatalytic properties of TiO 2 nanoparticles, revealing finite size effects that cannot be captured with standard model systems or functionals. This integrated development and application strategy ensures that SPARC-X remains a robust, efficient, and scalable software package for quantum simulations on current petascale and emerging exascale computing resources.

97 MATHEMATICS AND COMPUTING↗