Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “KS”

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

Accurate Kohn-Sham auxiliary system from the ground-state density of solids

The Kohn-Sham (KS) system is an auxiliary system whose effective potential is unknown in most cases. It is in principle determined by the ground-state density and it has been found numerically for some low-dimensional systems by inverting the KS equations starting from a given accurate density. For solids, only approximate results are available. In this work, we determine accurate exchange-correlation (xc) potentials for Si and NaCl using the ground-state densities obtained from auxiliary field quantum Monte Carlo calculations. We show that these xc potentials can be rationalized as an ensemble of a few local functions of the density, whose form depends on the specific environment and can be well characterized by the gradient of the density and the local kinetic energy density. Further, the KS band structure can be obtained with high accuracy. The true KS band gap turns out to be larger than the prediction of the local density approximation, but significantly smaller than the measurable photoemission gap, which confirms previous estimates. Finally, our findings show that the conjecture that very different xc potentials can lead to very similar densities and other KS observables is true also in solids, which questions the meaning of details of the potentials and, at the same time, confirms the stability of the KS system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Effects of pH on potassium sulfite based alkali pretreatment of corn stover under mild conditions

Sulfite salt-based alkali pretreatment is increasingly recognized as an effective method for biomass pretreatment under mild temperature conditions. However, the effects of pH on the pretreatment process are not yet fully understood. In this study, the pretreatment of corn stover using potassium sulfite (KS) combined with aqueous ammonia (AA) at 80 °C was investigated within a pH range of 9–14. Results demonstrated that high pH environment established by combining KS with AA is critical for enhancing sulfonation modifications, effectively removing lignin, facilitating deacetylation, and significantly improving the total sugar yield. For example, when comparing KS alone at pH 9.7, AA alone at pH 13, and KS combined with AA at pH 13, the delignification rates were 13.82 %, 67.71 %, and 71.27 %, respectively; the deacetylation rates were 25.84 %, 81.19 %, and 82.54 %, respectively; and the total sugar yields were 24.85 %, 78.96 %, and 92.12 %, respectively. The X-ray photoelectron spectroscopy (XPS) results confirmed that sulfonation reactions occurred under low-temperature and high pH conditions, leading to enhanced delignification efficiency. Results also showed that the cellulase adsorption, water retention, and the enzymatic digestibility of the pretreated solids increased along with rising pH. Moreover, nuclear magnetic resonance (NMR) analysis revealed that KS-based alkali pretreatment retained more β-O-4 native bonds in the pretreated lignin-rich hydrolysate than alkali pretreatment alone did, providing promising opportunities for the subsequent upgrading of lignin into high-value products.

09 BIOMASS FUELS↗
DOE OSTI2025

Engineering Polyketide Stereocenters with Ketoreductase Domain Exchanges

Polyketide synthases (PKSs) are versatile biosynthetic megasynthases capable of producing a diverse range of natural products with many applications, including in pharmaceuticals. The stereochemical precision of PKSs makes them a powerful tool for engineering tailored, unnatural polyketides; however, modifying the stereocenters of a PKS product while maintaining production levels remains a significant challenge. In this study, we systematically tested and evaluated strategies for ketoreductase (KR) domain exchanges, the domain responsible for setting stereocenters of polyketide products. After first optimizing the method for KR exchanges, we then performed 44 KR domain exchanges on three different PKSs to obtain high production of all four stereoisomers in vivo. By testing both one- and two-module PKS systems, we investigated how downstream modules process intermediates with altered stereochemistry and found that the configuration of the α-substituents was critical for gatekeeping by the ketosynthase (KS). To overcome this constraint, we investigated two different strategies for altering the KS domain, including introducing targeted mutations in the downstream KS, and exploring boundaries in exchanging the entire functional unit from the donor PKS. Both strategies successfully modified the KS stereocontrol with distinct trade-offs; the functional unit exchange resulted in higher titer improvements, though it was more likely to break the entire PKS. This study demonstrates a comprehensive approach to successfully engineering all four stereochemical configurations in multiple PKS systems, advancing our understanding of and ability to rationally modify polyketide stereochemistry through multiple engineering strategies.

Keiser, Leah S. [Joint BioEnergy Institute (JBEI),↗
DOE OSTI2022

Comparing correlation components and approximations in Hartree–Fock and Kohn–Sham theories via an analytical test case study

The asymmetric Hubbard dimer is a model that allows for explicit expressions of the Hartree–Fock (HF) and Kohn–Sham (KS) states as analytical functions of the external potential, Δv, and of the interaction strength, U. We use this unique circumstance to establish a rigorous comparison between the individual contributions to the correlation energies stemming from the two theories in the {U, Δv} parameter space. Within this analysis of the Hubbard dimer, we observe a change in the sign of the HF kinetic correlation energy, compare the indirect repulsion energies, and derive an expression for the “traditional” correlation energy, i.e., the one that corrects the HF estimate, in a pure site-occupation function theory spirit [Eq. (45)]. Next, we test the performances of the Liu–Burke and the Seidl–Perdew–Levy functionals, which model the correlation energy based on its weak- and strong-interaction limit expansions and can be used for both the traditional and the KS correlation energies. Further, our results show that, in the Hubbard dimer setting, they typically work better for the HF reference, despite having been originally devised for KS. These conclusions are somewhat in line with prior assessments of these functionals on various chemical datasets. However, the Hubbard dimer model allows us to show the extent of the error that may occur in using the strong-interaction ingredient for the KS reference in place of the one for the HF reference, as has been carried out in most of the prior assessments.

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

Dynamical structure factors of warm dense matter from time-dependent orbital-free and mixed-stochastic-deterministic density functional theory

Abstract We present the first calculations of the inelastic part of the dynamical structure factor (DSF) for warm dense matter (WDM) using time-dependent orbital-free density functional theory (TD-OF-DFT) and mixed-stochastic-deterministic (mixed) Kohn Sham TD-DFT (KS TD-DFT). WDM is an intermediate phase of matter found in planetary cores and laser-driven experiments, where the accurate calculation of the DSF is critical for interpreting x-ray Thomson scattering measurements. Traditional TD-DFT methods, while highly accurate, are computationally expensive, motivating the exploration of TD-OF-DFT and mixed TD-KS-DFT as more efficient alternatives. We applied these methods to experimentally measured WDM systems, including solid-density aluminum and beryllium, compressed beryllium, and carbon–hydrogen mixtures. Our results show that TD-OF-DFT requires a dynamical kinetic energy potential in order to qualitatively capture the plasmon response. Additionally, it struggles with capturing bound electron contributions. In contrast, mixed TD-KS-DFT offers greater accuracy in distinguishing bound and free electron effects, aligning well with experimental data, though at a higher computational cost. This study highlights the trade-offs between computational efficiency and accuracy, demonstrating that TD-OF-DFT remains a valuable tool for rapid scans of parameter space, while mixed TD-KS-DFT should be preferred for high-fidelity simulations. Our findings provide insight into the future development of DFT methods for WDM and suggest potential improvements for TD-OF-DFT.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Advances and Challenges of SCAN and r 2 SCAN Density Functionals in Transition‐Metal Compounds

ABSTRACT Transition‐metal compounds (TMCs) with open‐shell d ‐electrons are characterized by a complex interplay of lattice, charge, orbital, and spin degrees of freedom, giving rise to various fascinating applications. Often exhibiting exotic properties, these compounds are commonly classified as correlated systems due to strong inter‐electronic interactions called Hubbard U . This inherent complexity presents significant challenges to Kohn‐Sham density functional theory (KS‐DFT), the most widely used electronic structure method in condensed matter physics and materials science. While KS‐DFT is, in principle, exact for the ground‐state total energy, its exchange‐correlation energy must be approximated in practice. The mean‐field nature of KS implementations, combined with the limitations of current exchange‐correlation density functional approximations, has led to the perception that DFT is inadequate for correlated systems, particularly TMCs. Consequently, a common workaround involves augmenting DFT with an on‐site Hubbard‐like U correction. In recent years, the strongly constrained and appropriately normed (SCAN) density functional, along with its refined variant r 2 SCAN, has achieved remarkable progress in accurately describing the structural, energetic, electronic, magnetic, and vibrational properties of TMCs, challenging the traditional perception of DFT's limitations. This review explores the design principles of SCAN and r 2 SCAN, highlights their key advancements in studying TMCs, explains the mechanisms driving these improvements, and addresses the remaining challenges in this evolving field.

Zhang, Yubo [Minjiang Collaborative Center for The↗
DOE OSTI2020

Electrode-molecule energy level offsets in a gold-benzene diamine-gold single molecule tunnel junction

One means for describing electron transport across single molecule tunnel junctions (MTJs) is to use density functional theory (DFT) in conjunction with a nonequilibrium Green’s function formalism. This description relies on interpreting solutions to the Kohn–Sham (KS) equations used to solve the DFT problem as quasiparticle (QP) states. Many practical DFT implementations suffer from electron self-interaction errors and an inability to treat charge image potentials for molecules near metal surfaces. For MTJs, the overall effect of these errors is typically manifested as an overestimation of electronic currents. Correcting KS energies for self-interaction and image potential errors results in MTJ current–voltage characteristics in close agreement with measured currents. An alternative transport approach foregoes a QP picture and solves for a many-electron wavefunction on the MTJ subject to open system boundary conditions. It is demonstrated that this many-electron method provides similar results to the corrected QP picture for electronic current. The analysis of these two distinct approaches is related through corrections to a junction’s electronic structure beyond the KS energies for the case of a benzene diamine molecule bonded between two gold electrodes.

Szepieniec, Mark S.↗
DOE OSTI2021

Efficacy of the radial pair potential approximation for molecular dynamics simulations of dense plasmas

Macroscopic simulations of dense plasmas rely on detailed microscopic information that can be computationally expensive and is difficult to verify experimentally. In this work, we delineate the accuracy boundary between microscale simulation methods by comparing Kohn–Sham density functional theory molecular dynamics (KS-MD) and radial pair potential molecular dynamics (RPP-MD) for a range of elements, temperature, and density. By extracting the optimal RPP from KS-MD data using force matching, we constrain its functional form and dismiss classes of potentials that assume a constant power law for small interparticle distances. Our results show excellent agreement between RPP-MD and KS-MD for multiple metrics of accuracy at temperatures of only a few electron volts. The use of RPPs offers orders of magnitude decrease in computational cost and indicates that three-body potentials are not required beyond temperatures of a few eV. Due to its efficiency, the validated RPP-MD provides an avenue for reducing errors due to finite-size effects that can be on the order of ∼20%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

String theory and grand unification suggest a submicroelectronvolt QCD axion

Axions, grand unification, and string theory are each compelling extensions of the Standard Model. We show that combining these frameworks imposes strong constraints on the QCD axion mass. Using perturbative unitarity arguments and explicit string compactifications—such as those from the Kreuzer-Skarke (KS) type IIB ensemble—we find that the axion mass is favored to lie within the range 10 −11 eV ≲ 𝑚 𝑎 ≲ 10 −8 eV. This range is directly relevant for near-future axion dark matter searches, including ABRACADABRA/DMRadio and CASPEr. We argue that grand unification and the absence of proton decay suggest a compactification volume that keeps the string scale above the unification scale ( ∼ 10 16 GeV), which in turn limits how heavy the axion can be. The same requirements limit the KS axiverse to have at most ∼ 47 axions. As an additional application of our methodology, we search for axions in the KS axiverse that could explain the recent Dark Energy Spectroscopic Instrument hints of evolving dark energy but find none with high enough decay constant (𝑓 𝑎 ≳ 2.5 ×10 17 GeV); we comment on why such high decay constants and low axion masses are difficult to obtain in string compactifications more broadly.

Axions↗
DOE OSTI2020

Fast and Universal Kohn-Sham Density Functional Theory Algorithm for Warm Dense Matter to Hot Dense Plasma

Understanding many processes, e.g., fusion experiments, planetary interiors, and dwarf stars, depends strongly on microscopic physics modeling of warm dense matter and hot dense plasma. This complex state of matter consists of a transient mixture of degenerate and nearly free electrons, molecules, and ions. This regime challenges both experiment and analytical modeling, necessitating predictive ab initio atomistic computation, typically based on quantum mechanical Kohn-Sham density functional theory (KS-DFT). However, cubic computational scaling with temperature and system size prohibits the use of DFT through much of the warm dense matter regime. A recently developed stochastic approach to KS-DFT can be used at high temperatures, with the exact same accuracy as the deterministic approach, but the stochastic error can converge slowly and it remains expensive for intermediate temperatures (< 50 eV). Here we have developed a universal mixed stochastic-deterministic algorithm for DFT at any temperature. This approach leverages the physics of KS-DFT to seamlessly integrate the best aspects of these different approaches. We demonstrate that this method significantly accelerated self-consistent field calculations for temperatures from 3 to 50 eV, while producing stable molecular dynamics and accurate diffusion coefficients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Automated Identification of Characteristic Droplet Size Distributions in Stratocumulus Clouds Utilizing a Data Clustering Algorithm

Abstract Droplet-level interactions in clouds are often parameterized by a modified gamma fitted to a “global” droplet size distribution. Do “local” droplet size distributions of relevance to microphysical processes look like these average distributions? This paper describes an algorithm to search and classify characteristic size distributions within a cloud. The approach combines hypothesis testing, specifically, the Kolmogorov–Smirnov (KS) test, and a widely used class of machine learning algorithms for identifying clusters of samples with similar properties: density-based spatial clustering of applications with noise (DBSCAN) is used as the specific example for illustration. The two-sample KS test does not presume any specific distribution, is parameter free, and avoids biases from binning. Importantly, the number of clusters is not an input parameter of the DBSCAN-type algorithms but is independently determined in an unsupervised fashion. As implemented, it works on an abstract space from the KS test results, and hence spatial correlation is not required for a cluster. The method is explored using data obtained from the Holographic Detector for Clouds (HOLODEC) deployed during the Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA) field campaign. The algorithm identifies evidence of the existence of clusters of nearly identical local size distributions. It is found that cloud segments have as few as one and as many as seven characteristic size distributions. To validate the algorithm’s robustness, it is tested on a synthetic dataset and successfully identifies the predefined distributions at plausible noise levels. The algorithm is general and is expected to be useful in other applications, such as remote sensing of cloud and rain properties. Significance Statement A typical cloud can have billions of drops spread over tens or hundreds of kilometers in space. Keeping track of the sizes, positions, and interactions of all of these droplets is impractical, and, as such, information about the relative abundance of large and small drops is typically quantified with a “size distribution.” Droplets in a cloud interact locally, however, so this work is motivated by the question of whether the cloud droplet size distribution is different in different parts of a cloud. A new method, based on hypothesis testing and machine learning, determines how many different size distributions are contained in a given cloud. This is important because the size distribution describes processes such as cloud droplet growth and light transmission through clouds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Production of ent-kaurene from lignocellulosic hydrolysate in Rhodosporidium toruloides

Abstract Background Rhodosporidium toruloides has emerged as a promising host for the production of bioproducts from lignocellulose, in part due to its ability to grow on lignocellulosic feedstocks, tolerate growth inhibitors, and co-utilize sugars and lignin-derived monomers. Ent -kaurene derivatives have a diverse range of potential applications from therapeutics to novel resin-based materials. Results The Design, Build, Test, and Learn (DBTL) approach was employed to engineer production of the non-native diterpene ent -kaurene in R. toruloides . Following expression of kaurene synthase (KS) in R. toruloides in the first DBTL cycle, a key limitation appeared to be the availability of the diterpene precursor, geranylgeranyl diphosphate (GGPP). Further DBTL cycles were carried out to select an optimal GGPP synthase and to balance its expression with KS, requiring two of the strongest promoters in R. toruloides , ANT (adenine nucleotide translocase) and TEF1 (translational elongation factor 1) to drive expression of the KS from Gibberella fujikuroi and a mutant version of an FPP synthase from Gallus gallus that produces GGPP. Scale-up of cultivation in a 2 L bioreactor using a corn stover hydrolysate resulted in an ent -kaurene titer of 1.4 g/L. Conclusion This study builds upon previous work demonstrating the potential of R. toruloides as a robust and versatile host for the production of both mono- and sesquiterpenes, and is the first demonstration of the production of a non-native diterpene in this organism.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

GLBRC Soil Yearlong Incubation 13C-SIP-Lipidomics

Data package for Lipids represent a dynamic, yet stable pool of microbially-derived soil carbon This data is published under a CC0 license. The authors encourage data reuse and request attribution by referencing the below citations for the data packages and associated manuscript. Please cite as: Rempfert KR, Bell SL, Kasanke CP, Kyle JE, Hofmockel KS. 2025. GLBRC Soil Yearlong Incubation 13C-SIP-Lipidomics. [Data Set] PNNL DataHub. doi: Rempfert KR, Bell SL, Kasanke CP, Kyle JE, Hofmockel KS. 2025. MSV000097435: GLBRC soil yearlong incubation 13C-SIP-Lipidomics [Data Set] MassIVE. doi:10.25345/C57659T3K Rempfert KR, Bell SL, Kasanke CP, Kyle JE, Hofmockel KS. 2025. Lipids represent a dynamic, yet stable pool of microbially-derived soil carbon. In Prep This data package consists of compound-specific 13C SIP-lipidomics data from a yearlong tracer incubation experiment designed to investigate microbial lipid persistence in switchgrass bioenergy crop soils. In order to explore how lipid structure may modulate the persistence of C in soil lipids, we leveraged soils from two sites (Michigan - sandy texture, Wisconsin - silty texture) operated by the U.S. Department of Energy-funded Great Lakes Bioenergy Research Center (GLBRC). These sites had comparable climates, identical management practices, but contrasting soil textures, allowing us to assess the variability of lipid accrual or degradation in soils as well as provide insight regarding the degree to which edaphic properties may regulate the retention of soil lipids. Untargeted lipidomics analyses were performed to identify 13C-labeled lipids in the soil microbiome after long-term incubation. Soils were supplemented with 100 micrograms glucose per gram dry soil (99 atom % 13C or natural abundance for paired control) and incubated; samples were collected two months and one year after glucose addition. Lipid extracts (MPLEx) were analyzed by LC-MS/MS and identified using LIQUID. Calculation of isotopic enrichment of lipids was performed by targeted approach using TarMet to quantify lipid isotopologues and IsoCorrectoR to correct for natural abundance isotopes. Contents: Data package contents reported here are the first version and contain downstream analysis files for the raw LC-MS mass spectrometry files (.mzXML) deposited at the MassIVE database repository under accession MSV000097435 (80 experimental runs; 5.85 GB) | MassIVE DOI: 10.25345/C57659T3K. Support files include the additional data download 'Read Me' file containing data descriptor information. Reported data download contents are structured for compliance with project data sharing guidelines, community standards initiatives, and sponsor stakeholder policies supporting FAIR data principles. Data processing software, analysis tools, and data workflows are listed below corresponding to the host repository long-term location. Available Data Downloads (0.3 GB): "GLBRC soil yearlong incubation 13C-SIP-Lipidomics_readme.txt" - 'Read Me' data package content file (txt) "GLBRC_DataPackage_analysis files" - Data processing files (Rmd) and saved intermediate data processing outputs (rds, csv, xlsx) "GLBRC_13C_lipidomics_dataset.xlsx" - processed data in tabular format (xlsx) Linked Software: LIQUID LC-MS Analysis Software | 10.5281/zenodo.6459462 Lipid Mini-On Software Tools | 10.5281/zenodo.1492803 pmartR Omics Statistical Software | 10.5281/zenodo.6108667 xcms (v4.3.3) TarMet (v1.1.1) IsoCorrectoR (1.24.0) Funding Acknowledgments: This research was supported by an Early Career Research Program award funded by the U.S. Department of Energy, Office of Science, Office of Biological and Environmental Research (OBER) Genomic Science program under FWP 68292, FWP 07880 and EMSL Exploratory Research Project 51095. A portion of this work was performed in the William R. Wiley Environmental Molecular Sciences Laboratory, a national scientific user facility sponsored by OBER and located at Pacific Northwest National Laboratory (PNNL). PNNL is a multi-program national laboratory operated by Battelle for the DOE under Contract DE-AC05-76RLO1830.

Rempfert, Kaitlin R [Pacific Northwest National La↗
DOE OSTI2026

Deep X-Ray Observation of NGC 3221: Everything Everywhere All at Once

We present a comprehensive analysis of 475 ks (438 ks unpublished and 37 ks archival) XMM-Newton/EPIC-pn observations of a nearby, highly inclined, star-forming, luminous infrared galaxy NGC 3221 through spatial, temporal, and spectral information. We confirm the presence of a low-luminosity (presumably Compton-thick) active galactic nucleus (AGN). The 0.4–12 keV luminosity and the hardness ratio of the six ultraluminous X-ray sources previously identified in Chandra data exhibit diverse variability on day scales. The collective emission from unresolved sources exhibits a different day-scale variability. We have also discovered two new predominantly soft (<1 keV) sources. One of these has an enigmatic spectral shape featuring a soft component, which we interpret as a superbubble in NGC 3221, and a variable hard component from a compact object, unresolved from the superbubble. We do not confidently detect any X-ray emission from SN 1961L. The hot gas in the interstellar medium (ISM, out to ±6 kpc from the disk plane) and that in the extraplanar region (6–12 kpc) both require two thermal phases at ∼0.15 keV and ∼0.55 keV. The ∼0.55 keV component is fainter in the ISM than the ∼0.15 keV component, but the emission from the latter falls off more steeply with disk height than the emission from the former. This makes the extraplanar region hotter and less dense than the ISM. The proximity of NGC 3221 and the occurrence of the underluminous AGN offer a unique observing opportunity to study the hot diffuse medium along with nuclear and diskwide point sources.

Das, Sanskriti [Stanford Univ., CA (United States)↗
DOE OSTI2024

Measurement of the Photoproduction Cross-Section of f1(1285) in the Exclusive Reactions p¿ p'K±KSp± AT 7.5 < E¿ < 11.5 GeV with GlueX at Jefferson Lab

The f1(1285) meson is commonly understood to belong to the axial-vector nonet as the low-mass isoscalar partner to the f1(1420) but has been suggested as a candidate for a KK?+c.c molecule. A nearly mass-degenerate 0?+ state has been observed in ? p scattering that calls into question the established branching ratios of the f1(1285) decays. Recently, the differential cross-section for the photoproduction of f1(1285) was measured by the CLAS experiment and the results disagreed with theoretical predictions. Additionally, the CLAS experiment did not observe a mass-degenerate 0?+ state. We present the results of the first photoproduction cross-section measurement d?/dt of f1(1285) at 7.5 < E? < 11.5 GeV from the GlueX experiment at Thomas Jefferson National Accelerator Facility in the reactions ? p ? p?K?Ks?+ and ? p ? p?K+Ks??. We find that the production mechanism of f1(1285) is consistent with t-channel exchange. We also observe a difference in the differential cross-section of f1(1285) as a function of ?t as measured in the charge conjugated K?Ks?± decay modes that we hypothesize to be due to baryonic N? interference. The observed distributions of M(KK) and M(K?) are consistent with both a q ¯ q state and predictions for a KK?+c.c molecular state. Additionally, we find that the cross-section at low M(KK?) is dominated by the 1++ state.

Viducic, Tyler↗
DOE OSTI2023

$\mathrm{DRAGON}$: A multi-GPU orbital-free density functional theory molecular dynamics simulation package for modeling of warm dense matter

As progress in electronic structure theoretical methods is made, ab initio molecular dynamics (MD) based on orbital-free density functional theory (OF-DFT) is becoming increasingly more successful at substituting the traditional, very accurate but computationally costly Kohn–Sham (KS) approach for simulations of matter at the challenging warm dense matter (WDM) regime. However, despite the significant cost alleviation of eliminating the dependence on the KS orbitals, OF-DFT MD runs require ~10 2 to 10 3 CPU cores running for days, or even weeks, for simulations of systems comprised of 10 2 to 10 3 atoms, depending on thermodynamic conditions. Here we present DRAGON, a multi-GPU OF-DFT MD code for fast and efficient simulations of WDM. With a relatively small allocation of resources (4 to 8 GPU devices) it can provide an order of magnitude speedup for simulations containing $\mathscr{O}$(10 4 ) atoms and target systems composed of $\mathscr{O}$(10 5 ) atoms at conditions within the WDM regime, which is currently outside the capabilities of CPU codes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Achieving performance portability in Gaussian basis set density functional theory on accelerator based architectures in NWChemEx

The numerical integration of the exchange–correlation (XC) potential is one of the primary computational bottlenecks in Gaussian basis set Kohn–Sham density functional theory (KS-DFT). To achieve optimal performance and accuracy, care must be taken in this numerical integration to preserve local sparsity as to allow for near linear weak scaling with system size. This leads to an integration scheme with several performance critical kernels which must be hand optimized for each architecture of interest. As the set of available accelerator hardware goes more diverse, a key challenge for developers of KS-DFT software is to maintain performance portability across a wide range of computational architectures. In this article, we examine a modular software design pattern which decouples the implementation details of performance critical kernels from the expression of high-level algorithmic workflows in a device-agnostic language such as C++; thus allowing for developers to target existing and emerging accelerator hardware within a single code base. We consider the efficacy of such a design pattern in the numerical integration of the XC potential by demonstrating its ability to achieve performance portability across a set of accelerator architectures which are representative of those on current and future U.S. Department of Energy Leadership Computing Facilities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Real-Space Pseudopotential Method for the Calculation of Third-Row Elements X-ray Photoelectron Spectroscopic Signatures

X-ray photoelectron spectroscopy (XPS) is a powerful characterization technique that unveils subtle chemical environment differences via core– electron binding energy (CEBE) analysis. We extend the development of realspace pseudopotential methods to calculating 1s, 2s, and 2p 3/2 CEBEs of third-row elements (S, P, and Si) within the framework of Kohn–Sham density-functional theory (KS-DFT). The new approach systematically prevents variational collapse and simplifies core-excited orbital selection within dense energy level distributions. However, careful error cancellation analysis is required to achieve accuracy comparable to all-electron methods and experiments. Combined with real-space KS-DFT implementation, this development enables large-scale simulations with both Dirichlet boundary conditions and periodic boundary conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗