Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Electronic Structure Software”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Integrated fluorescence light microscopy-guided cryo-focused ion beam-milling for in situ montage cryo-ET

Cryogenic-electron tomography (cryo-ET) permits the in situ visualization of biological macromolecules at the molecular level. Owing to the variable thickness of cells, tissues and organisms, frozen specimens may need to be thinned by cryo-focused ion beam (FIB) milling to produce thin (<500 nm) cryo-lamellae suitable for cryo-ET. Locating regions of interest remains a challenge because untargeted milling can lead to inadvertent ablation and removal of regions of interest. Correlative light and electron microscopy, combined with cryo-FIB milling, can guide the identification of labeled targets in the cellular milieu. Multiple transfers between cryo-imaging instruments, cumbersome correlation algorithms, limited accuracy and low throughput have hindered the routine adoption of cryo-FIB milling within a multimodal correlative workflow for in situ structural biology. Here, in this study, we present a workflow for 3D correlative cryo-fluorescence light microscopy-FIB-ET that streamlines fluorescence light microscopy-guided FIB milling, improving throughput while preserving both structural and contextual information. The complete integration of hardware and software described here minimizes sample contamination from cross-platform exchanges and greatly enhances the efficiency of 3D targeting in cryo-milling. We then describe procedures for implementing montage parallel array cryo-ET (MPACT), which can be easily adapted to any modern life-science transmission electron microscope. MPACT supports high-throughput cryo-ET acquisitions (10 tilt series in 1.5 h) for structure determination and comprehensive contextual understanding of macromolecules within their native surroundings. A complete session from sample preparation to MPACT data processing takes 5−7 d for an individual experienced in both cryo-EM and cryo-FIB milling.

Yang, Jie E. [Univ. of Wisconsin, Madison, WI (Uni↗
DOE OSTI2024

2D Polyhedral Template Matching for Atomic Resolution Microscopy

SAND2024-13879O The 2D Polyhedral Template Matching for Atomic Resolution Microscopy is a suite of functions for analyzing atomic resolution electron microscopy images using the 2D polyhedral template matching (2D-PTM) method. This software analyzes atomic resolution microscopy data from electron microscopic imaging or scanning probe microscopies. The primary application is for identifying different crystal phases, crystal orientations, and defect structures obtained in such atomic resolution images. Written in MATLAB and starting from an atomic resolution image, the code identifies the positions of atomic intensity peaks. It then matches predefined structural templates to the local atomic environments. Outputs include the local structural identification, the template scaling factor and rotation angle, root-mean-squared deviation (RMSD), and centrosymmetry parameter.

Medlin, Douglas↗
DOE OSTI2025

Fragme∩t: An Open‐Source Framework for Multiscale Quantum Chemistry Based on Fragmentation

Fragment-based quantum chemistry offers a means to circumvent the nonlinear computational scaling of conventional electronic structure calculations, by partitioning a large calculation into smaller subsystems then considering the many-body interactions between them. Variants of this approach have been used to parameterize classical force fields and machine learning potentials, applications that benefit from interoperability between quantum chemistry codes. However, there is a dearth of software that provides interoperability yet is purpose-built to handle the combinatorial complexity of fragment-based calculations. To fill this void we introduce “Fragme∩t”, an open-source software application that provides a tool for community validation of fragment-based methods, a platform for developing new approximations, and a framework for analyzing many-body interactions. Fragme∩t includes algorithms for automatic fragment generation and structure modification, and for distance- and energy-based screening of the requisite subsystems. Checkpointing, database management, and parallelization are handled internally and results are archived in a portable database. Interfaces to various quantum chemistry engines are easy to write and exist already for Q-Chem, PySCF, xTB, Orca, CP2K, MRCC, Psi4, NWChem, GAMESS, and MOPAC. Applications reported here demonstrate parallel efficiencies around 96% on more than 1000 processors but also showcase that the code can handle large-scale protein fragmentation using only workstation hardware, all with a codebase that is designed to be usable by non-experts. Fragme∩t conforms to modern software engineering best practices and is built upon well established technologies including Python, SQLite, and Ray. The source code is available under the Apache 2.0 license.

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

Implementation of Relativistic Coupled Cluster Theory for Massively Parallel GPU-Accelerated Computing Architectures

In this paper, we report reimplementation of the core algorithms of relativistic coupled cluster theory aimed at modern heterogeneous high-performance computational infrastructures. The code is designed for parallel execution on many compute nodes with optional GPU coprocessing, accomplished via the new ExaTENSOR back end. The resulting ExaCorr module is primarily intended for calculations of molecules with one or more heavy elements, as relativistic effects on the electronic structure are included from the outset. In the current work, we thereby focus on exact two-component methods and demonstrate the accuracy and performance of the software. The module can be used as a stand-alone program requiring a set of molecular orbital coefficients as the starting point, but it is also interfaced to the DIRAC program that can be used to generate these. We therefore also briefly discuss an improvement of the parallel computing aspects of the relativistic self-consistent field algorithm of the DIRAC program.

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

Cartesian message passing neural networks for directional properties: Fast and transferable atomic multipoles

The message passing neural network (MPNN) framework is a promising tool for modeling atomic properties but is, until recently, incompatible with directional properties, such as Cartesian tensors. We propose a modified Cartesian MPNN (CMPNN) suitable for predicting atom-centered multipoles, an essential component of ab initio force fields. The efficacy of this model is demonstrated on a newly developed dataset consisting of 46 623 chemical structures and corresponding high-quality atomic multipoles, which was deposited into the publicly available Molecular Sciences Software Institute QCArchive server. We show that the CMPNN accurately predicts atom-centered charges, dipoles, and quadrupoles and that errors in the predicted atomic multipoles have a negligible effect on multipole–multipole electrostatic energies. The CMPNN is accurate enough to model conformational dependencies of a molecule’s electronic structure. This opens up the possibility of recomputing atomic multipoles on the fly throughout a simulation in which they might exhibit strong conformational dependence.

Glick, Zachary L. (ORCID:0000000309002849)↗
DOE OSTI2025

CASM Monte Carlo: Calculations of the thermodynamic and kinetic properties of complex multicomponent crystals

Monte Carlo techniques play a central role in statistical mechanics approaches that connect macroscopic thermodynamic and kinetic properties to the electronic structure of a material. This paper describes the implementation of Monte Carlo techniques for the study of multicomponent crystalline materials within the Clusters Approach to Statistical Mechanics (CASM) software suite, and demonstrates their use in model systems to calculate free energies and kinetic coefficients, study phase transitions, and construct phase diagrams from first principles. Many crystal structures are complex, with multiple sublattices occupied by differing sets of chemical species, along with the presence of vacancies or interstitial species. This imposes constraints on concentration variables, the form of thermodynamic potentials, and the values of kinetic transport coefficients. The framework used by CASM to formulate thermodynamic potentials and kinetic transport coefficients accounting for arbitrarily complex crystal structures is presented and demonstrated with examples of increasing complexity. Additionally, an overview of the capabilities of the CASM software specific to Monte Carlo methods is given, and a new CASM software package is introduced, casm-flow, which helps automate the setup, submission, management, and analysis of Monte Carlo simulations.

Cluster expansion↗
DOE OSTI2025

Nuclear Materials Packaging, Transportation, and Systems Analysis Group Software Quality Assurance Plan: ANSYS Mechanical Finite Element Analysis Software Version 2023R1

ANSYS Inc. develops and markets engineering simulation software and services used in the aerospace, automotive, manufacturing, electronics, biomedical, energy, defense, and many other industries. ANSYS is dedicated to engineering simulation and is the world’s leading software provider. ANSYS was founded in 1970 and is headquartered in Canonsburg, Pennsylvania. ANSYS provides an engineering analysis tool combining structural, thermal, computational fluid dynamics, acoustic, and electromagnetic simulation capabilities. ANSYS has two main programs, which use the same solvers: (1) Mechanical APDL (ANSYS Design Parametric Language), a Fortran-based coding platform, and (2) ANSYS Workbench, which uses a graphical user interface to aid in finite element analysis implementation. This plan covers both APDL and Workbench. The ANSYS computer program is a large-scale, multipurpose finite element program that can be used to solve several classes of engineering analyses. The analysis capabilities of ANSYS include the ability to solve static and dynamic structural analyses, steady-state and transient heat transfer problems, mode-frequency and buckling eigenvalue problems, static or time-varying magnetic analyses, and various types of field and coupled-field applications. The program contains many special features that allow nonlinearities or secondary effects such as plasticity, large strain, hyperelasticity, creep, swelling, large deflections, contact, stress stiffening, temperature dependency, material anisotropy, and radiation to be included in the solution. As ANSYS has been developed, other special capabilities such as substructuring, submodeling, random vibration, kinetostatics, kinetodynamics, free convection fluid analysis, acoustics, magnetics, piezoelectrics, coupled-field analysis, and design optimization have been added to the program. These capabilities contribute further to making ANSYS a multipurpose analysis tool for varied engineering disciplines. The ANSYS program has been in commercial use for over 50 years and has been used extensively in the aerospace, automotive, construction, electronic, energy services, manufacturing, nuclear, plastics, oil, and steel industries. Additionally, many consulting firms and hundreds of universities have used ANSYS for analysis, research, and educational purposes. ANSYS is recognized worldwide as one of the most widely used and capable programs of its type. Ansys design analysis software is the first created within a quality system with ISO 9001 certification, the internationally accepted quality standard. Product development, testing, maintenance and support processes also meet the United States Nuclear Regulatory Commission's quality requirements, as they have for nearly four decades. The Quality Assurance Service Agreement is suitable for the customers working in the nuclear industry who need to meet specific federal regulations including 10CRF50 Appendix B and provisions of 10CFR21. ANSYS has retained its original International Organization for Standardization (ISO) 9001 accreditation certificate since1995-05-04, It’s current certificate is valid until 2027-05-29.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Electron Microscopy Studies of Soft Nanomaterials

This review highlights recent efforts on applying electron microscopy (EM) to soft (including biological) nanomaterials. We will show how developments of both the hardware and software of EM have enabled new insights into the formation, assembly, and functioning (e.g., energy conversion and storage, phonon/photon modulation) of these materials by providing shape, size, phase, structural, and chemical information at the nanometer or higher spatial resolution. Specifically, we first discuss standard real-space two-dimensional imaging and analytical techniques which are offered conveniently by microscopes without special holders or advanced beam technology. The discussion is then extended to recent advancements, including visualizing three-dimensional morphology of soft nanomaterials using electron tomography and its variations, identifying local structure and strain by electron diffraction, and recording motions and transformation by in situ EM. On these advancements, we cover state-of-the-art technologies designed for overcoming the technical barriers for EM to characterize soft materials as well as representative application examples. Here, the even more recent integration of machine learning and its impacts on EM are also discussed in detail. With our perspectives of future opportunities offered at the end, we expect this review to inspire and stimulate more efforts in developing and utilizing EM-based characterization methods for soft nanomaterials at the atomic to nanometer length scales in academic research and industrial applications.

Imaging↗
DOE OSTI2024

GDB-9-Ex_TD-DFT-PBE0: Dataset containing Time Dependent Density Functional Theory (TDDFT) calculations for organic molecules of the GDB-9-Ex dataset.

This dataset contains data-intensive quantum chemical electronic structure calculations for 96,766 organic molecules of the GDB-9-Ex dataset. Calculations were performed using the Time Dependent Density Functional Theory (TDDFT) first principles method using the ORCA software. It provides UV-vis spectra calculations of molecules with a high level of accuracy. The optical spectra behavior was collected based on the optimized molecular geometries in the DFTB method with 3ob parameters. All calculations utilized the def2-TZVP basis sets with the auxiliary def2/J and def2-TZVP/C basis sets. The time-dependent density-functional theory (TDDFT) approach with the PBE0 exchange-correlation functional and ORCAs default integration grid was employed. For the excitation energy calculations, the lowest 50 excitation states were calculated.

AI dataset↗
DOE OSTI2022

libwfa: Wavefunction analysis tools for excited and open‐shell electronic states

Abstract An open‐source software library for wavefunction analysis, libwfa, provides a comprehensive and flexible toolbox for post‐processing excited‐state calculations, featuring a hierarchy of interconnected visual and quantitative analysis methods. These tools afford compact graphical representations of various excited‐state processes, provide detailed insight into electronic structure, and are suitable for automated processing of large data sets. The analysis is based on reduced quantities, such as state and transition density matrices (DMs), and allows one to distill simple molecular orbital pictures of physical phenomena from intricate correlated wavefunctions. The implemented descriptors provide a rigorous link between many‐body wavefunctions and intuitive physical and chemical models, for example, exciton binding, double excitations, orbital relaxation, and polyradical character. A broad range of quantum‐chemical methods is interfaced with libwfa via a uniform interface layer in the form of DMs. This contribution reviews the structure of libwfa and highlights its capabilities by several representative use cases. This article is categorized under: Software > Quantum Chemistry Theoretical and Physical Chemistry > Spectroscopy

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

Ground and excited state gradients with end-to-end differentiable semiempirical quantum chemistry

Accurate and efficient gradients of molecular energy with respect to nuclear degrees of freedom are essential for geometry optimization and molecular dynamics, including simulations that go beyond the Born–Oppenheimer regime. A common approach involves deriving analytical formulas for new electronic structure methods, which is often conceptually difficult and requires tedious coding. Here, we implement analytical, semi-numerical, and automatic differentiation (AD)-based gradient pathways for semiempirical Hamiltonian models in the PYSEQM software package, leveraging both graphics processing unit (GPU) and central processing unit (CPU) architectures. We further extend these capabilities to excited states calculated using the configuration interaction singles and time-dependent Hartree–Fock ansätze. We benchmark wall time, peak memory usage, and accuracy across three molecular families of varying chemical complexity, including systems of up to a thousand atoms. For ground-state simulations, analytical and AD gradients achieve near-identical GPU runtimes, while semi-numerical gradients are slower on GPU but remain competitive on CPU. For excited states, both analytical and custom AD approaches using implicit differentiation show similar performance and low memory requirements, whereas gradients with full AD are memory-limited. AD gradients match analytical ones in accuracy across all tested systems, aided by a quaternion-based diatomic frame rotation for two-center quantities that ensures smooth energy surfaces. Overall, automatic differentiation emerges as a practical alternative to analytical gradients in semiempirical quantum chemistry, offering high accuracy while allowing seamless integration in AI-driven workflows and popular packages, such as PyTorch and JAX. Our results provide actionable guidance for selecting optimal gradient strategies in large-scale ground- and excited-state molecular dynamics simulations.

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

GDB-9-Ex_EOM-CCSD: Dataset containing Equation of Motion Coupled Cluster (EOM-CCSD) calculations for organic molecules of the GDB-9-Ex dataset.

This dataset contains data-intensive quantum chemical electronic structure calculations for 80,593 organic molecules of the GDB-9-Ex dataset. Calculations were performed using the Equation of Motion Coupled Cluster (EOM-CCSD) first principles method using the ORCA software. It provides UV-vis spectra calculations of molecules with a high level of accuracy. The optical spectra behavior was collected based on the optimized molecular geometries in the DFTB method with 3ob parameters. All calculations utilized the def2-TZVP basis sets with the auxiliary def2/J and def2-TZVP/C basis sets. The similarity-transformed EOM-CCSD method that used domain-based local pair natural orbitals (DLPNO) approximation which constitutes the STEOM-DLPNO-CCSD method was used. This method is based on the STEOM approach and was found to make accurate predictions of transition energies for organic molecules. For the excitation energy calculations, the lowest 50 excitation states were calculated.

AI dataset↗
DOE OSTI2025

2025 ASMS Investigation of the Collision-Induced Dissociation Mechanism of Protonated TODGA with IRIS

Title (20 words): Investigation of the Collision-Induced Dissociation Mechanism of Protonated TODGA with IRIS Introduction (120 words): One of the challenges facing wide-spread adoption of nuclear power is the development of efficient separation processes for used nuclear fuel. The molecules in separation processes are subjected to an extreme environment due to the high radiation fields from the used fuel and highly acidic media used for fuel dissolution, which results in significant molecular degradation, leading to reduced process efficiency. These degradation products must be identified and studied so mitigation strategies can be developed to maintain process efficiency. However, complex systems can have many degradation products, complicating identification. Untargeted analysis tools could be used to understand radiation chemistry in complex systems. However, this would necessitate improved understanding of the gas-phase fragmentation mechanisms of fuel cycle molecules like tetraoctyldiglycolamide (TODGA). Methods (120 words): The gas-phase fragmentation of protonated TODGA was investigated using collision-induced dissociation (CID), resonance ejection, and infrared ion spectroscopy (IRIS). CID and resonance ejection experiments were conducted using a Bruker Daltonics (Bremen, Gemany) SolariX XR fourier transform ion cyclotron resonance (FT-ICR) mass spectrometer. IRIS spectra of protonated TODGA and its two CID fragmentation products were measured using a modified Bruker amaZon Speed ETD 3D quadrupole ion trap mass spectrometer coupled to the Free Electron Lasers for Infrared eXperiments (FELIX) free electron laser. Measured spectra were compared with density functional theory (DFT) calculations using the Gaussian 16, Revision C.02 software package with the ?B97X-D functional and def2-TZVPP basis sets. Candidate structures were generated using the CREST 3.0 conformational sampling software tool. Preliminary Data (300 words): Collision-induced dissociation of protonated TODGA ([C36H73N2O3]+, m/z=581.562) results two fragment ions, one at m/z=340.285 assigned as [C20H38NO3]+ and the other at m/z=312.290, assigned as [C19H38NO2]+. Based on the assigned formula and the structure of protonated TODGA, the fragment at m/z=340.285 is likely formed from elimination of neutral dioctylamine. Comparison of the IRIS spectrum of m/z=340.285 with DFT predictions suggests it contains a ring structure, and is assigned as N-octyl-N-(6-oxo-1,4-dioxan-2-ylidene)octan-1-aminium. Based on this structure and the structure of protonated TODGA, we hypothesize this fragment formed from elimination of neutral dioctylamine followed by a ring closure mechanism. Comparison of the IRIS spectrum of the fragment at m/z=312.290 with DFT predictions also indicated the presence of a ring structure, assigned as N-(1,3-dioxolan-4-ylidene)-N-octyloctan-1-aminium. This product could be formed from elimination of carbon monoxide from the ring of m/z=340.285 as a sequential fragmentation or formed directly from protonated TODGA via elimination of neutral N,N-dioctylformamide followed by a ring closure. Resonance ejection experiments where m/z=340 was continuously ejected from the IRC cell showed no decrease in intensity of m/z=312.290 across several collision energies, suggesting that the later, direct formation mechanism, dominates. The location of the ionizing proton in protonated TODGA is important for modeling the fragmentation mechanisms. DFT calculations suggested that the position of bands involving the coupled vibrations of the amide C—N and C=O bonds in TODGA are the most sensitive to proton location. Evaluation of the IRIS spectrum of protonated TODGA suggests that the ionizing proton is located between the two amid oxygens. This protonation location was calculated to lie approximately 30 kJ/mol lower in energy than the next lowest energy location, with the proton located solely on one of the amide oxygens. Novel aspect (20 words): Infrared ion spectroscopy combined with resonance ejection experiments and density functional theory to probe the collision-induced dissociation mechanism of tetraoctyldiglycolamide.

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

Nonadiabatic Excited-State Molecular Dynamics for Open-Shell Systems

Nonadiabatic Molecular Dynamics (NAMD) of excited states has been widely used in the simulation of photoinduced phenomena. However, the inability to treat bond breaking and forming processes with single-reference electronic structure methods limits their application in photochemistry for extended molecular systems. In this work, the extension of excited-state NAMD for open-shell systems is developed and implemented in the NEXMD software. We present the spin-unrestricted CIS and TDSCF formalism for the ground and excited states, analytical derivatives, and nonadiabatic derivative couplings for the respective potential energy surfaces. Next, this methodology is employed to study the photochemical reaction of three model molecules. The results demonstrate the advantage of the open-shell approach in modeling photochemical reactions, especially involving bond breaking processes. We find that the open-shell method lowers the reaction barrier at the bond-breaking limits resulting in larger calculated photochemical quantum yields compared to the respective closed-shell results. We also address problems related to spin contamination in the open-shell method, especially when molecular geometries are far from equilibrium.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Electron-phonon coupling from GW perturbation theory: Practical workflow combining BerkeleyGW, ABINIT, and EPW

Here we present a workflow of practical calculations of electron-phonon (e-ph) coupling with many-electron correlation effects included using the GW perturbation theory (GWPT). This workflow combines BerkeleyGW, ABINIT, and EPW software packages to enable accurate e-ph calculations at the GW self-energy level, going beyond standard calculations based on density functional theory (DFT) and density-functional perturbation theory (DFPT). This workflow begins with DFT and DFPT calculations (ABINIT) as starting point, followed by GW and GWPT calculations (BerkeleyGW) for the quasiparticle band structures and e-ph matrix elements on coarse electron k- and phonon q-grids, which are then interpolated to finer grids through Wannier interpolation (EPW) for computations of various e-ph coupling determined physical quantities such as the electron self-energies or solutions of anisotropic Eliashberg equations, among others. A gauge-recovering symmetry unfolding technique is developed to reduce the computational cost of GWPT (as well as DFPT) while fulfilling the gauge consistency requirement for Wannier interpolation.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Spinor $GW$ Bethe-Salpeter calculations in BerkeleyGW: Implementation, symmetries, benchmarking, and performance

Computing the GW quasiparticle band structure and Bethe-Salpeter equation (BSE) absorption spectra for materials with spin-orbit coupling have commonly been done by treating GW corrections and spin-orbit coupling (SOC) as separate perturbations to density-functional theory. However, accurate treatment of materials with strong spin-orbit coupling (such as many topological materials of recent interest, and thermoelectrics) often requires a nonperturbative approach using spinor wave functions in the Kohn-Sham equation and GW/BSE. Such calculations have only recently become available, in particular for the BSE. Here, we have implemented this approach in the plane-wave pseudopotential GW/BSE code BerkeleyGW, which is highly parallelized and widely used in the electronic-structure community. We present reference results for quasiparticle band structures and optical absorption spectra of solids with different strengths of spin-orbit coupling, including Si, Ge, GaAs, GaSb, CdSe, Au, and Bi 2 Se 3 . The calculated quasiparticle band gaps of these systems are found to agree with experiment to within a few tens of meV. SOC splittings are found to be generally in better agreement with experiment, including quasiparticle corrections to band energies. The absorption spectrum of GaAs is not significantly impacted by the inclusion of spin-orbit coupling due to its relatively small value (0.2 eV) in the Λ direction, while the absorption spectrum of GaSb calculated with the spinor GW/BSE captures the large spin-orbit splitting of peaks in the spectrum. For the prototypical topological insulator Bi 2 Se 3 , we find a drastic change in the low-energy band structure compared to that of DFT, with the spinorial treatment of the GW approximation correctly capturing the parabolic nature of the valence and conduction bands after including off-diagonal self-energy matrix elements. We present the detailed methodology, approach to spatial symmetries for spinors, comparison against other codes, and performance compared to spinless GW/BSE calculations and perturbative approaches to SOC. This work aims to spur further development of spinor GW/BSE methodology in excited-state research software and enables a more accurate and detailed exploration of electronic and optical properties of materials containing elements with large atomic numbers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

JARVIS-Leaderboard: a large scale benchmark of materials design methods

Abstract Lack of rigorous reproducibility and validation are significant hurdles for scientific development across many fields. Materials science, in particular, encompasses a variety of experimental and theoretical approaches that require careful benchmarking. Leaderboard efforts have been developed previously to mitigate these issues. However, a comprehensive comparison and benchmarking on an integrated platform with multiple data modalities with perfect and defect materials data is still lacking. This work introduces JARVIS-Leaderboard, an open-source and community-driven platform that facilitates benchmarking and enhances reproducibility. The platform allows users to set up benchmarks with custom tasks and enables contributions in the form of dataset, code, and meta-data submissions. We cover the following materials design categories: Artificial Intelligence (AI), Electronic Structure (ES), Force-fields (FF), Quantum Computation (QC), and Experiments (EXP). For AI, we cover several types of input data, including atomic structures, atomistic images, spectra, and text. For ES, we consider multiple ES approaches, software packages, pseudopotentials, materials, and properties, comparing results to experiment. For FF, we compare multiple approaches for material property predictions. For QC, we benchmark Hamiltonian simulations using various quantum algorithms and circuits. Finally, for experiments, we use the inter-laboratory approach to establish benchmarks. There are 1281 contributions to 274 benchmarks using 152 methods with more than 8 million data points, and the leaderboard is continuously expanding. The JARVIS-Leaderboard is available at the website: https://pages.nist.gov/jarvis_leaderboard/

36 MATERIALS SCIENCE↗
DOE OSTI2023

Scalable generalized screening for high-order terms in the many-body expansion: Algorithm, open-source implementation, and demonstration

The many-body expansion lies at the heart of numerous fragment-based methods that are intended to sidestep the nonlinear scaling of ab initio quantum chemistry, making electronic structure calculations feasible in large systems. In principle, inclusion of higher-order n-body terms ought to improve the accuracy in a controllable way, but unfavorable combinatorics often defeats this in practice and applications with n ≥ 4 are rare. Here, we outline an algorithm to overcome this combinatorial bottleneck, based on a bottom-up approach to energy-based screening. This is implemented within a new open-source software application (“Fragme∩t”), which is integrated with a lightweight semi-empirical method that is used to cull subsystems, attenuating the combinatorial growth of higher-order terms in the graph that is used to manage the calculations. This facilitates applications of unprecedented size, and we report four-body calculations in (H2O)64 clusters that afford relative energies within 0.1 kcal/mol/monomer of the supersystem result using less than 10% of the unique subsystems. We also report n-body calculations in (H2O)20 clusters up to n = 8, at which point the expansion terminates naturally due to screening. These are the largest n-body calculations reported to date using ab initio electronic structure theory, and they confirm that high-order n-body terms are mostly artifacts of basis-set superposition error.

Chemistry↗