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

Device-Centric Ransomware Detection using Machine Learning-Based Memory Forensics for Smart Inverters

Ransomware attacks are the fastest-growing form of cyberattacks worldwide. Recently, ransomware attacks have targeted industrial control systems (ICSs), including power grids. Lessons learned from recent incidents in ICSs show that ransomware groups can deliver ransomware into not only the organization’s control servers, but also the operational technology (OT) devices such as smart inverters and smart grid devices. This paper proposes a machine learning (ML)- based memory forensics method enabling the detection of ransomware binaries stored in the memory of a commercial smart inverter. Device firmware binary files are extracted from a Serial Peripheral Interface (SPI) flash memory, and samples of both benign and ransomware binaries are generated by a binary manipulation method and a real-world ransomware encryption, separately. A deep transfer learning (DTL) method is used to retrain a convolutional neural network (CNN)-based ransomware detection algorithm using the generated samples. The experimental result validates that the proposed ML-based memory forensics method can accurately detect ransomware files.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Defect Energetics in Pseudo-Cubic Mixed Halide Lead Perovskites from First-Principles

Owing to the increasing popularity of lead-based hybrid perovskites for photovoltaic (PV) applications, it is crucial to understand their defect energetics and its influence on their optoelectronic properties. In this work, we simulate various point defects in pseudocubic structures of mixed iodide-bromide and bromide-chloride methylammonium lead perovskites with the general formula MAPbI 3-y Br y or MAPbBr 3-y Cl y (where y is between 0 and 3), and use first-principles based density functional theory computations to study their relative formation energies and charge transition levels. We identify vacancy defects and Pb on MA antisite defect as the lowest energy native defects in each perovskite. Here, we observe that while the low energy defects in all MAPbI 3-y Br y systems only create shallow transition levels, the Br or Cl vacancy defects in the Cl-containing pervoskites have low energy and form deep levels which become deeper for higher Cl content. We examine the structures and density of states of pure and defect-containing perovskite systems to obtain an understanding of the nature of defect levels. Further, we study extrinsic substitution by different elements at the Pb site in MAPbBr 3 , MAPbCl 3 , and the 50-50 mixed halide perovskite, MAPbBr 1.5 Cl 1.5 , and identify some transition metals that create lower energy defects than the dominant intrinsic defects and also create midgap charge transition levels.

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

Perspectives on the integration between first-principles and data-driven modeling

Efficiently embedding and/or integrating mechanistic information with data-driven models is essential if it is desired to simultaneously take advantage of both engineering principles and data-science. Further the opportunity for hybridization occurs in many scenarios, such as the development of a faster model of an accurate high-fidelity computer model; the correction of a mechanistic model that does not fully-capture the physical phenomena of the system; or the integration of a data-driven component approximating an unknown correlation within a mechanistic model. At the same time, different techniques have been proposed and applied in different literatures to achieve this hybridization, such as hybrid modeling, physics-informed Machine Learning (ML) and model calibration. In this paper we review the methods, challenges, applications and algorithms of these three research areas and discuss them in the context of the different hybridization scenarios. Moreover, we provide a comprehensive comparison of the hybridization techniques with respect to their differences and similarities, as well as advantages and limitations and future perspectives. Finally, we apply and illustrate hybrid modeling, physics-informed ML and model calibration via a chemical reactor case study.

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

Toward Enhancing Wastewater Treatment with Resource Recovery in Integrated Assessment and Computable General Equilibrium Models

Sustainable water management is essential to increasing water availability and decreasing water pollution. The wastewater sector is expanding globally and beginning to incorporate technologies that recover nutrients from wastewater. Nutrient recovery increases energy consumption but may reduce the demand for nutrients from virgin sources. We estimate the increase in annual global energy consumption (1,100 million GJ) and greenhouse gas emissions (84 million t CO2e) for wastewater treatment in the year 2030 compared to today’s levels to meet sustainable development goals. To capture these trends, integrated assessment and computable general equilibrium models that address the energy-water nexus must evolve. We reviewed 16 of these models to assess how well they capture wastewater treatment plant energy consumption and GHG emissions. Only three models include biogas production from the wastewater organic content. Four explicitly represent energy demand for wastewater treatment, and eight include explicit representation of wastewater treatment plant greenhouse gas emissions. Of those eight models, six models quantify methane emissions from treatment, five include representation of emissions of nitrous oxide, and two include representation of emissions of carbon dioxide. Our review concludes with proposals to improve these models to better capture the energy-water nexus associated with the evolving wastewater treatment sector.

42 ENGINEERING↗
DOE OSTI2021

Preferential Direction of Electron Transfers at a Dye–Metal Oxide Interface with an Insulating Fluorinated Self-Assembled Monolayer and MgO

Slowing non-desirable electron transfer reactions at metal oxide-dye interfaces is important for many technologies. Specifically, after an interfacial photoinduced charge separation event at a metal oxide-dye interface, it is critically important to limit the rate of electron transfer reactions back to sensitizers and to limit electron transfer reactions between the electrolyte and the metal oxide. Ruthenium-based dyes at metal oxide interfaces are widely used in many fields; however, these dyes often have poor surface insulation resulting in fast recombination kinetics with transition metal-based redox shuttles in an electrolyte. This work explores two semiconductor surface modification strategies designed to minimize recombination events of electrons in TiO 2 with oxidized redox shuttles using a fluorinated siloxane insulator (PFTS) and a metal oxide insulator (MgO) with a well-known Ru dye, B11. Additionally, the influence of these treatments on the rate and duration of photoinduced interfacial charge separation at the TiO 2 -dye interface was examined. The TiO 2 -dye-redox shuttle systems were studied via dye-sensitized solar cell current-voltage curve, incident photon-to-current conversion efficiency, small-modulated photovoltage transient, time-correlated single photon counting, and transient absorption spectroscopy measurements. MgO was found to decrease the rate of the electron transfer reaction from the metal oxide to the electrolyte, decrease the rate of the electron transfer reaction to the oxidized dye from TiO 2 , increase the electron transfer reaction rate from a reduced redox shuttle to an oxidized dye, and decrease the electron injection rate from the photoexcited dye to TiO2. Interestingly, PFTS was found to desirably improve on these rates relative to MgO or untreated TiO 2 . Here, a model based on electrostatic interactions is presented to explain the exceptional behavior of PFTS with DFT computational analysis of PFTS supporting this model.

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

Exploring the influence of transition metals on f-element bonding

The United States aims to triple its nuclear energy production by 2050, which will result in increased uranium usage and spent nuclear fuel generation. This highlights the need for a comprehensive understanding of actinide coordination chemistry, which is crucial for the extraction, processing, purification, and fabrication of uranium-based fuels. Moreover, it is vital for the reprocessing or safe disposal of nuclear waste and effective remediation efforts. Achieving this successfully requires an in-depth understanding of f-electron behavior, as “the role of 5f electrons in bond formation remains a fundamental topic in actinide chemistry”. Introducing a second metal into the system can increase structural dimensionality and diversify structural architecture. Heterometallic systems can also alter material properties, such as magnetic and spectroscopic characteristics, luminescence, and actinide mobility. Additionally, secondary transition metals, even when present only in the second coordination sphere and not directly coordinated, can influence the electron density at the actinyl metal center. In this study, uranium heterometallic single crystals were synthesized by incorporating transition metals such as iron(III), iron(II), nickel(II), manganese(II), copper(I), and cobalt(II). The crystals were formed using 2,6-pyridine dicarboxylic acid and other structurally similar ligands with varying functional groups. The synthesized crystals were analyzed using an extensive array of analytical and computational characterization techniques, including single crystal X-ray diffraction, Raman and infrared spectroscopy, and density functional theory calculations.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗
DOE OSTI2023

High-throughput ab initio reaction mechanism exploration in the cloud with automated multi-reference validation

Quantum chemical calculations on atomistic systems have evolved into a standard approach to studying molecular matter. These calculations often involve a significant amount of manual input and expertise, although most of this effort could be automated, which would alleviate the need for expertise in software and hardware accessibility. Here, we present the AutoRXN workflow, an automated workflow for exploratory high-throughput electronic structure calculations of molecular systems, in which (i) density functional theory methods are exploited to deliver minimum and transition-state structures and corresponding energies and properties, (ii) coupled cluster calculations are then launched for optimized structures to provide more accurate energy and property estimates, and (iii) multi-reference diagnostics are evaluated to back check the coupled cluster results and subject them to automated multi-configurational calculations for potential multi-configurational cases. All calculations are carried out in a cloud environment and support massive computational campaigns. Key features of all components of the AutoRXN workflow are autonomy, stability, and minimum operator interference. We highlight the AutoRXN workflow with the example of an autonomous reaction mechanism exploration of the mode of action of a homogeneous catalyst for the asymmetric reduction of ketones.

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

MACCS (MELCOR Accident Consequence Code System) User Guide Version 4.0, Revision 1

The MELCOR Accident Consequence Code System (MACCS) is used by Nuclear Regulatory Commission (NRC) and various national and international organizations for probabilistic consequence analysis of nuclear power accidents. This User Guide is intended to assist analysts in understanding the MACCS/WinMACCS model and to provide information regarding the code. This user guide version describes MACCS Version 4.0. Features that have been added to MACCS in subsequent versions are described in separate documentation. This User Guide provides a brief description of the model history, explains how to set up and execute a problem, and informs the user of the definition of various input parameters and any constraints placed on those parameters. This report is part of a series of reports documenting MACCS. Other reports include the MACCS Theory Manual, MACCS Verification Report, Technical Bases for Consequence Analyses Using MACCS, as well as documentation for preprocessor codes including SecPop, MelMACCS, and COMIDA2.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

Data from: Coupled machine learning-ecosystem ensemble models substantially improve predictions of nitrous oxide (N 2 O) fluxes from US croplands

Nitrous oxide (N₂O) is a potent and persistent greenhouse gas, with rising atmospheric concentrations driven in part by inefficient use of synthetic nitrogen (N) fertilizers in agriculture. Predicting soil N₂O emissions is challenging due to high spatial and temporal variability arising from complex soil biogeochemical processes. Process-based ecosystem models and standalone machine learning (ML) approaches without extensive site-specific calibration often miss high emission episodes. Here, we show how an Ensemble Modeling System (EMS) based on outputs from an ensemble of ecosystem models coupled to an ensemble of ML models can improve predictions and understanding of N2O fluxes from US cropland. Trained and validated on approximately 12,000 N2O chamber measurements at 17 U.S. Midwest sites (six crops, 35 management practices), the EMS accurately predicted daily fluxes of N2O at both training (R² = 0.84, RMSE = 16.4 g N ha⁻¹ d⁻¹) and held-out testing sites (R² = 0.84, RMSE = 6.2 g N ha⁻¹ d⁻¹). Analyses identified six dominant N₂O drivers: soil organic carbon (SOC), NH₄⁺, NO₃⁻, water-filled pore space (WFPS), soil temperature, and biomass production. Wet, warm soils produced large N₂O peaks only with sufficient SOC and mineral N; in low-SOC soils, fluxes remained low. Incorporating these drivers into process-based models might significantly improve their predictive capacity. The EMS demonstrates a strong potential to predict N₂O fluxes at unseen sites, enabling more reliable regional inventories, improved gap-filling where measurements are sparse, and enhanced understanding of mechanisms to advance targeted mitigation strategies in food, feed, and bioenergy crops.

agricultural sciences↗
DOE OSTI2023

Feasibility of a Novel Density Functional Method Outside the Kohn-Sham Framework for Modeling Global Potential Energy Surfaces of Molecular Chemical Reactions (Final Technical Report)

The project aimed towards the construction of computational methods capable of modeling the global potential energy surface (PES) of small molecules, molecular ions, and radicals — including the parts of the PES which correspond to chemical reactions, and the reaction paths connecting reactants, intermediates, and reaction products. This goal may seem humble at first glance, but for chemical systems with more than about six atoms in total, at the time the project was proposed, there were no established theoretical methods capable of simulating such systems reliably, not even for small molecules in the gas phase which are electronically benign. This restriction severely hampers our ability to model and control chemical processes under harsh conditions. With the goal of constructing a method capable of modeling such chemical systems, we proposed to pursue a novel approach towards a Multi-Configuration (MC) DFT outside the traditional frameworks of Kohn-Sham theory and other methods of coupling wave function theory with DFT. Rather than being a complete active space (CAS) method, the proposed DOCI-DFT would employ a special restricted form of the active space wave function, called Doubly-Occupied CI (DOCI)—this wave function form is sufficient to describe not only heterolytic, but also homolytic bond dissociation processes at the zeroth order (i.e., as active space wave function). There are no other standard mean field methods which can do so. In the original proposal, we also outlined strong formal and practical arguments speaking for this method. The project proposal was accepted by the Department of Energy and provided two years of seed funding for one graduate student, as well as two years of two weeks PI summer-salary for the PI. However, ultimately the project could not be effectively pursued due to severe interferences outside the context scientific problems and was cancelled. Apart from two weeks of PI salary (and associated fringe & overhead costs issues by the performing institution), no costs were charged against the award, with the entire rest of funding returned to the Department of Energy. This report reiterates the primary background information regarding the project, its original goals, and outlines the preliminary work performed during the two weeks of DOE funding charged.

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

Combining Computational Modeling with Reaction Kinetics Experiments for Elucidating the In Situ Nature of the Active Site in Catalysis

Microkinetic modeling based on density functional theory (DFT) derived energetics is important for addressing fundamental questions in catalysis. The quantitative fidelity of microkinetic models (MKMs), however, is often insufficient to conclusively infer the mechanistic details of a specific catalytic system. Here, this can be attributed to a number of factors such as an incorrect model of the active site for which DFT calculations are performed, deficiencies in the hypothesized reaction mechanism, inadequate consideration of the surface environment under reaction conditions, and intrinsic errors in the DFT exchange-correlation functional. Despite these limitations, we aim at developing a rigorous understanding of the reaction mechanism and of the nature of the active site for heterogeneous catalytic chemistries under reaction conditions. By achieving parity between experimental and modeling outcomes through robust parameter estimation and by ensuring coverage-consistency between DFT calculations and MKM predictions, it is possible to systematically refine the mechanistic model and, thereby, our understanding of the catalytic active site in situ.

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

Utilizing Essential Symmetry Breaking in Auxiliary-Field Quantum Monte Carlo: Application to the Spin Gaps of the C 36 Fullerene and an Iron Porphyrin Model Complex

We present three distinct examples where phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) can be reliably performed with a single-determinant trial wave function with essential symmetry breaking. Essential symmetry breaking was first introduced by Lee and Head-Gordon [ Phys. Chem. Chem. Phys. 2019, 21, 4763-4778, 10.1039/C8CP07613H]. We utilized essential complex and time-reversal symmetry breaking with ph-AFQMC to compute the triplet–singlet energy gap in the TS12 set. We found statistically better performance of ph-AFQMC with complex-restricted orbitals than with spin-unrestricted orbitals. We then showed the utilization of essential spin symmetry breaking when computing the singlet–triplet gap of a known biradicaloid, C 36 . ph-AFQMC with spin-unrestricted Hartree–Fock (ph-AFQMC+UHF) fails catastrophically even with spin-projection and predicts no biradicaloid character. With approximate Brueckner orbitals obtained from regularized orbital-optimized second-order Møller–Plesset perturbation theory (κ-OOMP2), ph-AFQMC quantitatively captures strong biradicaloid character of C 36 . Lastly, we applied ph-AFQMC to the computation of the quintet–triplet gap in a model iron porphyrin complex where brute-force methods with a small active space fail to capture the triplet ground state. We show unambiguously that neither triplet nor quintet is strongly correlated using UHF, κ-OOMP2, and coupled-cluster with singles and doubles (CCSD) performed on UHF and κ-OOMP2 orbitals. There is no essential symmetry breaking in this problem. By virtue of this, we were able to perform UHF+ph-AFQMC reliably with a cc-pVTZ basis set and predicted a triplet ground state for this model geometry. The largest ph-AFQMC in this work correlated 186 electrons in 956 orbitals. Our work highlights the utility, scalability, and accuracy of ph-AFQMC with a single-determinant trial wave function with essential symmetry breaking for systems mainly dominated by dynamical correlation with little static correlation.

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

Review—Micro-Fuel Cell Principal Biosensors for Monitoring Transdermal Volatile Organic Compounds in Humans

Knowledge of transduction mechanisms in biosensing applications paves the way for ultrasensitive and dynamic detection in living systems. Real-world biosensing applications where ultra-sensitivity and dynamic detection are paramount include monitoring the anesthetic agent concentration during surgery; the slightest variation in concentration can potentially result in a life-threatening overdose or, on the other end of the spectrum, the patient’s awareness during the procedure. We review the benefits and functions of the transcutaneous biosensor device compared with other current technology and discuss the sensor’s capability to accurately measure volatile anesthetic gas concentration in blood using fuel cell technology. We review fundamental concepts of fuel-cell technology for wearable bio-sensing applications. The fuel cell sensor can also continuously monitor other volatile organic compounds making it versatile with numerous potential applications.

42 ENGINEERING↗
DOE OSTI2025

Spectroscopic and Computational Studies of Spin-Orbit Coupling of Lanthanide Oxides

Spin-orbit (SO) coupling makes it possible for spin forbidden transitions or reactions within non-relativistic quantum theory feasible. Thus, quantification of such interactions has important implications in photophysics and chemical catalysis. This work aims to quantify SO interactions using lanthanide (Ln) oxides as target molecular systems. Through these systems, we examine the impact on the SO coupling by electron configurations and 4f orbital occupancies of Ln elements and sizes of the metal oxides. By examining these factors, we explore how SO coupling is affected by the number of Ln 4f electrons for a given size of molecules and if the Ln 4f orbitals remain atomic in nature in these small clusters. Ln oxides are produced in laser ablation molecular beams, identified with time-of-fight mass spectrometry, and characterized with laser spectroscopy and relativistic quantum chemical computations. Spectroscopic measurements include mass-analyzed threshold ionization (MATI), zero electron kinetic energy, and slow electron velocity-map imaging spectroscopies. Relativistic computations treat scalar relativistic corrections, electron correlations, and SO interactions. The main results are SO terms and energies of the neutral molecules and singly charged cations, ionization energies of the neutral species, metal-metal and metal-oxygen vibrational frequencies of the ions and, in some cases, neutrals as well, and charge effects on the bonding and structures.

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

Theoretical studies of chemical reactions related to the formation and growth of polycyclic aromatic hydrocarbons (PAH) and molecular properties of their key intermediates (Final Progress Report)

The formation mechanisms of polycyclic aromatic hydrocarbons, (PAHs) – organic molecules carrying fused benzene rings – are of great interest to scientists and engineers due to their importance in combustion chemistry and astrochemistry. On Earth, PAHs are largely produced in incomplete combustion of fossil fuel and are considered as critical precursors to unwanted soot particles leading to combustion inefficiency and causing air pollution along with detrimental health effects. Simple PAH molecules initially formed in the gas phase, are further involved in a build-up process in combustion flames leading to larger PAH, bowl-shaped nanostructures, fullerenes, and solid-phase species including carbonaceous dust, graphene particles, and soot. In deep space, PAH and their derivatives are potential key intermediates and nucleation sites leading eventually to carbonaceous nanoparticles (“interstellar grains”). Therefore, the understanding of the key processes in the synthesis of PAHs along with their precursors and their degradation mechanisms in combustion systems and in interstellar, circumstellar, and planetary atmospheric environments will provide critical insights into how complex aromatic structures, carbonaceous nanoparticles, and fullerenes are formed and destroyed. Achieving this understanding is an important step in the development of the efficient combustion processes and of the ecofriendly devices with reduced environmental pollution as well as technological strategies for the production of hydrogen and solid carbon through thermal or plasma-assisted pyrolysis of natural gas and biomass. Also, the understanding of the key processes of PAH and soot growth will help in our comprehension of chemical evolution in the universe. Detailed information on the mechanisms and reliable rate constants of the key elementary chemical reactions involved in PAH formation and destruction processes and in inception of soot particles is often missing, with the main deficiencies being the absence of temperature- and pressure-dependent rate constants for the broad range of conditions occurring in various terrestrial and interstellar processes and the lack of data on the reaction products and their branching ratios. Complementary to experimental studies, these gaps in knowledge can be filled by using quantum chemical calculations of reaction potential energy surfaces providing us with accurate energies of reaction products, intermediates, and transition states, revealing the reaction mechanism, and giving the molecular properties required to compute rate constants for relevant reaction steps and product branching ratios using the RRKM-Master Equation (ME) method. Molecular dynamics (MD) simulations can be used in cases when a reaction rate cannot be properly described by statistical theories. During the terminal renewal project period we employed these ab initio/RRKM-ME and MD approaches to complete our studies on several key reactions relevant to the formation/growth of PAH and inception of soot particles including (1) the reaction mechanism and kinetics of the resonance stabilized fulvenallenyl radical with propargyl and C 3 H 4 isomers; (2) the reaction mechanism and kinetics for the C + indene and C 2 + styrene reactions producing naphthyl or azulenyl radicals in low-temperature environments; (3) the MD study of non-equilibrium dimerization of acepyrene and coronene and its radical. The information derived from our theoretical calculations contributed to a better fundamental understanding of the reaction mechanisms and provide missing critical kinetic data to improve combustion models of hydrocarbon fuels and astrochemical models of the growth of carbonaceous molecules and particles in cold molecular clouds, circumstellar envelopes, and planetary atmospheres.

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

Using anion photoelectron spectroscopy of cluster models to gain insights into mechanisms of catalyst-mediated H 2 production from water

We report that metal oxide cluster models of catalyst materials offer a powerful platform for probing the molecular-scale features and interactions that govern catalysis. This perspective gives an overview of studies implementing the combination of anion photoelectron (PE) spectroscopy and density functional theory calculations toward exploring cluster models of metal oxides and metal-oxide supported Pt that catalytically drive the hydrogen evolution reaction (HER) or the water–gas shift reaction. The utility in the combination of these experimental and computational techniques lies in our ability to unambiguously determine electronic and molecular structures, which can then connect to results of reactivity studies. In particular, we focus on the activity of oxygen vacancies modeled by suboxide clusters, the critical mechanistic step of forming proximal metal hydride and hydroxide groups as a prerequisite for H2 production, and the structural features that lead to trapped dihydroxide groups. The pronounced asymmetric oxidation found in heterometallic group 6 oxides and near-neighbor group 5/group 6 results in higher activity toward water, while group 7/group 6 oxides form very specific stoichiometries that suggest facile regeneration. Studies on the trans-periodic combination of cerium oxide and platinum as a model for ceria supported Pt atoms and nanoparticles reveal striking negative charge accumulation by Pt, which, combined with the ionic conductivity of ceria, suggests a mechanism for the exceptionally high activity of this system towards the water–gas shift reaction.

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

Local self-interaction correction method with a simple scaling factor

A recently proposed local self-interaction correction (LSIC) method [Zope et al., J. Chem. Phys., 2019, 151, 214108] when applied to the simplest local density approximation provides a significant improvement over standard Perdew–Zunger SIC (PZSIC) for both equilibrium properties such as total or atomization energies as well as properties involving stretched bond such as barrier heights. The method uses an iso-orbital indicator to identify the single-electron regions. To demonstrate the LSIC method, Zope et al. used the ratio z σ of von Weizsäcker τ W σ and total kinetic energy densities τ σ , (z σ = τ W σ /τ σ ) as a scaling factor to scale the self-interaction correction. The present work further explores the LSIC method using a ratio of orbital and spin densities as a simpler scaling factor in place of the ratio of kinetic energy densities. We compute a wide array of both, equilibrium and non-equilibrium properties using LSIC and orbital scaling methods using this simple scaling factor and compare them with previously reported results. Our study shows that LSIC with the simple scaling factor performs better than PZSIC, with results comparable to those obtained by LSIC( z σ ) for most properties, but has slightly larger errors than LSIC( z σ ). Furthermore, we study the binding energies of small water clusters using both scaling factors. Our results show that LSIC with zσ has limitations in predicting the cluster binding energies of weakly bonded systems due to the inability of z σ to distinguish weakly bonded regions from slowly varying density regions. LSIC when used with the density ratio as a scaling factor, on the other hand, provides a good description of water cluster binding energies, thus highlighting the appropriate choice of the iso-orbital indicator.

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

Random Phase Approximation Correlation Energy Using Real-Space Density Functional Perturbation Theory

We present a real-space method for computing the random phase approximation (RPA) correlation energy within Kohn–Sham density functional theory, leveraging the low-rank nature of the frequency-dependent density response operator. In particular, we employ a cubic-scaling formalism based on density functional perturbation theory that circumvents the calculation of the response function matrix, instead relying on the ability to compute its product with a vector through the solution of the associated Sternheimer linear systems. We develop a large-scale parallel implementation of this formalism using the subspace iteration method in conjunction with the spectral quadrature method while employing the Kronecker product-based method for the application of the Coulomb operator and the conjugate orthogonal conjugate gradient method for the solution of the linear systems. We demonstrate convergence with respect to key parameters and verify the method’s accuracy by comparing with plane-wave results. We show that the framework achieves good strong scaling to many thousands of processors, reducing the time to solution for a lithium hydride system with 128 electrons to around 150 s on 4608 processors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗