Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamic inference”

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

Excitons in strained and suspended monolayer WSe 2

Here, we study suspended membranes of atomically thin WSe 2 as hosts of excitons. We perform optical reflectance measurements to probe the exciton physics and obtain the peak energies for the 1$s$, 2$s$, and 3$s$ states of the $A$ exciton in suspended WSe 2 and consider supported membranes as a reference. We find that elimination of the influence of the dielectric environment enables a strong electron–hole interaction and a concomitant increase in the exciton binding energy in suspended monolayer (1L) WSe 2 . Based on the experimental results, we calculate the excitonic binding energies by employing the recently developed quantum electrostatic heterostructure model and the commonly employed Rytova–Keldysh potential model. We see that the binding energy of the ground state $A$ exciton increases from about 0.3 eV (on a substrate) to above 0.4 eV (suspended). We also exploit the tunability of the excitons in suspended samples via mechanical strain. By applying external gas pressure of 2.72 atm to a 1L suspended over a circular hole of 8 μm diameter, we strain the WSe 2 and obtain a reversible 0.15 eV redshift in the exciton resonance. The linewidth of the $A$ exciton decreases by more than half, from about 50 to 20 meV under 1.5% biaxial strain at room temperature. This line narrowing is due to the suppression of intervalley exciton–phonon scattering. By making use of the observed strain-dependent optical signatures, we infer the two-dimensional (2D) elastic moduli of 1L and 2L WSe 2 . Our results exemplify the use of suspended 2D materials as novel systems for fundamental studies, as well as for strong and dynamic tuning of their optical properties.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Empirical Comparison of Machine Learning Approaches for Black-Box Modeling of Power Conversion System Dynamics

Inverter-based resources are key components in modern power systems, but accurately modeling their complex behavior can be challenging. Standard, generic converter models often oversimplify inverter dynamics, leading to significant errors in predicting performance. In this work, we compare several data-driven machine learning (ML) approaches for inverter modeling, performing experiments on power conversion systems, systematically varying input conditions, and recording the resulting voltages and currents. The ML models were then trained on this measured data to capture the inverter's dynamic response and to predict the inverter's output current. A performance comparison between the four ML models under study is conducted, laying the foundation for future work on hardware implementation for real-time inference.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2023

Dynamics and Power Balance of Near Unity Target Gain Inertial Confinement Fusion Implosions

The change in the power balance, temporal dynamics, emission weighted size, temperature, mass, and areal density of inertially confined fusion plasmas have been quantified for experiments that reach target gains up to 0.72. It is observed that as the target gain rises, increased rates of self-heating initially overcome expansion power losses. This leads to reacting plasmas that reach peak fusion production at later times with increased size, temperature, mass and with lower emission weighted areal densities. Analytic models are consistent with the observations and inferences for how these quantities evolve as the rate of fusion self-heating, fusion yield, and target gain increase. At peak fusion production, it is found that as temperatures and target gains rise, the expansion power loss increases to a near constant ratio of the fusion self-heating power. Finally, this is consistent with models that indicate that the expansion losses dominate the dynamics in this regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Learning physics-based reduced-order models from data using nonlinear manifolds

Here we present a novel method for learning reduced-order models of dynamical systems using nonlinear manifolds. First, we learn the manifold by identifying nonlinear structure in the data through a general representation learning problem. The proposed approach is driven by embeddings of low-order polynomial form. A projection onto the nonlinear manifold reveals the algebraic structure of the reduced-space system that governs the problem of interest. The matrix operators of the reduced-order model are then inferred from the data using operator inference. Numerical experiments on a number of nonlinear problems demonstrate the generalizability of the methodology and the increase in accuracy that can be obtained over reduced-order modeling methods that employ a linear subspace approximation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Quantification and prediction of solidification textures under additive manufacturing conditions

Crystallographic textures are a major determinant of the macroscale anisotropic properties of polycrystalline metallic alloys produced in a wide range of additive manufacturing (AM) processes. Here, we introduce a statistical method that can accurately quantify the degree of orientational order of textures despite the large random fluctuations in the orientation of individual grains inherent in AM processes. The method, demonstrated for laser and resolidification of AlSi thin films, extends Z-scoring to a dynamical regime to assess the statistical significance of observed textures compared to randomly generated ones at different stages of solidification. We further show that, combined with phase-field modeling, this method can be used to infer fundamental anisotropic properties of the solid-liquid interface that are essential for texture prediction, and are compared here to the results of atomistic simulations. In addition, phase-field modeling reveals that, even at rapid AM solidification rates, the observed 〈110〉-dominated textures in the AlSi thin films are controlled predominantly by the anisotropy of the interface free-energy and sheds light on the physical mechanism of grain competition. These results significantly enhance both the existing tools for the quantification and prediction of AM crystallographic textures and our basic understanding of their formation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Development of Whole System Digital Twins for Advanced Reactors: Leveraging Graph Neural Networks and SAM Simulations

Here, in this work, we introduce a novel method to develop whole system digital twins (DTs) for advanced nuclear reactors. This method treats a complex reactor system as a heterogeneous graph: with the system components as different types of graph nodes and their physical interconnections as edges. Based on the heterogeneous graph, a graph neural network combining graph convolution and temporal node attention is developed as the DT, facilitating a comprehensive understanding of the system's dynamic behavior. By utilizing the System Analysis Module (SAM) code for simulating various operational transients, we develop a graph-based database that trains the DT. This DT is characterized by two primary functions: It can infer the entire system's status using sparse node information, and it can predict the progress of transients based on current and historical system information. Our approach is validated through case studies on the Experimental Breeder Reactor II (EBR-II) system and a generic Fluoride-salt-cooled High-temperature Reactor (gFHR), demonstrating the DT's accuracy in forecasting operational transients. The DT's rapid computation capabilities enhance its potential for supporting advanced reactor operations, offering benefits in intelligent simulation, autonomous control, and anomaly detection, paving the way for improved safety analysis and intelligent component health management for advanced reactor systems and reducing their operations and maintenance cost.

EBR-II↗
DOE OSTI2020

Predicting and Managing Risk to Bats at Commercial Wind Farms using Acoustics

Bat populations in North America face novel threats from white-nose syndrome and widespread turbine-related mortality related to the rapidly expanding wind power industry in addition to long-standing pressures from habitat loss and degradation. Bats, unlike most small mammals, are long-lived and slow to reproduce, highlighting the importance of understanding and managing anthropogenic sources of mortality. My dissertation research used acoustic bat detectors to measure bat activity at commercial wind projects, predict patterns in risk, and design strategic measures to reduce fatality rates by curtailing turbine operation during periods when bats are most active. Bats collide with wind turbines only when their rotors are spinning, and risk of turbine-related fatality is therefore a dynamic factor that can be manipulated by curtailing turbine operation when bats are active. We first measured inter-detector variation in metrics of acoustic bat activity to understand how the acoustic detection process may affect inferences related to spatial and temporal variation in bat activity. Using acoustic detectors mounted on top of wind turbines at two commercial wind farms in West Virginia, we then demonstrated that the amount of bat activity recorded when turbines were operating aligned closely with bat fatality rates on multiple scales. Accordingly, the metric of bat activity exposed to turbine operation provides a meaningful, quantitative indicator of turbine-related bat fatality risk. Further, bats responded consistently to changing wind speed and temperature at turbines in both wind farms across multiple years, enabling exposed bat activity to be predicted accurately among turbines and years. Building on these results, we simulated exposure of bats to turbine operation and energy loss for curtailment strategies recommended by state and federal agencies in the United States and Canada. By adjusting parameters such as cut-in wind speeds and temperature thresholds, we demonstrated the ability to design strategic curtailment programs that achieve equivalent or greater predicted reductions in bat activity exposure for substantially less energy-production loss. Characterizing fatality risk on a finer scale using acoustics will help regulatory agencies and the wind industry alike reduce risks of population-level impacts to vulnerable bat species while continuing to expand large-scale renewable energy generation.

17 WIND ENERGY↗
DOE OSTI2021

Dynamic flow stress of pure polycrystalline aluminum: Pressure-shear plate impact experiments and extension of dislocation-based modeling to large strains

The dynamic thermo-mechanical behavior of pure aluminum has attracted renewed interest lately due to experimental observations of an anomalous increase in Hugoniot Elastic Limit (HEL) at incipient plasticity and elevated temperatures in polycrystalline pure metals. In context of current dislocation-mediated plasticity models for metals, this increase in dynamic strength is indicative of a transition in the rate controlling mechanism for dislocation glide from being thermally assisted to phonon-drag restricted due to increase in phonon viscosity at elevated temperatures. Though these studies have helped to shed light on these important mechanisms operative in FCC metals at incipient plasticity, the extent to which they contribute to flow stress, particularly at larger strains, remains unclear. Here, we address these questions through a combined experimental and modeling effort focused on investigating the evolution of dynamic flow stress in polycrystalline aluminum using experimental data gathered from a series of combined pressure-and-shear plate impact (PSPI) experiments designed to reveal the flow stress of pure aluminum at strain rates ~ 10 5 /s, plastic strains of up to 40% and temperatures ranging from room to 866 K. In all cases, the flow stress of aluminum, as inferred from the measured transverse particle velocity histories at the free surface of an elastic tungsten carbide target plate, reveals saturation with increasing plastic strains at stress levels that decrease with increasing test temperatures. Numerical simulations are performed to correlate the experimentally observed temperature and strain rate dependence of flow stress at small and large plastic strains using the Austin-McDowell dislocation-mediated plasticity model parametrized to normal plate impact experiments conducted in an earlier study by Zaretsky and Kanel (2012). Extensions to the model are made to better represent dynamic behavior of pure aluminum at larger plastic strains as observed in elevated temperature split Hopkinson Pressure bar (SPHB) experiments of Samanta (1971) and Lindholm and Yeakly (1965), and the combined pressure-and-shear plate-impact experiments conducted in the present study. The main theoretical extension to the Austin-McDowell model made in this paper is the introduction of a new rate- and temperature- dependent dynamic recovery function, which can potentially allow for an effective reduction in the rate of accumulation of dislocations at large plastic strains. The numerical predictions of the revised and re-calibrated plasticity model are brought into agreement with the experimental observations, correlating sufficiently well with dynamic yield stress at incipient plasticity and the flow stress levels at larger plastic strains, plastic strain rates in the range 10 3 – 10 6 /s, and elevated temperatures up to near melt. The model, in concert with the experimental measurements, suggests that at incipient plasticity the rate governing mechanism for plastic flow is phonon-drag restricted dislocation glide, whereas, at higher magnitudes of plastic strain, it transitions to stress-assisted thermally activated glide. The transition between these two rate governing mechanisms is controlled by the evolution of dislocations throughout the deformation process.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A single-cell Arabidopsis root atlas reveals developmental trajectories in wild-type and cell identity mutants

In all multicellular organisms, transcriptional networks orchestrate organ development. The Arabidopsis root, with its simple structure and indeterminate growth, is an ideal model to investigate the spatiotemporal transcriptional signatures underlying developmental trajectories. To map gene expression dynamics across root cell types and developmental time, we built a comprehensive, organ-scale atlas at single cell resolution. In addition to estimating developmental progressions in pseudotime, we employed the mathematical concept of optimal transport to infer developmental trajectories and identify their underlying regulators. To demonstrate the utility of the atlas to interpret new datasets, we profiled mutants for two key transcriptional regulators at single cell resolution, shortroot and scarecrow. We report transcriptomic and in vivo evidence for tissue trans-differentiation underlying a mixed cell identity phenotype in scarecrow. Our results support the atlas as a rich community resource for unraveling the transcriptional programs that specify and maintain cell identity to regulate spatiotemporal organ development.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Scanning SQUID study of ferromagnetism and superconductivity in infinite-layer nickelates

Infinite-layer nickelates R 1-x Sr x NiO 2 (R = La, Pr, Nd) are a class of superconductors with structural similarities to cuprates. Although long-range antiferromagnetic order has not been observed for these materials, magnetic effects such as antiferromagnetic spin fluctuations and spin-glass behavior have been reported. Different experiments have drawn different conclusions about whether the pairing symmetry is s or d wave. In this paper, we applied a scanning superconducting quantum interference device (SQUID) to probe the magnetic behavior of film samples of three infinite-layer nickelates (La 0.85 Sr 0.15 NiO 2 , Pr 0.8 Sr 0.2 NiO 2 , and Nd 0.775 Sr 0.225 NiO 2 ) grown on SrTiO 3 (STO), each with a nominal thickness of 20 unit cells. In all three films, we observed a ferromagnetic background. We also measured the magnetic susceptibility above the superconducting critical temperature in Pr 0.8 Sr 0.2 NiO 2 and La 0.85 Sr 0.15 NiO 2 and identified a non-Curie-Weiss dynamic susceptibility. Both magnetic features are likely due to NiO x nanoparticles. We investigated superconductivity in Pr 0.8 Sr 0.2 NiO 2 and Nd 0.775 Sr 0.225 NiO 2 , which exhibited inhomogeneous diamagnetic screening. The superfluid density inferred from the diamagnetic susceptibility in relatively homogeneous regions shows T-linear behavior in both samples. Finally, we observed superconducting vortices in Nd 0.775 Sr 0.225 NiO 2 . We determined a Pearl length of 330 µm for Nd 0.775 Sr 0.225 NiO 2 at 300 mK, both from the strength of the diamagnetism and from the size and shape of the vortices. These results highlight the importance of considering NiO x particles when interpreting experimental results for these films.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Using Community Science to Reveal the Global Chemogeography of River Metabolomes

River corridor metabolomes reflect organic matter (OM) processing that drives aquatic biogeochemical cycles. Recent work highlights the power of ultrahigh-resolution mass spectrometry for understanding metabolome composition and river corridor metabolism. However, there have been no studies on the global chemogeography of surface water and sediment metabolomes using ultrahigh-resolution techniques. Here, we describe a community science effort from the Worldwide Hydrobiogeochemistry Observation Network for Dynamic River Systems (WHONDRS) consortium to characterize global metabolomes in surface water and sediment that span multiple stream orders and biomes. We describe the distribution of key aspects of metabolomes including elemental groups, chemical classes, indices, and inferred biochemical transformations. We show that metabolomes significantly differ across surface water and sediment and that surface water metabolomes are more rich and variable. We also use inferred biochemical transformations to identify core metabolic processes shared among surface water and sediment. Finally, we observe significant spatial variation in sediment metabolites between rivers in the eastern and western portions of the contiguous United States. Our work not only provides a basis for understanding global patterns in river corridor biogeochemical cycles but also demonstrates that community science endeavors can enable global research projects that are unfeasible with traditional research models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Enhanced accuracy through ensembling of randomly initialized auto-regressive models for dynamical systems

Computational mechanics simulations using traditional finite element methods (FEM) require prohibitively expensive computational resources for real-time engineering applications, design optimization, and digital twin implementations. While machine learning (ML) surrogate models offer significant computational speedups, autoregressive ML models for time-dependent mechanical systems suffer from error accumulation that compromises long-term prediction reliability - a critical concern for engineering applications where accuracy over extended time horizons is essential for safety and performance assessments. Here, we propose a deep ensemble framework specifically designed to address this challenge in computational mechanics applications, where multiple ML surrogate models with random weight initializations are trained in parallel and their predictions aggregated during inference. This approach leverages statistical diversity to maximize information gain from a fixed set of training data and to mitigate error propagation, while maintaining the computational efficiency that makes ML surrogates attractive for engineering practice. We validate the framework on three representative problems spanning critical areas of computational mechanics: stress field evolution in heterogeneous microstructures under complex loading (relevant to advanced materials design and composite analysis), planetary-scale shallow water dynamics (applicable to environmental and geotechnical engineering), and Gray-Scott reaction-diffusion systems (relevant to mass transport and chemical process engineering). Across all test cases, the ensemble approach demonstrates consistent error reduction of 15-33% compared to individual models. The codes for this work are available on GitHub (https://github.com/Graham-Brady-Research-Group/AutoregressiveEnsemble_SpatioTemporal_Evolution).

autoregressive prediction↗
DOE OSTI2020

Effects of social structure and management on risk of disease establishment in wild pigs

Contact heterogeneity among hosts determines invasion and spreading dynamics of infectious disease, thus its characterization is essential for identifying effective disease control strategies. Yet, little is known about the factors shaping contact networks in many 28 wildlife species and how wildlife management actions might affect contact networks. Wild pigs in North America are an invasive, socially-structured species that pose a health concern for domestic swine given their ability to transmit numerous devastating diseases 2 such as African swine fever (ASF). Using proximity loggers and GPS data from 48 wild pigs in Florida and South Carolina, USA, we employed a probabilistic framework to estimate weighted contact networks. We determined the effects of sex, social group, and spatial distribution (monthly home range overlap and distance) on wild pig contact. We also estimated the impacts of management-induced perturbations on contact and inferred their effects on ASF establishment in wild pigs with simulation. Social group membership was the primary factor influencing contacts. Between-group contacts depended primarily on space use characteristics, with fewer contacts among groups separated by >2 km and no contacts among groups >4 km apart within a month. Modeling ASF dynamics on the contact network demonstrated that indirect contacts resulting from baiting (a typical method of attracting wild pigs or game species to a site to enhance recreational hunting) increased the risk of disease establishment by ~33% relative to direct contact. Low-intensity population reduction (<5.9% of the population) had no detectable impact on contact structure but reduced predicted ASF establishment risk relative to no population reduction. Here, we demonstrate an approach for understanding the relative role of spatial, social, and individual-level characteristics in shaping contact networks and predicting their effects on disease establishment risk, thus providing insight for optimizing disease control in spatially- and socially-structured wildlife species.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Origins of Conformational Heterogeneity in Peptoid Helices Formed by Chiral N -1-Phenylethyl Sidechains

N-substituted glycines (polypeptoids) containing chiral hydrophobic sidechains are known to fold into biomimetic alpha helices. These helix formers often produce conformationally heterogeneous structures and are difficult to characterize at a sub-nanometer resolution. Previously, peptoid N-1-phenylethyl (S)-enantiomer sidechains (Nspe) were inferred from various experiments to form right-handed helices and (R)-enantiomers (Nrpe), left-handed helices. Prior computational work for N(s/r)pe oligomers has struggled to reproduce this trend. Herein, quantum mechanics calculations and molecular dynamics simulations are used to understand the origins of this discrepancy. Results from DFT and molecular mechanics calculations on a variety of Nspe and Nrpe oligomers as a function of chain length are in agreement, showing that Nspe and Nrpe prefer left- and right-handed helices, respectively. Additional metadynamics simulations are used to study Nrpe and Nspe oligomers folding in water. These results show that the free-energy driving forces for assembly into a helical backbone configuration are very small (within ~k B T). Lastly, we compare DFT calculations for other experimentally characterized peptoid sidechains, N(r/s)sb, N(r/s)tbe, and N(r/s)npe. In this analysis, we show that peptoid sidechains determined to be more robust experimentally (tbe and npe) have helical preferences opposite the trend seen in less robust assemblies formed by N(r/s)pe and N(r/s)sb chemistries. Furthermore, the more robust tbe and nnpe favor the (S)-enantiomer to right-handed and the (R)-enantiomers to left-handed helices.

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

Enhanced Frequency Support Scheme of Generic Inverter-Based Resource Models for Renewable-Dominated Power Grids

The frequency response of SG-dominated power grids is predictable ahead of an occurrence of a frequency event because the frequency response of SGs is consistent, and it can be inferred from the swing equation [1]. However, increasing the portion of IBRs in an SG-dominated power grid might make the characteristics of the conventional power grids no longer valid because this changing resource mix affects grid dynamics and controls [2]. Thus, maintaining these characteristics greatly benefits the control and operation of the power grids with high penetration of IBRs. To maintain these characteristics in IBR-dominated power grids, IBRs should have frequency response capability similar to that of an SG. The WECC modeling validation subcommittee has developed generic IBR models for large system planning [3]-[5]. These models can represent various vendors' dynamic behavior for WTG, PV, and ESS [5]. The current generic IBR models approved by WECC can provide frequency response only from droop control loops in REPC models [6], [7]. The contribution of the loops is proportional to the frequency deviation from the nominal frequency. Thus, it presents an insufficient contribution to arrest frequency variation compared to the frequency response of SGs because it allows a high ROCOF in the early stage of frequency events. This shortfall will become greater as the PL of IBRs increases in power grids. Controller enhancement for the generic IBR models is required to secure the frequency stability under high PL of IBRs as in the SG-dominated power grids. This paper proposes a control extension for the generic IBR models to enhance the frequency support capabilities and discusses the classification of frequency support for the different types of IBR considering their operating constraints. An inertial control scheme is implemented in the REPC and REEC models of the generic IBR models to achieve these objectives. The inertial control scheme includes the following stages: Control area data acquisition, inertia time constant estimation, IBR-related constraint check, IBR contribution determination, and inertial response provision. In the scheme, a REPC acquires control area data from a system operator and estimates a total inertia time constant for the control area the applicable IBR power plant belongs. Then, the estimated inertial time constant is transferred to each IBR controller—REEC—within the power plant. Each REEC checks the availability of applicable IBR for inertial response participation. If the IBR is available, the REEC amplifies the estimated inertial time constant to utilize it for inertial response provision. In this way, the proposed inertial response scheme extends the functionality of the generic IBR models to provide SG-like frequency response within their constraints. Various scenarios considering different IBR types, IBR penetration levels, and frequency control schemes were simulated and compared in an IEEE 39-bus system using PSCAD simulator to verify the effectiveness of the proposed scheme.

Kim, Jinho↗
DOE OSTI2025

Joint Modeling of Quasar Variability and Accretion Disk Reprocessing Using Latent Stochastic Differential Equations

Quasars are bright active galactic nuclei powered by the accretion of matter around supermassive black holes at the center of galaxies. Their stochastic brightness variability depends on the physical properties of the accretion disk and black hole. The upcoming Rubin Observatory Legacy Survey of Space and Time (LSST) is expected to observe tens of millions of quasars, so there is a need for efficient techniques like machine learning that can handle the large volume of data. Quasar variability is believed to be driven by an X-ray corona, which is reprocessed by the accretion disk and emitted as UV/optical variability. We are the first to introduce an auto-differentiable simulation of the accretion disk and reprocessing. We use the simulation as a direct component of our neural network to jointly model the driving variability and reprocessing, trained with supervised learning on simulated LSST-like 10 yr quasar light curves. We encode the light curves using a transformer encoder, and the driving variability is reconstructed using latent stochastic differential equations, a physically motivated generative deep learning method that can model continuous-time stochastic dynamics. By embedding the physical processes of the driving signal and reprocessing into our network, we achieve a model that is more robust and interpretable. We demonstrate that our model outperforms a Gaussian process regression baseline and can infer accretion disk parameters and time delays between wave bands, even for out-of-distribution driving signals. Our approach provides a powerful framework that can be adapted to solve other inverse problems in multivariate time series.

Fagin, Joshua [City Univ. of New York (CUNY), NY (↗
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↗