Search NASA⌕ Search

SEARCH · Search NASA

Results for “NWChem”

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 OSTI2025

SparcleQC: Automated Input File Creation for QM/MM Studies of Protein:Ligand Complexes

SparcleQC is a Python package that, given a protein:ligand complex in the Protein Data Bank (PDB) file format, can create quantum mechanics/molecular mechanics (QM/MM)-like input files for the electronic structure theory packages PSI4, QChem, and NWChem. The resulting input files include quantum mechanical representations of the ligand and a small section of the protein, surrounded by point charges that represent the rest of the protein. Creation of these QM/MM input files includes cutting and capping the QM subregion, obtaining point charges for the protein, and adjusting charges at the QM/MM boundary; and each of these tasks are automated by the software. In this article, we describe the details of SparcleQC’s procedure, show examples of the Python API, and explain additional features that are helpful in protein:ligand interaction studies. Finally, we show that SparcleQC enables automated preparation of input files for QM/MM calculations, which can return can return accurate interaction energies in minutes, while a fully quantum mechanical computation on the protein:ligand complex could take days, if it is even possible.

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

WFOT: A Wave Function Overlap Tool between Single- and Multi-Reference Electronic Structure Methods for Spectroscopy Simulation

We report the development of a novel diagnostic tool, named wave function overlap tool (WFOT), designed to evaluate the overlap between wave functions computed at single-reference [i.e., time-dependent density functional theory or configuration interaction singles (CIS)] and multireference (i.e., CASSCF/CASPT2) electronic structure levels of theory. It relies on truncating the single- and multireference WFs to CIS-like expansions spanning the same configurational space and maximizing the molecular orbital overlap by means of a unitary transformation. To demonstrate the functionality of the tool, we calculate the transient spectrum of acetylacetone by evaluating excited state absorption signals with multireference quality on top of single-reference on-the-fly dynamics simulations. Semiautomatic spectra generation is facilitated by interfacing the tool with the COBRAMM package, which also allows one to use WFOT with several quantum chemistry codes such as Gaussian, NWChem, and OpenMolcas. Finally, other exciting possibilities for the utilization of the code beyond the simulation of transient absorption spectroscopy are eventually discussed.

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

Quantum Chemistry Common Driver and Databases (QCDB) and Quantum Chemistry Engine (QCEngine): Automation and interoperability among computational chemistry programs

We report that community efforts in the computational molecular sciences (CMS) are evolving toward modular, open, and interoperable interfaces that work with existing community codes to provide more functionality and composability than could be achieved with a single program. The Quantum Chemistry Common Driver and Databases (QCDB) project provides such capability through an application programming interface (API) that facilitates interoperability across multiple quantum chemistry software packages. In tandem with the Molecular Sciences Software Institute and their Quantum Chemistry Archive ecosystem, the unique functionalities of several CMS programs are integrated, including CFOUR, GAMESS, NWChem, OpenMM, Psi4, Qcore, TeraChem, and Turbomole, to provide common computational functions, i.e., energy, gradient, and Hessian computations as well as molecular properties such as atomic charges and vibrational frequency analysis. Both standard users and power users benefit from adopting these APIs as they lower the language barrier of input styles and enable a standard layout of variables and data. These designs allow end-to-end interoperable programming of complex computations and provide best practices options by default.

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

Distributed-Memory Sparse Deep Neural Network Inference Using Global Arrays

Partitioned Global Address Space (PGAS) models exhibit tremendous promise in developing efficient and productive distributed-memory parallel applications. They have been used extensively in scientific computations due to conveniently offering a ``shared-memory''-like model and convenient interfaces that separate communication with synchronization. Traditionally, PGAS communication models have been applied to dense/contiguously distributed data, but most modern applications depict varied levels of sparsity. Existing PGAS models require certain adaptations to support distributed sparse computations, since associated computations often require matrix arithmetic, in addition to data movement. The Global Arrays toolkit from Pacific Northwest National Laboratory (PNNL) is one of the earliest PGAS models to combine one-sided data communication and distributed matrix operations and is still used in the popular NWChem quantum chemistry suite. Recently, we have expanded the Global Arrays toolkit to support common sparse operations, like sparse matrix-dense matrix multiplies (SpMM), sparse matrix-sparse matrix multiplication (SpGEMM) and Sampled Dense-Dense Matrix Multiplication (SDDMM). As it turns out, these operations are the bedrock of sparse Deep Learning (DL); sparse deep neural networks and Graph Neural Networks (GNNs) have gained increasing attention recently in achieving speedups on training and inference with reduced memory footprints. Unlike scientific applications in High Performance Computing (HPC), modern (distributed-memory capable) DL toolkits often rely on non-standardized and closed-source vendor software optimizations, creating challenges in software-hardware co-design at scale. Our goal is to support a variety of distributed-memory sparse matrix operations and helper functions in the newly created Sparse Global Arrays (SGA), such that it is possible to build portable and productive Machine Learning scenarios for algorithm/software and hardware codesign purposes. Contemporary data-parallel schemes for training/inference are undergoing a major overhaul since model replication limits scalability and causes resource inefficiencies. As such, we have adopted tensor parallelism in decomposing the model and inputs, to mitigate memory issues. Current implementation is built on top of MPI and uses CPUs to maximize the portability across the platforms.

Distributed computing, machine learning↗
DOE OSTI2022

Gas‐Phase Ion‐Molecule Interactions in a Collision Reaction Cell with ICP‐MS / MS : Investigations with CO 2 as the Reaction Gas

Carbon dioxide (CO 2 ) was used as a reaction gas to investigate the gas‐phase ion‐molecule interactions using the Agilent 8900 ICP‐MS/MS. A solution containing forty‐five elements representative of the periodic table was used to supply the ions to react with CO 2 in the collision/reaction cell (CRC). The only significant product ions formed were monoxides. The general reactivity was shown to be consistent with density functional theory (DFT)‐predicted reaction enthalpies, such that all predicted exothermic reactions produced product ions at levels of at least 1% of the unreacted ion. Most endothermic reactions observed had sufficient kinetic energy in excess of the reaction enthalpies. Our results suggest that reaction enthalpy is a reasonable predictor of reactivity with CO 2 on the timescales of the interactions in non‐thermal ICP‐MS/MS systems. The ease and rapidity of data collection with the ICP‐MS/MS and DFT calculations using the NWChem suite has value given the scarcity of thermochemical data of CO 2 reactions in the literature. These studies are especially useful for the identification of targeted reaction chemistries to be leveraged for analytical method development, such as for the inline separation of isobaric interferences from analytes of interest.

36 MATERIALS SCIENCE↗
DOE OSTI2023

pnnl/PTMPSI

PTM-Psi is a Python 3 package that combines several capabilities to streamline the workflow to interrogate the impact of PTMs on proteins using well-established software packages. The workflow of the PTM-Psi software package includes input files and launch instances from standard packages such as AlphaFold, NWChem, GROMACS, and the Autodock Suite

Mejia-Rodriguez, Daniel↗
DOE OSTI2024

pnnl/emp-gnn

Efficient Graph Neural Network for Predicting Molecular Properties software can compare the prediction quality with ab initio DFT results reported by the high-performance state-of-theart NWChem quantum chemistry package [1] through Mean Absolute Error (MAE) obtained by the fitting between DFT and model predictions i.e, MPNN [2]. We also demonstrated the performance benefits by grouping large molecules by atom sizes, and executing GNN models on different types of resources. Since the training times depend on the number of atoms, we demonstrate the impact of distributing the workloads on two GPUs with varying capabilities (e.g., NVIDIA A100 vs. GeForce RTX 2080 Ti) to optimize the efficiency. We are at the precipice of broad adoption of GNNs for molecular property prediction tasks; hence, our work is timely by comparing model prediction against classical approaches with the intent of providing initial screening for specific classes of molecules.

Lee, Hyungro↗
DOE OSTI2020

An Atmospheric Aerosol Short Course for Early Career Scientists

The first Department of Energy (DOE) Environmental and Molecular Science Laboratory (EMSL)/Atmospheric Radiation Measurement (ARM) aerosol summer school was held July 15-19, 2019 at the EMSL user facility located at Pacific Northwest National Laboratory. The school set out to provide graduate students and early-career scientists an overview of atmospheric aerosol science and an introduction to the measurement techniques available through the EMSL and ARM user facilities to address current aerosol science issues. Instructors for the course were drawn from scientists from the EMSL (Laskin et al., 2016; China et al., 2018) and ARM facilities (McComiskey and Ferrare, 2016) as well as from the DOE Atmospheric System Research (ASR) program. Student support for the summer school was also jointly provided by EMSL, ARM, and ASR. The school spanned five days with each day covering a particular theme. 25 students from 21 institutions attended morning lectures on aerosol composition, properties, theory, measurements, and modeling. During afternoons, the students explored instruments that collect and analyze aerosol samples at EMSL and PNNL’s Atmospheric Measurements Laboratory. They also tested relevant computer modeling simulations and learned about NWChem, EMSL’s open-source computational chemistry software suite.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Theory and Simulation of Ultrafast Multidimensional Nonlinear X-ray Spectroscopy of Molecules (Final Report)

Emerging X-ray free electron laser (XFEL) beam sources offer new types of probes of matter with unprecedented spatial and temporal resolutions. These experimental advances must be met by robust theoretical and computational tools that provide predictive modeling capacity of the underlining electronic and structural dynamics. The latter will be essential for the design of sophisticated multi-pulse experiments and for their interpretation. The proposed research effort will focus on developing cutting- edge simulation tools for nonlinear multidimensional X-ray/optical spectroscopies and aims to address key questions in Priority Research Opportunities 1 (Probing and controlling electron motion within a molecule) and 3 (Capturing rare events and intermediate states in the transformation of matter) as mentioned in the BES Roundtable Report “Opportunities for Basic Research at the Frontiers of XFEL Ultrafast Science”. XFEL multidimensional nonlinear techniques, which combine sequences of X-ray and possibly optical pulses, provide a unique experimental toolbox for probing the dynamics of core and valence electronic excitations, as well as material structure. Predictive modeling of these dynamical processes requires the combination of analytical theory for nonlinear interactions of light and matter, robust quantum-chemical methodologies for the accurate description of electronic structure of various materials, and multiscale ab initio electron and nuclear dynamics techniques operating beyond Born-Oppenheimer approximation. These challenges will be addressed with three research thrusts (i) Develop and implement theoretical apparatus for modeling a broad range of multidimensional spectroscopic techniques enabled by present and upcoming XFEL facilities. This thrust also includes the incorporation of a computational module in the DOE supported open-source NWChem computational chemistry package as well as the development of other open-source codes ready for dissemination across a broad user base; (ii) Propose and design new multi-pulse experiments that make use of the capabilities of the incoming LCLS-II facility; (iii) Perform selected applications to specific molecular systems that can be carried out at LCLS-II and demonstrate how these X-ray sources may be used to study nonadiabatic dynamics through conical interactions, electronic correlations in multi-core excitons, and charge transfer/energy transfer processes. The proposed research will be carried out by a multi-disciplinary four-institution research team which combines academia and national laboratories and spans the broad and necessary expertise in theoretical spectroscopy, nonlinear optics, quantum chemistry, molecular non-adiabatic dynamics and code development. The work will be performed in a highly interactive team environment with junior researchers shared between institutions thus cementing cross- disciplinary interactions. The developed simulation tools will be immediately deployed for XFEL facility users, both experimentalists and theorists, via freely distributed codes and databases. Altogether, this project will facilitate establishing XFEL-based multidimensional spectroscopies as a novel diagnostic tool for monitoring electronic and structural dynamics in molecular materials.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

DOE STTR Phase I Final Report Report: Machine-learning Based Prediction of Thermal Limits for Conjugated Organic Materials

The SCANN-DT Phase I STTR project led by NLM Photonics and partnering with the National Renewable Energy Laboratory (NREL) sought to apply machine learning techniques based on graph neural networks (GNNs) towards the prediction of decomposition temperatures (Td) of organic semiconductors, based on prior work on bond dissociation energy (BDE) prediction as implemented in NREL’s ALFABET prediction tool. Using a curated set of experimental decomposition energies, the project examined GNN-based, classical quantitative structure-property relationship (QSPR) based on DFT calculations, and combinations of both methods to predict Td. While the best MAEs in Td achieved were near 40°C, below project targets, the project led to improvements in the ALFABET model, improvements in cloud-based implementations of NWChem software, and a substantial dataset of calculations on medium-sized conjugated organic molecules.

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

Cytochrome Nanowire XRD

Microbial metabolism plays an important role in the global cycling of carbon, nutrients, and metals. Recent discovery of microorganisms sharing energy using direct electrical connections might represent a prominent strategy by which anaerobic microorganisms interact with the diversity of environments. The purpose of the proposed research is to apply these newly developed principles of electron exchange in microbial communities to the study of biogeochemical cycling in climate and subsurface systems. Previous studies revealed that soil bacteria Geobacter sulfurreducens produce conductive pili nanofilaments that facilitate long-range electron transport to extracellular electron acceptors and other cells. This discovery has already transformed our understanding of the function of microbial communities in diverse environments, making basic contributions to microbial ecology and helping to improve practical applications such as increasing current output of microbial fuel cells and enhancing the conversion of organic waste to methane. Microbial electron exchange via pili may also contribute to methane production in terrestrial environments that are major sources of atmospheric methane. However, the mechanism of electron transport through pili is not well understood. Initial studies have suggested that pili exhibit metallic-like conductivity similar to conducting polymer polyaniline. However, these studies were limited to electrical measurements. X-ray diffraction studies have suggested that aromatic amino acids tyrosine and phenylalanine are closely packed in pili. But these studies used homology modeling to predict the pilus assembly and information about the number of monomers in each pilus assembly unit and their organization are still lacking. The central hypothesis of this proposal is aromatics in pili enable intermolecular electron delocalization due to pi-stacking that can give rise to metallic-like conductivity. By analyzing pili structure using infrared scattering-scanning near field optical microscope (IR s-SNOM) and combining with solution NMR, Focused Ion Beam (FIB and computational capabilities, this work will provide comprehensive understanding of the conduction mechanism from a structure-function perspective. The specific aims are: 1) With infrared scattering scanning near field optical microscope (IR s-SNOM), image the organization and density of pilin monomers and pH-induced conformational changes in a pilus filament. 2) With solution NMR, resolve pH-induced structural changes in the pilin monomer. 3) Using Focused Ion Beam (FIB), draw electrodes on single pili to measure their conductivity. 4) With computational modeling using the NWChem, simulate the effect of environment on the electronic properties of pili as a function of pH. The proposed work will gain the nanoscale insight into the structure of the pili, as the basis for understanding interspecies electrical communications. The improved understanding will be harnessed to accelerate bioenergy production and bioremediation using pili as well as for predictive modeling of carbon cycling in diverse environments to reach the central goals of DOE-BER to develop sustainable energy sources, regulate the contaminants in the subsurface and understand the effects of greenhouse gas emissions on biosphere. The proposed studies will provide design principles for engineering of electrical interactions among these communities by manipulating the amino acid composition of pili. Moreover, these studies will identify the forces that maintain the function of these communities in the face of environmental perturbation. Therefore, these studies will help to address the need "to incorporate process understanding of biogeochemical cycling into a scalable hierarchy of predictive capabilities" specified in EMSL's Science Themes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

On the Efficient Evaluation of the Exchange Correlation Potential on Graphics Processing Unit Clusters

The predominance of Kohn–Sham density functional theory (KS-DFT) for the theoretical treatment of large experimentally relevant systems in molecular chemistry and materials science relies primarily on the existence of efficient software implementations which are capable of leveraging the latest advances in modern high-performance computing (HPC). With recent trends in HPC leading toward increasing reliance on heterogeneous accelerator-based architectures such as graphics processing units (GPU), existing code bases must embrace these architectural advances to maintain the high levels of performance that have come to be expected for these methods. In this work, we purpose a three-level parallelism scheme for the distributed numerical integration of the exchange-correlation (XC) potential in the Gaussian basis set discretization of the Kohn–Sham equations on large computing clusters consisting of multiple GPUs per compute node. In addition, we purpose and demonstrate the efficacy of the use of batched kernels, including batched level-3 BLAS operations, in achieving high levels of performance on the GPU. We demonstrate the performance and scalability of the implementation of the purposed method in the NWChemEx software package by comparing to the existing scalable CPU XC integration in NWChem.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2005

Insulator Surface Charge as a Function of Pressure: Theory and Simulation

A two-phase equilibrium model was developed to explain the discontinuous surface charge decay versus atmospheric pressure of insulators that had been charged triboelectrically. The two-phase model is an electrostatic form of the Langmuir Isotherm for ions adsorbed on a surface in equilibrium with ions in the gas phase. In this paper, the model was extended to account for vibrational states of the adsorbed surface ions via the vibrational partition function. An analysis is performed that rules out Paschen discharge as the cause of the discharge observed. Also, a numerical simulation is performed using NWChem to calculate the adsorption energies of ions on insulator surfaces for comparison to curve fit adsorption energies developed from the model and experimental data.

Hogue, Michael D.↗
DOE OSTI2020

Scalable Heterogeneous Execution of a Coupled-Cluster Model with Perturbative Triples

The CCSD(T) coupled-cluster model with perturbative triples is considered a gold standard for computational modeling of the correlated behavior of electrons in molecular systems. A fundamental constraint is the relatively small global-memory capacity in GPUs compared to the main-memory capacity on host nodes, necessitating relatively smaller tile sizes for high-dimensional tensor contractions in NWChem's GPU-accelerated implementation of the CCSD(T) method. A coordinated redesign is described to address this limitation and associated data movement overheads, including a novel fused GPU kernel for a set of tensor contractions, along with inter-node communication optimization and data caching. The new implementation of GPU-accelerated CCSD(T) improves overall performance by 3.4x. Finally, we discuss the trade-offs in using this fused algorithm on current and future supercomputing platforms.

Kim, Jinsung↗
DOE OSTI2020

Electron transfer calculations between edge sharing octahedra in hematite, goethite, and annite

A key reaction underlying the charge transport in iron containing oxides, clays, micas is the Fe$^{2+}$-Fe$^{3+}$ exchange reaction between edge-sharing iron octahedra. These reactions facilitate conduction in these minerals by the thermally-activated hopping of small polarons across the lattice. Depending on the mineral and local charge state the small polaron can either encase an electron or hole. The probability for conduction of small polarons depends strongly on the height and adiabicity of the reaction barrier, with larger and more diabatic barriers yielding slow conduction associated with either weak coupling or a large prerequisite rearrangement of the lattice during charge transport. To model these reactions, a first principle electron transfer (ET) method was developed to model the small polaron hopping between the edge-sharing octahedra sites in hematite ($e^{-}$ polaron), goethite ($e^{-}$ polaron), and annite ($h^{+}$ polaron) bulk structures. The ET method is based on electronic structure methods (i.e., plane-wave Density Functional Theory) capable of performing calculations with periodic cells and large size systems efficiently while at the same time being accurate enough to be used in the estimation of the electron-transfer coupling matrix element, $V_{AB}$, and the electron transfer transmission factor, $\kappa_{el}$. Additionally, the calculations confirmed the existence of small polarons in all three minerals, and the reactions were predicted to be strongly adiabatic. It was found that transfer of a hole in the octahedral layer of annite had an adiabatic barrier of $0.311$ eV, and the transfer of an extra electron in hematite and goethite had adiabatic barriers of $0.242$ eV and $0.232$ eV respectively. The electronic coupling parameters,$V_{AB}$, were found to be $0.188$ eV, $0.196$ eV, and $0.102$ eV respectively for hematite, goethite, and annite. While similar bonding topologies pertain, the underlying basis for the differences is the subtle differences in local structures.

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

A Filon-like integration strategy for calculating exact exchange in periodic boundary conditions: a plane-wave DFT implementation

An efficient and accurate approach for calculating exact exchange and other two-electron integrals has been developed for periodic electronic structure methods. Traditional approaches used for integrating over the Brillouin zone in band structure calculations, e.g. trapezoidal or Monkhorst-Pack, are not accurate enough for two-electron integrals. This is because their integrands contain multiple singularities over the double integration of the Brillouin zone, which with simple integration methods lead to very inaccurate results. A common approach to this problem has been to replace the Coulomb interaction with a screened Coulomb interaction that removes singularities from the integrands in the two-electron integrals, albeit at the inelegance of having to introduce a screening factor which must precomputed or guessed. Instead of introducing screened Coulomb interactions in an ad hoc way, the method developed in this work derives an effective screened potential using a Filon-like integration approach that is based only on the lattice parameters. This approach overcomes the limitations of traditionally defined screened Coulomb interactions for calculating two-electron integrals, and makes chemistry many-body calculations tractable in periodic boundary conditions. This method has been applied to several systems for which conventional DFT methods do not work well, including the reaction pathways for the addition of H 2 to phenol and Au$^{–}_{20}$ nanoparticle, and the electron transfer of a charge trapped state in the Fe(II) containing mica, annite.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Quantum Solvers for Plane-Wave Hamiltonians: Abridging Virtual Spaces Through the Optimization of Pairwise Correlations

For many-body methods such as MCSCF and CASSCF, in which the number of one-electron orbitals is optimized and independent of the basis set used, there are no problems with using plane-wave basis sets. However, for methods currently used in quantum computing such as select configuration interaction (CI) and coupled cluster (CC) methods, it is necessary to have a virtual space that is able to capture a significant amount of electron-electron correlation in the system. The virtual orbitals in a pseudopotential plane-wave Hartree–Fock calculation, because of Coulomb repulsion, are often scattering states that interact very weakly with the filled orbitals. As a result, very little correlation energy is captured from them. The use of virtual spaces derived from the one-electron operators has also been tried, and while some correlations are captured, the amount is quite low. To overcome these limitations, we have been developing new classes of algorithms to define virtual spaces by optimizing orbitals from small pairwise CI Hamiltonians, which we term as correlation optimized virtual orbitals with the abbreviation COVOs. With these procedures, we have been able to derive virtual spaces, containing only a few orbitals, which are able to capture a significant amount of correlation. The focus in this manuscript is on using these derived basis sets to target full CI (FCI) quality results for H2 on near-term quantum computers. However, the initial results for this approach were promising. We were able to obtain good agreement with FCI/cc-pVTZ results for this system with just 4 virtual orbitals, using both FCI and quantum simulations. The quality of the results using COVOs suggests that it may be possible to use them in other many-body approaches, including coupled cluster and Møller–Plesset perturbation theories, and open up the door to many-body calculations for pseudopotential plane-wave basis set methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗