Search NASA⌕ Search

SEARCH · Search NASA

Results for “equilibrium complex concentrations”

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 OSTI2025

Modeling phase transformations in Mn-rich disordered rocksalt cathodes with machine-learning interatomic potentials

Mn-rich disordered rocksalt (DRX) cathode materials exhibit a phase transformation from a disordered to a partially disordered spinel-like structure (δ-phase) during electrochemical cycling. Here, in this computational study, we use charge-informed molecular dynamics with a fine-tuned CHGNet foundation potential to investigate the phase transformation in LixMn 0.8 Ti 0.1 O 1.9 F 0.1 . Our results indicate that transition metal migration occurs and reorders to form the spinel-like ordering in an FCC anion framework. The transformed structure contains a higher concentration of nontransition metal (0-TM) face-sharing channels, which are known to improve Li transport kinetics. Analysis of the Mn valence distribution suggests that the appearance of tetrahedral Mn 2+ is a consequence of spinel-like ordering, rather than the trigger for cation migration as previously suggested. Calculated equilibrium intercalation voltage profiles demonstrate that the δ-phase, unlike the ordered spinel, exhibits solid-solution signatures at low voltage. A higher Li capacity is obtained than in the DRX phase. This study provides atomic insights into solid-state phase transformation and its relation to experimental electrochemistry, highlighting the potential of machine-learning interatomic potentials for understanding complex oxide materials.

Zhong, Peichen [University of California, Berkeley↗
DOE OSTI2021

Progress Report on Performance of A709 and G91 Steels in Sodium

Specimens in six different processing and heat treatment conditions of A709 H58776 were exposed to sodium at 550, 600, and 650°C, respectively for various exposure times. G91 base metal was tested in sodium at 550, 600, and 650°C, and G91 weldment at 550 and 600°C, respectively. Parallel thermal aging experiments were conducted on G91 and A709 to obtain thermal aging data for comparison with sodium exposure results to separate the thermal and sodium effects. The corrosion data obtained on A709 steel continue to show low corrosion rates, which indicate good compatibility of A709 with sodium when oxygen content is controlled. G91 has also shown acceptable corrosion rates over the temperature range of 550-650°C in sodium. Thermal aging or sodium exposures reduced the tensile strength, uniform elongation, and total elongation for A709 ESR, AOD specimens. These specimens consistently showed higher yield stress and tensile strength and lower uniform and total elongations after sodium exposure than after thermal aging under comparable testing conditions. It implies that there was an additional effect of sodium exposure that gave rise to an increase in tensile strength of A709. Dynamic strain aginginduced flow serrations were completely removed by either thermal aging or sodium exposure at 650°C, which imply that carbon or nitrogen in solution was removed from the solution and formed precipitates during aging or sodium exposures. The tensile data ruled out the possibility of decarburization under the sodium exposure conditions. The HOMO specimens behaved somewhat differently from the ESR and AOD specimens in sodium. While thermal aging and sodium exposures at 550 and 600°C decreased the yield stress and the ultimate tensile strength of G91, there was virtually no additional effect resulting from sodium exposures at these two temperatures. In contrast, sodium exposures at 650°C had a drastic effect on the yield stress and the ultimate tensile strength of G91. It is suggested that G91 experienced decarburization in the 650°C sodium environments. Carbon concentrations in sodium in the SMT-1 and SMT-2 loops were determined by a foil equilibration method. The estimated carbon concentration was in the range of 0.8-1.2 ppm in the SMT-1 loop and 0.3-0.7 ppm in the SMT-2 loop. The carbon activity in sodium in the SMT-2 loop was estimated to be 0.03-0.08 at 600°C, and 0.08-0.2 at 550°C. The carbon activity in sodium at 650°C in the SMT-1 loop was 0.04-0.07. Equilibrium simulation of the carburization – decarburization processes was conducted for A709 and G91 steels exposed in sodium environments at temperatures of 550-700°C. The carbon activity-concentration relationship for G91 was re-evaluated by considering four phases in G91, i.e. bcc ferrite, M 23 C 6 , NbC and VC carbides. It was found that the carburization-decarburization process in G91 steel was dictated by M 23 C 6 carbides at high carbon activities, while NbC carbides dominated the process at low carbon activities. Formation of NbC increases the decarburization resistance of G91 steel. It remains to be understood whether the decarburization resistance provided by MC carbides in G91 can be maintained during long-term operations of SFRs. The carbon concentration-activity relationship for A709 was calculated based on the equilibrium of fcc-austenite and M 23 C 6 carbide phase. The calculations showed that A709 decarburizes at 650°C in the SMT-1 loop environments, which is different from the experimental findings. The carbon concentration-activity relationship in A709 was further evaluated by considering four phases, including fcc-austenite, M 23 C 6 , NbC, and TiC carbide phases. TiC is the most stable carbides among the three carbide phases. The carburization-decarburization process in A709 is dictated by M 23 C 6 carbides at high carbon activities and by TiC at low carbon activities. Formation of TiC increases the decarburization resistance of A709 steel. Because of the complexity of the precipitation process in A709 during sodium exposures, detailed characterization of precipitates in sodium-exposed specimens and improved thermodynamic models are needed to understand the carburization-decarburization behavior of A709. Kinetic analysis of carbon transfer will be investigated for G91 and A709 in future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Insights into Sorption and Molecular Transport of Aqueous Glucose into Zeolite Nanopores

Liquid-phase heterogeneous catalysis using zeolites is important for biomass conversion to fuels and chemicals. There is a substantial body of work on gas-phase sorption in zeolites with different topologies; however, studies investigating the diffusion of complex molecules in liquid medium into zeolitic nanopores are scarce. Here, we present a molecular dynamics study to understand the sorption and diffusion of aqueous β-d-glucose into β-zeolite silicate at T = 395 K and P = 1 bar. Through 2-μs-long molecular dynamics trajectories, we reveal the role of the solvent, the kinetics of the pore filling, and the effect of the water model on these properties. We find that the glucose and water loading is a function of the initial glucose concentration. Although the glucose concentration increases monotonically with the initial glucose concentration, the water loading exhibits a nonmonotonic behavior. At the highest initial concentration (~20 wt %), we find that the equilibrium loading of glucose is approximately five molecules per unit cell and displays a weak dependence on the water model. Glucose molecules follow a single-file diffusion in the nanopores due to confinement. Here, the dynamics of glucose and water molecules slows significantly at the interface. The average residence time for glucose molecules is an order of magnitude larger than that in the bulk solution, while it is about twice as large for the water molecules. Our simulations reveal critical molecular details of the glucose molecule’s local environment inside the zeolite pore relevant to catalytic conversion of biomass to valuable chemicals.

09 BIOMASS FUELS↗
DOE OSTI2020

Collaborative Research: Natural Organic Matter and Microbial Controls on Mobilization/Immobilization of I and Pu in Soils and Waters Affected by Radionuclide Releases in USA and Japan

In this project, the relationship between natural organic matter (NOM) and two radioactive elements that are relevant to nuclear waste disposal were studied: Plutonium (Pu) and Iodine. The human and environmental risks associated with Pu stem mainly from the very long half-lives of several of its isotopes ( 238 Pu, 88 yr; 239 Pu, 24,100 yr; 240 Pu, 6560 yr) and its radiotoxicity. Understanding Pu biogeochemical behavior in both near-field (>10-11M) and far-field scales (<10-11M) is imperative to the development of approaches for reprocessing Pu, remediation of Pu contamination and accurate assessment of risks posed by disposal practices for Pu-bearing wastes. The environmental mobility of Pu can be affected by redox potential, pH, adsorption, precipitation, complexation, colloid formation, and microbial activity, of which the first characteristic has the most profound influence. Numerous studies have shown high affinity of Pu towards NOM, as well as to mineral phases. NOM is ubiquitous in the environment, e.g., both fulvic and humic acids are able to reduce Pu(V,VI) to Pu(IV) and the redox potential of NOM is positively related to the abundance of phenolic/acidic OH groups. NOM can either facilitate or limit actinide migration, depending on specific biogeochemical conditions including pH, mineral and organic matter characteristics, etc. The other radionuclide of interest is radioiodine ( 129 I). 129I is a major by-product of nuclear fission and of serious concern to the Department of Energy (DOE) as it is among the top risk drivers at existing and potential radiowaste-contaminated sites. The risk of 129 I stems largely from its high bioconcentration factor (90% of the body’s iodine is accumulated in the thyroid), a high inventory at source terms, a very long-half life (16M years), and rapid mobility in the subsurface environment. As a consequence, 129 I has the lowest drinking water standard (1 pCi/L) among all radionuclides in the Federal Register. With a novel and sensitive gas chromatography-mass spectrometry (GC-MS) method developed in our lab, it is possible to quickly and simultaneously determine the distribution of 129 I and stable 127 I forms in environments, as low as 2 pCi/L for 129 I. This method was subsequently validated using accelerator mass spectroscopy, AMS. IO 3 - and organo-I were determined as major species in the groundwater of SRS and the Hanford Site, contrary to thermodynamic predictions that I- should be the dominant species at these sites. Mobility of 129 I was also demonstrated to depend greatly on the I species and its concentration, sediment pH, and redox state, with times to achieve equilibrium taking up to 12 weeks. Along the groundwater pathway in the F-Area of SRS, 129 I- supplied from the seepage basins was transformed to 129 IO 3 - and organo- 129 I with increasing iodine sediment sorption, causing the lower total 127 I and 129 I concentrations along the gradient transect of the waste plume. By contrast, groundwater 129 I concentrations in the wetlands (as high as 1617.3 pCi/L) were greatly elevated with respect to the source term (159.3 pCi/L). While the NOM promoted the uptake of 129 I to the wetland sediment, it also promoted the formation of soluble organic fraction. A small fraction of NOM that is bound to iodine can behave as a mobile organo-I source. Iodide was enzymatically incorporated into NOM, whereas both iodide and iodate were abiotically bound to NOM, under certain conditions. Iodate removal from the mobile aqueous phase can also occur through incorporation into carbonate (e.g., at the Hanford Site, USA). Thus immobilization and re-mobilization of iodine species were influenced by pH, Eh and the presence of NOM and metal oxides, which adds to the complexity of site remediation action. A ground-breaking result was to elucidate the products (i.e. organo-iodine moieties formed via enzymatic and non-enzymatic processes) at the molecular level by nuclear magnetic resonance (NMR) and electrospray ionization Fourier transform ion cyclotron resonance mass spectrometry (ESI-FTICRMS). We found that iodine-NOM interactions may be influenced by NOM hydrophobic aliphatic moieties. From the perspective of ESI-FTICRMS, organo-iodine formulas were ascribed to the groups of unsaturated hydrocarbons, lignins and proteins. Iodate is likely abiotically reduced to reactive iodine species by lignin- and tannin-like compounds or carboxylic-rich alicyclic molecules (CRAM). We also investigated microbial mechanisms in iodine incorporation into NOM. We established that soil bacteria isolated from F-Area of SRS did not accumulate significant amounts of I- (0.2-2%). Intracellular uptake of I- decreases with increasing pH when pH ranged from 4 to 6. In contrast, 44 out of 84 strains isolated from the F-Area of SRS can transform I- to IO 3 - and organo-iodine. In some cases, oxidation was facilitated in the presence of H 2 O 2 . Microbes can also excrete organic acids that enhance I- oxidation by lowering the ambient pH and reacting with H 2 O 2 to form peroxy carboxylic acids. At lower pH values (≤5), H 2 O 2 hydrolysis was the driving force for iodide-oxidation; whereas, at pH ≥ 6, spontaneous decomposition of peroxy carboxylic acids, originating from H 2 O 2 and organic acids were the primary cause of iodide oxidation. Lastly, it was determined that microbial processes involved in Mn (II) are capable of directly oxidizing I- via enzymatic catalysis (i.e., multicopper oxidases), or indirectly through the formation of reactive oxygen species (ROS) and/or biogenic manganese oxides. ROS-mediated oxidation of I- was found to predominate at pH >5, whereas the enzymatic and Mn oxide pathways were more active at pH < 5. Together, this project has resulted in 9 publications in high-impact journals, and the training of 1 Ph.D and 4 undergraduate students.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Spatially Patterned Architectures to Modulate CO 2 Reduction Cascade Catalysis Kinetics

Electrochemical CO 2 reduction using renewable sources of electrical energy holds promise for converting CO 2 into fuels and chemicals. The complex interactions among chemical/electrochemical reactions and mass transport make it difficult to analyze the effect of an individual process on electrode performance based only on experimental methods. Here, we developed a generalized steady-state simulation to describe an electrode surface in which sequential cascade catalysts are patterned in a periodic trench design. If appropriately constructed, this trench geometry is hypothesized to be able to yield a higher net current density for a CO 2 reduction (CO 2 R) cascade reaction. We have used realistic experimental reaction kinetics to investigate the role of trench geometry in mass transport, local microenvironments, and selectivity for a model CO 2 R cascade reaction. The model considers local concentration gradients of bicarbonate species at quasi-equilibrium and catalytic surface reactions based on concentration-dependent Butler–Volmer kinetics. Our results suggest that varying the spatial distribution of active sites plays a significant role in facilitating effective mass transport between active sites, modulating selectivity for the cascade reaction, and enhancing the yield of desirable cascade products. Moreover, we observe that this trench geometry significantly alters the cascade reaction rate by affecting the local pH, which can cause inadvertent depletion of available aqueous CO 2 to limit the CO 2 R cascade kinetics and modest suppression of the hydrogen evolution reaction (HER). The results highlight the trade-offs between mass transport, pH, and reaction kinetics that become apparent only when considering the coupled physics of all processes at the electrode surface. Here, this model can thus serve as a primary tool to build more selective and efficient patterned architectures for the CO 2 R cascade catalysis.

CO2 reduction↗
DOE OSTI2022

Master equation study of hydrogen abstraction from HCHO by OH via a chemically activated intermediate

The abstraction reaction of hydrogen from formaldehyde by OH radical plays an important role in formaldehyde oxidation. The reaction involves a bimolecular association to form a chemically activated hydrogen-bonded reaction complex followed by a unimolecular reaction of the complex to generate the products. The reaction rate is usually considered to be pressure-independent by assuming equilibrium between the reactants and the complex. However, our nonequilibrium calculations based on the chemically significant eigenmode of the master equation, carried out with our recently developed TUMME program, indicate that the reaction complex makes the rate constant dependent on pressure at low temperatures (T < 200 K). The calculations include anharmonicity, variational effects, and multi-dimensional tunneling. We find that the reaction rate constant reaches a low-pressure limit at pressures below 10 Torr over the whole investigated temperature range (20–1800 K), which explains why the available low-temperature experiments, which are for pressures below 2 Torr, did not observe the pressure dependence. A new extension of the TUMME master-equation program is used to explore the time evolutions of the concentrations of the OH radical and the complex under pseudo-first-order conditions. The time-dependent evolution of the concentrations of the complex at a low temperature provide direct evidence for the stabilization of the reaction complex at high pressures, and it shows the negligible role of the stabilized reaction complex at low pressures. Finally, the picture that emerges is qualitatively consistent with our previous study of the reaction of methanol with OH in that the tunneling in the unimolecular step from the complex to the products affects the phenomenological reaction rate constants differently at high and low pressures and leads to a significant pressure effect.

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

Modeling diffusion and types I-V sorption of water vapor in heterogeneous systems

Here, a sorption and diffusion model is developed using a mobile-immobile decomposition of the vapor concentration. Henry’s mode represents the mobile species which diffuses through the solid, while Langmuir and pooling modes represent possible immobilization processes that cause nonlinearities in sorption capacity. The model can simulate the five classical sorption types and can incorporate additional dynamics compared to models based on thermodynamic equilibrium. Additionally, the framework allows for discontinuities in material properties encountered in multi-material systems. The resulting non-linear coupled equations are solved by employing a finite element method to discretize in space, and a backward Euler method to discretize in time. The discretized system of equations is solved via a Newton-type iteration scheme at each time-step. Four different materials are parameterized for the model, which are then used in examples to demonstrate the model’s ability to capture complex sorption processes in multi-material systems.

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

Bubble Mass Transport Measurement in the Large-Scale Test Loops at Oak Ridge National Laboratory

Molten salts are complex fluids that incorporate multi-phase behavior depending on the chemistry and the physical properties of the entrained components. These components include the carrier salt, the actinide fuel, and fission, activation, and corrosion products, the concentrations of which depends on the burnup history of the salt. Radionuclide transport from the salt into the cover gas / off-gas system depends on volatility as predicted by thermochemistry, but data from the Molten Salt Reactor Experiment (MSRE) conducted in the late 1960s suggest that bubble formation and transport are also important. Anomalously high fractions of noble metals were found in the off-gas system and were attributed to transportation with parent salt aerosols and their association with rising noble gas bubbles. The prediction of such phenomena requires coupled neutronic, thermal hydraulic, and chemical equilibrium calculations, the framework of which is being developed within the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program supported by the US Department of Energy (DOE) Office of Nuclear Energy (NE). However, separate effects tests and models of experiments using tools such as SAM and Thermochimica can guide model developers through important processes. Conversely, model development guides the choice of experiments and systems to provide data that are relevant for validation. This report describes several experiments that tracked gas transport in molten salts, ranging from small-scale systems to large-scale loops. Gas transportation in a molten salt, LiCl-KCl, has been studied using the shadowgraph technique. Sensors such as residual gas analysis, Raman spectroscopy, and laser-induced breakdown spectroscopy have been tested for off-gas measurements. These data were used to interpret how the gases move through the upstream salt / cover gas and the interface between them. Differential pressure measurements were able to detect individual gas bubbles as they popped at the liquid–gas interface. Salt aerosols were collected on a cascade impactor. Their formation was also observed directly via the shadowgraph method, and most of these aerosols were launched ballistically into the cover gas. Fine mists could also be observed. Convection currents through the salt were visualized and can be used to calculate the thermophysical properties of the salt itself. The apparatuses described in this report and in a previous work (McFarlane et al. 2023) have been commissioned and are available for use in making further measurements of salt/surrogate fission product behavior. Novel approaches using neutron imaging are planned for the study of fluoride salts, which cannot be contained in quartz, so shadowgraph visualization is not available. Bubble transport in convective flow and in a slow-moving salt column are planned. The mobile laser-induced breakdown spectroscopy (LIBS) system is available for several applications, including iodine capture in a molten hydroxide scrubber, H2 transport though molten salts to complement Raman analysis, and online tracking of salt aerosol generation, transport, and deposition. This report summarizes the findings from FY24 and the plans for FY25. The work completes milestone M2AT-24OR0702013 of the DOE-NE Advanced Reactor Technology, Molten Salt Reactor Campaign, DOE-NE-5.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Unravelling Magnetic Nanochain Formation in Dispersion for In Vivo Applications

Self-assembly of iron oxide nanoparticles (IONPs) into 1D chains is appealing, because of their biocompatibility and higher mobility compared to 2D/3D assemblies while traversing the circulatory passages and blood vessels for in vivo biomedical applications. In this work, parameters such as size, concentration, composition, and magnetic field, responsible for chain formation of IONPs in a dispersion as opposed to spatially confining substrates, are examined. In particular, the monodisperse 27 nm IONPs synthesized by an extended LaMer mechanism are shown to form chains at 4 mT, which are lengthened with applied field reaching 270 nm at 2.2 T. The chain lengths are completely reversible in field. Using a combination of scattering methods and reverse Monte Carlo simulations the formation of chains is directly visualized. The visualization of real-space IONPs assemblies formed in dispersions presents a novel tool for biomedical researchers. This allows for rapid exploration of the behavior of IONPs in solution in a broad parameter space and unambiguous extraction of the parameters of the equilibrium structures. Additionally, it can be extended to study novel assemblies formed by more complex geometries of IONPs.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Comprehensive characterization of irradiation induced defects in ceria: Impact of point defects on vibrational and optical properties

Validation of multiscale microstructure evolution models can be improved when standard microstructure characterization tools are coupled with methods sensitive to individual point defects. We demonstrate how electronic and vibrational properties of defects revealed by optical absorption and Raman spectroscopies can be used to compliment transmission electron microscopy (TEM) and x-ray diffraction (XRD) in the characterization of microstructure evolution in ceria under non-equilibrium conditions. Experimental manifestation of non-equilibrium conditions was realized by exposing cerium dioxide (CeO2) to energetic protons at elevated temperature. Two sintered polycrystalline CeO2 samples were bombarded with protons accelerated to a few MeVs. These irradiation conditions produced a microstructure with resolvable extended defects and a significant concentration of point defects. A rate theory (RT) model was parametrized using the results of TEM, XRD, and thermal conductivity measurements to infer point defect concentrations. An abundance of cerium sublattice defects suggested by the RT model is supported by Raman spectroscopy measurements, which show peak shift and broadening of the intrinsic T2g peak and emergence of new defect peaks. Additionally, spectroscopic ellipsometry measurements performed in lieu of optical absorption reveals the presence of Ce3+ ions associated with oxygen vacancies. This work lays the foundation for a coupled approach that considers a multimodal characterization of microstructures to guide and validate complex defect evolution models.

Chauhan, Vinay S.↗
DOE OSTI2023

Preclinical studies of a PARP targeted, Meitner-Auger emitting, theranostic radiopharmaceutical for metastatic ovarian cancer

Advanced ovarian cancer currently has few therapeutic options. Poly(ADP-ribose) polymerase (PARP) inhibitors bind to nuclear PARP and trap the protein-inhibitor complex to DNA. This work investigates a theranostic PARP inhibitor for targeted radiopharmaceutical therapy of ovarian cancer in vitro and PET imaging of healthy mice in vivo. Methods: [ 77 Br]RD1 was synthesized and assessed for pharmacokinetics and cytotoxicity in human and murine ovarian cancer cell lines. [ 76 Br]RD1 biodistribution and organ uptake in healthy mice were quantified through longitudinal PET/CT imaging and ex vivo radioactivity measurements. Organ-level dosimetry following [ 76/77 Br]RD1 administration was calculated using RAPID, an in-house platform for absorbed dose in mice, and OLINDA for equivalent and effective dose in human. Results: The maximum specific binding (B max ), equilibrium dissociation constant (K d ), and nonspecific binding slope (NS) were calculated for each cell line. These values were used to calculate the cell specific activity uptake for cell viability studies. The half maximal effective concentration (EC 50 ) was measured as 0.17 (95 % CI: 0.13–0.24) nM and 0.46 (0.13–0.24) nM for PARP(+) and PARP(–) expressing cell lines, respectively. The EC 50 was 0.27 (0.21–0.36) nM and 0.30 (0.22–0.41) nM for BRCA1(–) and BRCA1(+) expressing cell lines, respectively. When measuring the EC 50 as a function of cellular activity uptake and nuclear dose, the EC 50 ranges from 0.020 to 0.039 Bq/cell and 3.3–9.2 Gy, respectively. Excretion through the hepatobiliary and renal pathways were observed in mice, with liver uptake of 2.3 ± 0.4 %ID/g after 48 h, contributing to estimated absorbed dose values in mice of 19.3 ± 0.3 mGy/MBq and 290 ± 10 mGy/MBq for [ 77 Br]RD1 and [ 76 Br]RD1, respectively. Conclusion: [ 77 Br]RD1 cytotoxicity was dependent on PARP expression and independent of BRCA1 status. Finally, the in vitro results suggest that [ 77 Br]RD1 cytotoxicity is driven by the targeted Meitner-Auger electron (MAe) radiotherapeutic effect of the agent. Further studies investigating the theranostic potential, organ dose, and tumor uptake of [ 76/77 Br]RD1 are warranted.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2022

Modified Energy Span Analysis of Catalytic Parallel Pathways and Selectivity

Mechanistic modeling provides vital insights into catalytic reactions. To analyze complex reaction networks with parallel pathways, we leverage the graph theory approach of the Energy Span Model (ESM) to develop a modified energy span analysis (MESA). A new method of cycle plots is proposed to perform reaction pathways analysis visually. We demonstrate this method on two published models: one describing carbon monoxide oxidation and the other simulating ethylene conversion to propanal via hydroformylation or ethane via hydrogenation. Fundamental insights explain kinetic observables, such as a reactant’s negative reaction order. General principles are revealed, such as rate-determining surface species being outside the primary reaction flux cycle and pathway selectivity being a purely kinetic property when reaction conditions are not near equilibrium. Lastly, we demonstrate MESA’s consistency with published microkinetic modeling results, highlighting this technique’s extension of the ESM to heterogeneous catalysts using collision theory to describe adsorption steps and concentration effects.

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

Structure of the electrical double layer at the ice–water interface

The surface of ice in contact with water contains sites that undergo deprotonation and protonation and can act as adsorption sites for aqueous ions. Therefore, an electrical double layer should form at this interface and existing models for describing the electrical double layer at metal oxide–water interfaces should be able to be modified to describe the surface charge, surface potential, and ionic occupancy at the ice–water interface. I used a surface complexation model along with literature measurements of the zeta potential of ice in brines of various strength and pH to constrain equilibrium constants. I then made predictions of ion site occupancy, surface charge density, and partitioning of counterions between the Stern and diffuse layers. The equilibrium constant for cation adsorption is more than 5 orders of magnitude larger than the other constants, indicating that this reaction dominates even at low salinity. Deprotonated OH sites are predicted to be slightly more abundant than dangling O sites, consistent with previous work. Surface charge densities are on the order of ±0.001 C/m2 and are always negative at the moderate pH values of interest to atmospheric and geophysical applications (6–9). In this pH range, over 99% of the counterions are contained in the Stern layer. This suggests that diffuse layer polarization will not occur because the ionic concentrations in the diffuse layer are nearly identical to those in the bulk electrolyte and that electrical conduction and polarization in the Stern layer will be negligible due to reduced ion mobility.

Daigle, Hugh (ORCID:0000000260628321)↗
DOE OSTI2024

Generative Physics-Informed Neural Network Solving Multi-Scale and Multi-Phase Plasma Chemical Flow Field

Low-temperature plasmas (LTPs) are non-equilibrium systems with near-room-temperature gas and highly energetic electrons. This makes them ideal for delicate applications in biomedicine and semiconductor manufacturing, enabling processes like wound healing, sterilization, etching, and plasma-enhanced chemical vapor deposition without thermal damage. However, LTPs involve complex chemistries, with hundreds of species and thousands of reactions, complicating their diagnosis, prediction, and control. Conventional diagnostics, such as Fourier-transform infrared spectroscopy (FTIR), laser-induced fluorescence (LIF), and optical emission spectroscopy (OES), offer limited species detection, while mass spectrometry (MS) struggles with low-sensitivity species. Additionally, LTP simulations face multi-scale challenges, as macroscopic fluid dynamics and microscopic particle collisions operate on vastly different timescales. To address these issues, we developed an artificial intelligence (AI) based diagnostic system: a generative physics-informed neural network (PINN-Gen) that can predict spatially resolved species concentrations and temperatures in LTPs by integrating experimental data from planar LIF with microscopic plasma chemical kinetics and macroscopic fluid mechanics, including plasma-liquid interactions at the interface between two phases. PINN-Gen solves no equations but checks the errors of physical laws by substituting the output from neural network, and the comparison with the experimental results. Thus, it naturally avoids the multi-scale difficulty of numerical simulations and predicts the results of conventionally unsolvable multi-scale and multi-phase problems. The real-time prediction will be robust due to the physical information used in the training of such a neural network, and only very limited input of condition required due to its generative feature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Assemblies of Polyacrylonitrile-Derived Photoactive Polymers as Blue and Green Light Photo-Cocatalysts for Cu-Catalyzed ATRP in Water and Organic Solvents

Photoluminescent nanosized quasi-spherical polymeric assemblies prepared by the hydrothermal reaction of polyacrylonitrile (PAN), ht-PLP PAN , were demonstrated to have the ability to photo-induce atom transfer radical polymerization (ATRP) catalyzed by low, parts per million concentrations of Cu II complex with tris(2-pyridylmethyl)amine (TPMA). Such photo induced ATRP reactions of acrylate and methacrylate monomers were performed in water or organic solvents, using ht-PLP PAN as the photo-cocatalyst under blue or green light irradiation. Mechanistic studies indicate that ht-PLP PAN helps to sustain the polymerization by facilitating the activation of alkyl bromide species by two modes: 1) green or blue light-driven photoreduction of the Cu II catalyst to the activating Cu I form, and 2) direct activation of dormant alkyl bromide species which occurs only under blue light. The photoreduction of the Cu II complex by ht-PLP PAN was confirmed by linear sweep voltammetry performed under illumination. Analysis of the polymerization kinetics in aqueous media indicated even though Cu I complexes comprised only 1–1.4% of all Cu species at equilibrium, they exhibited high activation rate constant and activated the alkyl bromide initiators five to six orders of magnitude faster than ht-PLP PAN .

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

Resolving atomic transport through iron oxide under irradiation using isotopic tracers

Structural materials in nuclear environments are exposed to complex and harsh conditions, including high temperatures, corrosive solutions, and irradiation. Material stability and likewise degradation are fundamentally controlled by atomic transport mechanisms. Here, we investigate the influence of radiation-induced, non-equilibrium point defect populations on mass transport within model oxide systems using isotopic tracers and atom probe tomography. An isotopically labelled hematite (a-Fe2O3) film was synthesized with 16O- and 18O-rich layers (>30nm and 10 nm, respectively) sequentially deposited via molecular beam epitaxy. One region of the film was proton irradiated to ~0.1 dpa at 450°C, and regions outside the beam spot provided a non-irradiated control of thermal annealing only. 3D APT was used to directly visualize the elemental and O-isotopic redistribution following both thermal and combined thermal and irradiation conditions. Complementary scanning transmission electron microscopy (STEM) was also employed to study the microstructural evolution. Error function fits of the 18O-concentration profiles enabled quantification of diffusion coefficients.

isotopic tracers, diffusion, proton irradiation, h↗
DOE OSTI2023

Prediction of alkaline earth metal ion adsorption on goethite for various background electrolytes with the CD-MUSIC model

As water scarcity drives the use of more saline water sources, contaminant fate and transport models must capture the impact of high concentrations of alkaline earth metal ions (AEMs) and background electrolytes in these more complex waters. By utilizing macroscopic adsorption data from various electrolyte systems, a Charge Distribution – Multisite Complexation (CD-MUSIC) model, capable of incorporating electrolyte adsorption, was able to accurately simulate the adsorption behavior of alkaline earth metal ions onto goethite. The modeling effort was guided by previous spectroscopic and surface complexation modeling of alkaline earth metal adsorption and built on previous CD-MUSIC modeling that accounted for changes in crystal face contributions to the surface site density as a function of specific surface area. Further, the model was constrained to consider only two dominant surface complex species for each metal ion adsorption reaction. These two species were selected from 44 possible species through objective curve fitting of single-solute macroscopic adsorption data. While most of the alkaline earth metal surface complexes formed outer-sphere complexes at the goethite surface, an inner-sphere species was utilized for Mg 2+ . With the surface complex species and equilibrium constants obtained from this study, the calibrated model successfully predicted alkaline earth metal ion adsorption over a wide range of solution and surface conditions; the model predictions encompassed a wide range of pH (5–11), solute/solid ratio (1.37 × 10 -5 – 8.33 × 10 -4 mol -solute /g -solid ), ionic strengths (0.01 M – 0.7 M), and background electrolytes (Na + , Cs + , Rb + , Cl - , and NO 3 - ) using the same crystal face contribution methodology for site density, capacitance values, and surface acidity constants adopted for proton and cadmium adsorption in previous work (Han and Katz, 2019). Model simulations for a range of background water chemistries demonstrated the potential for Mg 2+ to reduce Cd 2+ adsorption to goethite in model seawater and oil- and gas-produced waters.

42 ENGINEERING↗
DOE OSTI2022

Toward Better Understanding of Microphysical Processes and Resulting Precipitation Physics: A Merger Of Observations and Cloud Models (Final Report)

A large database of global disdrometer observations and a diverse set of simulations from the Regional Atmospheric Modeling System (RAMS) were compared using principal component analysis (PCA) in order to better understand warm and ice-based precipitation processes. The analysis demonstrated that six distinct precipitation groups (PGs) with common characteristics were revealed in both the observations and model simulations. These PGs were defined by similar co-variability of rain parameters. However, the model showed large concentrations of drops with a given size compared to the observations. Detailed investigations determined the parameterization of drop breakup was forcing the drops to an equilibrium size with a stronger constraint than was indicated by the observations. A series of sensitivity studies permutating the rain shape parameter in the RAMS 2-moment bulk scheme demonstrated how the assumed shape of the rain distribution influences the microphysics as well as precipitation characteristics. Results were compared to the same simulations performed with a bin microphysics scheme where the rain shape parameter is freely evolving. In the bin scheme, a wide range of shape parameters were found, with horizontal and vertical variability, in addition to being a function of storm lifetime. The results of this study highlight the limitations of a fixed assumed rain distribution in 2-moment microphysics schemes. Model simulations were used to probe the microphysical origins of the six PGs. Rain budgets from a supercell case showed that several of the groups hypothesized to be associated with strong ice microphysical processes had significant contributions from melting hail, but the complexity of ice and warm-rain processes in the supercell precluded definitive determination of microphysical origins for many of the PGs. A case study from the Mid-latitude Continental Convection and Clouds Experiment (MC3E) compared the model PGs to disdrometer observations, and demonstrated the spatial and temporal variability of the PGs which is not possible from disdrometer point measurements alone. The analyses performed during this project demonstrated how statistical analysis can provide a robust method for comparing model simulations and observations, which ultimately resulted identification of limitations in some current microphysics parameterizations as well as insights into precipitation variability which ultimately benefits both observations and modeling efforts.

54 ENVIRONMENTAL SCIENCES↗