Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computer systems organization”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Tsdat: An Open-Source Data Standardization Framework for Marine Energy and Beyond: Preprint

Many organizations are tasked with the collection and processing of large quantities of data from various measurement devices. Data reported from these sources are often not interoperable with datasets and software used by analysts and other organizations in the same field, introducing barriers for collaboration on large-scale projects. This poses a particular problem for cross-device comparisons and machine learning applications. To address these challenges, the open source Time-Series Data Pipelines (Tsdat) software was developed by a joint collaboration between Pacific Northwest National Laboratory, the National Renewable Energy Laboratory, and Sandia National Laboratories to facilitate collaboration and accelerate advancements in the Marine Energy domain through the development of an open-source ecosystem of tools. This paper will describe the Tsdat software and the data standards within which the framework operates. A beta version of the framework has been released and is currently being used by several projects in marine energy, wind energy, and building energy systems.

big data↗
DOE OSTI2021

Gas-phase ion-molecule interactions in a collision reaction cell with triple quadrupole-inductively coupled plasma mass spectrometry: Investigations with N 2 O as the reaction gas

Nitrous oxide (N 2 O) was used as a reaction gas to investigate the gas phase ion-molecule interactions using the Agilent 8900 QQQ-ICP-MS. A multi-element standard containing 45 elements with masses ranging from 9 to 208 u was measured in the presence and absence of N 2 O. The main product ion species observed were oxides and nitrides. Comparison of the N 2 O reaction results with similar measurements conducted with O 2 revealed that N 2 O was more effective at forming oxides in general: the elements Cd and Pb were shown to produce oxides with N 2 O where the reaction did not occur with O 2 . Nitrous oxide was also shown to produce a significant amount of nitride species in a few cases. The general reactivity was shown to be consistent with density functional theory (DFT)-predicted reaction enthalpies, such that all predicted exothermic reactions produced product ions at levels at least 1% of the unreacted ion. Our results show that reaction enthalpy is a reasonable predictor of reactivity with N2O on the timescales of the interactions in non-thermal ICP-MS/MS systems. Our work demonstrates the utility of two relatively new platforms (commercial elemental ICP-MS/MS and EMSL Arrows interface to the NWChem program suite), which allows for the study of a large number of elements within a short period. While DFT with the basis sets utilized here is not the most accurate computational method, it is also not computationally expensive and is shown to be suitable for predicting gas phase reactivity in the QQQ-ICP-MS for the majority of ions studied. Here, the ease and rapidity of data collection and DFT calculations has the potential to be very impactful for the identification of targeted reaction chemistries to be leveraged for analytical method development, such as for the inline separation of isobaric interferences from analytes of interest.

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

HIPPO – A Software Platform for Electricity Market Research and Development

The goal of this project is to provide Regional transmission organizations (RTOs) and independent system operators (ISOs) a market design and prototyping software, High-Performance Power-Grid Optimization (HIPPO), that they can evaluate electricity market design options, calculate market planning strategies and operational performance. With the high standards and strict reliability requirements for operating power systems, impacts of new technologies need to be fully investigated prior to any consideration for adoption. A market design and prototyping software tool which can be used to prototype electricity market design options, to calculate market planning strategies and operational performance with high precision, and to investigate the impacts for integrating future power grid technologies will be valuable to RTOs/ISOs who operate power systems, to vendors like GE and ABB who provide the market solvers, and to market participants and researchers who are actively doing market research. HIPPO is a such tool that can be used to improve the current market operations and provide capabilities for rigorous forward-looking design and prototyping of next-generation energy markets. HIPPO has a high-resolution model for the day-ahead SCUC, which was validated with MISO and GE-Grid Solutions. HIPPO is built with parallel and distributed computing capabilities and can be executed in both multi-thread and high-performance computing (HPC) settings. This capability provides fast solution speed necessary to handle the larger and more complex SCUC problems of real-world cases and the potentially growing size and complexity of future scenarios. In addition, HIPPO has a concurrent optimizer (CO) which manages multiple algorithm executions simultaneously and leverages the advantages from different algorithms. This structure provides flexibility to better benchmark competing approaches. Highly accurate market model, fast solution technologies and flexible model and algorithm control are the features which will make HIPPO an extensible platform for developing and testing multiple approaches to meet a wide range of future market needs.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Newton-Okounkov bodies of chemical reaction systems

Despite their noted potential in polynomial-system solving, there are few concrete examples of Newton-Okounkov bodies arising from applications. Accordingly, in this paper, we introduce a new application of Newton-Okounkov body theory to the study of chemical reaction networks and compute several examples. An important invariant of a chemical reaction network is its maximum number of positive steady states Here, we introduce a new upper bound on this number, namely the ‘Newton-Okounkov body bound’ of a chemical reaction network. Further, through explicit examples, we show that the Newton-Okounkov body bound of a network gives a good upper bound on its maximum number of positive steady states. We also compare this Newton-Okounkov body bound to a related upper bound, namely the mixed volume of a chemical reaction network, and find that it often achieves better bounds.

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

Enhanced Twist-Averaging Technique for Magnetic Metals: Applications Using Quantum Monte Carlo

In this study, we propose an improved twist-averaging (TA) scheme for quantum Monte Carlo methods that use converged Kohn–Sham or Hartree–Fock orbitals as the reference. This TA technique is tailored to sample the Brillouin zone of magnetic metals, although it naturally extends to nonmagnetic (NM) conducting systems. The proposed scheme aims to reproduce the reference magnetization and achieves charge neutrality by construction, thus avoiding the large energy fluctuations and the postprocessing needed to correct the energies. It shows the most robust convergence of total energy and magnetism to the thermodynamic limit (TDL) when compared to four other TA schemes. Diffusion Monte Carlo applications are shown on NM Al and ferromagnetic α-Fe. The cohesive energy of Al in the TDL shows an excellent agreement with the experimental result. Furthermore, the magnetic moments in α-Fe exhibit rapid convergence with an increasing number of twists.

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

Forces from Stochastic Density Functional Theory under Nonorthogonal Atom-Centered Basis Sets

We develop a formalism for calculating forces on the nuclei within the linear-scaling stochastic density functional theory (sDFT) in a nonorthogonal atom-centered basis set representation (Fabian et al. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2019, 9, e1412, 10.1002/wcms.1412) and apply it to the Tryptophan Zipper 2 (Trp-zip2) peptide solvated in water. We use an embedded-fragment approach to reduce the statistical errors (fluctuation and systematic bias), where the entire peptide is the main fragment and the remaining 425 water molecules are grouped into small fragments. We analyze the magnitude of the statistical errors in the forces and find that the systematic bias is of the order of 0.065 eV/Å (~1.2 × 10 -3 E h / a 0 ) when 120 stochastic orbitals are used, independently of system size. This magnitude of bias is sufficiently small to ensure that the bond lengths estimated by stochastic DFT (within a Langevin molecular dynamics simulation) will deviate by less than 1% from those predicted by a deterministic calculation.

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

Field-Driven Simulations to Probe the Impact of Ionic Correlations on Solution Transport Coefficients in Binary, Ternary, and Reciprocal Quaternary Aqueous Electrolytes

Ionic correlations play a critical role in governing transport properties of mixed-salt aqueous electrolyte solutions, yet their quantitative characterization remains challenging, particularly for multicomponent aqueous systems. Here, we develop a general nonequilibrium molecular dynamics framework to efficiently compute Onsager transport coefficients and ionic correlations in mixed-salt aqueous solutions. Using field-driven simulations, we obtain accurate Onsager matrices for LiCl/ KCl, KCl/KBr, and LiBr/KCl electrolyte solutions with significantly reduced computational cost relative to equilibrium Green−Kubo methods. The framework enables direct assessment of how attractive cation−anion and repulsive like-ion interactions contribute to conductivity and salt diffusivity across compositions. While static ion pairing exhibits strong composition dependence, dynamic ion correlations remain nearly invariant, leading to constant deviations from Nernst−Einstein predictions. These results highlight the disconnect between static ion association and dynamic transport correlations, and they establish a transferable approach for analyzing ion transport in complex electrolyte environments relevant to separation processes and electrochemical systems.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Unraveling Local Structure of Molten Salts via X-ray Scattering, Raman Spectroscopy, and Ab Initio Molecular Dynamics

In this work, we resolve a long-standing issue concerning the local structure of molten MgCl 2 by employing a multimodal approach, including X-ray scattering and Raman spectroscopy, along with the theoretical modeling of the experimental spectra based on ab initio molecular dynamics (AIMD) simulations utilizing several density functional theory (DFT) methods. We demonstrate the reliability of AIMD simulations in achieving excellent agreement between the experimental and simulated spectra for MgCl 2 and 50 mol % MgCl 2 + 50 mol % KCl, and ZnCl 2 , thus allowing structural insights not directly available from experiment alone. A thorough computational analysis using five DFT methods provides a convergent view that octahedrally coordinated magnesium in pure MgCl 2 upon melting preferentially coordinates with five chloride anions to form distorted square pyramidal polyhedra that are connected via corners and to a lesser degree via edges. This is contrasted with the results for ZnCl 2 , which does not change its tetrahedral coordination on melting. Although the five-coordinate MgCl 5 3– complex was not considered in the early literature, together with an increasing tendency to form a tetrahedrally coordinated complex with decreasing the MgCl 2 content in the mixture with alkali metal chloride systems, current work reconciles the results of most previous seemingly contradictory experimental studies.

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

Competitive Valerate Binding Enables RuO 2 -Mediated Butene Electrosynthesis in Water

The (non)-Kolbe oxidation of valeric acid, sourced from a hydrolysis product of cellulose, provides a sustainable synthetic route to access value-added products, such as butene. An essential mechanistic step preceding product formation involves the oxidative and decarboxylative cleavage of a C–C bond. Yet, the role of the electrode surface in mediating this oxidative step remains an open question: the electron transfer can occur either via an inner-sphere or outer-sphere mechanism. Here, we report the electrochemical, in situ spectroscopic, computational, and reactivity studies of RuO 2 -mediated oxidative decarboxylation of valeric acid to butene in aqueous electrolytes. We find that carboxylates bind to RuO 2 anode surfaces at potential values where decarboxylation products are observed. Our results are consistent with a reaction scheme where the competitive and catalytic oxygen evolution reaction (OER) is impeded by these bound carboxylate species while these species are inert toward butene formation. Our results implicate an outer-sphere electron transfer mechanism for decarboxylation where the surface chemistry of the RuO 2 electrode serves to enable higher non-Kolbe reaction selectivity by suppressing the parasitic OER. Furthermore, our findings delineate interfacial design principles for selective electrochemical systems that utilize water as the ultimate oxidant for sustainable decarboxylation.

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

TAMM: Tensor algebra for many-body methods

Tensor algebra operations such as contractions in computational chemistry consume a significant fraction of the computing time on large-scale computing platforms. The widespread use of tensor contractions between large multi-dimensional tensors in describing electronic structure theory has motivated the development of multiple tensor algebra frameworks targeting heterogeneous computing platforms. In this paper, we present Tensor Algebra for Many-body Methods (TAMM), a framework for productive and performance-portable development of scalable computational chemistry methods. TAMM decouples the specification of the computation from the execution of these operations on available high-performance computing systems. With this design choice, the scientific application developers (domain scientists) can focus on the algorithmic requirements using the tensor algebra interface provided by TAMM, whereas high-performance computing developers can direct their attention to various optimizations on the underlying constructs, such as efficient data distribution, optimized scheduling algorithms, and efficient use of intra-node resources (e.g., graphics processing units). The modular structure of TAMM allows it to support different hardware architectures and incorporate new algorithmic advances. We describe the TAMM framework and our approach to the sustainable development of scalable ground- and excited-state electronic structure methods. We present case studies highlighting the ease of use, including the performance and productivity gains compared to other frameworks.

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

Computational Methods in Heterogeneous Catalysis

The unprecedented ability of computations to probe atomic-level details of catalytic systems holds immense promise for the fundamentals-based bottom-up design of novel heterogeneous catalysts, which are at the heart of the chemical and energy sectors of industry. Here, we critically analyze recent advances in computational heterogeneous catalysis. First, we will survey the progress in electronic structure methods and atomistic catalyst models employed, which have enabled the catalysis community to build increasingly intricate, realistic, and accurate models of the active sites of supported transition-metal catalysts. We then review developments in microkinetic modeling, specifically mean-field microkinetic models and kinetic Monte Carlo simulations, which bridge the gap between nanoscale computational insights and macroscale experimental kinetics data with increasing fidelity. Here, we finally review the advancements in theoretical methods for accelerating catalyst design and discovery. Throughout the review, we provide ample examples of applications, discuss remaining challenges, and provide our outlook for the near future.

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

How the Hydrophobic Interface between a Perfluorosulfonic Acid Polymer and Water Vapor Controls Membrane Hydration

Stable hydration in perfluorinated polyelectrolyte membranes such as Nafion is essential to maintain good ion conductivity and manage permeation, especially in vapor-fed devices where water content depends on relative humidity in a gas stream. Extensive studies in the literature have shown that Nafion hydration in water vapor is controlled by its interfacial transport resistance. Nafion forms a fluorine-rich layer at the polymer-gas interface, and it has been proposed that this layer blocks water transport due to its hydrophobicity. To develop a molecular-level description of the physics underlying transport resistance in this system, we have performed a computational reaction-diffusion kinetics study of water evaporation from Nafion. Two distinct models are examined, one mimicking the blocking function proposed in the literature and the other assuming that there is no blocking, treating instead water evaporation as a dynamic balance between uptake from the gas and desorption from the polymer surface. Simulation results are compared to time-dependent infrared data over a range of 100-0% relative humidity from the literature. Only the dynamic model successfully reproduces experimental observations. This indicates that the physical nature of interfacial transport resistance is not slow diffusion across an interfacial layer; rather, it is due to the competition between dehydration and rehydration. The simulation data provide details on the accompanying water distributions throughout the membrane and on interfacial kinetics, showing that they are characterized by strong fluctuations.

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

Solid–liquid separation of lignocellulosic sugars from biomass by rotating ceramic disc filtration

In many biomass conversion processes, the separation of cellulosic sugars from residual, lignin-rich solids is a critical step, and achieving high recovery yields of sugars by conventional tangential crossflow and vacuum filtration is challenged by the presence of biomass solids, which rapidly foul filters, resulting in decreased throughput. Considering the performance limitations of traditional filtration methods, dynamic filtration, which generates high shear at the membrane surface to decrease fouling, is emerging as a viable alternative for demanding solid–liquid separations. For high solids separations, there is little available information regarding the performance, limitations, and energy consumption of dynamic filtration. To that end, here we characterized the performance of a dynamic filtration module, specifically a rotating ceramic disc (RCD) filter, for the aseptic recovery of cellulosic sugars from biomass solids following pretreatment and enzymatic hydrolysis. We show how RCD rotational velocity and percent biomass solids impact the filter throughput. Additionally, we used computational fluid dynamics (CFD) simulations to estimate the shear rate at the membrane surface and to visualize hydrodynamic profiles within the module. With the combined CFD simulations and experimental results, we estimated the energy demand and operating expenses for a viable dynamic filtration system operating with a lignocellulosic feed slurry. Our results indicate that an RCD filter can achieve ≥95% recovery of sugars and produce a retentate slurry containing 12 wt% insoluble solids with low energy consumption (a 2.2-fold improvement over cross-flow filtration) and low operating costs ($\$$0.06 per kg sugars). These results show a viable path towards operationally reliable, energy efficient, and cost-effective separations of sterilized cellulosic sugars from biomass solids and highlight the potential of dynamic filtration systems for challenging solid–liquid separations.

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

University Data Management Pilot Utilizing the Nuclear Research Data System

Background In 2022, the Office of Science and Technology Policy (OSTP) issued a memo that significantly reshaped the landscape of access to federally funded research. The memo mandated that all taxpayer-funded research be made available to the public without delay upon publication, without an embargo period, superseding the 2013 OSTP public access policy. This public access policy promotes transparency and the democratization of knowledge, ensuring that the fruits of scientific endeavors funded by federal agencies could be immediately accessed and built upon by scientists, educators, students, and the public at large. To implement the requirements of the OSTP guidance and DOE Public Access Plan, the Office of Nuclear Energy (NE) has implemented public access plan guidance and has identified several areas where better data management practices would further expand public access to important nuclear energy related scientific data, reports, and other technical products. Significant NE supported efforts are already underway for data management and public access to important nuclear energy related data.1 2 To address gaps in data management practices, and improve retention and accessibility of data, NE is actively exploring enhanced data management options utilizing its high-performance computing resources administered by its Nuclear Scientific User Facility Program. A newly piloted system, the Nuclear Research Data System (NRDS) acts as a portal for data collection and dissemination. Nuclear Energy University Program Research and Development Portfolio According to Web of Science, NEUP has produced 2,345 journal publication that have been cited more than 61,000 times3 and countless conference proceedings. These publications are publicly available through OSTI.gov and in the open literature. Additional scientific and technical products including project milestones that are not publications and NEUP project final reports are vetted through OSTI.gov and released once reviewed and approved by DOE. Since 2009, NEUP has awarded close to 1,000 different R&D projects in technical areas across the NE research programs. As of June 2023, 512 NEUP reports are publicly available on OSTI. The underlying data for projects is still held at universities, and data transfer, co-location, and dissemination has not occurred in a systematic way. NEUP data is currently accessible through myriad university-based data repositories, or through direct requests to PIs. The program identified this patchwork of repositories, or often lack of publicly available data, as a significant barrier to an organized, accessible, and comprehensive solution to sharing data with the larger nuclear energy community. Approach The goal of this pilot project is to establish a pathway to a consolidated long-term repository for NEUP project data. To accomplish this goal, the pilot strives to accomplish the following objectives: Establish data collection standards, including a standard set of required supplementary information to contextualize and support raw data files. Work with the HPC group collect and upload information and to modify the NRDS system, as needed, to support a standardized approach. Resolve potential barriers to successful roll out of an expanded data collection strategy, including modifying data management plan guidelines and establishing a document and data release process that accounts for potential intellectual property and/or export control concerns. Results Overall, the pilot was successful in collecting 8,982 raw and processes data files, 220 reports, 56 calibration files, and 5,931 other supplementary documents. Supplementary documents included experimental plans, methods, journal publications and conference proceedings, milestone reports, and final reports. Figure 2 shows the number of data sets and supplementary project information provided by each project. Projects has significantly different input, depending on experimental data produced and completeness of the datasets provided.

Data collection↗
DOE OSTI2026

Model Parameter Development for Complex Materials: Species-Specific Diffusion Barriers in 316 Stainless Steel from Systematic DFT Calculations

Vacancy-mediated diffusion barriers in 316 stainless steel have been systematically calculated using density functional theory to provide essential parameters for mesoscale microstructure evolution models. A statistical sampling approach employing 210 nudged elastic band calculations across multiple special quasi-random structures captures the effects of local chemical environments in this concentrated alloy. The computational methodology addresses challenges specific to chemically disordered systems, including proper magnetic treatment throughout multi-step calculations and validation against experimental structural properties. The calculated activation barriers reveal clear species-dependent diffusion behavior with the hierarchy Ni >> Fe ˜ Cr >> Mo. Nickel exhibits the highest barriers (0.74–1.31 eV, mean 1.045 eV), confirming its role as the slowest-diffusing major component. Iron and chromium show similar moderate barriers averaging 0.587 eV and 0.522 eV, respectively. Remarkably, molybdenum demonstrates exceptionally low barriers (0.12–0.28 eV, mean 0.194 eV), suggesting much higher mobility than previously recognized and potentially significant implications for precipitation kinetics and microstructure evolution. The barrier ranges remain consistent across different 316 SS compositions, supporting parameter transferability for modeling applications. The overall mean barrier of 0.64 eV provides a practical approximation for phase field simulations, while species-specific values enable detailed treatments of diffusion-controlled processes. This systematic approach establishes a validated framework for generating diffusion parameters in other concentrated alloys where experimental data are limited, while providing the first systematic set of species-specific barriers for predictive modeling of 316 stainless steel microstructure evolution.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effect of Methyl Groups on Formation of Ordered or Layered Graphitic Materials from Aromatic Molecules

Abstract Developing functionally complex carbon materials from small aromatic molecules requires an understanding of how the chemistry and structure of its constituent molecules evolve and crosslink, to achieve a tailorable set of functional properties. Here, molecular dynamics (MD) simulations are used to isolate the effect of methyl groups on condensation reactions during the oxidative process and evaluate the impact on elastic modulus by considering three monodisperse pyrene‐based systems with increasing methyl group fraction. A parameter to quantify the reaction progression is designed by computing the number of new covalent bonds formed. Utilizing the previously developed MD framework, it is found that increasing methylation leads to an almost doubling of bond formation, a larger fraction of the new bonds oriented in the direction of tensile stress, and a higher basal plane alignment of the precursor molecules along the direction of tensile stress, resulting in enhanced tensile modulus. Additionally, via experiments, it is demonstrated that precursors with a higher fraction of methyl groups result in a higher alignment of molecules. Moreover, increased methylation results in the lower spread of single molecule alignment which may lead to smaller variations in tensile modulus and more consistent properties in carbon materials derived from methyl‐rich precursors.

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

Electronic Structure and Bonding of US, SUO, and US 2

Anion photoelectron spectra of US – and US 2 – were recorded using the third (355 nm) and fourth (266 nm) harmonics of an Nd:YAG laser, which yielded vertical detachment energies (VDEs) of 1.71 and 2.02 eV, respectively. The experimental results are supported by extensive relativistic ab initio calculations, primarily at the coupled cluster level of theory, with systematic sequences of correlation consistent basis sets. Calculations include the closely related SUO and SUO – molecules, as well as the oxide congeners UO/UO – and UO 2 /UO 2 – which are well-known experimentally and provide benchmark systems for the sulfide calculations. Adiabatic electron detachment energies (ADEs) are computed for UO – , UO 2 – , US – , SUO – , and US 2 – using the Feller–Peterson–Dixon (FPD) composite approach. Additionally, ADEs are determined for UO – and US – using a spinor-based coupled cluster approach where spin–orbit coupling is included at the orbital level. VDEs are derived from the ab initio results from Franck–Condon simulations of the photoelectron spectra.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based Ab Initio Molecular Dynamics. 1. Theory, Algorithm, and Performance

By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal density functional theory (DFT) and thereby furnish a more accurate and reliable description of the underlying electronic structure in systems throughout biology, chemistry, physics, and materials science. However, the high computational cost associated with the evaluation of all required EXX quantities has limited the applicability of hybrid DFT in the treatment of large molecules and complex condensed-phase materials. To overcome this limitation, we describe a linear-scaling approach that utilizes a local representation of the occupied orbitals (e.g., maximally localized Wannier functions (MLWFs)) to exploit the sparsity in the real-space evaluation of the quantum mechanical exchange interaction in finite-gap systems. In this work, we present a detailed description of the theoretical and algorithmic advances required to perform MLWF-based ab initio molecular dynamics (AIMD) simulations of large-scale condensed-phase systems of interest at the hybrid DFT level. We focus our theoretical discussion on the integration of this approach into the framework of Car–Parrinello AIMD, and highlight the central role played by the MLWF-product potential (i.e., the solution of Poisson’s equation for each corresponding MLWF-product density) in the evaluation of the EXX energy and wave function forces. We then provide a comprehensive description of the exx algorithm implemented in the open-source Quantum ESPRESSO program, which employs a hybrid MPI/OpenMP parallelization scheme to efficiently utilize the high-performance computing (HPC) resources available on current- and next-generation supercomputer architectures. Furthermore, this is followed by a critical assessment of the accuracy and parallel performance (e.g., strong and weak scaling) of this approach when AIMD simulations of liquid water are performed in the canonical (NVT) ensemble. With access to HPC resources, we demonstrate that exx enables hybrid DFT-based AIMD simulations of condensed-phase systems containing 500–1000 atoms (e.g., (H₂O)₂₅₆) with a wall time cost that is comparable to that of semilocal DFT. In doing so, exx takes us one step closer to routinely performing AIMD simulations of complex and large-scale condensed-phase systems for sufficiently long time scales at the hybrid DFT level of theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗