Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “REACTIVITY”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Numerical Studies of the On-Line Reactivity Monitoring Techniques in Subcritical Systems: Revisit of the Area-Ratio Method, the Neutron Source-Modulation Method, and the Flux-to-Current Ratio Method

In the subcritical systems, the subcriticality level must be monitored continuously to ensure the safe operation of the systems. Techniques such as the neutron source-modulation method and the flux-to-current ratio method have been proposed for on-line reactivity measurements of the subcritical systems. However, the accuracy of these methods is unknown in the actual applications. In this report, the fundamental concepts of these techniques are visited, with numerical simulation techniques developed to validate these on-line reactivity monitoring methodologies. Specifically, the Monte Carlo numerical simulation results show that the neutron source-modulation method is not accurate in determining the reactivity of the subcritical systems. The flux-to-current ratio method can provide on-line reactivity measurements, and its spatial effects have been studied and well understood through the numerical simulations of the fuel loading of the neutron source facility of Ukraine. In addition, a modified version of the area-ratio method is proposed for the on-line measurements of the absolute reactivity of the subcritical systems. Monte Carlo numerical simulations showed that this method can provide the absolute reactivity of the subcritical systems within ~5% error, if the neutron detectors are placed close to the fuel assembly interface with the reflector zone or in the reflector zones away from the neutron source.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Exploration of Fuel Property Impacts on the Combustion of Late Post Injections Using Binary Blends and High-Reactivity Ether Bioblendstocks

In this study, the impacts of fuel volatility and reactivity on combustion stability and emissions were studied in a light-duty single-cylinder research engine for a three-injection catalyst heating operation strategy with late post-injections. Here, N-heptane and blends of farnesane/2,2,4,4,6,8,8-heptamethylnonane were used to study the impacts of volatility and reactivity. The effect of increased chemical reactivity was also analysed by comparing the baseline #2 diesel operation with a pure blend of mono-ether components (CN > 100) representative of potential high cetane oxygenated bioblendstocks and a 25 vol.% blend of the mono-ether blend and #2 diesel with a cetane number (CN) of 55. At constant reactivity, little to no variation in combustion performance was observed due to differences in volatility, whereas increased reactivity improved combustion stability and efficiency at late injection timings. Fuels with higher reactivity were found to reduce engine-out hydrocarbon and carbon monoxide emissions while also achieving stable combustion at post-injection timings later than those achievable with #2 diesel fuel. The pure ether blend had the latest achievable post-injection timing of +30.5 CAD while still maintaining stable combustion (coefficient of variation of gross-indicated mean effective pressure < 5%). With post-injection timing adjusted to achieve a matched exhaust temperature of 300 °C, the ether-diesel CN 55 blend was observed to have slightly higher thermal efficiency in comparison to the baseline #2 diesel fuel. The results also indicate that cetane number may serve as a good indicator of combustion characteristics at late injection timings used for aftertreatment thermal management operation.

09 BIOMASS FUELS↗
DOE OSTI2023

Molecular determinants of cross-reactivity and potency by VH3-33 antibodies against the Plasmodium falciparum circumsporozoite protein

IGHV3-33 -encoded antibodies are prevalent in the human humoral response against the Plasmodium falciparum circumsporozoite protein (PfCSP). Among VH3-33 antibodies, cross-reactivity between PfCSP major repeat (NANP), minor (NVDP), and junctional (NPDP) motifs is associated with high affinity and potent parasite inhibition. However, the molecular basis of antibody cross-reactivity and the relationship with efficacy remain unresolved. Here, we perform an extensive structure-function characterization of 12 VH3-33 anti-PfCSP monoclonal antibodies (mAbs) with varying degrees of cross-reactivity induced by immunization of mice expressing a human immunoglobulin gene repertoire. We identify residues in the antibody paratope that mediate cross-reactive binding and delineate four distinct epitope conformations induced by antibody binding, with one consistently associated with high protective efficacy and another that confers comparably potent inhibition of parasite liver invasion. Our data show a link between molecular features of cross-reactive VH3-33 mAb binding to PfCSP and mAb potency, relevant for the development of antibody-based interventions against malaria.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Influence of physical and chemical hydrology on bioremediation of a U-contaminated aquifer informed by reactive transport modeling incorporating 238 U/ 235 U ratios

Microbially-catalyzed reductive immobilization of aqueous uranium (U) as a solid phase has been proposed as a U remediation technique. Both laboratory and field experiments have demonstrated that this reduction reaction alters the 238 U/ 235 U ratio, producing a 238 U-enriched U(IV) solid. In contrast, other major U reactive transport processes fractionate these isotopes much less. This suggests the potential to quantify the extent of bioreduction occurring in groundwater containing U using the 238 U/ 235 U ratio, which would substantially improve upon current practices largely relying on U concentration measurements alone. Current reactive transport models for uranium dynamics include only concentration measurements, which are strongly influenced by highly coupled reactive transformation and aqueous transport processes. The complex physical and chemical behavior of U potentially compromises such quantitative analysis of U storage and release. Here we report the first numerical reactive transport model which explicitly incorporates variations in the 238 U/ 235 U ratio of U and demonstrates improved interpretation of the principal chemical reactions and groundwater transport processes affecting the subsurface mobility and distribution of this widespread contaminant. Recent U bioreduction studies performed in a contaminated aquifer in Rifle, Colorado, USA applied Rayleigh distillation models to interpret U stable isotope fractionation observed as a result of acetate amendment. These simplified models were unable to resolve the spatiotemporal pattern of U isotope fractionation recorded in the aqueous solutes. Here, we employ the multi-component, isotope-enabled CrunchTope reactive transport software to interpret these measured U isotope ratios, and demonstrate accurate reproduction of observed trends in both geochemistry and 238 U/ 235 U ratios for two consecutive years of field experiments. Overall, our results indicate that accurately modeling both the U concentration and isotope ratio distributions greatly constrains the parameter space of the model. We find that the transport properties of U in the Rifle aquifer are governed by the presence of low-permeability regions, which the isotopes are uniquely sensitive to. When U reduction is spatially constrained by these low-permeability regions, the shift in the 238 U/ 235 U ratio becomes more muted. Accurate modeling of observed U isotope ratios thus provides a powerful means to better understand bioremediation, and the current study serves to advance the application of this novel method.

238U/235U↗
DOE OSTI2022

Impacts of engine lubrication oil-derived ash on soot oxidative reactivity on a catalytic gasoline particulate filter

Characterizing soot oxidation kinetics is crucial for understanding how gasoline particulate filters (GPFs) perform both in terms of filtration efficiency and pressure drop. The most common method for measuring soot oxidative reactivity is thermogravimetric analysis (TGA). Because TGA is an offline method, it may inaccurately predict how soot oxidizes on a GPF, especially if a catalytic washcoat is present that may enhance oxidation rates through surface interactions as the soot loads on the filter. In this work, a novel in-situ soot oxidative reactivity measurement method was developed. The method involved loading a catalytic GPF under high temperature conditions conducive to soot oxidation and measuring the filtration efficiency of 100 nm particles. A correlation was made between filtration efficiency and loaded soot mass on the filter to allow calculation of the soot oxidation rate. The method was evaluated in experiments using a 2.0 L gasoline direct injection (GDI) engine with three oils of varying additive packages, including an oil with a high zinc dialkyldithiophosphate (ZDDP) concentration, a non-additive, pure poly-alpha olefin (PAO) oil, and an oil with a high concentration of calcium sulfonate. Calculated specific soot oxidation rates ranged from ~0.06–0.6 min –1 and showed strong dependence on GPF inlet gas temperature and lubrication oil additive type. The results clearly demonstrated a catalytic effect of calcium-containing ash particles. Reactivity of soot produced by the engine running with the high calcium containing lubrication oil was increased. Similarly, the results indicated that the catalytic washcoat on the filter increased soot oxidation rates on-filter, especially at low soot loading. Conversely, the ZDDP oil additive exhibited a reactivity inhibiting effect, resulting in lower soot oxidation rates. Finally, this work represents the first known in-situ soot oxidative reactivity measurements on a GPF and elucidates the effect of lubrication oil additives and catalytic washcoat on oxidation rate.

42 ENGINEERING↗
DOE OSTI2025

A quantitative risk assessment framework for fault reactivation in underground hydrogen storage: Coupled simulation and deep learning approach

Underground hydrogen storage (UHS) is emerging as a critical solution for large-scale energy storage. However, like all subsurface fluid injection activities, UHS poses the risk of injection-induced fault reactivation. Accurate risk assessment is essential to ensuring the safety and efficiency of UHS operations. This study presents the development of deep-learning surrogate models for fault reactivation prediction in UHS, trained on a comprehensive database of fully coupled fluid flow-geomechanics simulations. Our findings reveal that analytical models often yield unreliable estimates, with errors up to 54% in the allowable injection pressure, potentially leading to a 40% reduction in UHS operational capacity. The developed surrogate models were incorporated into a quantitative risk assessment (QRA) framework, enabling probabilistic evaluation of fault reactivation risk while accounting for uncertainties in the input variables. Site-specific features, such as horizontal stress gradients, fault’s dip and strike angles, and operational parameters like bottom-hole injection pressure and well-fault distance, were identified as the primary drivers of fault reactivation across various stress regimes. Whereas other hydraulic, geological, and poroelastic reservoir properties were found to have a secondary impact. Notably, we observed that the risk of fault reactivation for a critically oriented fault with a static friction coefficient greater than 0.55 remains below 10% in a normal faulting stress regime. However, the risk significantly increases as the stress regime transitions from normal to strike-slip and ultimately to reverse faulting conditions. These findings underscore the importance of rigorous site characterization and comprehensive QRA evaluations to optimize UHS performance and minimize geomechanical risks.

25 ENERGY STORAGE↗
DOE OSTI2021

Physics-based, neural network force fields for reactive molecular dynamics: Investigation of carbene formation from [EMIM + ][OAc - ]

Reactive molecular dynamics simulations enable detailed understanding of solvent effects on chemical reaction mechanisms and reaction rates. While classical molecular dynamics using reactive force fields allows significantly longer simulation time scales and larger system sizes compared with ab initio molecular dynamics, constructing reactive force fields is a difficult and complex task. In this work, we describe a general approach following the Empirical Valence Bond (EVB) framework for constructing ab initio reactive force fields for condensed phase simulations by combining physics-based methods with neural networks (PB/NN). The physics-based terms ensure correct asymptotic behavior of electrostatic, polarization, and dispersion interactions, and are compatible with existing solvent force fields. Neural networks are utilized for versatile description of short-range orbital interactions within the transition state region, and accurate rendering of vibrational motion of the reacting complex. Herein, we demonstrate our methodology for a simple deprotonation reaction of the 1-ethyl-3-methylimidazolium (EMIM+) cation with acetate to form 1-ethyl- 3-methylimidazol-2-ylidene and acetic acid. Our PB/NN force field exhibits ~ 1 kJ/mol MAE accuracy within the transition state region for the gas-phase complex. To characterize solvent modulation of the reaction profile, we compute potentials of mean force (PMFs) for the gas-phase reaction as well as the reaction within a four ion cluster, and benchmark against ab initio molecular dynamics simulations. We find that the surrounding ionic environment significantly destabilizes formation of the carbene product, and we show that this effect is accurately captured by the reactive force field. By construction, the PB/NN potential may be directly employed for simulations of other solvents/chemical environments without additional parameterization.

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

Structural features of xylan dictate reactivity and functionalization potential for bio-based materials

Plant-based materials have the potential to replace some petroleum-based products, offering compostability and biodegradability as critical advantages. Xylan-rich biomass sources are gaining recognition due to their abundance and underutilization in current industrial applications. Research of potential xylan applications has been complicated by the complex and heterogeneous structure that varies for different xylan feedstocks. Acylation is a broadly used reaction in functionalization of polysaccharides at an industrial scale. However, the efficiency of this reaction varies with the xylan source. To optimize xylan valorization, a systematic understanding of structure–reactivity relationships is essential. This study explores, characterizes, and compares various xylan feedstocks in the acylation process. Xylan feedstocks were analyzed for their chemical composition, degree of polymerization, branching, solubility, and presence of impurities. These features were correlated with xylan glycotypes’ reactivity toward functionalization with succinic anhydride in an optimized DMSO/KOH condition, achieving carboxyl contents of up to 1.46. We used principal component analysis and hierarchical clustering to identify key structural features of xylan that promote its reactivity. Our findings reveal that xylans with higher xylose content and lower degrees of branching exhibit enhanced reactivity, achieving higher carboxyl content and yields. Structural analyses confirmed successful modification, and light scattering analyses showed dramatic changes in the solution properties. Succinylation improves the solubility and film-forming properties of native xylans. This study shows key structure–reactivity relationships in xylan succinylation, establishing that low branching, high xylose content, and reduced lignin impurity enhance chemical functionalization. The results offer a framework for selecting optimal biomass feedstocks and support future efforts in genetic and synthetic biology to design plants with tunable xylan architectures. These findings advance the hemicellulose valorization for applications in coatings and packaging.

Acylation↗
DOE OSTI2022

The Contribution of Proton-Donor pKa on Reactivity Profiles of [FeFe]-hydrogenases

The [FeFe]-hydrogenases are enzymes that catalyze the reversible activation of H 2 coupled to the reduction–oxidation of electron carriers. Members of the different taxonomic groups of [FeFe]-hydrogenases display a wide range of preference, or bias, for H 2 oxidation or H 2 production reactions, despite sharing a common catalytic cofactor, or H-cluster. Identifying the properties that control reactivity remains an active area of investigation, and models have emerged that include diversity in the catalytic site coordination environments and compositions of electron transfer chains. The kinetics of proton-coupled electron transfer at the H-cluster might be expected to be a point of control of reactivity. To test this hypothesis, systematic changes were made to the conserved cysteine residue that functions in proton exchange with the H-cluster in the three model enzymes: CaI, CpII, and CrHydA1. CaI and CpII both employ electron transfer accessory clusters but differ in bias, whereas CrHydA1 lacks accessory clusters having only the H-cluster. Changing from cysteine to either serine (more basic) or aspartate (more acidic) modifies the sidechain pKa and thus the barrier for the proton exchange step. The reaction rates for H 2 oxidation or H 2 evolution were surveyed and measured for model [FeFe]-hydrogenases, and the results show that the initial proton-transfer step in [FeFe]-hydrogenase is tightly coupled to the control of reactivity; a change from cysteine to more basic serine favored H 2 oxidation in all enzymes, whereas a change to more acidic aspartate caused a shift in preference toward H 2 evolution. Overall, the changes in reactivity profiles were profound, spanning 10 5 in ratio of the H 2 oxidation-to-H 2 evolution rates. The fact that the change in reactivity follows a common trend implies that the effect of changing the proton-transfer residue pKa may also be framed as an effect on the scaling relationship between the H-cluster di(thiolmethyl)amine (DTMA) ligand pKa and E m values of the H-cluster. Experimental observations that support this relationship, and how it relates to catalytic function in [FeFe]-hydrogenases, are discussed.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Single-Well Push–Pull Tracer Test Analyses to Determine Aquifer Reactive Transport Parameters at a Former Uranium Mill Site (Grand Junction, Colorado)

At a former uranium mill site where tailings have been removed, prior work has determined several potential ongoing secondary uranium sources. These include locations with uranium sorbed to organic carbon, uranium in the unsaturated zone, and uranium associated with the presence of gypsum. To better understand uranium mobility controls at the site, four single-well push–pull tests (with a drift phase) were completed with the goal of deriving aquifer flow and contaminant transport parameters for inclusion in a future sitewide reactive transport model. This goes beyond the traditional use of a constant sorption distribution coefficient (K d ) and allows for the evaluation of alternative remedial injection fluids, which can produce variable K d values. Dispersion was first removed from the resulting data to determine possible reactions before conducting reactive transport simulations. These initial analyses indicated the potential need to include cation exchange, uranium sorption, and gypsum dissolution. A reactive transport model using multiple layers to account for partially penetrating wells was completed using the PHT-USG reactive transport modeling code and calibrated using PEST. The model results quantify the hydraulic conductivity and dispersion parameters using the injected tracer concentrations. Uranium sorption, cation exchange, and gypsum dissolution parameters were quantified by comparing the simulated versus observed geochemistry. All simulations required some cation exchange and calcite equilibrium, and one simulation required gypsum dissolution to improve the model fit for calcium and sulfate. Uranium sorption parameters were not strongly influenced by the other parameter values but were highly influenced by uranium concentrations during the drift phase, with possible kinetic rate limitations. Thus, a future recommendation for such push–pull tests is to collect more geochemical data during the drift phase. The final uranium sorption parameters were within the range of values determined from prior column testing. The flow and transport parameters derived from these single-well push–pull tests will provide initial parameters for any future sitewide reactive transport model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Modifying the Reactivity of Single Pd Sites in a Trimetallic Sn‐Pd‐Ag Surface Alloy: Tuning CO Binding Strength

Abstract Improving control over active‐site reactivity is a grand challenge in catalysis. Single‐atom alloys (SAAs) consisting of a reactive component doped as single atoms into a more inert host metal feature localized and well‐defined active sites, but fine tuning their properties is challenging. Here, a framework is developed for tuning single‐atom site reactivity by alloying in an additional inert metal, which this work terms an alloy‐host SAA. Specifically, this work creates about 5% Pd single‐atom sites in a Pd 33 Ag 67 (111) single crystal surface, and then identifies Sn based on computational screening as a suitable third metal to introduce. Subsequent experimental studies show that introducing Sn indeed modifies the electronic structure and chemical reactivity (measured by CO desorption energies) of the Pd sites. The modifications to both the electronic structure and the CO adsorption energies are in close agreement with the calculations. These results indicate that the use of an alloy host environment to modify the reactivity of single‐atom sites can allow fine‐tuning of catalytic performance and boost resistance against strong‐binding adsorbates such as CO.

Mohrhusen, Lars↗
DOE OSTI2025

Corrosion Resistance of an AlCeMg/Stainless-Steel Reactive Bond

A major issue for metal components in many industries is corrosion as it can substantially reduce their lifetime. This issue is especially problematic for materials used in heat exchanger applications. Al–Ce–Mg alloys, which exhibit corrosion resistance and can reactively bond with other metals, may be a viable solution to this problem. This investigation studied the corrosion behavior of Al–2Ce–6Mg (atomic percent)/stainless-steel (SS) reactive bond interfaces after full immersion in nitric, sulfuric, formic, and mixed acid for 267 h. This particular Al–Ce–Mg alloy was chosen due to its good castability. The results of scanning electron microscope characterization showed that reactive bond formations repeatedly occurred throughout the length of the casting in the as-cut samples and that these formations maintained a secure bonding between the alloy and the stainless-steel tubes. Transmission electron microscopy results showed that there was a clear compositional and microstructural transition across the reactive bond. The results of the immersion tests indicated that the nitric, sulfuric, and the mixed acid did not have an observably negative effect on the reactive bond structure. As for the sample exposed to formic acid only, noticable changes were seen in both the microstructural appearance and the elemental profile across the bond, suggesting that oxide formation occurred.

Brechtl, Jamieson [ORNL] (ORCID:0000000217394283)↗
DOE OSTI2021

Role of surface diffusion in formation of unique reactivity for graphite oxidation: Time-resolved measurements in a pulsed diffusion reactor

Quantification of oxidation kinetics is essential to develop graphitic materials for diverse applications: from refractories found in gas-cooled nuclear reactors to catalysts needed for chemical manufacturing. In this work, using well-defined highly oriented pyrolytic graphite, low-pressure isotopic transient experiments combined with controlled annealing periods, we resolve the role of surface diffusion and quantify oxidation kinetics with nanomole-precision. We observe an unexpected increase in reactivity following annealing which is explained by the role of surface diffusion increasing the probability for trapping mobile oxygen at more reactive edge sites. Here, the locus of adsorption and spillover to the basal plane is distinct from the trapping location creating a more active oxygen species. Isotopic products reflect the population dynamics of oxygen added at the edge and surface diffusion that relocates basal plane oxygen to more reactive edge sites. Since this process proceeds in parallel with direct oxidation reactions, it is not likely to be observed using steady-state or conventional ‘bulk’ characterization techniques. Our unique time-resolved non-equilibrium measurement in a well-defined transport regime, enables observation of three distinct behaviors: short-term deactivation due to the balance of rates in oxygen supply/product formation, reactivity increases due to surface diffusion and longer-term reactivity increase with oxygen accumulation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Pozzolanic reactivity of natural pozzolans, ground glasses and coal bottom ashes and implication of their incorporation on the chloride permeability of concrete

It is well understood that the partial replacement of portland cement with pozzolans in sufficient proportions results in improved long-term performance of concrete. However, dwindling availability of the most widely used pozzolan in the construction industry, fly ash, is forecast due to the termination of coal-fired electricity power stations around the world, necessitating the exploration of alternative pozzolans. Experimental investigations are conducted to study the performance of traditional pozzolans such as fly ash, silica fume and natural pozzolans (metakaolin, pumice, perlite and lassenite) together with novel pozzolans, including ground glasses (high-alkali and low-alkali), and ground bottom ash. Twenty-one concrete mixtures are examined for compressive strength, electrical conductivity, electrical resistivity, chloride permeability, and chloride migration coefficient. The reactivity of pozzolans studied is also determined in terms of compressive strength in lime-pozzolan mortar using a modified lime-reactivity test method. Many of the materials tested demonstrate the potential to be used in place of fly ash. However, a wide variation in the performance of these materials is evident which highlights the need for a reliable test to determine the level of reactivity of a pozzolan. The results indicate that the 7-day strength in the modified lime-reactivity test provides a good indication of the pozzolanic reactivity of the material and of how the material can be expected to contribute to the strength and permeability of concrete. The use of electrical resistivity as an indicator of the performance of pozzolans in terms of strength development and chloride penetration resistance when used to partially replace portland cement in concrete is also discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effect of the ..beta..-Hydroxy Group on Ester Reactivity: Combustion Kinetics of Methyl Hexanoate and Methyl 3-Hydroxyhexanoate

One of the major biofuels used today is biodiesel, composed of fatty acid methyl esters. The fats and oils feedstocks used to make biodiesel are in relatively short supply. To meet increasing global demand, engineered microorganisms have been developed to catalyze conversion of sugar to fatty acid esters that contain unique ..beta..-hydroxy-esters. This study investigated the effect of a hydroxyl group on the combustion characteristics of methyl esters by comparing the chemical behavior of methyl hexanoate (MHx) and methyl 3-hydroxyhexanoate (M3OHHx) - used as surrogates for diesel-boiling-range esters with longer fatty acid chains. The oxidation of these esters was studied experimentally in a flow reactor at 0.84 and 10 bar, 600 to 1,100 K, and stoichiometric conditions; and in a constant volume combustion chamber at 5 and 10 bar, 600 to 900 K, and equivalence ratios of 0.3 and 0.6. MHx was more reactive in the constant volume chamber at temperatures below 800 K at 10 bar and equivalence ratio of 0.6. MHx also exhibited a higher indicated cetane number (16.4) than M3OHHx (8.1). We investigated this reactivity trend using an updated MHx kinetic model and M3OHHx model developed as part of this work. The kinetic models predicted that radicals formed from MHx were consumed by low-temperature chemistry reactions, whereas the M3OHHx reactivity was significantly governed by the ..beta..-radical (on the same carbon as the OH group) chemistry which mainly terminated via the less reactive chain propagation pathway to methyl-3-oxohexanoate + HO2. This study extends our understanding of structural effects on ester reactivity, which will allow for accurate surrogate formulation and performance simulations.

AFIDA↗
DOE OSTI2025

Improving Bond Dissociations of Reactive Machine Learning Potentials through Physics-Constrained Data Augmentation

In the field of computational chemistry, predicting bond dissociation energies (BDEs) presents well-known challenges, particularly due to the multireference character of reactive systems. Many chemical reactions involve configurations where single-reference methods fall short, as the electronic structure can significantly change during bond breaking. As generating training data for partially broken bonds is a challenging task, even state-of-the-art reactive machine learning interatomic potentials (MLIPs) often fail to predict reliable BDEs and smooth dissociation curves. By contrast, simple and inexpensive physics-based models, such as the well-established Morse potential, do not suffer from any such limitations. This work leverages the Morse potential to improve reactive MLIPs by augmenting the training data set with inexpensive Morse data along the dissociation pathways. Further, this physics-constrained data augmentation (PCDA) approach results in MLIPs with smooth bond dissociation curves as well as near coupled-cluster level BDEs, all without requiring any expensive multireference quantum mechanical calculations. A case study for methane combustion demonstrates how the PCDA approach can improve an existing reactive MLIP, namely, ANI-1xnr. In conclusion, not only are the BDEs and bond dissociation curves for all radicals and molecules significantly improved compared to ANI-1xnr but the PCDA-trained MLIP retains the reliability of ANI-1xnr when performing reactive molecular dynamics simulations.

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

Exploring the frontiers of condensed-phase chemistry with a general reactive machine learning potential

Abstract Atomistic simulation has a broad range of applications from drug design to materials discovery. Machine learning interatomic potentials (MLIPs) have become an efficient alternative to computationally expensive ab initio simulations. For this reason, chemistry and materials science would greatly benefit from a general reactive MLIP, that is, an MLIP that is applicable to a broad range of reactive chemistry without the need for refitting. Here we develop a general reactive MLIP (ANI-1xnr) through automated sampling of condensed-phase reactions. ANI-1xnr is then applied to study five distinct systems: carbon solid-phase nucleation, graphene ring formation from acetylene, biofuel additives, combustion of methane and the spontaneous formation of glycine from early earth small molecules. In all studies, ANI-1xnr closely matches experiment (when available) and/or previous studies using traditional model chemistry methods. As such, ANI-1xnr proves to be a highly general reactive MLIP for C, H, N and O elements in the condensed phase, enabling high-throughput in silico reactive chemistry experimentation.

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

Simulation of the Fast Reactor Fuel Assembly Duct-Bowing Reactivity Effect Using Monte Carlo Neutron Transport and Finite Element Analysis

This paper discusses a new method of simulating the fuel assembly duct-bowing reactivity coefficient for EBR-II run 138B. Quantification of the fuel assembly duct-bowing reactivity effect in liquid metal–cooled fast reactors has been a persistent problem since they were first designed and operated. Simulation of the duct-bowing reactivity effect is difficult because the level of detail required to simulate the effect has exceeded most modeling capabilities. The new method outlined in this paper utilizes the finite element analysis code ANSYS to analyze the thermal and structural components. Here, the displacement of the fuel assembly duct due to thermal expansion and mechanical interaction was calculated by ANSYS using recorded EBR-II run 138B temperature and power boundary value data. The displacement values were incorporated into to a Monte Carlo model of EBR-II run 138B and keff was calculated. Multiple Monte Carlo calculations were performed with duct displacement values corresponding to different reactor temperatures. Using the calculated keff values associated with the different duct displacement results allowed calculation of the duct-bowing reactivity coefficient. The duct-bowing reactivity coefficient was calculated to be –14.5 × 10 –4 $/°C/ ± 4.4%.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗