Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex mixtures”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

In Situ Determination of Speciation and Local Structure of NaCl–SrCl 2 and LiF–ZrF 4 Molten Salts

Understanding the local environment of the metal atoms in salt melts is important for modeling the properties of melts and predicting their behavior and thus helping enable the development of technologies such as molten salt reactors and solar-thermal power systems and new approaches to recycling rare-earth metals. Toward that end, we have developed an in situ approach for measuring the coordination of metals in molten salt coupling X-ray absorption spectroscopy (XAS) and Raman spectroscopy. Our approach was demonstrated for two salt mixtures (1.9 and 5 mol % SrCl 2 in NaCl, 0.8 and 5 mol % ZrF 4 in LiF) at up to 1100 °C. Near-edge (X-ray absorption near-edge structure, XANES) and extended X-ray absorption fine structure (EXAFS) spectra were measured. The EXAFS response was modeled using ab initio FEFF calculations. Strontium’s first shell is observed to be coordinated with chlorine (Sr 2+ –Cl – ) and zirconium’s first shell is coordinated by fluorine (Zr 4+ –F – ), both having coordination numbers that decrease with increasing temperature. Multiple zirconium complexes are believed to be present in the melt, which may interfere and distort the EXAFS spectra and result in an anomalously low zirconium first shell coordination number. Finally, the use of boron nitride (BN) powder as a salt diluent for XAFS measurements was found to not interfere with measurements and thus can be used for investigations of such systems.

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

Reductive catalytic fractionation of cotton stalks: catalytic strategy for tuning the selectivity of phenolic monomers

Pretreatment of lignocellulosic biomass is a primary step to delink lignin from the lignin-carbohydrate complex for bioethanol production and other value-added products. Presently, these processes yield lignin (technical) with significant structural changes compared to native lignin and make it difficult to valorize. Reductive catalytic fractionation (RCF) is an advanced pretreatment process to valorize native lignin to selective phenolic monomers before capitalizing on carbohydrates. Herein, cotton stalks (CS) were pretreated with ethanol: water mixture catalyzed by bimetallic, Ru-Ni/HY (RNY) and Ru-Fe/HY (RFY), and trimetallic, Ru-Fe-Ni/HY (RFNY), catalysts using the endogenously produced H 2 . The favorable catalytic properties of the RFNY catalyst, including high reducibility, acidity, metal dispersion, the synergistic effect of all the metals, and appropriate pore and particle size resulted in the highest catalytic activity. Maximum 92% delignification efficiency was achieved with 91% holocellulose retention. The extracted lignin after depolymerization in the same step produced a 19 wt.% yield of phenolic monomers. Increasing the duration of the reaction from 2 to 6 h increased the selectivity to saturated alkyl chain phenolics (1 to 16 wt.%) at the expense of hemicellulose (50 to 9%). The fate of carbohydrates left behind after the RCF (CS-RCF) process was analyzed through enzymatic saccharification and compared with raw CS and acid-pretreated CS. In conclusion, CS-RCF yielded maximum glucose yield (45 g/L) within an incubation period of 72 h.

09 BIOMASS FUELS↗
DOE OSTI2025

Bioconversion of self-neutralized chemically depolymerized lignin streams into polyhydroxyalkanoates

Lignin is a chemically complex, diverse, and abundant plant polymer mainly composed of aromatic monomers. These aromatic monomers make lignin a potential source of aromatics and a viable substitute for petrochemically-derived aromatics. However, the structural recalcitrance of lignin requires harsh reagents from the chemical process and expensive catalysts for effective depolymerization. This chemical process often results in poor yields of chemical intermediate mixtures of varying bioavailability and/or toxicity. Furthermore, the cost of additional reagents required to separate or detoxify these intermediates makes processing lignin impractical. We report progress towards the use of such chemically depolymerized lignin streams by employing bacterial strains to produce polyhydroxyalkanoates (PHA). PHAs are a group of biodegradable microbial polyesters that have potential as a replacement for petroleum-based plastics. In this study we utilized two distinct lignin streams obtained after chemical depolymerization of lignin under alkaline and acidic pH in the presence of catalysts. We mixed the alkali-treated depolymerized stream with the acid-treated depolymerized stream to create a solution of neutral-pH chemically depolymerized lignin (CDL). We found moderate to substantial growth of both native and non-native PHA producers on the mixture of CDL as well as its aliphatic and aromatic components. PHA was detected by Sudan Black B staining in C. necator H16, P. putida KT2440, and E. coli LSBJ STQKAB grown on mixed CDL as the sole carbon source. For C. necator H16 and E. coli LSBJ STQKAB we found that PHA content was greater when grown on mixed CDL when compared to their preferred carbon source by GC-FID quantification. Our study provided progress towards a cost-competitive, sustainable, and industrially relevant use for lignin.

Winkler, Gordon L. W. [Univ. of Wisconsin, Madison↗
DOE OSTI2020

Sequential ensemble transform for Bayesian inverse problems

In this work, we present the Sequential Ensemble Transform (SET) method, an approach for generating approximate samples from a Bayesian posterior distribution. The method explores the posterior distribution by solving a sequence of discrete optimal transport problems to produce a series of transport plans which map prior samples to posterior samples. We prove that the sequence of Dirac mixture distributions produced by the SET method converges weakly to the true posterior as the sample size approaches infinity. Furthermore, our numerical results indicate that, when compared to standard Sequential Monte Carlo (SMC) methods, the SET approach is more robust to the choice of Markov mutation kernels and requires less computational efforts to reach a similar accuracy when used to explore complex posterior distributions. Finally, we describe adaptive schemes that allow to completely automate the use of the SET method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Interplay between Mixed and Pure Exciton States Controls Singlet Fission in Rubrene Single Crystals

Singlet fission (SF) is a multielectron process in which one singlet exciton S converts into a pair of separated triplet excitons T. SF is widely studied as it may help overcome the Shockley−Queisser efficiency limit for semiconductor photovoltaic cells. To elucidate and control the SF mechanism, great attention has been given to the identification of intermediate states in SF materials, which often appear elusive due to the complexity and fast time scales of the SF process. Here, we apply 14 fs-1 ms transient absorption techniques to high-purity rubrene single crystals to disentangle the intrinsic fission dynamics from the effects of defects and grain boundaries and to identify reliably the fission intermediates. Our data demonstrates that above-gap excitation directly generates a hybrid vibronically assisted mixture of singlet state and triplet-pair multiexciton [S/TT], which rapidly (<100 fs) and coherently branches into pure singlet or triplet excitations. The relaxation of [S/TT] to S is followed by a relatively slow and temperature-activated (48 meV activation energy) incoherent fission process. The SF competing pathways and intermediates revealed here unify the observations and models presented in previous studies of SF in rubrene and offer alternative strategies for the development of SF-enhanced photovoltaic materials.

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

Machine learning prediction of the mechanical properties of refractory multicomponent alloys based on a dataset of phase and first principles simulation

In this work, a dataset including structural and mechanical properties of refractory multicomponent alloys was developed by fusing computations of phase diagram (CALPHAD) and density functional theory (DFT). The refractory multicomponent alloys, also named refractory complex concentrated alloys (CCAs) which contain 2–5 types of refractory elements were constructed based on Special Quasi-random Structure (SQS). The phase of alloys was predicted using CALPHAD and the mechanical property of alloys with stable and single body-centered cubic (BCC) at high temperature (over 1,500°C) was investigated using DFT-based simulation. As a result, a dataset with 393 refractory alloys and 12 features, including volume, melting temperature, density, energy, elastic constants, mechanical moduli, and hardness, were produced. To test the capability of the dataset on supporting machine learning (ML) study to investigate the property of CCAs, CALPHAD, and DFT calculations were compared with principal components analysis (PCA) technique and rule of mixture (ROM), respectively. It is demonstrated that the CALPHAD and DFT results are more in line with experimental observations for the alloy phase, structural and mechanical properties. Furthermore, the data were utilized to train a verity of ML models to predict the performance of certain CCAs with advanced mechanical properties, highlighting the usefulness of the dataset for ML technique on CCA property prediction.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Exploring the Nature of f-Element Soft Donor Interactions Using Electronically Tunable Azolate Ionic Liquids

This project was undertaken to advance the understanding of how f-elements interact with moderately soft donors, a heavily investigated yet open question which is of prime importance in spent nuclear fuel processing and fundamental inorganic chemistry. During the course of the project, based on exciting results, a stretch goal was developed to study the hydrolysis products of transuranic actinide metals, a somewhat understudied field even with its significance in nuclear fuel cycle and impacts in environmental chemistry. The stretch goal was to take our serendipitous discovery of an easy route to isolation of crystalline multinuclear ƒ-element hydrolysis products, and apply it to gaining a mechanistic understanding of Pu(III/IV) colloid formation. The simplicity of our techniques should lend themselves to the remote handling required for study of many transuranic elements. We developed several methodologies using azolium azolate chemistry to overcome ƒ-element hydrolysis problems that prohibit the isolation of ƒ-element soft donor complexes and to build a crystallographic library of ƒ-element N-donor complexes as a means to understand the fundamental differences between actinide and lanthanide interactions with moderately soft donor ligands. Our next major endeavor will be to transfer this chemistry from 4ƒ elements to transuranic elements, particularly in the study of hydrolysis of Pu(III/IV). While our work is fundamental in nature, applications of the knowledge we are generating should be felt in such diverse fields as catalysis, separations in general, nuclear waste remediation specifically, and many other applications in f-element magnetic and luminescent properties. The potential ramifications of the consistent and robust formation of hydrolysis controlled hexanuclear lanthanide structures are enormous, with future uses being catalyst formation, higher-nuclearity structure synthesis using our hexanuclear motif as a fundamental building block, Pu waste remediation, separations, and many other potential applications resulting from characteristic magnetic and luminescent properties of lanthanide polynuclear structures. Three synthetic methodologies (direct mixing with variable stoichiometries, use of volatile solvent, metathesis) were developed starting with 7 acidic and 6 basic azoles to obtain ionic liquids suitable for f-element coordination. Proton transfer by acidic/basic azole combination led to suitable low melting salts and two cocrystals. Acid/base reaction of azoles with soft-donor permanent cations of ([X 4444 ][OH] (where [X 4444 ] + = tetrabutylammonium [N 4444 ] + or tetrabutylphosphonium [P 4444 ] + ) with weakly acidic azoles including imidazole, 1,2,3-triazole, 1,2,4-triazole, 5-aminotetrazole, 4,5-dicyanoimidazole, and 2-amino-4,5-dicyanoimidazole) revealed several suitable low-melting salts. Metathesis reactions of Na(azolate) were conducted by first using weakly acidic azoles including 4,5-dicyanoimidazole, 2-amino-4,5-dicyanoimidazole, 5-aminotetrazole, and 1,2,4-triazole to form sodium or lithium salts using group(I) hydroxides in methanolic solutions. The best results were obtained by reacting the basic and acidic azoles directly in 1:1 or 3:1 ratios at elevated temperatures. Twenty-two azole mixtures were identified which are either low melting solids or room temperature liquids. Each of the low melting solids was confirmed to be a new solid phase, each of which is being further characterized. The liquids and solids are anticipated to be ILs, eutectics, or partially ionized systems, all of which will be suitable for the dissolution of f-element salts. Five new synthetic methodologies were developed to finding suitable crystallization conditions for f-element complexation with the goal of finding simple one pot reaction syntheses and crystallization strategies that could be used under the demanding conditions of transuranic chemistries. These synthetic methods yield many new crystalline phases which were studied by single crystal X-ray diffraction.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

AmeriFlux US-xTA NEON Talladega National Forest (TALL)

This is the AmeriFlux version of the carbon flux data for the site US-xTA NEON Talladega National Forest (TALL). Site Description - The TALL terrestrial field site covers 5300 hectares within the larger Oakmulgee District of the Talladega National Forest (475,000 hectares) in west-central Alabama. TALL is located within the Gulf Coastal Plain which is the dominant physiographic province (52%) of NEON's Domain 8: Ozarks Complex. It is also where the upper coastal plain gives rise to the Appalachian foothills. The lands of this region are rich a mosaic of forest types and habitats. Steep ridges to rolling hills fading to hardwood bottoms make up the diversity of forest conditions found in Oakmulgee. While known for its longleaf forests, over 40% of Oakmulgee is covered with a mixture of hardwoods and wetlands.

Network), NEON (National Ecological Observatory↗
DOE OSTI2023

AmeriFlux FLUXNET-1F US-xTA NEON Talladega National Forest (TALL)

This is the AmeriFlux Management Project (AMP) created FLUXNET-1F version of the carbon flux data for the site US-xTA NEON Talladega National Forest (TALL). This is the FLUXNET version of the carbon flux data for the site US-xTA NEON Talladega National Forest (TALL) produced by applying the standard ONEFlux (1F) software. Site Description - The TALL terrestrial field site covers 5300 hectares within the larger Oakmulgee District of the Talladega National Forest (475,000 hectares) in west-central Alabama. TALL is located within the Gulf Coastal Plain which is the dominant physiographic province (52%) of NEON's Domain 8: Ozarks Complex. It is also where the upper coastal plain gives rise to the Appalachian foothills. The lands of this region are rich a mosaic of forest types and habitats. Steep ridges to rolling hills fading to hardwood bottoms make up the diversity of forest conditions found in Oakmulgee. While known for its longleaf forests, over 40% of Oakmulgee is covered with a mixture of hardwoods and wetlands.

Network), NEON (National Ecological Observatory↗
DOE OSTI2026

Stable Co-valorization of Carbon Dioxide and Methane via Dynamic Reconstruction of a Metal Oxide Solid Solution Catalyst

Dry reforming of methane (DRM) is a process that converts two greenhouse gases (methane and carbon dioxide) into syngas, a mixture of H 2 and CO, that can lead to a variety of value-added chemicals. Owing to its endothermic nature, high reaction temperatures up to 800 °C are typically required and the grand challenge lies in developing robust catalysts without sintering and coking-induced deactivation during the long-term on-stream operation. Towards this aim, herein, a robust complex oxide-supported NiCu alloy catalyst was generated in situ during DRM. By leveraging the configurational stability of a solid oxide solution precursor, tightly anchored NiCu bimetallic nanoparticles were in situ exsoluted and acted as the active sites in DRM. The as-afforded catalyst exhibited stable performance for DRM due to the ability to repel coke off the surface as the reaction proceeds. Kinetic experiments along with top surface characterization detail the reconstruction behavior of the solid oxide solution under DRM reaction conditions. In conclusion, the fundamental insights from this work provide guidance on generating resistant and flexible catalysts via in situ active sites formation from easily synthesized metal oxide solid solutions.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Design of Multi-Stage Solvent Extraction Process for Separation of Rare Earth Elements

Flowsheet design and stage determination for the separation of rare earth elements (REEs) using solvent extraction (SX) is a challenging task because of the chemical similarity of the REEs. Low separation factors between the elements and complex equilibrium chemistry provide unique challenges to designing an efficient flowsheet for the separation of elements. The multi-stage nature of the SX process adds further complexity, making the assessment of products for a proposed design and stage combination difficult. Therefore, to develop a SX flowsheet, it is essential to quantify the performance for various design and separation conditions. This paper attempts to address the challenge by utilizing an equilibrium and process modeling approach. Results from a bench-scale study performed on a 10 g/L rare earth salt mixture were used in studying the extraction/stripping behavior and developing equilibrium models. DEHPA with TBP as a phase modifier was used as an extractant, while hydrochloric acid was utilized as a stripping agent. The results obtained were used in developing extraction/stripping models, which were integrated into a process framework of a SX train in a Matlab/Simulink environment. The models were programmed as a function block routine and used for developing a flowsheet, which was simulated for differing separation and design conditions. To identify optimum stage combinations, a particle swarm optimization (PSO) routine was developed and implemented for each SX train. Recovery and purity of elements of interest were used as objective function criteria. The stage combination leading to the minimization of the objective function was used to identify the optimum stage combination for a series of SX trains to attempt a balance of purity and recovery. The models and optimization method were implemented to separate a feed mixture containing REEs, which indicated that 99.52 and 85.41 percent purity is achievable for Yttrium and Lanthanum separation using 8-12-3 and 10-3-5 stage combination for loading, scrubbing, and striping. The model also indicated difficult separability between neodymium, praseodymium, and cerium.

Srivastava, Vaibhav (ORCID:0000000212645987)↗
DOE OSTI2024

Credible Criticality Safety Margin in the 30B Package with LEU+ UF 6 and Hypothetical Water Ingress

The commercial nuclear industry is pursuing advancements in fuel and reactor design that increase the uranium enrichment above 5 wt. % 235 U. These advancements will necessitate the ability to transport bulk quantities of UF 6 at increased enrichments. Currently, the 30B cylinder is the primary container used by the industry for UF 6 storage and transportation and has a long history of successful shipments. This container can support up to 2,277 kg of UF 6 at a maximum enrichment of 5 wt. % 235 U. Previous evaluations have assessed the potential impact of criticality safety for 30B transport at higher enrichments but assumed moderator intrusion would not require evaluation. Although current regulations allow for the exception of moderator intrusion for UF 6 packages through the design and quality control of the package content, this exception is limited to enrichments up to 5 wt. % 235 U. Thus, an investigation of moderator intrusion into a 30B cylinder should be performed. Moderator intrusion into a 30B cylinder is a unique condition for criticality safety evaluation in transportation because of the violent chemical reactions that occur between UF 6 and H 2 O. The resulting intrusion is strongly dependent on the breach size, breach location, breach interface solid/ullage of the UF 6 content, the UF 6 distribution (which is temperature dependent), and temperature/pressure conditions which are dynamic in accident conditions. Additionally, the complexity of the HF-UO 2 F 2 -H 2 O interface that occurs during the event can influence the potential solubility of uranium in the system, as well as influence the amount of UF 6 reaction with H 2 O. With the operating experience from Orano Federal Services LLC, the UF 6 chemical expertise from Oak Ridge National Laboratory (ORNL), and the criticality safety expertise from ORNL, this paper evaluates the neutronic conditions (i.e., k eff ) that apply the understood chemistry and experimental conditions that occur during moderator intrusion of a breached 30B cylinder under postulated accidents. This report examines the historical evaluations of UF 6 transport and expands these evaluations for the enrichments expected for nuclear industry advancement. These simulations primarily feature a homogeneous mixture of UF 6 and H 2 O as an infinite media system and as a sphere with water reflection given an impurity limit of 0.5 wt. % UF 6 . This report demonstrates that this H/U limit is valid up to 8 wt. % enrichment for both H 2 O and HF as the moderating mixture. The water-reflected homogeneous mixed spheres evaluated demonstrate the amount of safety margin applied by restricting the 30B cylinder impurity limit to 0.5 wt. %. Additionally, this report evaluates moderator ingress scenarios of a 30B cylinder. The 30B cylinder simulations vary the cylinder orientation, the mass of UF 6 in the system, the mass of H 2 O in the system, and how much H 2 O has reacted with UF 6 . The moderator ingress of a 30B cylinder was evaluated for the water-reflected homogeneous mixed spheres to demonstrate the amount of safety margin applied by restricting the 30B cylinder impurity limit to 0.5 wt. %. These simulations are standard practice and are independent of the UF 6 -H 2 O reaction. This work also explored more complex simulations that layer the H 2 O over the UF 6 , which incorporates the amount of reacted UF 6 as a separate layer and considers solubility limits of UO 2 F 2 in H 2 O as it is produced. This approach is intended to simulate the postulated event of a large hairline crack occurring on a submerged cylinder. Finally, this work considers the potential of H 2 O mixing into the UF 6 solid heterogeneously using a sponge-like model. The water ingress into the UF 6 is treated as random size spheres, and the reaction products form into a layer over the spheres.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Bayesian sparse learning with preconditioned stochastic gradient MCMC and its applications

Deep neural networks have been successfully employed in an extensive variety of research areas, including solving partial differential equations. Despite its significant success, there are some challenges in effectively training DNN, such as avoiding overfitting in over-parameterized DNNs and accelerating the optimization in DNNs with pathological curvature. Here, we propose a Bayesian type sparse deep learning algorithm. The algorithm utilizes a set of spike-and-slab priors for the parameters in the deep neural network. The hierarchical Bayesian mixture will be trained using an adaptive empirical method. That is, one will alternatively sample from the posterior using preconditioned stochastic gradient Langevin Dynamics (PSGLD), and optimize the latent variables via stochastic approximation. The sparsity of the network is achieved while optimizing the hyperparameters with adaptive searching and penalizing. A popular SG-MCMC approach is Stochastic gradient Langevin dynamics (SGLD). However, considering the complex geometry in the model parameter space in nonconvex learning, updating parameters using a universal step size in each component as in SGLD may cause slow mixing. To address this issue, we apply a computationally manageable preconditioner in the updating rule, which provides a step-size parameter to adapt to local geometric properties. Moreover, by smoothly optimizing the hyperparameter in the preconditioning matrix, our proposed algorithm ensures a decreasing bias, which is introduced by ignoring the correction term in the preconditioned SGLD. According to the existing theoretical framework, we show that the proposed algorithm can asymptotically converge to the correct distribution with a controllable bias under mild conditions. Numerical tests are performed on both synthetic regression problems and learning solutions of elliptic PDE, which demonstrate the accuracy and efficiency of the present work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Deep Eutectic Solvents: A Review of Fundamentals and Applications

Deep eutectic solvents (DESs) are an emerging class of mixtures characterized by significant depressions in melting points compared to those of the neat constituent components. These green materials are promising for applications as inexpensive ’designer’ solvents exhibiting a host of tunable physicochemical properties. Here, a detailed review of the current literature reveals the lack of predictive understanding of the microscopic mechanisms that govern the structure-property relationships in this class of solvents. Complex hydrogen bonding behavior is postulated as the root cause of their melt-ing point depressions and physicochemical properties; to understand these hydrogen bonded networks, it is imperative to study these systems as dynamic entities using both simulations and experiments. This review emphasizes recent research efforts in order to elucidate the next steps needed to develop a fundamental framework for understanding DESs. It covers recent developments in DES research, frames outstanding scientific questions, and identifies gainful research thrusts that would advance the field towards predictive models and fundamental understanding of these solvents.

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

Building molecular model series from heterogeneous CryoEM structures using Gaussian mixture models and deep neural networks

Cryogenic electron microscopy (CryoEM) produces structures of macromolecules at near-atomic resolution. However, building molecular models with good stereochemical geometry from those structures can be challenging and time-consuming, especially when many structures are obtained from datasets with conformational heterogeneity. Here we present a model refinement protocol that automatically generates series of molecular models from CryoEM datasets, which describe the dynamics of the macromolecular system and have near-perfect geometry scores. This method makes it easier to interpret the movement of the protein complex from heterogeneity analysis and to compare the structural dynamics observed from CryoEM data with results from other experimental and simulation techniques.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Local Bayesian Dirichlet mixing of imperfect models

Abstract To improve the predictability of complex computational models in the experimentally-unknown domains, we propose a Bayesian statistical machine learning framework utilizing the Dirichlet distribution that combines results of several imperfect models. This framework can be viewed as an extension of Bayesian stacking. To illustrate the method, we study the ability of Bayesian model averaging and mixing techniques to mine nuclear masses. We show that the global and local mixtures of models reach excellent performance on both prediction accuracy and uncertainty quantification and are preferable to classical Bayesian model averaging. Additionally, our statistical analysis indicates that improving model predictions through mixing rather than mixing of corrected models leads to more robust extrapolations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Atomic cluster expansion potential for large scale simulations of hydrocarbons under shock compression

We present an Atomic Cluster Expansion (ACE) machine learned potential developed for high-fidelity atomistic simulations of hydrocarbons, targeting pressures and temperatures near and above supercritical fluid regimes for molecular fluids. A diverse set of stoichiometries were covered in training, including 1:0 (pure carbon), 1:4 (methane), and 1:1 (benzene), and rich bonding environments sampled at supercritical temperatures, hydrogen rich, reactive mixtures where metastable stoichiometries arise, including 1:2 (ethylene) and 1:3 (ethane). A high-fidelity training database was constructed by performing large-scale quantum molecular dynamic simulations [density functional theory (DFT) MD] of diamond, graphite, methane, and benzene. A novel approach to selecting structures from DFT MD is also presented, which allows for the rapid selection of unique DFT MD frames from complex trajectories. Comparisons to DFT and experimental data demonstrate that the presented ACE potential accurately reproduces isotherms, carbon melting curves, radial distribution functions, and shock Hugoniots for carbon and hydrocarbon systems for pressures up to 100 GPa and temperatures up to 6000 K for hydrocarbon systems and up to 9000 K for pure carbon systems. This work delivers a potential that can be used for accurate, large-scale simulations of shocked hydrocarbons and demonstrates a methodology for fitting and validating machine learning interatomic potentials to complex molecular environments, which can be applied to energetic materials in future works.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗