Search NASA⌕ Search

SEARCH · Search NASA

Results for “atomic & molecular clusters”

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 OSTI2022

Homogeneous ice nucleation in an ab initio machine-learning model of water

Molecular simulations have provided valuable insight into the microscopic mechanisms underlying homogeneous ice nucleation. While empirical models have been used extensively to study this phenomenon, simulations based on first-principles calculations have so far proven prohibitively expensive. Here, we circumvent this difficulty by using an efficient machine-learning model trained on density-functional theory energies and forces. We compute nucleation rates at atmospheric pressure, over a broad range of supercoolings, using the seeding technique and systems of up to hundreds of thousands of atoms simulated with ab initio accuracy. The key quantity provided by the seeding technique is the size of the critical cluster (i.e., a size such that the cluster has equal probabilities of growing or melting at the given supersaturation), which is used together with the equations of classical nucleation theory to compute nucleation rates. We find that nucleation rates for our model at moderate supercoolings are in good agreement with experimental measurements within the error of our calculation. We also study the impact of properties such as the thermodynamic driving force, interfacial free energy, and stacking disorder on the calculated rates.

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

Trends in Siting of Metals in Heterometallic Nd–Yb Metal–Organic Frameworks and Molecular Crystals

Several studies suggest that metal ordering within metal–organic frameworks (MOFs) is important for understanding how MOFs behave in relevant applications; however, these siting trends can be difficult to determine experimentally. To garner insight into the energetic driving forces that may lead to nonrandom ordering within heterometallic MOFs, we employ density functional theory (DFT) calculations on several bimetallic metal–organic crystals composed of Nd and Yb metal atoms. We also investigate the metal siting trends for a newly synthesized MOF. Our DFT-based energy of mixing results suggest that Nd will likely occupy sites with greater access to electronegative atoms and that local homometallic domains within a mixed-metal Nd–Yb system are favored. We also explore the use of less computationally extensive methods such as classical force fields and cluster expansion models to understand their feasibility for large system sizes. Here, this study highlights the impact of metal ordering on the energetic stability of heterometallic MOFs and crystal structures.

energy of mixing↗
DOE OSTI2023

Impact of environmental oxygen on nanoparticle formation and agglomeration in aluminum laser ablation plumes

Here, the role of ambient oxygen gas (O 2 ) on molecular and nanoparticle formation and agglomeration was studied in laser ablation plumes. As a lab-scale surrogate to a high explosion detonation event, nanosecond laser ablation of an aluminum alloy (AA6061) target was performed in atmospheric pressure conditions. Optical emission spectroscopy and two mass spectrometry techniques were used to monitor the early to late stages of plasma generation to track the evolution of atoms, molecules, clusters, nanoparticles, and agglomerates. The experiments were performed under atmospheric pressure air, atmospheric pressure nitrogen, and 20% and 5% O 2 (balance N 2 ), the latter specifically with in situ mass spectrometry. Electron microscopy was performed ex situ to identify crystal structure and elemental distributions in individual nanoparticles. We find that the presence of ≈20% O 2 leads to strong AlO emission, whereas in a flowing N 2 environment (with trace O 2 ), AlN and strong, unreacted Al emissions are present. In situ mass spectrometry reveals that as O 2 availability increases, Al oxide cluster size increases. Nanoparticle agglomerates formed in air are found to be larger than those formed under N 2 gas. High-resolution transmission electron microscopy demonstrates that Al 2 O 3 and AlN nanoparticle agglomerates are formed in both environments; indicating that the presence of trace O 2 can lead to Al 2 O 3 nanoparticle formation. The present results highlight that the availability of O 2 in the ambient gas significantly impacts spectral signatures, cluster size, and nanoparticle agglomeration behavior. These results are relevant to understanding debris formation in an explosion event, and interpreting data from forensic investigations.

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

Hydrogen Bonding Inside Anionic Polymeric Brush Layer: Machine Learning-Driven Exploration of the Relative Roles of the Polymer Steric Effect, Charging, and Type of Screening Counterions

This paper employs a combination of all-atom molecular dynamics (MD) simulations and unsupervised machine learning (ML) for studying the water-water hydrogen bonds (HBs) inside the anionic poly-acrylic acid (PAA) brushes modeled using all-atom MD simulations. PAA brush layer with different charge fraction (f), namely f=0, f=0.25, and f=1, is considered. Water-water interactions, both inside and outside the brush layer, are represented through distinct clusters of tupules of variables representing distances associated with the interacting water molecules. While clusters representing the HBs are present for water inside and outside the brushes, several clusters representing the long-range water-water interactions are missing for the water molecules inside the highly charged (f=1) PAA brushes. More importantly, inside highly charged brushes, the edge of the clusters representing the water-water HBs is progressively shortened, as compared to that in the bulk. Both these results stem from the presence of the PAA brushes imparting the steric effect and the charge effect, or the effect associated with enhanced interactions of water molecules with PE charges and counterions, thereby disrupting the water connectivity. This water-charged-species interaction also increases the water-water HB angle, i.e., makes the water-water HBs less stable inside the highly charged PAA brush layer. The narrowing of the clusters representing the HBs and the alteration of the angle characterizing the HBs confirm that the conditions defining the water-water HBs change inside the PAA brush layer as a function of the charges on the PAA brush layer. Furthermore, we show that the use of the generic definition of HBs, as compared to using our simulation-motivated modified definition of water-water HBs, overpredict the number of water-water HBs inside the PAA brush layer. Finally, we employ this all-atom-MD-ML framework to quantify the effect of other types of screening counterions (Li + , Ca 2+ , and Y 3+ ions) in determining the water-water interactions and water-water HB properties inside the PAA brush layer. Furthermore, the findings of the present study, confirming the weakening of water-water HBs inside the PAA brush layer, points to the possibility that the water molecules will be more available for hydrating the brush layer and counterions, thereby leading to a more pronounced wetting of the PAA brush layer.

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

Effect of Helium Flux on Near-surface Helium Accumulation in Plasma-Exposed Tungsten

We report results of object kinetic Monte Carlo (OKMC) simulations aimed at understanding the effect of helium flux on the near-surface helium accumulation in plasma-facing tungsten, which is initially defect-free and has a W(100) surface orientation. These OKMC simulations are performed at 933 K for fluxes ranging from 10 22 to 4 × 10 25 He/m 2 s, with 100 eV helium atoms impinging on a W(100) surface up to a maximum fluence of 4×10 19 He/m 2 . In the near-surface region, helium clusters interact elastically with the free surface. The interaction is attractive and results in the drift of mobile helium clusters towards the surface as well as increased trap mutation rates. The associated kinetics and energetics of the above-mentioned processes obtained from molecular dynamics simulations are also considered. The OKMC simulations indicate that as the flux decreases, the retention of implanted helium decreases, and its depth distribution shifts to deeper below the surface in initially pristine tungsten. Furthermore, the fraction of retained helium diffusing into the bulk increases as well, so much so that for 10 22 He/m 2 s, almost all of the retained helium diffused into the bulk with minimal/negligible near-surface helium accumulation. At a given flux, with increasing fluence, the fraction of retained helium initially decreases and then starts to increase after reaching a minimum. The occurrence of the retention minimum shifts to higher fluences as the flux decreases. Although the near-surface helium accumulation spreads deeper into the material with decreasing flux and increasing fluence, the spread appears to saturate at depths between 80 and 100 nm. Finally, we present a detailed analysis of the influence of helium flux on the size and depth distribution of total helium and helium bubbles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

pnnl/mlprotein

The molecular graphics program, PyMol, was employed to produce all images in the training/validation/testing datasets. A PyMol script identified the clusters of interest from a list of all proteins containing 1) Fe (III) and Fe (II) ion and 2) FeS, respectively. A PyMol script and a Python webscraper were also used to classify [2Fe-2S] cofactors into three classes: standard, Rieske, and other and to identify iron atoms with four cysteine residues. We also used the pyMol script to rotate amino acids for training

Oostrom, Marjolein↗
DOE OSTI2023

Probing the Potential Energy Profile of the I + (H 2 O) 3 → HI + (H 2 O) 2 OH Forward and Reverse Reactions: High Level CCSD(T) Studies with Spin-Orbit Coupling Included

Three different pathways for the atomic iodine plus water trimer reaction I + (H 2 O) 3 → HI + (H 2 O) 2 OH were preliminarily examined by the DFT-MPW1K method. Related to previous predictions for the F/Cl/Br + (H 2 O) 3 reactions, three pathways for the I + (H 2 O) 3 reaction are linked in terms of geometry and energetics. To legitimize the results, the “gold standard” CCSD(T) method was employed to investigate the lowest-lying pathway with the correlation-consistent polarized valence basis set up to cc-pVQZ(-PP). According to the CCSD(T)/cc-pVQZ(-PP)//CCSD(T)/cc-pVTZ(-PP) results, the I + (H 2 O) 3 → HI + (H 2 O) 2 OH reaction is predicted to be endothermic by 47.0 kcal mol -1 . The submerged transition state is predicted to lie 43.7 kcal mol -1 above the separated reactants. The I···(H 2 O) 3 entrance complex lies below the separated reactants by 4.1 kcal mol -1 , and spin-orbit coupling has a significant impact on this dissociation energy. The HI···(H 2 O) 2 OH exit complex is bound by 4.3 kcal mol -1 in relation to the separated products. Compared with simpler I + (H 2 O) 2 and I + H 2 O reactions, the I + (H 2 O) 3 reaction is energetically between them in general. It is speculated that the reaction between the iodine atom and the larger water clusters may be energetically analogous to the I + (H 2 O) 3 reaction. The iodine reaction I + (H 2 O) 3 is connected with the analogous valence isoelectronic bromine/chlorine reactions Br/Cl + (H 2 O) 3 but much different from the F + (H 2 O) 3 reaction. Significant difference with other halogen systems, especially for barrier heights, are seen for the iodine systems.

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

Revealing the Hidden Third Dimension of Point Defects in Two-Dimensional MXenes

Point defects govern many important functional properties of two-dimensional (2D) materials. However, resolving the three-dimensional (3D) arrangement of these defects in multi-layer 2D materials remains a fundamental challenge, hindering rational defect engineering. Here, we overcome this limitation using an artificial intelligence-guided electron microscopy workflow to map the 3D topology and clustering of atomic vacancies in Ti3C2TX MXene. Our approach reconstructs the 3D coordinates of vacancies across hundreds of thousands of lattice sites, generating robust statistical insight into their distribution that can be correlated with specific synthesis pathways. This large-scale data enables us to classify a hierarchy of defect structures-from isolated vacancies to nanopores-revealing their preferred formation and interaction mechanisms, as corroborated by molecular dynamics simulations. This work provides a generalizable framework for understanding and ultimately controlling point defects across large volumes, paving the way for the rational design of defect-engineered functional 2D materials.

2D materials↗
DOE OSTI2025

Ab initio informed 20 Ne(p, pα) 16 O reaction elucidates the emergence of alpha clustering from chiral potentials

We report on the first ab initio informed α knock-out reaction in the intermediate-mass region, with the aim to probe the underlying chiral potential and its impact on the emergence of alpha clustering in this mass region. The theoretical predictions of the α+ 16 O clustering in the 20 Ne ground state, based on the ab initio symmetry-adapted no-core shell model with continuum, yield a triple differential cross section for 20 Ne(p, pα) 16 O that is in a remarkable agreement with the data. This allows us to examine predictions of surface and in-medium α-cluster features that emerge from the underlying realistic nucleon-nucleon interaction with no parameters fitted to nuclear data beyond the two-body system, and to compare these to the successful antisymmetrized molecular dynamics approach.

Sargsyan, G. H. [Michigan State Univ., East Lansin↗
DOE OSTI2022

Attosecond spectroscopy of size-resolved water clusters

Electron dynamics in water are of fundamental importance for a broad range of phenomena, but their real-time study faces numerous conceptual and methodological challenges. Here, in this paper, we introduce attosecond size-resolved cluster spectroscopy and build up a molecular-level understanding of the attosecond electron dynamics in water. We measure the effect that the addition of single water molecules has on the photoionization time delays of water clusters. We find a continuous increase of the delay for clusters containing up to four to five molecules and little change towards larger clusters. We show that these delays are proportional to the spatial extension of the created electron hole, which first increases with cluster size and then partially localizes through the onset of structural disorder that is characteristic of large clusters and bulk liquid water. These results indicate a previously unknown sensitivity of photoionization delays to electron-hole delocalization and indicate a direct link between electronic structure and attosecond photoionization dynamics. Our results offer new perspectives for studying electron-hole delocalization and its attosecond dynamics.

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

Quantification of Brønsted Acid Sites in Zeolites by Water Desorption Thermogravimetry

The many studies that have been conducted on water adsorption on zeolites were carried out in “adsorption” mode, from adsorption isotherms or single-point measurements, with the samples being first degassed at high temperature, and then exposed to water vapor at variable partial pressure. Here, we report Thermogravimetric Analysis (TGA) studies of HZSM-5 saturated with water at room temperature, thus studying water adsorption in “desorption” mode. In this work, TGA was used to ascertain the amount of water remaining as a function of the pretreatment temperature, and to identify the temperature at which one water molecule per framework Al atom persists. Solid-state 1H NMR spectroscopy, combined with molecular modeling, provided evidence for water molecules hydrogen bonded to Brønsted acid sites and formation of oxonium species for distinct water per framework Al atom. TGA is shown to be a simple and expeditious technique to measure the number of Brønsted acid sites in aluminosilicate zeolites, with advantages over the traditional methods of IR or TPD of dangerous adsorbed ammonia or pyridine. TGA of adsorbed water also provides information on the water clusters formed at the acid sites of HZSM-5, which are proportional to the number of acid sites.

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

Emerging Electrochemical Techniques for Probing Site Behavior in Single-Atom Electrocatalysts

Single-atom catalysts (SACs) have aroused tremendous interest over the past decade, particularly in the community of energy and environment-related electrocatalysis. A rapidly growing number of recent publications have recognized it as a promising candidate with maximum atomic utilization, distinct activity, and selectivity in comparison to bulk catalysts and nanocatalysts. However, the complexity of localized coordination environments and the dispersion of isolated sites lead to significant difficulties when it comes to gaining insight into the intrinsic behavior of electrocatalytic reactions. Furthermore, the low metal loadings of most SACs make conventional ensemble measurements less likely to be accurate on the subnanoscale. Thus, it remains challenging to probe the activity and properties of individual atomic sites by available commercial instruments and analytical methods. In spite of this, continuing efforts have lately focused on the development of advanced measurement methodologies, which are very useful to the fundamental understanding of SACs. There have recently been a number of in situ/operando techniques applied to SACs, such as electron microscopy, spectroscopy, and other analysis methods, which support relevant functions to identify the active sites and reaction intermediates and to investigate the dynamic behavior of localized structures of the catalytic sites. This Account aims to present recent electrochemical probing techniques which can be used to identify single-atomic catalytic sites within solid supports. First, we describe the basic principles of molecular probe methods for the study and analysis of electrocatalytic site behavior. In particular, the in situ probing technique enabled by surface interrogation scanning electrochemical microscopy (SISECM) can measure the active site density and kinetic rate with high resolution. An alternative electrochemical probing technique is further demonstrated on the basis of single-entity electrochemistry, which allows the unique electrochemical imaging of the size and catalytic rate of single atoms, molecules, and clusters. Next, the merits and limitations of different electrochemical techniques are then discussed, along with perspectives for future prospects. Apart from this, we further showcase the powerful capability of emerging electrochemical probing techniques for determining significant effects and properties of SACs for various electrocatalytic reactions, including oxygen reduction and evolution, hydrogen evolution, and nitrate reduction. Overall, electrochemical techniques with atomic resolution have greatly increased opportunities for observing, measuring, and understanding the surface and interface chemistry during energy conversion. In the future, it is anticipated that the development of electrochemical probing techniques will be advanced with innovative perspectives on the behavior and features of SACs. We hope that this Account can contribute in several ways to promoting the fundamental knowledge and technical progress of emerging electrochemical measurements for studying SACs.

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

Blue phosphorene reactivity on the Au(111) surface

Abstract The synthesis of blue phosphorene by molecular beam epitaxy (MBE) has recently come under the spotlight due to its potential applications in electronic and optoelectronic devices. However, this synthesis remains a significant challenge. The surface reactivity between the P atoms and the Au atoms should be considered for the P/Au(111) system. In the MBE process, the temperature of the substrate is a key parameter for the growth of blue phosphorene. During the initial growth stage, irregularly shaped Phosphorus clusters grow on top of Au(111) surface at room temperature. When the substrate temperature is increased, these clusters transform into a phosphorene-like structure with a honeycomb lattice. An atom exchange reaction is observed between the P and first layer Au atoms under thermal activation at higher temperature, where the P atoms replace Au atoms to form a blue phosphorene structure within the top Au layer and at the step edges.

Materials Science↗
DOE OSTI2025

Resolving the Solvation Structure and Transport Properties of Aqueous Zinc Electrolytes from Salt-in-Water to Water-in-Salt Using Neural Network Potential

Zn Cl 2 solutions are promising electrolytes for aqueous zinc-ion batteries. Here, we report a joint computational and experimental study of the structural and dynamic properties of aqueous Zn Cl 2 electrolytes with concentrations ranging from salt-in-water to water-in-salt (WIS). By developing a neural network potential (NNP) model, we perform molecular dynamics (MD) simulations with accuracy but at much larger lengths and longer timescales. The NNP predicted structures are validated by the structure factors measured by X-ray total scattering experiments. The MD trajectories provide a comprehensive and quantitative picture of the Zn 2 + solvation shell structures. Additionally, we find that the O − H covalent bonds in water are strengthened with increasing salt concentration, thus expanding the electrochemical stability window of aqueous electrolytes. In terms of dynamic properties, the calculated and experimentally measured conductivities are in good agreement. Through the analysis of the calculated cation transference number, we propose a three-stage charge carrier transport mechanism with increasing concentration: independent ion transport, strongly correlated ion transport, and small positive charge carrier diffusion through negatively charged polymeric clusters. Our study provides fundamental atomic scale insights into the structure and transport properties of the Zn Cl 2 electrolyte that can aid the optimization and development of WIS electrolytes. Published by the American Physical Society 2025

25 ENERGY STORAGE↗
DOE OSTI2022

Transferability of data-driven, many-body models for CO 2 simulations in the vapor and liquid phases

Here, extending on the previous work by Riera et al. [J. Chem. Theory Comput. 16, 2246–2257 (2020)], we introduce a second generation family of data-driven many-body MB-nrg models for CO 2 and systematically assess how the strength and anisotropy of the CO 2 –CO 2 interactions affect the models’ ability to predict vapor, liquid, and vapor–liquid equilibrium properties. Building upon the many-body expansion formalism, we construct a series of MB-nrg models by fitting one-body and two-body reference energies calculated at the coupled cluster level of theory for large monomer and dimer training sets. Advancing from the first generation models, we employ the charge model 5 scheme to determine the atomic charges and systematically scale the two-body energies to obtain more accurate descriptions of vapor, liquid, and vapor–liquid equilibrium properties. Challenges in model construction arise due to the anisotropic nature and small magnitude of the interaction energies in CO 2 , calling for the necessity of highly accurate descriptions of the multidimensional energy landscape of liquid CO 2 . These findings emphasize the key role played by the training set quality in the development of transferable, data-driven models, which, accurately representing high-dimensional many-body effects, can enable predictive computer simulations of molecular fluids across the entire phase diagram.

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

Mechanisms of Retention and Transport of Fission Products in Virgin and Irradiated Nuclear Graphite

To fulfill the design and licensing requirements for advanced reactors, the U.S. Department of Energy (DOE) is specifically interested in targeted research to quantify the mechanisms of fission product transport and retention in graphitic grades that will likely be used in next-generation very/high temperature reactors, fluoride salt-cooled high temperature reactors and molten salt reactors. We have examined an integrated experimental-computational approach to determine the diffusivities of selected fission products – ruthenium (Ru) and silver (Ag) in five nuclear graphite grades – POCO AXF-5Q, POCO ZXF-5Q, PCEA, IG110, and NBG-18. Experimental investigation for diffusion of other fission products (iodine, cesium, and krypton) in graphite could not be carried out due to extreme evaporation due to high vapor pressure, despite several attempts incorporating various experimental schemes. Experimental investigation with irradiated graphite could not be carried out because necessary approval to transport irradiated graphite could not be obtained due to procedural reasons at ORNL and national policy change in UK. An effort was made to “produce” by ion sputtering, but it was discontinued because the outcome did not adequately represent irradiated graphite. Thin film and “cup-and-cap” methods were employed to deposit fission products, diffusion annealed were carried out in an argon atmosphere, and concentration profiles were determined by dynamic secondary ion mass spectroscopy depth profiling technique. Regardless of graphite type, anomalous diffusion behavior of Ru (i.e., lower than expected based on Arrhenius temperature dependence) at an intermediate temperature, around 600 to 700 °C, was observed. Diffusion behavior of silver on the other hand agreed well to the Arrhenius temperature dependence although the excessive evaporation of silver hampered the accurate determination of diffusivities at temperature higher than 900 °C. In general, silver was observed to diffuse faster than Ru in the temperature range compared from 500 to 900 °C, had consistently had higher pre-exponential factor corresponding to higher jump frequency. To understand the abnormal diffusional behavior of Ru in graphite in intermediate temperature range, a detailed microstructural analysis was carried out. Clustering of Ru into particles on the surface of graphite was observed around 700 °C corresponding to de-wetting and spheroidization. This change would influence the boundary condition of diffusion, although features of clustering was much smaller (<1 mm) than the SIMS raster size (120 x 120 mm). Raman spectroscopy also demonstrated that Raman peak intensity for Ru decreased significantly for samples annealed higher than 700 °C and presence of Ru in graphite altered the crystallinity of graphite. To further elucidate the abnormal diffusion behavior observed in intermediate temperature investigated, the vibrational properties of graphene, bilayer graphene and Kr atoms in bilayer graphene were investigated with molecular dynamics. Krypton, which is a fission product, was introduced between a bilayer graphene system to evaluate (i) whether the phonon dispersion curves change with high temperature in bilayer graphene (ii) whether the addition of Kr atoms between the bilayer introduces any significant change in the vibrational properties, and (iii) whether any change in the vibrational properties can be correlated to the diffusive behavior of Kr atoms. The results show that the Kr atoms indeed exhibit an anomalous diffusive behavior at 550°C. Phonon analysis shows that a few phonon frequencies are enhanced as much as 20 to 30 % with the addition of just 0.2 at.% Kr. So, this may explain the unexpected self-diffusivity of fission products.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

On the Nature of Extra-Framework Aluminum Species and Improved Catalytic Properties in Steamed Zeolites

Steamed zeolites exhibit improved catalytic properties for hydrocarbon activation (alkane cracking and dehydrogenation). The nature of this practically important phenomenon has remained a mystery for the last six decades and was suggested to be related to the increased strength of zeolitic Bronsted acid sites after dealumination. We now utilize state-of-the-art infrared spectroscopy measurements and prove that during steaming, aluminum oxide clusters evolve (due to hydrolysis of Al out of framework positions with the following clustering) in the zeolitic micropores with properties very similar to (nano) facets of hydroxylated transition alumina surfaces. The Bronsted acidity of the zeolite does not increase and the total number of Bronsted acid sites decreases during steaming. O 5 Al(VI)-OH surface sites of alumina clusters dehydroxylate at elevated temperatures to form penta-coordinate Al 1 O 5 sites that are capable of initiating alkane cracking by breaking the first C-H bond very effectively with much lower barriers (at lower temperatures) than for protolytic C-H bond activation, with the following reaction steps catalyzed by nearby zeolitic Bronsted acid sites. This explains the underlying mechanism behind the improved alkane cracking and alkane dehydrogenation activity of steamed zeolites: heterolytic C-H bond breaking occurs on Al-O sites of aluminum oxide clusters confined in zeolitic pores. Our findings explain the origin of enhanced activity of steamed zeolites at the molecular level and provide the missing understanding of the nature of extra-framework Al species formed in steamed/dealuminated zeolites.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Uncertainty-aware molecular dynamics from Bayesian active learning for phase transformations and thermal transport in SiC

Abstract Machine learning interatomic force fields are promising for combining high computational efficiency and accuracy in modeling quantum interactions and simulating atomistic dynamics. Active learning methods have been recently developed to train force fields efficiently and automatically. Among them, Bayesian active learning utilizes principled uncertainty quantification to make data acquisition decisions. In this work, we present a general Bayesian active learning workflow, where the force field is constructed from a sparse Gaussian process regression model based on atomic cluster expansion descriptors. To circumvent the high computational cost of the sparse Gaussian process uncertainty calculation, we formulate a high-performance approximate mapping of the uncertainty and demonstrate a speedup of several orders of magnitude. We demonstrate the autonomous active learning workflow by training a Bayesian force field model for silicon carbide (SiC) polymorphs in only a few days of computer time and show that pressure-induced phase transformations are accurately captured. The resulting model exhibits close agreement with both ab initio calculations and experimental measurements, and outperforms existing empirical models on vibrational and thermal properties. The active learning workflow readily generalizes to a wide range of material systems and accelerates their computational understanding.

36 MATERIALS SCIENCE↗