Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model interpretability”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Downstream Separation of Formic Acid with Anion-Exchange Resin from Electrocatalytic Carbon Dioxide (CO 2 ) Conversion: Adsorption, Kinetics, and Equilibrium Modeling

The electrocatalytic carbon dioxide reduction reaction (CO 2 RR) that produces liquid formic acid (FA) is one potential route to minimize the CO 2 concentration in the atmosphere by utilizing CO 2 as a feedstock. However, one of the main disadvantages of this route is the high energy demand required for purification through conventional distillation. This method of purification has been shown to result in a large carbon footprint in life-cycle analysis (LCA). Therefore, finding an alternative purification technique that can integrate with an electrochemical CO 2 RR process with the potential to lower the carbon footprint would be of great interest. Under the experimental conditions of the CO 2 RR, FA remains in its basic form as the formate anion in the catholyte. Anion exchange resins with their basic properties can separate the formate anion from the catholyte mixture as an alternative to energy intensive distillation. In this work, three anion exchange resins with increasing basicity order: Amberlite IRA-96 < Amberlite IRA-910 < Ambersep 900 are tested to separate FA in the presence of aqueous solution of potassium bicarbonate catholyte with varying levels of resin and FA. Kinetics and equilibrium studies data for the FA adsorption are interpreted using several kinetics and isotherm models. The kinetics data fit well with a pseudo-first-order model at high initial FA concentrations and a pseudo-second-order model at low initial FA concentrations. The experimental data can be best explained with the Freundlich isotherm model. The fitted results show that strongly basic Ambersep 900 is more effective at separating FA, with a maximum adsorption capacity of 336.7 mg/g compared to Amberlite IRA-96 (275.2 mg/g) and Amberlite IRA-910 (209.2 mg/g) in neutral water. However, weakly basic Amberlite IRA-96 (110.8 mg/g) can separate FA more efficiently from potassium bicarbonate catholyte than Amberlite IRA-910 (99.9 mg/g) and Ambersep 900 (40.0 mg/g). As a result, competitive adsorption occurs in the presence of potassium bicarbonate catholyte, but the formic acid adsorption capacity is still sufficient for adequate separation.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2023

An Interpretable Index of Social Vulnerability to Environmental Hazards

Index-based measures of social vulnerability to environmental hazards are commonly modeled from composites of population-level risk factors. These models overlook individual context in communities' experiences of environmental hazards, producing metrics that may hinder spatial decision support for mitigating and responding to hazards. This paper introduces an interpretable, high-resolution model for generating an individual-oriented social vulnerability index (IOSVI) for the United States built on synthetic populations that couples individual and social determinants of vulnerability. The IOSVI combines an individual vulnerability index (IVI) that ranks individuals in an area’s synthetic population based on intersecting risk factors, with a social vulnerability index (SVI) based on the population’s cumulative distribution of IVI scores. Interpretability of the IOSVI procedure is demonstrated through examples of national, metropolitan, and neighborhood (census tract) level spatial variation in index scores and IVI themes, as well as an exploratory analysis examining risk factors affecting a specific sub-population (military veterans) in areas of high social and environmental vulnerability.

Tuccillo, Joe↗
DOE OSTI2023

Development of fracture diagnostic methods for fluid distribution based on quantitative interpretation of distributed acoustic sensing and distributed temperature sensing

Multistage hydraulic fracturing design on horizontal wells has significantly evolved with larger fluid volume, more fracturing stages, and tighter perforation cluster spacing to efficiently stimulate unconventional reservoirs. From the published field observations, the recent fracturing design results in complex fracture networks or swarm of fractures. Fracture treatment evaluation is extremely challenging in such a case because of the large amount of variables in well completion and stimulation design. Combined measurements from different technologies can help in fracture diagnosis. Fluid distribution, either during fracture injection or during production, directly relates to the stimulation efficiency at the cluster level and at the stage level. Because it is unlikely in the real world to distribute the injected fluid uniformly among all the clusters, we need diagnostic techniques to generate the flow profile along a lateral. Fiber-optic measurements, such as distributed acoustic sensing (DAS) and distributed temperature sensing (DTS), are currently used to diagnose downhole flow conditions. This technology allows us to qualitatively confirm the fluid flow profile and other issues occurring downhole during fracturing such as leakage through plugs. For optimizing a fracturing design, we also need to understand how the design parameters are correlated with the stimulation efficiency. In this study, we combine two sets of models of DAS and DTS data interpretation for injected fluid volume distribution. The DAS is interpreted based on an empirical correlation between fluid flow rates and frequency band energy from the acoustic signals. The DTS is interpreted by performing temperature history match-based thermal energy conservation. Because of the completely different physics behind the interpretations, the confirmation of two interpretations provides confidence in fluid distribution.

Geochemistry & Geophysics↗
DOE OSTI2023

A data-driven framework for permeability prediction of natural porous rocks via microstructural characterization and pore-scale simulation

Understanding the microstructure–property relationships of porous media is of great practical significance, based on which macroscopic physical properties can be directly derived from measurable microstructural informatics. However, establishing reliable microstructure–property mappings in an explicit manner is difficult, due to the intricacy, stochasticity, and heterogeneity of porous microstructures. In this paper, a data-driven computational framework is presented to investigate the inherent microstructure–permeability linkage for natural porous rocks, where multiple techniques are integrated together, including microscopy imaging, stochastic reconstruction, microstructural characterization, pore-scale simulation, feature selection, and data-driven modeling. A large number of 3D digital rocks with a wide porosity range are acquired from microscopy imaging and stochastic reconstruction techniques. A broad variety of morphological descriptors are used to quantitatively characterize pore microstructures from different perspectives, and they compose the raw feature pool for feature selection. Here high-fidelity lattice Boltzmann simulations are conducted to resolve fluid flow passing through porous media, from which reliable permeability references are obtained. The optimal feature set that best represents permeability is identified through a performance-oriented feature selection process, upon which a cost-effective surrogate model is rapidly fitted to approximate the microstructure-permeability mapping via data-driven modeling. This surrogate model exhibits great advantages over empirical/analytical formulas in terms of prediction accuracy and generalization capacity, which can predict reliable permeability values spanning four orders of magnitude. Besides, feature selection also greatly enhances the interpretability of the data-driven prediction model, from which new insights into the mechanism of how microstructural characteristics determine intrinsic permeability are obtained.

58 GEOSCIENCES↗
DOE OSTI2020

Measurement-induced criticality in random quantum circuits

We investigate the critical behavior of the entanglement transition induced by projective measurements in (Haar) random unitary quantum circuits. Using a replica approach, we map the calculation of the entanglement entropies in such circuits onto a two-dimensional statistical-mechanics model. In this language, the area- to volume-law entanglement transition can be interpreted as an ordering transition in the statistical-mechanics model. We derive the general scaling properties of the entanglement entropies and mutual information near the transition using conformal invariance. We analyze in detail the limit of infinite on-site Hilbert space dimension in which the statistical-mechanics model maps onto percolation. In particular, we compute the exact value of the universal coefficient of the logarithm of subsystem size in the n th Rényi entropies for n ≥ 1 in this limit using relatively recent results for the conformal field theory describing the critical theory of two-dimensional (2D) percolation, and we discuss how to access the generic transition at finite on-site Hilbert space dimension from this limit, which is in a universality class different from 2D percolation. Here, we also comment on the relation to the entanglement transition in random tensor networks, studied previously in Vasseur et al. [Phys. Rev. B 100, 134203 (2019)].

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

A Probabilistic Autoencoder for Type Ia Supernova Spectral Time Series

We construct a physically parameterized probabilistic autoencoder (PAE) to learn the intrinsic diversity of Type Ia supernovae (SNe Ia) from a sparse set of spectral time series. The PAE is a two-stage generative model, composed of an autoencoder that is interpreted probabilistically after training using a normalizing flow. We demonstrate that the PAE learns a low-dimensional latent space that captures the nonlinear range of features that exists within the population and can accurately model the spectral evolution of SNe Ia across the full range of wavelength and observation times directly from the data. By introducing a correlation penalty term and multistage training setup alongside our physically parameterized network, we show that intrinsic and extrinsic modes of variability can be separated during training, removing the need for the additional models to perform magnitude standardization. We then use our PAE in a number of downstream tasks on SNe Ia for increasingly precise cosmological analyses, including the automatic detection of SN outliers, the generation of samples consistent with the data distribution, and solving the inverse problem in the presence of noisy and incomplete data to constrain cosmological distance measurements. We find that the optimal number of intrinsic model parameters appears to be three, in line with previous studies, and show that we can standardize our test sample of SNe Ia with an rms of 0.091 ± 0.010 mag, which corresponds to 0.074 ± 0.010 mag if peculiar velocity contributions are removed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Structure–property relations of binary ferrite melts

Molten ferrite systems are used in the smelting and refining processes in steelmaking, to reduce the loss of metals in slags and to accelerate reaction rates. Here, high-energy x-ray diffraction experiments have been performed on aerodynamically levitated molten spheres of 43BaO–57FeO X and 43SrO–57FeO X at 1873 K using laser beam heating. The composition was varied within the range of x = 1–1.5 by changing the oxygen partial pressure of the levitation gas. The corresponding x-ray pair distribution functions have been interpreted using empirical potential structure refinement (EPSR) modeling. In oxygen-rich melts (x = 1.5), our EPSR models indicate very similar structures for the different alkaline-earth liquids, with both the Ba–O and Sr–O coordination numbers to be ∼8.4 and the total Fe–O coordination numbers ∼5.7. However, our models show that in reducing environments, the Fe 3+ and Fe 2+ ions exhibit very different behaviors in the Ba- and Sr-ferrite liquids. In the Ba-ferrite melt, the Fe 3+ –O coordination number decreases from 5.7 (at x = 1.5) to 5.2 (at x = 1.07), whereas Fe 2+ –O remains constant at ∼5.0 across the same compositional range. In the Sr melts, both the Fe 2+ –O and Fe 3+ –O coordination numbers rise from ∼5.7 (at x = 1.5) to 6.3 (at x = 1.07). All models show the structures to be heterogeneous with intertwined nanometer sized clusters or channels of Ba/Sr–O and Fe–O polyhedra that grow as oxygen content is reduced. Changes in the viscosity and electrical properties are interpreted in terms of the number of bridging and non-bridging oxygens associated with FeO 4 tetrahedra and concentration of charge carriers, respectively.

Benmore, Chris J. [Argonne National Laboratory (AN↗
DOE OSTI2020

Modeling Arctic boundary layer cloud streets at grey-zone resolutions

To better understand how model resolution affects the formation of Arctic boundary layer clouds, we investigated the influence of grid spacing on simulating cloud streets that occurred near Utqiagvik (formerly Barrow), Alaska on 2 May 2013 and were observed by the Moderate Resolution Imaging Spectroradiometer (MODIS). The Weather Research and Forecasting model was used to simulate the clouds using nested domains with increasingly fine resolution ranging from a horizontal grid spacing of 27 km in the boundary-layer-parameterized mesoscale domain to a grid spacing of 0.111 km in the large-eddy-permitting domain. We investigated the model-simulated mesoscale environment, horizontal and vertical cloud structures, boundary layer stability, and cloud properties, all of which were subsequently used to interpret the observed roll-cloud case. Increasing model resolution led to a transition from a more buoyant boundary layer to a more shear-driven turbulent boundary layer. The clouds were stratiform-like in the mesoscale domain, but as the model resolution increased, roll-like structures, aligned along the wind field, appeared with ever smaller wavelengths. A stronger vertical water vapor gradient occurred above the cloud layers with decreasing grid spacing. With fixed model grid spacing at 0.333 km, changing the model configuration from a boundary layer parameterization to a large-eddy-permitting scheme produced a more shear-driven and less unstable environment, a stronger vertical water vapor gradient below the cloud layers, and the wavelengths of the rolls decreased slightly. In this study, only the large-eddy-permitting simulation with gird spacing of 0.111 km was sufficient to model the observed roll clouds.

Lai, Hui-Wen↗
DOE OSTI2022

Uncertainty-aware mixed-variable machine learning for materials design

Abstract Data-driven design shows the promise of accelerating materials discovery but is challenging due to the prohibitive cost of searching the vast design space of chemistry, structure, and synthesis methods. Bayesian optimization (BO) employs uncertainty-aware machine learning models to select promising designs to evaluate, hence reducing the cost. However, BO with mixed numerical and categorical variables, which is of particular interest in materials design, has not been well studied. In this work, we survey frequentist and Bayesian approaches to uncertainty quantification of machine learning with mixed variables. We then conduct a systematic comparative study of their performances in BO using a popular representative model from each group, the random forest-based Lolo model (frequentist) and the latent variable Gaussian process model (Bayesian). We examine the efficacy of the two models in the optimization of mathematical functions, as well as properties of structural and functional materials, where we observe performance differences as related to problem dimensionality and complexity. By investigating the machine learning models’ predictive and uncertainty estimation capabilities, we provide interpretations of the observed performance differences. Our results provide practical guidance on choosing between frequentist and Bayesian uncertainty-aware machine learning models for mixed-variable BO in materials design.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The influence of physical and algorithmic factors on simulated far-field waveforms and source–time functions of underground explosions using unsupervised machine learning

SUMMARY Characterizing explosion sources and differentiating between earthquake and underground explosions using distributed seismic networks becomes non-trivial when explosions are detonated in cavities or heterogeneous ground material. Moreover, there is little understanding of how changes in subsurface physical properties affect the far-field waveforms we record and use to infer information about the source. Simulations of underground explosions and the resultant ground motions can be a powerful tool to systematically explore how different subsurface properties affect far-field waveform features, but there are added variables that arise from how we choose to model the explosions that can confound interpretation. To assess how both subsurface properties and algorithmic choices affect the seismic wavefield and the estimated source functions, we ran a series of 2-D axisymmetric non-linear numerical explosion experiments and wave propagation simulations that explore a wide array of parameters. We then inverted the synthetic far-field waveform data using a linear inversion scheme to estimate source–time functions (STFs) for each simulation case. We applied principal component analysis (PCA), an unsupervised machine learning method, to both the far-field waveforms and STFs to identify the most important factors that control variance in the waveform data and differences between cases. For the far-field waveforms, the largest variance occurs in the shallower radial receiver channels in the 0–50 Hz frequency band. For the STFs, both peak amplitude and rise times across different frequencies contribute to the variance. We find that the ground equation of state (i.e. lithology and rheology) and the explosion emplacement conditions (i.e. tamped versus cavity) have the greatest effect on the variance of the far-field waveforms and STFs, with the ground yield strength and fracture pressure being secondary factors. Differences in the PCA results between the far-field waveforms and STFs could possibly be due to near-field non-linearities of the source that are not accounted for in the estimation of STFs and could be associated with yield strength, fracture pressure, cavity radius and cavity shape parameters. Other algorithmic parameters are found to be less important and cause less variance in both the far-field waveforms and STFs, meaning algorithmic choices in how we model explosions are less important, which is encouraging for the further use of explosion simulations to study how physical Earth properties affect seismic waveform features and estimated STFs.

58 GEOSCIENCES↗
DOE OSTI2023

Search for heavy Higgs bosons with flavour-violating couplings in multi-lepton plus b-jets final states in pp collisions at 13 TeV with the ATLAS detector

A search for new heavy scalars with flavour-violating decays in final states with multiple leptons and b-tagged jets is presented. The results are interpreted in terms of a general two-Higgs-doublet model involving an additional scalar with couplings to the top-quark and the three up-type quarks (ρ tt , ρ tc , and ρ tu ). The targeted signals lead to final states with either a same-sign top-quark pair, three top-quarks, or four top-quarks. The search is based on a data sample of proton-proton collisions at $\sqrt{s}$ = 13 TeV recorded with the ATLAS detector during Run 2 of the Large Hadron Collider, corresponding to an integrated luminosity of 139 fb –1 . Events are categorised depending on the multiplicity of light charged leptons (electrons or muons), total lepton charge, and a deep-neural-network output to enhance the purity of each of the signals. Masses of an additional scalar boson m H between 200 – 630 GeV with couplings ρ tt = 0.4, ρ tc = 0.2, and ρ tu = 0.2 are excluded at 95% confidence level. Additional interpretations are provided in models of R-parity violating supersymmetry, motivated by the recent flavour and (g – 2) μ anomalies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Error field correction strategies in preparation to MAST-U operation

In magnetic fusion devices, unwanted non-axisymmetric magnetic field perturbations, known as error fields (EF), can have detrimental effects on plasma stability and confinement. Such EFs may originate from several sources, such as axi-symmetric coil misalignments, coil feeds, 3D structures in the wall surrounding the plasma, presence of ferromagnetic materials near the plasma surface, blanket materials in future devices. To minimize their impact on plasma performance and on the available operational space, it is important to identify the EF sources and develop EF control strategies. MAST Upgrade (MAST-U) is a spherical tokamak which will operate in the near future after a series of enhancements from the previous MAST experiment. To deliver a machine with EF amplitude low enough to allow a high quality experimental programme, systematic analysis of the intrinsic EF sources has been carried out for poloidal eld (P) and divertor (D) coils. To deliver a machine with EF amplitude low enough to allow a high quality experimental programme, systematic analysis of the intrinsic EF sources has been carried out for poloidal field (P) and divertor (D) coils, whose magnetic eld measurements were available when writing this paper. Such measurements reveal that P and D coils are 3D deformed and thus are responsible for intrinsic EFs, with mainly n=1 and n=2 toroidal mode numbers. In preparation to MAST-U operation, both passive and active EF control strategies have been adopted for n=1 and n=2 EFs compensation. Passive EF control consisted of installing finely each D and P coil within MAST-U device so as to minimize the intrinsic n=1 EF amplitude. The optimal coil alignment has been determined based on magnetic eld measurements and the corresponding 3D electro-magnetic modelling, and envisaged coil shift and tilt of 3.2 mm and 0.7 mrad, respectively, in the case of P coil named P4. Conversely, active control will be use during MAST-U operation to reduce the n=2 EF which is associated mainly with P4 and P5, as well. Since these coils have been re-used from the MAST device, studies attempting n=2 EF control, based on MAST plasmas, have been modelled utilizing the MARS-F code to interpret experimental results in MAST and to give hints for future n=2 EF control studies in MAST-U. Here, a model-based control set for n=2 EF control has been identified which would allow for the minimization of rotation braking, of the resonant magnetic field at the q=2 and of the plasma displacement in MAST-U, simultaneously.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Interpreting Tafel behavior of consecutive electrochemical reactions through combined thermodynamic and steady state microkinetic approaches

Assessing the reaction pathway of multi-electron-transfer reactions is an essential yet difficult task for the rational design of electrocatalysts. In this work, we develop a heuristic approach that combines thermodynamic adsorption energetics calculated through density functional theory with microkinetic modeling using the steady state approximation to interpret the potential-dependent Tafel behavior of consecutive electrochemical reactions. In doing so, we introduce a kinetic framework for ab initio calculations that ensures self-consistent adsorption energetics based on kinetically limited adsorbate coverages. The approach is applied to experimental results on CoO x (OH) 2-x single crystal electrocatalyst particles yielding coverage dependent mechanistic information and identification of the rate-limiting step with standard rate constants for the oxygen evolution reaction on the ($11\bar{2}0$) surfaces of the β-Co(OH) 2 , β-CoOOH, and CoO 2 bulk phases. Finally, this generalizable method enables catalyst benchmarking based on determining the active species involved and associated intrinsic reaction rate constants in consecutive multi-electron-transfer reactions.

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

All-Atom Modeling and Simulation of Biopolymer Interface: Dual Role of Antifouling Polymer Brushes

Antifouling polymer brushes are well-known for their exceptional resistance to unwanted protein adsorption. While experimental studies have extensively characterized protein–polymer brush interactions, computational investigations remain limited, largely due to the challenges in accurately modeling and integrating these complex, multicomponent systems without resorting to oversimplification. To address this challenge, this study presents one of the most comprehensive and realistic model systems to date, comprising the substrate, grafted polymers, and proteins. Our work interprets the interactions between polymer brush and protein based on realistic modeling without simplification. In particular, this study utilizes molecular modeling and simulation of polycationic and polyzwitterionic brushes─poly(dimethylaminoethyl methacrylate) (PDMAEMA), poly(2-(N-oxide-N,N-dimethylamino)ethyl methacrylate) (PNOMA), and poly(2-(N-3-sulfopropyl-N,N-dimethylammonium)ethyl methacrylate) (PSBMA)─grafted onto α-quartz substrates via polymerization initiator linkers. The brush models were developed to closely replicate experimentally synthesized samples and to provide detailed insights into structural and dynamical changes at the molecular level during protein adsorption. Using steered molecular dynamics simulations, we show that the PSBMA brush, due to its high local density, exhibits the greatest resistance to protein insertion. Cα root-mean-square deviation and interaction pattern analyses further reveal that the PSBMA brush also induces the most significant destabilization of lysozyme, while the PDMAEMA brush enhances protein stability through ion-mediated interactions. The PNOMA brush, while requiring the lowest force for protein adsorption, induces greater protein destabilization than the PDMAEMA brush, primarily due to electrostatic repulsion caused by a short carbon spacer length. Hydration analysis reveals that both the PSBMA brush and the lysozyme interacting with it exhibit the most rapid dehydration, attributed to the brush’s high local chain density, which results in the greatest lysozyme destabilization and the highest adsorption force. These findings highlight the dual role of antifouling polymer brushes: resisting protein adsorption and modulating protein structural dynamics. In conclusion, this study provides valuable insights for the rational design of next-generation antifouling materials and offers a framework for realistic model development in complex multicomponent systems.

Adsorption↗
DOE OSTI2021

Electrification of Pakistan's Transport System: Modeling Electric Vehicle Penetration and Energy Supply Chain Impacts

Initial analysis of the National Electric Vehicle Policy (NEVP) of Pakistan was undertaken by an integrated energy planning team from November 2019 through August 2020. This technical reference document seeks to present the background, inputs and assumptions, methodology, and results of the policy analysis in detail. The intended audience is the technical modeler, analyst, or reviewer who seeks to understand specifically the modeling effort within this initial analytical phase, with the end-goal of interpreting, re-producing, or modifying the simulations or extending the models to conduct follow-on studies subsequently. These studies may examine in greater detail the energy sector, impacts to emissions, charging station infrastructure needs, and overall benefits-cost trade of the policy, among other areas of interest.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Machine Learning a Simple Interpretable Short-Range Potential for Silica

A wide array of models, spanning from computationally expensive ab initio methods to a spectrum of force-field approaches, have been developed and employed to probe silica polymorphs and understand growth processes and atomic-level dynamical transitions in silica. However, the quest for a model capable of making accurate predictions with high computational efficiency for various silica polymorphs is still ongoing. Recent developments in short-range machine-learned models, such as GAP and NNPScan, have shown promise in providing reasonable descriptions of silica, but their computational cost remains high compared to force fields such as BKS which are based on simple interpretable functional forms. Here, in this study, we build on the recent success of our reinforcement learning (RL) workflow to derive a new set of optimal parameters for a promising short-range BKS-based model proposed by Soules. We use RL to navigate the eight-dimensional parameter space of the Soules potential using an experimental training data set that includes both local and global structural features from approximately 21 experimentally realized silica polymorphs, including high density phases and porous zeolites. We compare the performance of our machine-learned ML-Soules model with other high quality models including our recent machine-learned parametrization of BKS (ML-BKS), a machine-learned potential (GAP), as well as predictions of ab initio calculations with the highly fidelity SCAN functional. The ML-Soules accurately captures the relative energetic ordering of various polymorphs as well as their structural features at a significantly reduced computational expense. The ML-Soules model also reasonably captures the structure, density, and elastic constants of quartz, as well as metastable silica polymorphs. We further discuss the limitations of the Soules functional form and propose potential enhancements, including the incorporation of additional three-body terms and/or the utilization of different short-ranged functional forms to achieve greater accuracy for both global and local features in the modeling of silica while retaining low computational cost.

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

Finding Liebig’s law of the minimum

Abstract Liebig’s law of the minimum (LLM) is often used to interpret empirical biological growth data and model multiple substrates co‐limited growth. However, its mechanistic foundation is rarely discussed, even though its validity has been questioned since its introduction in the 1820s. Here we first show that LLM is a crude approximation of the law of mass action, the state of art theory of biochemical reactions, and the LLM model is less accurate than two other approximations of the law of mass action: the synthesizing unit model and the additive model. We corroborate this conclusion using empirical data sets of algae and plants grown under two co‐limiting substrates. Based on our analysis, we show that when growth is modeled directly as a function of substrate uptake, the LLM model improperly restricts the organism to be of fixed elemental stoichiometry, making it incapable of consistently resolving biological adaptation, ecological evolution, and community assembly. When growth is modeled as a function of the cellular nutrient quota, the LLM model may obtain good results at the risk of incorrect model parameters as compared to those inferred from the more accurate synthesizing unit model. However, biogeochemical models that implement these three formulations are needed to evaluate which formulation is acceptably accurate and their impacts on predicted long‐term ecosystem dynamics. In particular, studies are needed that explore the extent to which parameter calibration can rescue model performance when the mechanistic representation of a biogeochemical process is known to be deficient.

54 ENVIRONMENTAL SCIENCES↗