Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic structure factor”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Improving Signal to Noise Ratios in Ion Mobility Spectrometry and Structures for Lossless Ion Manipulations (SLIM) using a High Dynamic Range Analog-to-Digital Converter

Signal digitization is a commonly overlooked part of ion mobility-mass spectrometry (IMS-MS) workflows, yet it is a significant contributor for determining signal-to-noise ratios and MS resolution. Here we report on the integration of a 2 GS/s, 14-bit ADC with a structures for lossless ion manipulations (SLIM)-IMS-MS and compare the performance to a commonly used 8-bit ADC. The 14-bit ADC provided an effective reduction in digitized noise by factor of ~6, owing largely to the use of smaller bit sizes. The low baseline allowed the threshold voltage levels to be set very close to the MCP baseline voltage, allowing for as much signal to be acquired as possible without causing overloading or excessive digitization of MCP baseline noise. Analyses of Agilent tuning mixture ions and a complex mixture of heavy labeled phosphopeptides showed that the 14-bit ADC (compared to the 8-bit ADC) provided a modest signal-to-noise increase (~1.5 to 2-fold) for high intensity ions, such as the Agilent tuning mixture ions and the 2+ and 3+ charge states of many phosphopeptide constituents. However, signal enhancements were as much as 10-fold for low intensity ions, and the 14-bit ADC enabled discernable signal intensities otherwise lost using an 8-bit digitizer. Additionally, the 14-bit ADC required ~14-fold fewer mass spectra to be averaged to produce a mass spectrum with similar S/N as the 8-bit ADC under identical conditions, potentially providing an order of magnitude higher measurement throughput. The high resolution, low baseline, and fast speed of the new 14-bit ADC enables high performance digitization of MS, IMS-MS, and SLIM-IMS-MS spectra, and allows a much fuller picture of analyte profiles in complex mixtures to be acquired.

data acquisition, digitization, ion mobility spect↗
DOE OSTI2021

Exploring Transport Behavior in Hybrid Perovskites Solar Cells via Machine Learning Analysis of Environmental-Dependent Impedance Spectroscopy

Hybrid organic–inorganic perovskites are one of the promising candidates for the next-generation semiconductors due to their superlative optoelectronic properties. However, one of the limiting factors for potential applications is their chemical and structural instability in different environments. Herein, the stability of (FAPbI3)0.85(MAPbBr3)0.15 perovskite solar cell is explored in different atmospheres using impedance spectroscopy. An equivalent circuit model and distribution of relaxation times (DRTs) are used to effectively analyze impedance spectra. DRT is further analyzed via machine learning workflow based on the non-negative matrix factorization of reconstructed relaxation time spectra. This exploration provides the interplay of charge transport dynamics and recombination processes under environment stimuli and illumination. The results reveal that in the dark, oxygen atmosphere induces an increased hole concentration with less ionic character while ionic motion is dominant under ambient air. Under 1 Sun illumination, the environment-dependent impedance responses show a more striking effect compared with dark conditions. In this case, the increased transport resistance observed under oxygen atmosphere in equivalent circuit analysis arises due to interruption of photogenerated hole carriers. The results not only shed light on elucidating transport mechanisms of perovskite solar cells in different environments but also offer an effective interpretation of impedance responses.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Cloud-based implementation of white-box model predictive control for a GEOTABS office building: A field test demonstration

Model predictive control (MPC) has been proven in simulations and pilot case studies to be a superior control strategy for large buildings. MPC can utilize the weather and occupancy schedule forecasts, together with the system model, to predict the future thermal behavior of the building and minimize the overall energy use and maximize thermal comfort. However, these advantages come with the cost of increased modeling effort, computational demands, communication infrastructure, and commissioning efforts. Thus a typical approach is to, often rapidly, simplify the building modeling and MPC optimization problem while paying a price of not reaching the full performance potential. It has been shown that by employing accurate physics-based models, MPC performance can be notably increased closer to its theoretical performance bound. However, implementation of such high-fidelity MPC in real buildings remains a challenge, resulting in a lack of successful field test studies. This work presents the methodology and field test demonstration of a computationally efficient implementation of the white-box MPC in an office building in Belgium. The detailed model of the building is based on first-principle physical equations. The deployment and supervision of MPC operation in a practical setting are supported by an automated cloud-based communication infrastructure. The motivating factor behind the cloud-based architecture is its compatibility with a commercially appealing control as a service concept. The building is equipped with a ground source heat pump (GSHP) and thermally activated building structures (TABS), where the combination of both is also known as GEOTABS. From a control perspective, GEOTABS buildings are particularly challenging systems due to large scale, complex heating, ventilation and air conditioning (HVAC) system, and slow dynamics with time delays. On the other hand, there is an increased potential for energy savings due to the high thermal mass, which acts as thermal storage. The MPC operation is demonstrated during the challenging transient seasons (switching between heating and cooling), and its performance is compared to a traditional rule-based controller (RBC). We provide a proof of concept of real MPC operation for the most difficult seasons with notable GSHP energy use savings equal to 53.5% and thermal comfort improvement by 36.9%. Other MPC applications found in the literature describe tests for only cooling or only heating, and up to now only for a black-box or a grey-box approach.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Structural and spectral dynamics of single-crystalline Ruddlesden-Popper phase halide perovskite blue light-emitting diodes

Achieving perovskite-based high–color purity blue-emitting light-emitting diodes (LEDs) is still challenging. Here, we report successful synthesis of a series of blue-emissive two-dimensional Ruddlesden-Popper phase single crystals and their high–color purity blue-emitting LED demonstrations. Although this approach successfully achieves a series of bandgap emissions based on the different layer thicknesses, it still suffers from a conventional temperature-induced device degradation mechanism during high-voltage operations. To understand the underlying mechanism, we further elucidate temperature-induced device degradation by investigating the crystal structural and spectral evolution dynamics via in situ temperature-dependent single-crystal x-ray diffraction, photoluminescence (PL) characterization, and density functional theory calculation. The PL peak becomes asymmetrically broadened with a marked intensity decay, as temperature increases owing to [PbBr 6 ] 4- octahedra tilting and the organic chain disordering, which results in bandgap decrease. This study indicates that careful heat management under LED operation is a key factor to maintain the sharp and intense emission.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Multiple Environmental Influences on the Lightning of Cold-Based Continental Cumulonimbus Clouds. Part I: Description and Validation of Model

In this two-part paper, influences from environmental factors on lightning in a convective storm are assessed with a model. In Part I, an electrical component is described and applied in the Aerosol–Cloud model (AC). AC treats many types of secondary (e.g., breakup in ice–ice collisions, raindrop-freezing fragmentation, rime splintering) and primary (heterogeneous, homogeneous freezing) ice initiation. AC represents lightning flashes with a statistical treatment of branching from a fractal law constrained by video imagery. The storm simulated is from the Severe Thunderstorm Electrification and Precipitation Study (STEPS; 19/20 June 2000). The simulation was validated microphysically [e.g., ice/droplet concentrations and mean sizes, liquid water content (LWC), reflectivity, surface precipitation] and dynamically (e.g., ascent) in our 2017 paper. Predicted ice concentrations (~10 L -1 ) agreed—to within a factor of about 2—with aircraft data at flight levels (-10° to -15°C). Here, electrical statistics of the same simulation are compared with observations. Flash rates (to within a factor of 2), triggering altitudes and polarity of flashes, and electric fields, all agree with the coincident STEPS observations. The “normal” tripole of charge structure observed during an electrical balloon sounding is reproduced by AC. It is related to reversal of polarity of noninductive charging in ice–ice collisions seen in laboratory experiments when temperature or LWC are varied. Positively charged graupel and negatively charged snow at most midlevels, charged away from the fastest updrafts, is predicted to cause the normal tripole. Total charge separated in the simulated storm is dominated by collisions involving secondary ice from fragmentation in graupel–snow collisions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Preliminary modeling of triply periodic minimal surface (TPMS) structures using RELAP5-3D

With the United States Department of Energy (DOE)’s goal of quadrupling the nation’s nuclear energy supply by 2050, and with the Advanced Fuels Campaign pushing for new types of advanced reactor fuels and geometries, the need has arisen for new nuclear fuel designs. One such design is to swap out current nuclear fuel geometries in exchange for another type of geometry, called a Triply Periodic Minimal Surface (TPMS). TPMSs are self-supporting, infinitely repeating lattices—attributes that lend themselves well to additive manufacturing. These surfaces also possess enhanced heat transfer properties thanks to their internal area changes and large surface-area-to-volume ratios. Their drawback, however, is an increased pressure drop. Given the small amount of correlations and data (Reynolds numbers in the 2,000–8,000 range), and the minimal amount of experience so far obtained by modeling TPMS structures using 1D systems codes such as the Reactor Excursion and Leak Analysis Program (RELAP5-3D), further research into this topic was needed. Using data from the University of Wisconsin - Madison (UW), curve fits were created for both a Heat Transfer Coefficient (HTC) correlation and a Darcy friction factor empirical coefficient correlation. The curves’ coefficients and multipliers were then output and utilized in RELAP5-3D models of two upcoming experiments—Flow Loop for INFLUX Pressure drop (FLIP) and Microreactor Agile Non-nuclear Experimental Test (MAGNET)—aimed at increasing the available data for Reynolds numbers to the 16,000–36,000 range for TPMS structures. The models were run under the conditions utilized by a Computational Fluid Dynamics (CFD) analysis performed by another group at Idaho National Laboratory. Only CFD pressure drop values were obtained from the FLIP test, and those values showed that the RELAP5-3D models had a lower rate of pressure increase in comparison to the CFD values. In addition, there seemed to be a vertical shift upward in the pressure drop for both models whenever the TPMS porosity decreased, and the RELAP5-3D models showed a higher vertical shift in comparison to the CFD values. The MAGNET results did not correspond to any CFD or experimental results against which they could be compared, so they were instead compared against the proposed CFD input conditions. These values were then compared with each other to make sure the model seemed to be performing as expected, paving the way for future tests that can be run for the purpose of further analyses and comparisons. The pressure drop increased with temperature and mass flow rate independently. The temperature change would decrease with increasing mass flow rate and temperature, which was just as we expected based on the fact that the lower viscosity and decreased density would result in higher friction and churning losses. The last metric that was assessed was the enthalpy flow change, which increased with increasing mass flow rate and decreasing temperature.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

X-Ray Afterglows from the Gamma-Ray Burst “Large-angle” Emission

In this work, we derive basic analytical results for the timing and decay of the gamma-ray burst (GRB) counterpart and delayed afterglow light curves for a brief emission episode from a relativistic surface endowed with angular structure, consisting of a uniform core of size ${\theta }_{c}$ (Lorentz factor ${{\rm{\Gamma }}}_{c}$ and surface emissivity ${i}_{\nu ^{\prime} }^{{\prime} }$ are angle independent) and an axially symmetric power-law envelope (${\rm{\Gamma }}\sim {\theta }^{-g}$). In this "large-angle emission" model, radiation produced during the prompt emission phase (GRB) at angles $\theta \gt {\theta }_{c}$ arrives at the observer well after the burst (delayed emission). The dynamical time range of the very fast decaying GRB "tail" and of the flat afterglow "plateau" and the morphology of the GRB counterpart/afterglow are all determined by two parameters: the core's parameter ${{\rm{\Gamma }}}_{c}{\theta }_{c}$ and the envelope's Lorentz factor index g, leading to three types of light curves that display three post-GRB phases (type 1: tail, plateau/slow decay, post-plateau/normal decay), two post-GRB phases (type 2: tail and fast decay), or just one (type 3: normal decay). We show how X-ray light-curve features can be used to determine core and envelope dynamical and spectral parameters. Testing of the large-angle emission model is done using the Swift/XRT X-ray emission of two afterglows of type 1 (GRB 060607A, GRB 061121), one of type 2 (GRB 061110A), and one of type 3 (GRB 061007). We find that the X-ray afterglows with plateaus require an envelope Lorentz factor ${\rm{\Gamma }}\sim {\theta }^{-2}$ and a comoving-frame emissivity ${i}_{\nu ^{\prime} }^{{\prime} }\sim {\theta }^{2}$; thus, for a typical afterglow spectrum ${F}_{\nu }\sim {\nu }^{-1}$, the lab-frame energy release is uniform over the emitting surface.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Integration and validation of some modules for modelling of high-speed chemically reactive flows in two-phase gas-droplet mixtures

Three modules are integrated into the built-in OpenFOAM rhoCentralFoam solver towards accurate and efficient modelling of high-speed chemically reactive flows in two-phase gas-droplet mixtures within the OpenFOAM 10.0 framework. The first module is the mixture-averaged diffusion model. The second module is the built-in OpenFOAM Lagrangian solver coupled with optimised droplet drag coefficient and convective heat transfer coefficient sub-models. The last module is a sparse stiff chemistry solver based on dynamic adaptive hybrid integration (AHI-S). The optimised droplet sub-models are first verified in correct implementation for subsequent simulations in this work. Further, they show good accuracy against experimental and analytical data in the modelling of ammonia droplet acceleration and cooling in the flowing and/or low-temperature air. The accuracy and efficiency gains related to the mixture-averaged diffusion model and the AHI-S chemistry solver are examined by simulating 1-D detonation propagation in ammonia droplet-free/laden ammoniaoxygen mixtures. Numerical results of detonation propagation speed, gaseous temperature, density, and species distributions around the induction zone show good agreement with experimental data and analytical solutions. Compared to the built-in OpenFOAM diffusion model, the mixture-averaged diffusion model provides different numerical predictions of pulsating instabilities in detonation propagation. It shows better accuracy in depicting the detonation structure within the droplet-free section attributed to improved multi-component diffusion modelling. Compared to the built-in OpenFOAM solver EulerImplicit (backward Euler), the AHI-S chemistry solver reduces the computational cost by around 50%. It achieves satisfactory accuracy in calculating detonation propagation speed within the droplet-free section with the optimal efficiency when the safety factor, β, equals 0.5.

42 ENGINEERING↗
DOE OSTI2020

Supporting Bioproducts Industry Growth with a System Dynamics Decision-Support Tool

A bio-based economy requires chemical products as well as fuels to be produced from biomass. Although a variety of universities, government agencies, start-ups and established firms have engaged in bioproduct development, many projects have failed to reach the point of commercialization and commercialized bioproducts struggle to capture and maintain market share. To date, there has been no general research into the factors that contribute to bioproduct failure or success. This work presents the Bioproduct Transition Dynamics (BTD) system dynamics model, a decision-support tool that simulates the bioproduct development process from pre-piloting research through construction of the first commercial-scale plant, as well as the processes of obtaining funding from investors and government agencies. The core of the BTD is a feedback loop between bioproduct developers and funders, which relates development progress measured with indicators such as net present value to funders’ decisions to continue investing. External factors such as feedstock prices, market size and growth, and the existence of bioproduct consumers also influence a development project’s chances of receiving follow-on funding. Virtually any bioproduct can be represented with the BTD: direct replacements, performance-advantaged products, niche and commodity markets can all be modeled. The goal of the BTD project is to inform decisions made by developers in both established firms and start-ups, investors, government agencies, and other stakeholders interested in growing the nascent U.S. bioproducts industry. This talk will cover the general structure and functionality of the BTD, and present results from an analysis performed with the BTD to demonstrate its use as a decision-support tool and the insights it can provide. The goal of the analysis is to identify the most critical factors that lead to direct replacement and performance-advantaged bioproduct projects emerging successfully from the “Valley of Death”, and to determine if these factors differ between the two bioproduct types.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2025

Microbiome dynamics in the congregate environment of U.S. Army Infantry training

Within military training and operational environments, individuals from diverse backgrounds share common spaces, follow structured routines and diets, and engage in physically demanding tasks. While there has been interest in leveraging microbiome features to predict and improve military health and performance, the longitudinal convergence of microbiomes in such constrained environments has not been established. To assess the degree of microbiome convergence, we performed shotgun metagenomic sequencing on swab samples from a military trainee cohort. Samples were taken across four different body sites, three timepoints, and two spatially distinct platoons. We observed evidence of convergence in one platoon, whereby similarity in microbiome composition increased over time, with numerous differentially abundant species. We found no indication of strain transfer between individuals, suggesting that convergence was influenced by external environmental factors, diet, and lifestyle. Microbial shifts observed in the convergence process included a decrease in fungal species, such as Malassezia restricta in nasal cavities, and a decrease in Prevotella species at inguinal regions across time. Shifts in multiple Corynebacterium species were also observed with varying magnitudes depending on the body site. Overall, we provide preliminary evidence of convergence of host microbial communities in military-associated environments that were distinguishable using shotgun metagenomic sequencing approaches. The data presented here on microbiome convergence, dynamics, and stability may inform risk-based mitigation in congregate military settings facilitating development of targeted microbial, dietary, or other interventions to optimize health and performance of military populations.

Biological and medical sciences↗
DOE OSTI2023

Self-consistent and detailed opacities from a non-equilibrium average-atom model

Modern density functional theory (DFT) is a powerful tool for accurately predicting self-consistent material properties such as equations of state, transport coefficients and opacities in high energy density plasmas, but it is generally restricted to conditions of local thermodynamic equilibrium (LTE) and produces only averaged electronic states instead of detailed configurations. We propose a simple modification to the bound-state occupation factor of a DFT-based average-atom model that captures essential non-LTE effects in plasmas—including autoionization and dielectronic recombination—thus extending DFT-based models to new regimes. Finally, we then expand the self-consistent electronic orbitals of the non-LTE DFT-AA model to generate multi-configuration electronic structure and detailed opacity spectra. This article is part of the theme issue ‘Dynamic and transient processes in warm dense matter’.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Virus diversity and activity is driven by snowmelt and host dynamics in a high-altitude watershed soil ecosystem

Background: Viruses impact nearly all organisms on Earth, including microbial communities and their associated biogeochemical processes. In soils, highly diverse viral communities have been identified, with a global distribution seemingly driven by multiple biotic and abiotic factors, especially soil temperature and moisture. However, our current understanding of the stability of soil viral communities across time and their response to strong seasonal changes in environmental parameters remains limited. Here, we investigated the diversity and activity of environmental soil DNA and RNA viruses, focusing especially on bacteriophages, across dynamics’ seasonal changes in a snow-dominated mountainous watershed by examining paired metagenomes and metatranscriptomes. Results: We identified a large number of DNA and RNA viruses taxonomically divergent from existing environmental viruses, including a significant proportion of fungal RNA viruses, and a large and unsuspected diversity of positive single-stranded RNA phages ( Leviviricetes ), highlighting the under-characterization of the global soil virosphere. Among these, we were able to distinguish subsets of active DNA and RNA phages that changed across seasons, consistent with a “seed-bank” viral community structure in which new phage activity, for example, replication and host lysis, is sequentially triggered by changes in environmental conditions. At the population level, we further identified virus-host dynamics matching two existing ecological models: “Kill-The-Winner” which proposes that lytic phages are actively infecting abundant bacteria, and “Piggyback-The-Persistent” which argues that when the host is growing slowly, it is more beneficial to remain in a dormant state. The former was associated with summer months of high and rapid microbial activity, and the latter with winter months of limited and slow host growth. Conclusion: Taken together, these results suggest that the high diversity of viruses in soils is likely associated with a broad range of host interaction types each adapted to specific host ecological strategies and environmental conditions. As our understanding of how environmental and host factors drive viral activity in soil ecosystems progresses, integrating these viral impacts in complex natural microbiome models will be key to accurately predict ecosystem biogeochemistry.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A simple technique to classify diffraction data from dynamic proteins according to individual polymorphs

One often observes small but measurable differences in the diffraction data measured from different crystals of a single protein. These differences might reflect structural differences in the protein and may reveal the natural dynamism of the molecule in solution. Partitioning these mixed-state data into single-state clusters is a critical step that could extract information about the dynamic behavior of proteins from hundreds or thousands of single-crystal data sets. Mixed-state data can be obtained deliberately (through intentional perturbation) or inadvertently (while attempting to measure highly redundant single-crystal data). To the extent that different states adopt different molecular structures, one expects to observe differences in the crystals; each of the polystates will create a polymorph of the crystals. After mixed-state diffraction data have been measured, deliberately or inadvertently, the challenge is to sort the data into clusters that may represent relevant biological polystates. Here, this problem is addressed using a simple multi-factor clustering approach that classifies each data set using independent observables, thereby assigning each data set to the correct location in conformational space. This procedure is illustrated using two independent observables, unit-cell parameters and intensities, to cluster mixed-state data from chymotrypsinogen (ChTg) crystals. It is observed that the data populate an arc of the reaction trajectory as ChTg is converted into chymotrypsin.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Evidence for Enhanced Tracer Diffusion in Densely Packed Interfacial Assemblies of Hairy Nanoparticles

Nearly monodisperse nanoparticle (NP) spheres attached to a nonvolatile ionic liquid surface were tracked by in situ scanning electron microscopy to obtain the tracer diffusion coefficient D tr as a function of the areal fraction Φ. The in situ technique resolved both tracer (gold) and background (silica) particles for ~1–2 min, highlighting their mechanisms of diffusion, which were strongly dependent on Φ. Structure and dynamics at low and moderate Φ paralleled those reported for larger colloidal spheres, showing an increase in order and a decrease in D tr by over 4 orders of magnitude. However, ligand interactions were more important near jamming, leading to different caging and jamming dynamics for smaller NPs. Finally, the normalized D tr at ultrahigh Φ depended on particle diameter and ligand molecular weight. Increasing the PEG molecular weight by a factor of 4 increased D tr by 2 orders of magnitude at ultrahigh Φ, indicating stronger ligand lubrication for smaller particles.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

LANSCE CCL Performance Limits

This report summarizes the performance limits of the LANSCE Coupled-Cavity Linac (CCL). These results are captured or summarized directly from the references cited. This report was written in support of the LANSCE Modernization Project (LAMP). Many factors contribute to the performance limits of the CCL accelerator system. One performance limit is set by the CCL mechanical structure and available cooling. Heating of the structure is ultimately linked to the operating RF duty factor of the CCL. Another performance limit is set by the maximum RF and beam duty factors that can be supported by the 805-MHz klystrons and associated HVDC power supplies. The RF and HVDC systems also set limits for the maximum peak beam current that can be accelerated in the CCL. And finally, a maximum beam current limit is set by the beam dynamics determined by the details of the CCL physics design (number of cells per CCL tank, accelerating gradients, magnetic focusing lattice, beam losses, etc.). This limit can be informed by both simulation results and beam measurement data, if available. A detailed explanation of each performance limit is given in the sections below. The results are summarized in the table below. In all cases, it is assumed that the CCL and the RF system are operating in their nominal beam production configuration – nominal magnet set points and nominal cavity fields, unless otherwise specified.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Data-driven linear time advance operators for the acceleration of plasma physics simulation

In this study, we demonstrate the application of data-driven linear operator construction for time advance with a goal of accelerating plasma physics simulation. We apply dynamic mode decomposition (DMD) to data produced by the nonlinear SOLPS-ITER (Scrape-off Layer Plasma Simulator - International Thermonuclear Experimental Reactor) plasma boundary code suite in order to estimate a series of linear operators and monitor their predictive accuracy via online error analysis. We find that this approach defines when these dynamics can be represented by a sequence of approximate linear operators and is essential for providing consistent projections when compared to an unconstrained application. For linear diffusion and advection–diffusion fluid test problems, we construct and apply operators within explicit and implicit time advance schemes, demonstrating that stability can be robustly guaranteed in each case. We further investigate the use of the linear time advance operators within several integration methods including forward Euler, backward Euler, and the matrix exponential. The application of this method to simulation data from SOLPS-ITER, with varying levels of Markov chain Monte Carlo numerical noise, shows that constrained DMD operators yield a capability to identify, extract, and integrate a (slow) subset of the present timescales. Example applications show that for projected speedup factors of [Formula: see text], and [Formula: see text], a mean relative error of 3%, 5%, and 8% and maximum relative error less than 20% are achievable, which appears acceptable for typical SOLPS-ITER steady-state simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Modeling bacterial microcompartment architectures for enhanced cyanobacterial carbon fixation

The carboxysome is a bacterial microcompartment (BMC) which plays a central role in the cyanobacterial CO 2 -concentrating mechanism. These proteinaceous structures consist of an outer protein shell that partitions Rubisco and carbonic anhydrase from the rest of the cytosol, thereby providing a favorable microenvironment that enhances carbon fixation. The modular nature of carboxysomal architectures makes them attractive for a variety of biotechnological applications such as carbon capture and utilization. In silico approaches, such as molecular dynamics (MD) simulations, can support future carboxysome redesign efforts by providing new spatio-temporal insights on their structure and function beyond in vivo experimental limitations. However, specific computational studies on carboxysomes are limited. Fortunately, all BMC (including the carboxysome) are highly structurally conserved which allows for practical inferences to be made between classes. Here, we review simulations on BMC architectures which shed light on (1) permeation events through the shell and (2) assembly pathways. These models predict the biophysical properties surrounding the central pore in BMC-H shell subunits, which in turn dictate the efficiency of substrate diffusion. Meanwhile, simulations on BMC assembly demonstrate that assembly pathway is largely dictated kinetically by cargo interactions while final morphology is dependent on shell factors. Overall, these findings are contextualized within the wider experimental BMC literature and framed within the opportunities for carboxysome redesign for biomanufacturing and enhanced carbon fixation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Effects of Composition and Oxidation States on the Structures of Chromium-Containing Sodium Silicate Glasses: Molecular Dynamics Simulations using Machine Learning Interatomic Potentials

Chromium represents a significant challenge for the vitrification of high-level nuclear waste into silicate and borosilicate glasses due to its low solubility and variable oxidation states, which can limit the waste loading due to promotion of crystallization or phase separation during processing. In this study, we modeled chromium containing silicate glasses using molecular dynamics simulations with three machine learning interatomic potentials (MLIPs), MACE, CHGNet, and PFP were employed, to gain insights on glass composition and oxidation states on the structures of these glasses. One of the goals is to evaluate their ability of these MLIPs to accurately represent the general structure of silicate glasses and chromium local environments as a function of chromium oxidation states. Density Functional Theory (DFT) based calculations and experimental data such as neutron structure factors were used to validate the structural models. It was found that the foundation models of all three MLIPs are able to reproduce general structural features of the sodium silicate glass structure consistent with experimental and DFT data, but only CHGNet and PFP can accurately capture the oxidation states and local environment of chromium: tetrahedral for Cr6+ and octahedral for Cr3+. Furthermore, we studied the effect of varying Cr3+/ Cr6+ (Cr3+/Crtotal) ratio and total chromium content using PFP. Our results show that Cr6+ enhances network polymerization by reducing non-bridging oxygens through Na? charge compensation required due to the formation of chromate (CrO42-) species, while Cr³? acts as a network modifier that disrupts connectivity. System size effects on the structural characteristics and chromium environments were also tested using the PFP potential. This work highlights the importance of careful validation on the precision, transferability, and potential of MLIPs for modeling glasses containing transition metal elements that can exist in multiple oxidation states. It is also encouraging to see the foundational models are all three MLFFs are able to reproduce the basic sodium silicate glass structures, while suggesting additional training or refining is needed to improve the description of more complex systems containing transition metals.

Puga, Christina L.↗