Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parallel computational chemistry”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Valence : A Massively Parallel Implementation of the Variational Subspace Valence Bond Method

This work describes the software package, Valence, for the calculation of molecular energies using the variational subspace valence bond (VSVB) method. VSVB is an ab initio electronic structure method based on nonorthogonal orbitals. Important features of practical value include high parallel scalability, wave functions that can be constructed automatically by combining orbitals from previous calculations, and ground and excited states that can be modeled with a single configuration or determinant. Here, the open-source software package includes tools to generate wave functions, a database of generic orbitals, example input files, and a library build intended for integration with other packages. We also present the interface to an external software package, enabling the computation of optimized molecular geometries and vibrational frequencies.

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

tih_vqe [SWR-23-32]

This software supports the paper, "Exploring the scaling limitations of the variational quantum eigensolver with the bond dissociation of hydride diatomic molecules," published in the International Journal of Quantum Chemistry, whose abstract is as follows: Materials simulations involving strongly correlated electrons pose fundamental challenges to state-of-the-art electronic structure methods but are hypothesized to be the ideal use case for quantum computing. To date, no quantum computer has simulated a molecule of a size and complexity relevant to real-world applications, despite the fact that the variational quantum eigensolver (VQE) algorithm can predict chemically accurate total energies. Nevertheless, because of the many applications of moderately-sized, strongly correlated systems, such as molecular catalysts, the successful use of the VQE stands as an important waypoint in the advancement toward useful chemical modeling on near-term quantum processors. In this paper, we take a significant step in this direction. We lay out the steps, write, and run parallel code for an (emulated) quantum computer to compute the bond dissociation curves of the TiH, LiH, NaH, and KH diatomic hydride molecules using VQE. TiH was chosen as a relatively simple chemical system that incorporates d orbitals and strong electron correlation. Because current VQE implementations on existing quantum hardware are limited by qubit error rates, the number of qubits available, and the allowable gate depth, recent studies have focused on chemical systems involving s and p block elements. Through VQE + UCCSD calculations of TiH, we evaluate the near-term feasibility of modeling a molecule with d-orbitals on real quantum hardware. We demonstrate that the inclusion of d-orbitals and the use of the UCCSD ansatz, which are both necessary to capture the correct TiH physics, dramatically increase the cost of this problem. We estimate the approximate error rates necessary to model TiH on current quantum computing hardware using VQE+UCCSD and show them to likely be prohibitive until significant improvements in hardware and error correction algorithms are available.

Graf, Peter↗
DOE OSTI2016

Morphological Expression of the Coherence and Relative Phase of Optical Inputs to the Photoelectrodeposition of Nanopatterned Se–Te Films

Highly anisotropic and ordered nanoscale lamellar morphologies can be spontaneously generated over macroscopic areas, without the use of a photomask or any templating agents, via the photoelectrodeposition of Se–Te alloy films. To form such structures, the light source can be a single, linearly polarized light source that need not necessarily be highly coherent. In this work, the variation in the morphologies produced by this deposition process was evaluated in response to differences in the coherence and relative phase between multiple simultaneous linearly polarized illumination inputs. Specifically, the morphologies of photoelectrodeposits were evaluated when two tandem same-wavelength sources with discrete linear polarizations, both either mutually incoherent or mutually coherent (with defined phase differences), were used. Additionally, morphologies were simulated via computer modeling of the interfacial light scattering and absorption during the photoelectrochemical growth process. The morphologies that were generated using two coherent, in-phase sources were equivalent to those generated using only a single source. In contrast, the use of two coherent, out-of-phase sources produced a range of morphological patterns. For small out-of-phase addition of orthogonal polarization components, lamellar-type patterns were observed. When fully out-of-phase orthogonal sources (circular polarization) were used, an isotropic, mesh-type pattern was instead generated, similar to that observed when unpolarized illumination was utilized. In intermediate cases, anisotropic lamellar-type patterns were superimposed on the isotropic mesh-type patterns, and the relative height between the two structures scaled with the amount of out-of-phase addition of the orthogonal polarization components. Similar results were obtained when two incoherent sources were utilized. In every case, the long axis of the lamellar-type morphology component aligned parallel to the intensity-weighted average polarization orientation. The observations consistently agreed with computer simulations, indicating that the observed morphologies were fully determined by the nature of the illumination utilized during the growth process. Here, the collective data thus indicated that the photoelectrodeposition process exhibits sensitivity toward the coherency, relative phase, and polarization orientations of all optical inputs and that this sensitivity is physically expressed in the morphology of the deposit.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2013

An atomic layer deposition reactor with dose quantification for precursor adsorption and reactivity studies

An atomic layer deposition reactor has been constructed with quantitative, precision dose control for studying precursor adsorption characteristics and to relate dose quantity and exposure dynamics to fluid flow in both the viscous and molecular flow regimes. A fixed volume of gas, held at a controlled temperature and measured pressure, is dosed into the reaction chamber by computer-controlled pneumatic valves. Dual in situ quartz crystal microbalances provide parallel mass measurement onto two differently coated substrates, which allows adsorption coverage and relative sticking coefficients to be determined. Gas composition in the reaction chamber was analyzed in situ by a quadrupole mass spectrometer. Absolute reactant exposure is unambiguously calculated from the impingement flux, and is related to dose, surface area, and growth rates. A range of control over the dose amount is demonstrated and consequences for film growth control are demonstrated and proposed.

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

Performance of explicit and IMEX MRI multirate methods on complex reactive flow problems within modern parallel adaptive structured grid frameworks

Large-scale multiphysics simulations are computationally challenging due to the coupling of multiple processes with widely disparate time scales. The advent of exascale computing systems exacerbates these challenges since these systems enable ever-increasing size and complexity. In recent years, there has been renewed interest in developing multirate methods as a means to handle the large range of time scales, as these methods may afford greater accuracy and efficiency than more traditional approaches of using implicit-explicit (IMEX) and low-order operator splitting schemes. However, to date there have been few performance studies that compare different classes of multirate integrators on complex application problems. In this work, we study the performance of several newly developed multirate infinitesimal (MRI) methods, implemented in the SUNDIALS solver package, on two reacting flow model problems built on structured mesh frameworks. The first model revisits prior work on a compressible reacting flow problem with complex chemistry that is implemented using BoxLib but where we now include comparisons between a new explicit MRI scheme with the multirate spectral deferred correction (SDC) methods in the original paper. The second problem uses the same complex chemistry as the first problem, combined with a simplified flow model, but runs at a large spatial scale where explicit methods become infeasible due to stability constraints. Two recently developed IMEX MRI multirate methods are tested. These methods rely on advanced features of the AMReX framework on which the model is built, such as multilevel grids and multilevel preconditioners. The results from these two problems show that MRI multirate methods can offer significant performance benefits on complex multiphysics application problems and that these methods may be combined with advanced spatial discretization to compound the advantages of both.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2012

Runtime Techniques to Enable a Highly-Scalable Global Address Space Model for Petascale Computing

Over the past decade, the trajectory to the petascale has been built on increased complexity and scale of the underlying parallel architectures. Meanwhile, software developers have struggled to provide tools that maintain the productivity of computational science teams using these new systems. In this regard, Global Address Space (GAS) programming models provide a straightforward and easy to use addressing model, which can lead to improved productivity. However, the scalability of GAS depends directly on the design and implementation of the runtime system on the target petascale distributed-memory architecture. In this paper, we describe the design, implementation, and optimization of the Aggregate Remote Memory Copy Interface (ARMCI) runtime library on the Cray XT5 2.3 PetaFLOPs computer at Oak Ridge National Laboratory. We optimized our implementation with the flow intimation technique that we have introduced in this paper. Our optimized ARMCI implementation improves scalability of both the Global Arrays programming model and a real-world chemistry application—NWChem—from small jobs up through 180,000 cores.

Computer Science↗
DOE OSTI2012

A DISLOCATION-BASED CLEAVAGE INITIATION MODEL FOR PRESSURE VESSEL

Efforts are under way to develop a theoretical, multi-scale model for the prediction of fracture toughness of ferritic steels in the ductile-to-brittle transition temperature (DBTT) region that accounts for temperature, irradiation, strain rate, and material condition (chemistry and heat treatment) effects. This new model is intended to address difficulties associated with existing empirically-derived models of the DBTT region that cannot be extrapolated to conditions for which data are unavailable. Dislocation distribution equations, derived from the theories of Yokobori et al., are incorporated to account for the local stress state prior to and following initiation of a microcrack from a second-phase particle. The new model is the basis for the DISlocation-based FRACture (DISFRAC) computer code being developed at the Oak Ridge National Laboratory (ORNL). The purpose of this code is to permit fracture safety assessments of ferritic structures with only tensile properties required as input. The primary motivation for the code is to assist in the prediction of radiation effects on nuclear reactor pressure vessels, in parallel with the EURATOM PERFORM 60 project.

Cochran, Kristine B↗
DOE OSTI2022

Constant-potential molecular dynamics simulations of molten salt double layers for FLiBe and FLiNaK

We report the results of constant-potential molecular dynamics simulations of the double layer interface between molten 2LiF–BeF 2 (FLiBe) and 23LiF–6NaF–21KF (FLiNaK) fluoride mixtures and idealized solid electrodes. Employing methods similar to those used in studies of chloride double layers, we compute the structure and differential capacitance of molten fluoride electric double layers as a function of applied voltage. The role of molten salt structure is probed through comparisons between FLiBe and FLiNaK, which serve as models for strong and weak associate-forming salts, respectively. In FLiBe, screening involves changes in Be–F–Be angles and alignment of the oligomers parallel to the electrode, while in FLiNaK, the electric field is screened mainly by rearrangement of individual ions, predominantly the polarizable potassium cation.

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

A Computational Study of RNA Tetraloop Thermodynamics, Including Misfolded States

An important characteristic of RNA folding is the adoption of alternative configurations of similar stability, often referred to as misfolded configurations. These configurations are considered to compete with correctly folded configurations, although their rigorous thermodynamic and structural characterization remains elusive. Tetraloop motifs found in large ribozymes are ideal systems for an atomistically detailed computational quantification of folding free energy landscapes and the structural characterization of their constituent free energy basins, including nonnative states. In this work, we studied a group of closely related 10-mer tetraloops using a combined parallel tempering and metadynamics technique that allows a reliable sampling of the free energy landscapes, requiring only knowledge that the stem folds into a canonical A-RNA configuration. Here we isolated and analyzed unfolded, folded, and misfolded populations that correspond to different free energy basins. We identified a distinct misfolded state that has a stability very close to that of the correctly folded state. This misfolded state contains a predominant population that shares the same structural features across all tetraloops studied here and lacks the noncanonical A-G base pair in its loop portion. Further analysis performed with biased trajectories showed that although this competitive misfolded state is not an essential intermediate, it is visited in most of the transitions from unfolded to correctly folded states. Moreover, the tetraloops can transition from this misfolded state to the correctly folded state without requiring extensive unfolding.

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

Computational determination of the pigment binding motif in the chlorosome protein a of green sulfur bacteria

We present a molecular-scale model of Bacteriochlorophyll a (BChl a) binding to the chlorosome protein A (CsmA) of Chlorobaculum tepidum, and the aggregated pigment–protein dimer, as determined from protein–ligand docking and quantum chemistry calculations. Our calculations provide strong evidence that the BChl a molecule is coordinated to the His25 residue of CsmA, with the magnesium center of the bacteriochlorin ring situated <3 Å from the imidazole nitrogen atom of the histidine sidechain, and the phytyl tail aligned along the nonpolar residues of the α-helix of CsmA. We also confirm that the Q y band in the absorption spectra of BChl a experiences a large (+16 to +43 nm) redshift when aggregated with another BChl a molecule in the CsmA dimer, compared to the BChl a in solvent; this redshift has been previously established by experimental researchers. We propose that our model of the BChl a–CsmA binding motif, where the dimer contains parallel aligned N-terminal regions, serves as the smallest repeating unit in a larger model of the para-crystalline chlorosome baseplate protein.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

The hydrogen tunneling splitting in malonaldehyde: A full-dimensional time-independent quantum mechanical method

The accurate time-independent quantum dynamics calculations on the ground-state tunneling splitting of malonaldehyde in full dimensionality are reported for the first time. This is achieved with an efficient method developed by us. In our method, the basis functions are customized for the hydrogen transfer process which has the effect of greatly reducing the size of the final Hamiltonian matrix, and the Lanczos method and parallel strategy are used to further overcome the memory and central processing unit time bottlenecks. The obtained ground-state tunneling splitting of 24.5 cm{sup −1} is in excellent agreement with the benchmark value of 23.8 cm{sup −1} computed with the full-dimensional, multi-configurational time-dependent Hartree approach on the same potential energy surface, and we estimate that our reported value has an uncertainty of less than 0.5 cm{sup −1}. Moreover, the role of various vibrational modes strongly coupled to the hydrogen transfer process is revealed.

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

Fine-Tuning the Activity of Metal–Organic Framework Supported Cobalt Catalysts for the Oxidative Dehydrogenation of Propane

Few-atom cobalt-oxide clusters, when dispersed on a Zr-based metal–organic framework (MOF) NU-1000, have previously been shown to be active for the oxidative dehydrogenation (ODH) of propane at low temperatures (< 230 °C), affording a selective and stable propene production catalyst. In our current work, a series of promoter ions with varying Lewis acidity, including Ni(II), Zn(II), Al(III), Ti(IV) and Mo(VI), are anchored as metal-oxide,hydroxide clusters to NU-1000 via SIM (solvothermal deposition within MOFs–specifically the nodes) followed by incorporation of Co(II) ions via vapor-phase AIM (atomic layer deposition (ALD) in MOFs). This process yields a series of NU-1000-supported bimetallic-oxo,hydroxo,aqua clusters. Using difference envelope density (DED) analyses, the spatial locations of the promoter ions and catalytic cobalt ions are determined. For all samples the SIM-anchored promoter ions are sited between pairs of Zr 6 nodes along the MOF c-axis (channel-aligned axis) whereas the location of the AIM-anchored cobalt ions varies depending on the identity of promoter metal ion. With Ni(II)-, Al(III)-, or Ti(IV)-containing clusters as promoters, the oxy-cobalt species are sited atop the promoter sites; with Mo(VI) they grow exclusively on the MOF nodes sites (hexa-Zr(IV)- oxo,hydroxo,aqua units); with Zn(II) they grow on both the node and promoter. The NU-1000- supported bimetallic-oxide clusters are active for propane ODH after thermal activation under O 2 to open a cobalt coordination site and to oxidize Co(II) to Co(III), as evidenced by operando Xray absorption spectroscopy at the Co K-edge. The cobalt component is exclusively responsible for the observed catalysis. In accord with the decreasing Lewis acidity of the promoter ion, catalytic activity increases in the order: Mo(VI)<Ti(IV)<Al(III)<Zn(II)<Ni(II). The finding is attributed to increasing ease of formation Co(III)-O• species and stabilization of a cobalt(III)- oxyl/propane transition state as the Lewis acidity of the promoter ions decreases. (The promoterfree catalyst is sited on the oxy-Zr(IV) node of the MOF and is inferior catalytically to the strongly promoted versions of the catalyst.) The relative Lewis acidities of the promoters are inferred from the initial pKa values of aqua complexes of the corresponding single metal ions. Departures from the (ODH-rate)/(promoter-pKa) correlation are paralleled by differences in cobalt catalyst siting or structure. The combined results point to an increasing ability to fine-tune, with close to atomic specificity, the structure-dependent activity of MOF-supported heterogeneous catalysts for gas phase chemical reactions. When coupled with mechanistic studies–computational or experimental–this ability may translate into informed prediction of improved catalysts for propane ODH and other chemical reactions.

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

Graph-based linear scaling electronic structure theory

We show how graph theory can be combined with quantum theory to calculate the electronic structure of large complex systems. The graph formalism is general and applicable to a broad range of electronic structure methods and materials, including challenging systems such as biomolecules. The methodology combines well-controlled accuracy, low computational cost, and natural low-communication parallelism. This combination addresses substantial shortcomings of linear scaling electronic structure theory, in particular with respect to quantum-based molecular dynamics simulations.

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

Collaborative simulations and experiments for a novel yield model of coal devolatilization in oxy-coal combustion conditions

Oxy-coal combustion is an emerging low-cost “clean coal” technology for emissions reduction and Carbon Capture and Sequestration (CCS). The use of Computational Fluid Dynamics (CFD) tools is crucial for the development of cost-effective oxy-fuel technologies and the minimization of environmental concerns at industrial scale. The coupling of detailed chemistry models and CFD simulations is still challenging, especially for large-scale plants, because of the high computational efforts required. The development of scale-bridging models is therefore necessary, to find a good compromise between computational efforts and the physical-chemical modeling precision. This paper presents a procedure for scale-bridging modeling of coal devolatilization, in the presence of experimental error, that puts emphasis on the thermodynamic aspect of devolatilization, namely the final volatile yield of coal, rather than kinetics. The procedure consists of an engineering approach based on dataset consistency and Bayesian methodology including Gaussian-Process Regression (GPR). Experimental data from devolatilization tests carried out in an oxy-coal entrained flow reactor were considered and CFD simulations of the reactor were performed. Jointly evaluating experiments and simulations, a novel yield model was validated against the data via consistency analysis. In parallel, a Gaussian-Process Regression was performed, to improve the understanding of the uncertainty associated to the devolatilization, based on the experimental measurements. Potential model forms that could predict yield during devolatilization were obtained. The set of model forms obtained via GPR includes the yield model that was proven to be consistent with the data. Finally, the overall procedure has resulted in a novel yield model for coal devolatilization and in a valuable evaluation of uncertainty in the data, in the model form, and in the model parameters.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

PRISMS-Fatigue computational framework for fatigue analysis in polycrystalline metals and alloys

Abstract The PRISMS-Fatigue open-source framework for simulation-based analysis of microstructural influences on fatigue resistance for polycrystalline metals and alloys is presented here. The framework uses the crystal plasticity finite element method as its microstructure analysis tool and provides a highly efficient, scalable, flexible, and easy-to-use ICME community platform. The PRISMS-Fatigue framework is linked to different open-source software to instantiate microstructures, compute the material response, and assess fatigue indicator parameters. The performance of PRISMS-Fatigue is benchmarked against a similar framework implemented using ABAQUS. Results indicate that the multilevel parallelism scheme of PRISMS-Fatigue is more efficient and scalable than ABAQUS for large-scale fatigue simulations. The performance and flexibility of this framework is demonstrated with various examples that assess the driving force for fatigue crack formation of microstructures with different crystallographic textures, grain morphologies, and grain numbers, and under different multiaxial strain states, strain magnitudes, and boundary conditions.

Chemistry↗
DOE OSTI2022

Unified Language Frontend for Physic-Informed AI/ML

Artificial intelligence and machine learning (AI/ML) are becoming important tools for scientific modeling and simulation as in several other fields such as image analysis and natural language processing. ML techniques can leverage the computing power available in modern systems and reduce the human effort needed to configure experiments, interpret and visualize results, draw conclusions from huge quantities of raw data, and build surrogates for physics based models. Domain scientists in fields like fluid dynamics, microelectronics and chemistry can automate many of their most difficult and repetitive tasks or improve the design times by use of the faster ML-surrogates. However, modern ML and traditional scientific highperformance computing (HPC) tend to use completely different software ecosystems. While ML frameworks like PyTorch and TensorFlow provide Python APIs, most HPC applications and libraries are written in C++. Direct interoperability between the two languages is possible but is tedious and error-prone. In this work, we show that a compiler-based approach can bridge the gap between ML frameworks and scientific software with less developer effort and better efficiency. We use the MLIR (multi-level intermediate representation) ecosystem to compile a pre-trained convolutional neural network (CNN) in PyTorch to freestanding C++ source code in the Kokkos programming model. Kokkos is a programming model widely used in HPC to write portable, shared-memory parallel code that can natively target a variety of CPU and GPU architectures. Our compiler-generated source code can be directly integrated into any Kokkosbased application with no dependencies on Python or cross-language interfaces.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Composite Behavior of Nanopore Array Large Memristors

Synthetic nanopores were recently demonstrated with memristive and nonlinear voltage-current behaviors, akin to ion channels in a cell membrane. Such ionic devices are considered a promising candidate for the development of brain-inspired neuromorphic computing techniques. In this work, we show the composite behavior of nanopore-array large memristors, formed with different membrane materials, pore sizes, electrolytes, and device arrangements. Anodic aluminum oxide (AAO) membranes with 5 nm and 20 nm diameter pores and track-etched polycarbonate (PCTE) membranes with 10 nm diameter pores are tested and shown to demonstrate memristive and nonlinear behaviors with approximately 107–1010 pores in parallel when electrolyte concentration across the membranes is asymmetric. Ion diffusion through the large number of channels induces time-dependent electrolyte asymmetry that drives the system through different memristive states. The behaviors of series composite memristors with different configurations are also presented. In addition to helping understand fluidic devices and circuits for neuromorphic computing, the results also shed light on the development of field-assisted ion-selection-membrane filtration techniques as well as the investigations of large neurons and giant synapses. Further work is needed to de-embed parasitic components of the measurement setup to obtain intrinsic large memristor properties.

Chemistry↗
DOE OSTI2022

Exploring the Use of Novel Spatial Accelerators in Scientific Applications

Driven by the need to find alternative accelerators which can viably replace GPUs in next-generation Supercomputing systems, this paper proposes a methodology to enable agile application/hardware co-design. The application-first methodology provides the ability to come up with design of accelerators while working with real-world workloads, available accelerators, and system software. The iterative design process targets a set of kernels in a workload for performance estimates that can prune the design space for later phases of detailed architectural evaluations. To this effect, in this paper, a novel data-parallel device model is introduced that simulates the latency of performance-sensitive operations in an accelerator including data transfers and kernel computation using multi-core CPUs. The use of off-the-shelf simulators, such as pre-RTL simulator Aladdin or multiple tools available for exploring the design of deep neural network accelerators (e.g., Timeloop) is demonstrated for evaluation of various accelerator designs using applications with realistic inputs. Examples of multiple device configurations that are instantiable in a system are explored to evaluate the performance benefit of deploying novel accelerators. The proposed device is integrated with a programming model and system software to potentially explore the impacts of high-level programming languages/compilers and low-level effects such as task scheduling on multiple accelerators. We analyze our methodology for a set of applications that represent high-performance computing (HPC) and graph analytics. The applications include a computational chemistry kernel realized using tensor contractions, triangle counting, GraphSAGE and Breadth-first Search. These applications include kernels such as dense matrix-dense matrix multiplication, sparse matrix-spare matrix multiplication, and sparse matrix-dense vector multiplication. Our results indicate potential performance benefits and insights for system design by including accelerators that realize these kernels along-side general purpose accelerators.

AI, codesign, Accelerated Computing, Modeling and ↗