Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Kernel”

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 1,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Machine Learning-Guided Approach for Studying Solvation Environments

Molecular-level understanding and characterization of solvation environments are often needed across chemistry, biology, and engineering. Toward practical modeling of local solvation effects of any solute in any solvent, in this work, we report a static and all-quantum mechanics-based cluster-continuum approach for calculating single-ion solvation free energies. This approach uses a global optimization procedure to identify low-energy molecular clusters with different numbers of explicit solvent molecules and then employs the smooth overlap for atomic positions learning kernel to quantify the similarity between different low-energy solute environments. From these data, we use sketch maps, a nonlinear dimensionality reduction algorithm, to obtain a two-dimensional visual representation of the similarity between solute environments in differently sized microsolvated clusters. After testing this approach on different ions having charges 2+, 1+, 1–, and 2–, we find that the solvation environment around each ion can be seen to usually become more similar in hand with its calculated single-ion solvation free energy. Without needing either dynamics simulations or an a priori knowledge of local solvation structure of the ions, this approach can be used to calculate solvation free energies within 5% of experimental measurements for most cases, and it should be transferable for the study of other systems where dynamics simulations are not easily carried out.

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

Liquid–Vapor Phase Equilibrium in Molten Aluminum Chloride (AlCl 3 ) Enabled by Machine Learning Interatomic Potentials

Molten salts are promising candidates in numerous clean energy applications, where knowledge of thermophysical properties and vapor pressure across their operating temperature ranges is critical for safe operations. Due to challenges in evaluating these properties using experimental methods, fast and scalable molecular simulations are essential to complement the experimental data. In this study, we developed machine learning interatomic potentials (MLIP) to study the AlCl 3 molten salt across varied thermodynamic conditions (T = 473–613 K and P = 2.7–23.4 bar), which allowed us to predict temperature-surface tension correlations and liquid–vapor phase diagram from direct simulations of two-phase coexistence in this molten salt. Two MLIP architectures, a Kernel-based potential and neural network interatomic potential (NNIP), were considered to benchmark their performance for AlCl 3 molten salt using experimental structure and density values. The NNIP potential employed in two-phase equilibrium simulations yields the critical temperature and critical density of AlCl 3 that are within 10 K (∼3%) and 0.03 g/cm 3 (∼7%) of the reported experimental values. An accurate correlation between temperature and viscosities is obtained as well. In doing so, we report that the inclusion of low-density configurations in their training is critical to more accurately represent the AlCl 3 system across a wide phase-space. The MLIP trained using PBE-D3 functional in the ab initio molecular dynamics (AIMD) simulations (120 atoms) also showed close agreement with experimentally determined molten salt structure comprising Al 2 Cl 6 dimers, as validated using Raman spectra and neutron structure factor. Furthermore, the PBE-D3 as well as its trained MLIP showed better liquid density and temperature correlation for AlCl 3 system when compared to several other density functionals explored in this work. Overall, the demonstrated approach to predict temperature correlations for liquid and vapor densities in this study can be employed to screen nuclear reactors-relevant compositions, helping to mitigate safety concerns.

Ab initio molecular dynamics↗
DOE OSTI2025

Physics-Informed Gaussian Process Inference of Liquid Structure from Scattering Data

We present a nonparametric Bayesian framework to infer radial distribution functions from experimental scattering measurements with uncertainty quantification using nonstationary Gaussian processes. The Gaussian process prior mean and kernel functions are designed to mitigate well-known numerical challenges with the Fourier transform, including discrete measurement binning and detector windowing, while encoding fundamental yet minimal physical knowledge of the liquid structure. We demonstrate uncertainty propagation of the Gaussian process posterior to unmeasured quantities of interest. Experimental radial distribution functions of liquid argon and water with uncertainty quantification are provided as both a proof of principle for the method and a benchmark for molecular models.

Chemical structure↗
DOE OSTI2016

Energetic and dynamic analysis of transport of Na + and K + through a cyclic peptide nanotube in water and in lipid bilayers

Potential of mean force (PMF) profiles and position-dependent diffusion coefficients of Na + and K + are calculated to elucidate the translocation of ions through a cyclic peptide nanotube, composed of 8 × cyclo[-(D-Leu-Trp) 4- ] rings, in water and in hydrated DMPC bilayers. The PMF profiles and PMF decomposition analysis for the monovalent cations show that favorable interactions of the cations with the CPN as well as the lipid bilayer and dehydration free energy penalties are two major competing factors which determine the free energy surface for ion transport through CPNs both in water and lipid bilayers, and that the selectivity of CPNs to cations mainly arises from favorable interaction energies of cations with CPNs and lipid bilayers that are more dominant than the dehydration penalties. Calculations of the position-dependent diffusion coefficients and dynamic friction kernels of the cations indicate that the dehydration process along with the molecular rearrangements occurring outside the channel and the coupling of the ion motions with the chain-structured water movements inside the channel lead to decrease of the diffusion coefficients far away from the channel entrance and also reduced coefficients inside the channel. Here the PMF and diffusivity profiles for Na + and K + reveal that the energetics of ion transport through the CPN are governed by global interactions of ions with all the components in the system while the diffusivity of ions through the channel is mostly determined by local interactions of ions with the confined water molecules inside the channel. Comparison of Na + and K + ion distributions based on overdamped Brownian dynamics simulations based on the PMF and diffusivity profiles with the corresponding results from molecular dynamics shows good agreement, indicating accuracy of the Bayesian inference method for determining diffusion coefficients in this application. In addition this work shows that position-dependent diffusion coefficients of ions are required to explain the dynamics and conductance of ions through the CPN properly.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Comparative Study on the Machine Learning-Based Prediction of Adsorption Energies for Ring and Chain Species on Metal Catalyst Surfaces

Computation of adsorption and transition state energies for a large number of surface intermediates for numerous active site models pose significant computational overhead in computational screening of catalysts. Machine learning (ML) techniques can be used to predict part of these energies. To predict the energies, ML models need to be fed appropriate metal and species descriptors. For complex surface chemistries, the structures of the intermediate species can vary greatly. In this paper, working with the hydrodeoxygenation of succinic acid on six different metal surfaces, we have studied the effect of linear and non-linear ML models used along with pen-and-paper based species descriptors and two categories of metal descriptors on two different categories of intermediate species: chain and ring. More specifically, our computations include the prediction of chain species when trained on only chain species and also when trained on both chain and ring species. Similar computations were performed for predictions of ring species. In each case, results of linear ML models were compared with kernel based non-linear models. Our results indicate that ring species data does not improve the prediction of chain species. Similarly, chain species data does not improve the prediction of ring species. The use of non-linear ML models, however, did help to minimize the prediction errors compared to the linear models. Furthermore, the study also shows that electronic or adsorption energy based metal descriptors along with bond count based species fingerprints can achieve a mean absolute error (MAE) of less than 0.2 eV for complex chain molecules when used with an appropriate machine learning model.

Adsorption↗
DOE OSTI2018

Prediction of Adsorption Energies for Chemical Species on Metal Catalyst Surfaces Using Machine Learning

Computational catalyst screening has the potential to significantly accelerate heterogeneous catalyst discovery. Typically, this involves developing microkinetic reactor models that are based on parameters obtained from density functional theory and transition-state theory. To reduce the large computational cost involved in computing various adsorption and transition-state energies of all possible surface states on a large number of catalyst models, linear scaling relations for surface intermediates and transition states have been developed that only depend on a few, typically one or two descriptors, such as the carbon atom adsorption energy. As a result, only the descriptor values have to be computed for various active site models to generate volcano curves in activity or selectivity. Unfortunately, for more complex chemistries the predictability of linear scaling relations is unknown. Also, the selection of descriptors is essentially a trial and error process. Here, using a database of adsorption energies of the surface species involved in the decarboxylation and decarbonylation of propionic acid over eight monometalic transition-metal catalyst surfaces (Ni, Pt, Pd, Ru, Rh, Re, Cu, Ag), we tested if nonlinear machine learning (ML) models can outperform the linear scaling relations in prediction accuracy when predicting the adsorption energy for various species on a metal surface based on data from the rest of the metal surfaces. We found linear scaling relations to hold well for predictions across metals with a mean-absolute error of 0.12 eV, and ML methods being unable to outperform linear scaling relations when the training dataset contains a complete set of energies for all of the species on various metal surfaces. Only when the training dataset is incomplete, namely, contains a random subset of species’ energies for each metal, a currently unlikely scenario for catalyst screening, do kernel-based ML models significantly outperform linear scaling relations. We also found that simple coordinate-free species descriptors, such as bond counts, achieve as good results as sophisticated coordinate-based descriptors. Finally, we propose an approach for automatic discovery of appropriate metal descriptors using principal component analysis.

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

Curing the Divergence in Time-Dependent Density Functional Quadratic Response Theory

While time-dependent density functional theory has emerged as a method of choice for computing electronic spectra and response of molecules and materials, its reliability is hindered by the adiabatic approximation that is commonly made. In this work, we address one problematic aspect that arises from this approximation: an incorrect pole structure in the quadratic response function, leading to unphysical divergences in excited state-to-state transition probabilities and hyperpolarizabilties. We find the form of the exact quadratic response kernel and derive a practical and accurate approximation that cures the divergence. Here, we demonstrate our results on excited state-to-state transition probabilities of a model system and of the LiH molecule.

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

Electron-Beam-Induced Synthesis of Hexagonal 1 H -MoSe 2 from Square β-FeSe Decorated with Mo Adatoms

Two-dimensional (2D) materials have generated interest in the scientific community because of the advanced electronic applications they might offer. Powerful electron beam microscopes have been used not only to evaluate the structures of these materials but also to manipulate them by forming vacancies, nanofragments, and nanowires or joining nanoislands together. In this work, we show that the electron beam in a scanning transmission electron microscope (STEM) can be used in yet another way: to mediate the synthesis of 2D 1H-MoSe 2 from Mo-decorated 2D β-FeSe and simultaneously image the process on the atomic scale. This is quite remarkable given the different crystal structures of the reactant (square lattice β-FeSe) and the product (hexagonal lattice 1H-MoSe 2 ). Herein, the feasibility of the transformation was first explored by theoretical calculations that predicted that the reaction is exothermic. Furthermore, a theoretical reaction path to forming a stable 1H-MoSe 2 nucleation kernel within pure β-FeSe was found, demonstrating that the pertinent energy barriers are smaller than the energy supplied by the STEM electron beam.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2019

Machine Learning the Voltage of Electrode Materials in Metal-Ion Batteries

Machine-learning (ML) techniques have rapidly found applications in many domains of materials chemistry and physics where large data sets are available. Aiming to accelerate the discovery of materials for battery applications, in this work, we develop a tool based on ML models to predict voltages of electrode materials for metal-ion batteries. To this end, we use deep neural network, support vector machine, and kernel ridge regression as ML algorithms in combination with data taken from the Materials Project database, as well as feature vectors from properties of chemical compounds and elemental properties of their constituents. We show that our ML models have predictive capabilities for different reference test sets and, as an example, we utilize them to generate a voltage profile diagram and compare it to density functional theory calculations. In addition, using our models, we propose nearly 5000 candidate electrode materials for Na- and K-ion batteries. We also make available a web-accessible tool that, within a minute, can be used to estimate the voltage of any bulk electrode material for a number of metal ions. Here, these results show that ML is a promising alternative for computationally demanding calculations as a first screening tool of novel materials for battery applications.

25 ENERGY STORAGE↗
DOE OSTI2019

Machine Learning Accelerates the Discovery of Design Rules and Exceptions in Stable Metal–Oxo Intermediate Formation

Metal–oxo moieties are important catalytic intermediates in the selective partial oxidation of hydrocarbons and in water splitting. Stable metal–oxo species have reactive properties that vary depending on the spin state of the metal, complicating the development of structure–property relationships. To overcome these challenges, we train machine-learning (ML) models capable of predicting metal–oxo formation energies across a range of first-row metals, oxidation states, and spin states. Using connectivity-only features tailored for inorganic chemistry as inputs to kernel ridge regression or artificial neural network (ANN) ML models, we achieve good mean absolute errors (4–5 kcal/mol) on set-aside test data across a range of ligand orientations. Analysis of feature importance for oxo formation energy prediction reveals the dominance of nonlocal, electronic ligand properties in contrast to other transition metal complex properties (e.g., spin-state or ionization potential). We enumerate the theoretical catalyst space with an ANN, revealing expected trends in oxo formation energetics, such as destabilization of the metal–oxo species with increasing d-filling, as well as exceptions, such as weak correlations with indicators of oxidative stability of the metal in the resting state or unexpected spin-state dependence in reactivity. We carry out uncertainty-aware evolutionary optimization using the ANN to explore a >37 000 candidate catalyst space. New metal and oxidation state combinations are uncovered and validated with density functional theory (DFT), including counterintuitive oxo formation energies for oxidatively stable complexes. This approach doubles the density of confirmed DFT leads in originally sparsely populated regions of property space, highlighting the potential of ML-model-driven discovery to uncover catalyst design rules and exceptions.

transition metal catalysis↗
DOE OSTI2026

Accelerating Structure–Property Relationship Discovery with Multimodal Machine Learning and Self-Driving Microscopy

Microscopy combined with local spectroscopy is widely used to correlate nanoscale structure with functional properties in materials, but conventional measurements rely heavily on human-selected sampling locations and predefined targets, limiting data set diversity and the potential for discovery. Here, we present a framework that integrates autonomous microscopy with dual-novelty deep kernel learning (DN-DKL) for adaptive data acquisition and a dual variational autoencoder (VAE) for representation learning. DN-DKL actively guides the microscopy toward structurally and spectroscopically novel regions, enabling efficient collection of large spectral data sets. Dual-VAE embeds local structures and spectroscopic responses into a shared latent manifold that serves as a structure–property relationship map. We applied this framework for the investigation of halide perovskite films by using conductive atomic force microscopy. The results reveal distinct hysteresis behaviors that are linked to specific nanoscale structural motifs, including grain boundary junction points that show hysteresis under different bias conditions and asymmetric grain boundaries that suppress the charge transport. This framework establishes a general strategy that leverages the complementary strengths of self-driving microscopy, machine learning, and human expertise to accelerate scientific discovery in functional materials.

atomic force microscopy↗
DOE OSTI2020

Lignin Monomers from beyond the Canonical Monolignol Biosynthetic Pathway: Another Brick in the Wall

Lignin is conventionally defined as being formed by the oxidative polymerization of three main monolignols, p-coumaryl alcohol, coniferyl alcohol, and sinapyl alcohol, that are derived from the general phenylpropanoid biosynthetic pathway. Many other phenolic compounds that are also derived from the phenylpropanoid pathway are also known to perform as genuine lignin monomers in many plants, as is the case of the monolignol ester conjugates, phenolic compounds arising from the truncated biosynthesis of monolignols, or ferulate esters. Recent investigations, however, have indicated that phenolic compounds arising from beyond the canonical phenylpropanoid pathway, namely flavonoids, hydroxystilbenes, and hydroxycinnamic amides, may also behave as authentic lignin monomers and are incorporated into the lignin in some plants, further challenging the traditional definition of lignin. This is the case of the flavone tricin that is incorporated into the lignin of grasses and other monocots, the hydroxystilbene piceatannol (together with resveratrol and isorhapontigenin, at lower levels) that has been found in the lignins of palm fruit shells, their respective O-glucosides (astringin, piceid, and isorhapontin) that are present in the lignin of Norway spruce bark, or the ferulic amides feruloyltyramine, incorporated into the lignin of tobacco and potato tubers, and diferuloylputrescine, which appears to be incorporated into maize kernel lignin. These valuable compounds are potentially available in high amounts and at low cost and may be obtained from the waste products from the processing of agricultural or forest biomass.

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

Hydride-Containing Eight-Electron Pt/Ag Superatoms: Structure, Bonding, and Multi-NMR Studies

Recent reports on hydride-doped noble metal nanoclusters strongly suggest that the encapsulated hydride is a part of the superatom core, but no accurate location of the hydride could be experimentally proved, so far. Here, we report herein a hydride-doped eight-electron platinum/silver alloy nanocluster in which the position of four-coordinated hydride was determined by neutron diffraction for the first time. X-ray structures of [PtHAg 19 (dtp/desp) 12 ] (dtp = S 2 P(O n Pr) 2 , 1; dsep = Se 2 P(O i Pr) 2 , 2) describe a central platinum hydride (PtH) unit encapsulated within a distorted Ag 12 icosahedron, the resulting (PtH)@Ag 12 core being stabilized by an outer sphere made up of 7 capping silver atoms and 12 dichalcogenolates. Solid-state structures of 1 and 2 differ somewhat in the spatial configuration of their outer spheres, resulting in overall different symmetries, C 1 and C 3 , respectively. Whereas the multi-NMR spectra of 2 in solution at 173 K reveal that the structure of C 3 symmetry is the predominant one, 1 H and 195 Pt NMR spectra of 1 at the same temperature disclose the presence of isomers of both C 1 and C 3 symmetry. DFT calculations found both isomers to be very close in energy, supporting the fact that they co-exist in solution. They also show that the [PtH@Ag 12 ] 5+ kernel can be viewed as a closed-shell superatomic core, the μ4-hydride electron contributing to its eight-electron count. On the other hand, the 1s(H) orbital contributes only moderately to the superatomic orbitals, being mainly involved in the building of a Pt–H bonding electron pair with the 5d z 2 (Pt) orbital.

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

Overdestabilization vs Overstabilization in the Theoretical Analysis of f-Orbital Covalency

The complex nature of the f-orbital electronic structures and their interaction with the chemical environment pose significant computational challenges. Advanced computational techniques that variationally include scalar relativities and spin–orbit coupling directly at the molecular orbital level have been developed to address this complexity. Among these, variational relativistic multiconfigurational multireference methods stand out for their high accuracy and systematic improvement in studies of f-block complexes. Additionally, these advanced methods offer the potential for calibrating low-scaling electronic structure methods such as density functional theory. However, studies on the Cl K-edge X-ray absorption spectra of the [Ce(III)Cl 6 ] 3– and [Ce(IV)Cl 6 ] 2– complexes show that time-dependent density functional theory with approximate exchange–correlation kernels can lead to inaccuracies, resulting in an overstabilization of 4f orbitals and incorrect assessments of covalency. In contrast, approaches utilizing small active space wave function methods may understate the stability of these orbitals. The results herein demonstrate the need for large active space, multireference, and variational relativistic methods in studying f-block complexes.

chemical calculations↗
DOE OSTI2024

Using Data Science Tools to Reveal and Understand Subtle Relationships of Inhibitor Structure in Frontal Ring-Opening Metathesis Polymerization

The rate of frontal ring-opening metathesis polymerization (FROMP) using the Grubbs generation II catalyst is impacted by both the concentration and choice of monomers and inhibitors, usually organophosphorus derivatives. Herein we report a data-science-driven workflow to evaluate how these factors impact both the rate of FROMP and how long the formulation of the mixture is stable (pot life). Using this workflow, we built a classification model using a single-node decision tree to determine how a simple phosphine structural descriptor (V bur-near ) can bin long versus short pot life. Additionally, we applied a nonlinear kernel ridge regression model to predict how the inhibitor and selection/concentration of comonomers impact the FROMP rate. Furthermore, the analysis provides selection criteria for material network structures that span from highly cross-linked thermosets to non-cross-linked thermoplastics as well as degradable and nondegradable materials.

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

Real Time Monitoring of the Dynamic Intracluster Diffusion of Single Gold Atoms into Silver Nanoclusters

Alloying metal materials with heterometal atoms is an efficient way to diversify the function of materials, but in-depth understanding of the dynamic heterometallic diffusion inside the alloying materials is rather limited, especially at the atomic level. Here, we report the real-time monitoring of the dynamic diffusion process of a single gold (Au) atom into an atomically precise silver nanocluster (Ag NC), Ag 25 (MHA) 18 (MHA = 6-mercaptohexanoic acid), by using in situ UV–vis absorption spectroscopy in combination with mass and tandem mass spectrometry. We found that the Au heteroatom first replaces the Ag atom at the surface Ag 2 (MHA) 3 motifs of Ag 25 (MHA) 18 . After that, the Au atom diffuses into the surface layer of the icosahedral Ag 13 kernel and finally occupies the center of the alloy NCs to form the thermodynamically stable [email protected] 24 (MHA) 18 product. Density functional theory (DFT) calculations reveal that the key thermodynamic driving force is the preference of the Au heteroatom for the central site of alloy NCs. Furthermore, the real-time monitoring approach developed in this study could also be extended to other metal alloy systems to reveal the reaction dynamics of intracluster diffusion of heteroatoms, as well as the formation mechanisms of metal alloy nanomaterials.

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

Evaluating cloud feedbacks and rapid responses in the ACCESS model

A cloud feedback diagnostic package is implemented in the Australian Community Climate and Earth-System Simulator General Circulation Model, based on the methodology of “cloud radiative kernels.” Using separate increased sea surface temperature and CO 2 experiments, both the “rapid response” cloud contribution to forcing and temperature-mediated cloud feedbacks are analyzed. Under increased temperature and CO 2 changes, temperature-mediated cloud radiative feedback dominates over the rapid response in the final radiative response. Cloud feedback is positive in both long and short wave, with short wave dominating global values. Contributing most to this are low to mid-level clouds, of medium-to-high optical thickness. As a means of illustration of the methodology, a number of key parameters related to clouds, precipitation, and convection that are typically used in “tuning” in the model are modified. These changes result in substantial impacts on the model's current climate, but only modest changes to rapid response and feedbacks occur globally, regionally, and as a function of cloud optical thickness and height. This limited set of experiments shows that cloud adjustments and feedbacks in this model are robust under these changes, lending confidence that both model climate change projections and the conclusions of attribution studies are not overly sensitive to such parameterization tuning. Of course, a considerably larger set of experiments would be needed to demonstrate that feedbacks and rapid response are robust under the wider set of tuning adjustments commonly undertaken.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Computational Benefit of GPU Optimization for the Atmospheric Chemistry Modeling

Abstract Global chemistry‐climate models are computationally burdened as the chemical mechanisms become more complex and realistic. Optimization for graphics processing units (GPU) may make longer global simulation with regional detail possible, but limited study has been done to explore the potential benefit for the atmospheric chemistry modeling. Hence, in this study, the second‐order Rosenbrock solver of the chemistry module of CAM4‐Chem is ported to the GPU to gauge potential speed‐up. We find that on the CPU, the fastest performance is achieved using the Intel compiler with a block interleaved memory layout. Different combinations of compiler and memory layout lead to ~11.02× difference in the computational time. In contrast, the GPU version performs the best when using a combination of fully interleaved memory layout with block size equal to the warp size, CUDA streams for independent kernels, and constant memory. Moreover, the most efficient data transfer between CPU and GPU is gained by allocating the memory contiguously during the data initialization on the GPU. Compared to one CPU core, the speed‐up of using one GPU alone reaches a factor of ~11.7× for the computation alone and ~3.82× when the data transfer between CPU and GPU is considered. Using one GPU alone is also generally faster than the multithreaded implementation for 16 CPU cores in a compute node and the single‐source solution (OpenACC). The best performance is achieved by the implementation of the hybrid CPU/GPU version, but rescheduling the workload among the CPU cores is required before the practical CAM4‐Chem simulation.

54 ENVIRONMENTAL SCIENCES↗