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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Structural, electronic and optical properties of novel carbonate fluorides ABCO{sub 3}F (A=K, Rb, Cs; B=Ca, Sr)

In the present study, we performed first principles calculations on the electronic and optical properties of layered alkali–alkaline earth carbonate fluorides which attract attention in the domain of nonlinear optics. The calculated lattice parameters and volumes with and without inclusion of van der Waals (vdW) correction methods to standard density functional methods were compared with experiments. We observed that vdW interactions are predominant in RbCaCO{sub 3}F and CsCaCO{sub 3}F as compared with other computed compounds. The calculated bulk modulus from single crystal elastic constants reveals that these materials are all relatively harder than the KH{sub 2}PO{sub 4} (KDP) crystal. We also found that these materials are wide band gap insulators as obtained from Tran–Blaha modified Becke–Johnson potential. The linear optical properties such as dielectric function, refractive indices, birefringence and absorption spectra are presented. Finally, the calculated birefringence values indicate that these crystals could be promising for producing phase matching in the deep ultra-violet region. - Graphical abstract: The co- and anti-parallel alignment of CO{sub 3} groups leads to larger and smaller SHG coefficients in (a) KCaCO{sub 3}F, (b) KSrCO{sub 3}F, (c) RbSrCO{sub 3}F, (d) RbCaCO{sub 3}F, and (e) CsCaCO{sub 3}F. - Highlights: • Effect of van der Waals interactions on structural properties. • Single-crystalline elastic constants and mechanical stability. • Electronic properties. • Linear optical properties. • Comparison of present computed compounds with other well known materials in non-linear optical materials.

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

SCEPTRE 2.2 Quick Start Guide

This report provides a summary of notes for building and running the Sandia Computational Engine for Particle Transport for Radiation Effects (SCEPTRE) code. SCEPTRE is a general- purpose C++ code for solving the li near Boltzmann transport equation in serial or parallel using unstructured spatial finite elements, multigroup energy treatment, and a variety of angular treatments including discrete ordinates and spherical harmonics. Either the first-order form of the Boltzmann equation or one of the second-order forms may be solved. SCEPTRE requires a small number of open-source Third Part y Libraries (TPL) to be available, and example scripts for building these TPLs are provided. The TPLs needed by SCEPTRE are Trilinos, boost, and netcdf. SCEPTRE uses an autotools build system , and a sample configure script is provided. Running the SCEPTRE code requires that the user provide a spatial finite-elements mesh in Exodus format and a cross section library in a format that will be described. SCEPTRE uses an xml-based input, and several examples will be provided.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2017

Probing the Energetics of Molecule–Material Interactions at Interfaces and in Nanopores

During the past decades, advances in interfacial chemistries at the molecular level are shaping our world by playing crucial roles in balancing global scale energy crisis and critical environmental concerns. However, systematic investigations into the binding energies, site distribution and their correlation with the molecular-level surface assemblages and structures at interfaces and in nanopores are rarely documented. In this review, we summarize a set of systematic calorimetric studies on surface energetics we performed during the last decade. These studies demonstrate how thermochemistry can reveal crucial energetic insights into a series of molecule – material interactions relevant to a number of applications, including carbon capture and sequestration, energy production, sustainable chemical processing, catalysis, and nanogeoscience. Calorimetric methodologies developed and applied include direct gas adsorption calorimetry, nearroom temperature solvent immersion/solution calorimetry and high temperature oxide melt solution calorimetry. Using these highly unique techniques, we reveal the thermodynamic complexity of carbon dioxide capture on metal – organic framework (MOF) sorbents with built-in and grafted nucleophilic functional groups (-OH and -NH2). These studies suggest that carbon dioxide adsorption on functionalized MOFs is a complex process involving multiple thermodynamic factors, as reflected by changes in surface phase and structure, chemical bonding and degree of disorder with varying temperature and gas loading. The fundamental insights obtained may help optimize the design, synthesis and application of MOF-based carbon dioxide sorbents for carbon capture and sequestration. In parallel, we also explore the energetics of interaction and competition between small molecules (water, carbon dioxide, methane, simple and complex organics) and inorganic materials (calcite, silica, zirconia, zeolites, mesoporous frameworks, alumina, and uranium), at interfaces and in nanopores. Combined with spectroscopic, diffraction, electron microscopic and computational techniques, the energetics of gas/liquid – solid interactions can be correlated with specific bonds, molecular configurations and nanostructures. Although the energetics evolves continuously from weak association to strong bonding to classical capping, distinct regions of rapidly changing stepwise energetics often separate the different regimes. These phenomena are closely related to the properties of inorganic material surfaces (hydrophobicity and acidity/basicity), the framework architectures, and the chemical nature of adsorbate molecules. These direct thermodynamic insights reinforce our understanding of complex small molecule – inorganic material interactions important to multiple disciplines of chemical engineering, materials science, nanogeoscience and environmental technology, including heterogeneous catalysis, molecular separation, material design and synthesis, biomineralization, contaminant and nutrient transport, carbonate formation, and water – organic competitions on material/mineral surfaces.

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

Enabling Modular Autonomous Feedback‐Loops in Materials Science through Hierarchical Experimental Laboratory Automation and Orchestration

Abstract Materials acceleration platforms (MAPs) operate on the paradigm of integrating combinatorial synthesis, high‐throughput characterization, automatic analysis, and machine learning. Within a MAP, one or multiple autonomous feedback loops may aim to optimize materials for certain functional properties or to generate new insights. The scope of a given experiment campaign is defined by the range of experiment and analysis actions that are integrated into the experiment framework. Herein, the authors present a method for integrating many actions within a hierarchical experimental laboratory automation and orchestration (HELAO) framework. They demonstrate the capability of orchestrating distributed research instruments that can incorporate data from experiments, simulations, and databases. HELAO interfaces laboratory hardware and software distributed across several computers and operating systems for executing experiments, data analysis, provenance tracking, and autonomous planning. Parallelization is an effective approach for accelerating knowledge generation provided that multiple instruments can be effectively coordinated, which the authors demonstrate with parallel electrochemistry experiments orchestrated by HELAO. Efficient implementation of autonomous research strategies requires device sharing, asynchronous multithreading, and full integration of data management in experimental orchestration, which to the best of the authors’ knowledge, is demonstrated for the first time herein.

36 MATERIALS SCIENCE↗
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↗
NASA NTRS1997

Computation of Reacting Flows in Combustion Processes

The main objective of this research was to develop an efficient three-dimensional computer code for chemically reacting flows. The main computer code developed is ALLSPD-3D. The ALLSPD-3D computer program is developed for the calculation of three-dimensional, chemically reacting flows with sprays. The ALL-SPD code employs a coupled, strongly implicit solution procedure for turbulent spray combustion flows. A stochastic droplet model and an efficient method for treatment of the spray source terms in the gas-phase equations are used to calculate the evaporating liquid sprays. The chemistry treatment in the code is general enough that an arbitrary number of reaction and species can be defined by the users. Also, it is written in generalized curvilinear coordinates with both multi-block and flexible internal blockage capabilities to handle complex geometries. In addition, for general industrial combustion applications, the code provides both dilution and transpiration cooling capabilities. The ALLSPD algorithm, which employs the preconditioning and eigenvalue rescaling techniques, is capable of providing efficient solution for flows with a wide range of Mach numbers. Although written for three-dimensional flows in general, the code can be used for two-dimensional and axisymmetric flow computations as well. The code is written in such a way that it can be run in various computer platforms (supercomputers, workstations and parallel processors) and the GUI (Graphical User Interface) should provide a user-friendly tool in setting up and running the code.

Keith, Theo G., Jr.↗
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↗
NASA NTRS1984

Interactive computer modeling of combustion chemistry and coalescence-dispersion modeling of turbulent combustion

An interactive computer code for simulation of a high-intensity turbulent combustor as a single point inhomogeneous stirred reactor was developed from an existing batch processing computer code CDPSR. The interactive CDPSR code was used as a guide for interpretation and direction of DOE-sponsored companion experiments utilizing Xenon tracer with optical laser diagnostic techniques to experimentally determine the appropriate mixing frequency, and for validation of CDPSR as a mixing-chemistry model for a laboratory jet-stirred reactor. The coalescence-dispersion model for finite rate mixing was incorporated into an existing interactive code AVCO-MARK I, to enable simulation of a combustor as a modular array of stirred flow and plug flow elements, each having a prescribed finite mixing frequency, or axial distribution of mixing frequency, as appropriate. Further increase the speed and reliability of the batch kinetics integrator code CREKID was increased by rewriting in vectorized form for execution on a vector or parallel processor, and by incorporating numerical techniques which enhance execution speed by permitting specification of a very low accuracy tolerance.

Pratt, D. T.↗
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↗
NASA NTRS1987

Molecular dynamics of the water liquid-vapor interface

The results of molecular dynamics calculations on the equilibrium interface between liquid water and its vapor at 325 K are presented. For the TIP4P model of water intermolecular pair potentials, the average surface dipole density points from the vapor to the liquid. The most common orientations of water molecules have the C2 nu molecular axis roughly parallel to the interface. The distributions are quite broad and therefore compatible with the intermolecular correlations characteristic of bulk liquid water. All near-neighbor pairs in the outermost interfacial layers are hydrogen bonded according to the common definition adopted here. The orientational preferences of water molecules near a free surface differ from those near rigidly planar walls which can be interpreted in terms of patterns found in hexagonal ice 1. The mean electric field in the interfacial region is parallel to the mean polarization which indicates that attention cannot be limited to dipolar charge distributions in macroscopic descriptions of the electrical properties of this interface. The value of the surface tension obtained is 132 +/- 46 dyn/cm, significantly different from the value for experimental water of 68 dyn/cm at 325 K.

NASA Center ARC↗
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↗