Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variability analysis”

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 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Machine learning inversion from scattering for mechanically driven polymers

A machine learning inversion method is developed for analyzing scattering functions of mechanically driven polymers and extracting the corresponding feature parameters, which include energy parameters and conformation variables. The polymer is modeled as a chain of fixed-length bonds constrained by bending energy, and it is subject to external forces such as stretching and shear. We generate a data set consisting of random combinations of energy parameters, including bending modulus, stretching and shear force, along with Monte Carlo-calculated scattering functions and conformation variables such as end-to-end distance, radius of gyration and off-diagonal component of the gyration tensor. The effects of the energy parameters on the polymer are captured by the scattering function, and principal component analysis ensures the feasibility of the machine learning inversion. Finally, we train a Gaussian process regressor using part of the data set as a training set and validate the trained regressor for inversion using the rest of the data. The regressor successfully extracts the feature parameters.

Gaussian process regressors↗
DOE OSTI2024

Advanced Signal Decomposition Analysis and Anomaly Detection in Photovoltaic Systems

With the rapid expansion of large-scale photovoltaic (PV) plants, it is paramount for solar stakeholders to understand the reliability and efficiency of their plants to inform maintenance decisions, increase production, and understand the design factors that impact performance. Diagnosing underperformance in PV plants is challenging due to the relatively few monitoring points with respect to the large geographic footprint of the plant. This work introduces a cutting-edge method that transforms the analysis and management of key factors influencing PV plant performance, including performance loss rate (PLR), recoverable soiling, and major system changes. Identifying these factors is critical for deriving actionable insights. Leveraging advanced analytical techniques such as wavelet transformation, robust regression, and extreme point analysis, this approach provides a nuanced understanding of these factors. This method has been tested across two synthetic datasets and one real dataset, consistently surpassing existing benchmarks by achieving a lower median mean absolute error and reduced error variability across all comparable components.

14 SOLAR ENERGY↗
DOE OSTI2020

Multicarrier Energy Systems: Shaping Our Energy Future

Multicarrier energy systems (MCESs) are characterized by strong coordination in operation and planning across multiple energy vectors and/or sectors to deliver reliable, cost-effective energy services to end users/customers with minimal impact on the environment. They have efficiency and flexibility benefits and are deployed in large and small scales on the supply and demand sides and at the network level but are more complex to control and manage. Here, MCESs are reviewed in the context of future low carbon energy systems based on electrification and very high variable renewable energy penetrations. Fully exploiting these systems requires some cost reductions, more sophisticated operations enabled by standardized communications and control capabilities detailed planning paradigms, and addressing their corresponding economic challenges. All these point toward the direction of analysis, markets, and technology research and development coupled with better policy and regulatory frameworks. One futuristic vision of a very low carbon energy system is proposed that illustrates potential pathways to an MCES-dominated energy future.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Stable isotopes of carbon dioxide (CO2), dissolved inorganic carbon (DIC), water and concentrations of greenhouse gases in a mountainous watershed (Crested Butte, Colorado)

Temperature and precipitation variability in response to climate change affects water cycling of watersheds and can potentially impact water quality, water availability, elemental and molecular fluxes, and biogeochemical processes. Here we report the measurements of light stable isotopes and greenhouse gases analysis to a large multidisciplinary project addressing watershed function. Water hydrogen (H) and oxygen (O) stable isotopes of snow and water saturated soil were measured to constrain source contributions to groundwater and streamflow. In an effort to better constrain existing nitrogen and carbon cycling models in high-altitude ecosystems, carbon dioxide (CO2), nitrous oxide (N2O), methane (CH4) greenhouse gases concentrations and stable isotope signatures of CO2, dissolved inorganic carbon (DIC), and water were measured. This dataset includes one CSV file and one Excel file containing data including stable isotopes of carbon dioxide (CO2), dissolved inorganic carbon (DIC), water and concentrations of greenhouse gases. The CSV and Excel files contain the same information/data but have included both versions for accessibility.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Interpretable Deep Learning for the Earth System with Fractal Nets

Focal Area 3: Explainable AI Our confidence in the projections made by Earth System Models (ESMs) depends on understanding them to be, in some important respects, faithful representations of the Earth system. Here we present an “explainable Artificial Intelligence (AI)” method that allows us to uncover the dynamical structure of the observed and modeled Earth system, discover hidden links across wide spatiotemporal scales, target model development efforts at poorly-represented dynamics, and optimize observed or modeled data collection to maximize predictive information. Science Challenge: Dynamical system science for the Earth system poses unique challenges given the large degree of internal climate variability. Thus, tools that help us understand how ESMs succeed and fail at representing these dynamics are crucial, particularly in relation to the observed system. Furthermore, the computational and memory constraints on ESM data output motivate in situ analysis of ESM dynamics, including automatic detection of dynamical shifts. Also, of key importance are procedures that leverage ESMs to optimize observational campaigns for improving process representation, reducing structural uncertainty and improving model skill.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Autonomous Synthesis of Metastable Materials Using a Modular Mixed-Flow Reactor

Understanding and controlling atomic-level processes at solid-liquid interfaces is key to advancing technologies in energy storage, carbon capture, critical element recovery, and materials synthesis. Many of these processes are dominated by the formation of short-lived intermediate precipitates that determine the final properties of synthesized materials. However, studying these intermediates is challenging due to their sensitivity and the reliance on trial-and-error methods. To address this, we developed an automated variable-volume mixed-flow reactor (MFR) to optimize metastable material synthesis and investigate rapid kinetic processes. This state-of-the-art MFR system, paired with an automated modeling framework, enables efficient synthesis and real-time analysis of transient phases. Benchmarking with advanced capabilities, such as wide-/small-angle X-ray scattering, allows us to resolve fast nucleation and growth dynamics that were previously inaccessible. By combining automation, ML-guided optimization, and tailored kinetic modeling, this approach provides a robust platform for improving material design and achieving precise control over solid-liquid reactions.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Genomic-based biosurveillance for avian influenza: whole genome sequencing from wild mallards sampled during autumn migration in 2022–2023 reveals a high co-infection rate on migration stopover site in Georgia

The Caucasus region, including Georgia, is an important intersection for migratory waterbirds, offering potential for avian influenza virus (AIV) transmission between populations from different geographic areas. In 2022 and 2023, wild ducks were sampled during autumn migration events in Georgia to study the genetic relationships and molecular characteristics of influenza strains. Sequencing and phylogenetic analysis were used to compare the sampled strains to reference sequences from Africa, Asia, and Europe, allowing assessment of genetic relationships and virus transmission between migratory birds. Protein language modeling identified potential co-infections. Of 225 duck samples, 128 tested positive for the influenza M gene. 55 influenza-positive samples underwent whole-genome sequencing, revealing significant diversity. Analysis of the hemagglutinin (HA) segment showed notable differences among subtypes. Most samples were H6N1 and H6N6, but co-infections with combinations like H6H3, N8N1, N6H9, N2N6, and H9H6/N1N2 were also identified. These findings demonstrate the high variability of influenza viruses in migratory waterbirds in Georgia, including a notable rate of co-infections. Some samples exhibited uncommon genetic characteristics compared to other strains from the same year, suggesting Georgia’s role as a mixing vessel for influenza viruses. This facilitates reassortment during co-infections and contributes to the genetic diversity observed across flyways.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

n-Propylbenzene Droplet Combustion: Experiments and Numerical Modeling

n-Propylbenzene (PB, C9H12, 432K boiling point) is important as a constituent in surrogates for some transportation fuels but little is understood about its droplet burning characteristics. We present an experimental and numerical study of a fuel droplet burning under conditions where the sole means of gas transport is fuel evaporation at the droplet surface thereby creating a one-dimensional gas transport that is amenable to detailed numerical modeling. The one-dimensional gas transport that results is amenable to detailed numerical modeling. Initial droplet diameters (Do) ranged from 0.3 mm to about 5 mm to better understand the mechanisms that control burning over such a wide range of Do. For Do < 1 mm droplets were burned while anchored to 14 m SiC fibers within a sealed container in free-fall during which time the droplets were ignited by spark discharge and the burning process recorded by video cameras. Larger and free-floating droplets were studied in the spaced-based platform of the orbiting International Space Station for Do > 1.0 mm which provided unlimited times for experimental observation. Measurements of droplet, flame and soot shell diameters were made and are discussed in the presentation. The experimental results showed formation of optically thick soot clouds developing during burning that required a robust image analysis algorithm to enable measurement of droplet diameter. The measured droplet diameters were well predicted by a numerical model that included soot formation, a PB kinetic mechanism comprising 576 species and 27369 reactions, unsteady gas and liquid transport, variable properties and radiative transport. Flame diameters showed more variability compared to the simulations, depending on how the flame was defined (peak gas temperature and OH concentration). Droplets burned slower as Do was increased and the PB sooting propensity increased with Do. Flame extinction was observed for Do > 3.3 mm indicating a strong effect of radiation as Do was increased compared to the smaller droplet diameters examined where extinction was not experimentally observed. For the largest droplets examined the possibility that the observed extinction is indicative of a transition to cool flame burning is still under investigation.

Guo, Songtao↗
DOE OSTI2026

HARMONY: Large-Scale Architecture Search for Efficient Hybrid Language Models

As large language models scale to trillions of parameters, their computational and memory requirements present critical challenges for efficient training and deployment. While Mixture of Experts (MoE) architectures enable efficient scaling through sparse parameter activation, and state-space models like Mamba offer linear-time complexity, principled methods for combining these paradigms remain undeveloped. We introduce HARMONY (Hybrid Architecture Research for Mamba, Optimized with Neural efficiencY), a multi-objective evolutionary neural architecture search framework for discovering efficient hybrid language models that integrate Transformer attention mechanisms, Mixture-of-Experts routing, and Mamba state-space components. Through large-scale distributed search using 16,384 MI250X GPUs on the Frontier supercomputer, HARMONY explores a comprehensive design space encompassing six attention variants (MHA, MQA, GQA, MLA, SWA, and Mamba-2), variable MoE configurations with both routed and shared experts, and extensive Mamba hyperparameters. Our framework discovers heterogeneous architectures that balance training performance with computational efficiency through multi-objective optimization incorporating latency penalties and fitness-based selection. Analysis of discovered architectures reveals that optimal hybrid designs favor heterogeneous component mixing rather than homogeneous patterns, with Mamba-2 and Multi-Head Latent Attention (MLA) emerging as preferred mechanisms. Discovered architectures demonstrate superior training efficiency: our best configuration achieves a final perplexity of 1.0874 with 2.38B parameters while processing 4,320 tokens/second, outperforming significantly larger manually designed models. Full-scale evaluation shows HARMONY's top architectures achieve better loss trajectories than equivalently-sized models using state-of-the-art configurations including Mixtral, Jamba, and Samba. Additionally, we demonstrate 91% weak scaling efficiency when training discovered 36B-parameter models across 1,024 GPUs. HARMONY is released as an open framework with comprehensive tools for building and training hybrid models using expert-data-pipeline parallelism, democratizing access to automated architecture design for next-generation language models.

Herron, Emily [ORNL] (ORCID:0000000273008172)↗
DOE OSTI2020

A new phenomenological model to describe root-soil interactions based on percolation theory

In his paper on net primary productivity of terrestrial communities predicted from climatological data, Rosenzweig (1968) argued that variability in productivity is well accounted for by (evapo)-transpiration, and that water from transpiration is, on global scales, the most variable component in the photosynthesis reaction. The goal of this paper is to investigate whether variability in plant growth on local scales and within species is primarily related to transpiration under several scenarios including different terrain curvature, slope aspect, soil characteristics, and climate ranges. Here, we test the hypothesis that this relationship exists because root growth into the surface soil layers (0–2 m) tends to follow paths with minima in resistance, which in turn maximizes water flow and nutrient delivery rates that regulate growth. The set of all connected paths with individual pore-to-pore flow resistances less than a critical, percolating, value forms a cluster with mass fractal dimensionality, d f . We propose that roots follow paths through the 2D percolation cluster, defining the set of all optimal flow paths, making the 2D value of d f from percolation relevant to root fractal dimensionality. The tortuosity of such optimal paths as defined in percolation theory should then relate root length to root radial extent, linking the parameters of root tortuosity and plant productivity. Our analysis of large data sets across species implies that root radial extent and tree height are both proportional to cumulative transpiration until trees approached maximum height, and their growth rates are proportional to the transpiration rate, not to the moisture content. Local variations in tree height as functions of the variables investigated appear generally consistent with deduced variations in transpiration. Here this correlation is investigated more closely in the context of studies addressing individual tree species.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Network analysis of earthquake ground motion spatial correlation: a case study with the San Jacinto seismic nodal array

SUMMARY The spatial correlation of earthquake ground motion intensity can be measured from strong motion data; however, the data used in past studies is sparsely sampled in space, and only the interstation distance was considered as a correlation variable. These limitations mean that we have only weak constraints on the true correlation structure of ground motion and that potentially important aspects of spatial correlation are unconstrained. In this study, we combine a large-N seismic array and graph analytics to explore this issue at a local scale using small local and regional earthquakes. Our result suggests site conditions, and how they interact with the incident seismic wavefield, strongly condition the spatial correlation of ground motion. Future progress in characterizing ground motion spatial variability will require dense wavefield measurements, either through nodal deployments, or perhaps distributed acoustic sensing measurements, of seismic wavefields. Aftershock sequences of major earthquakes would provide particularly data-rich targets of opportunity.

58 GEOSCIENCES↗
DOE OSTI2023

Voltage-Based Strategies for Preventing Battery Degradation under Diverse Fast-Charging Conditions

Maintaining safe operating conditions is a key challenge for high-performance lithium-ion battery applications. The lithium-plating reaction remains a risk during charging, but limited studies consider the highly variable charging conditions possible in commercial cells. Here we combine pseudo-2D electrochemical modeling with data visualization methods to reveal important relationships between the measurable cell voltage and difficult-to-predict Li-plating onset criteria. An extensively validated model is used to compute Li plating for thousands of multistep charging conditions spanning diverse rates, temperatures, states-of-charge, and cell aging. Here we observe an empirical cell operating voltage limit below which plating does not occur across all conditions, and this limit varies with the battery state-of-charge and aging. A model sensitivity analysis also indicates that, when comparing two charging voltage profiles, the capacity difference at 4.0 V correlates well with the difference in the plating onset capacity. These results encourage simple strategies for Li-plating prevention that are complementary to existing battery controls.

25 ENERGY STORAGE↗
DOE OSTI2023

BeyondPlanck: I. Global Bayesian analysis of the Planck Low Frequency Instrument data

We describe the BEYONDPLANCK project in terms of our motivation, methodology, and main products, and provide a guide to a set of companion papers that describe each result in more detail. Building directly on experience from ESA’s Planck mission, we implemented a complete end-to-end Bayesian analysis framework for the Planck Low Frequency Instrument (LFI) observations. The primary product is a full joint posterior distribution P(ω | d), where ω represents the set of all free instrumental (gain, correlated noise, bandpass, etc.), astrophysical (synchrotron, free-free, thermal dust emission, etc.), and cosmological (cosmic microwave background – CMB – map, power spectrum, etc.) parameters. Some notable advantages of this approach compared to a traditional pipeline procedure are seamless end-to-end propagation of uncertainties; accurate modeling of both astrophysical and instrumental effects in the most natural basis for each uncertain quantity; optimized computational costs with little or no need for intermediate human interaction between various analysis steps; and a complete overview of the entire analysis process within one single framework. As a practical demonstration of this framework, we focus in particular on low-ℓ CMB polarization reconstruction with Planck LFI. In this process, we identify several important new effects that have not been accounted for in previous pipelines, including gain over-smoothing and time-variable and non-1/f correlated noise in the 30 and 44 GHz channels. Modeling and mitigating both previously known and newly discovered systematic effects, we find that all results are consistent with the ΛCDM model, and we constrained the reionization optical depth to τ = 0.066 ± 0.013, with a low-resolution CMB-based χ 2 probability to exceed of 32%. This uncertainty is about 30% larger than the official pipelines, arising from taking a more complete instrumental model into account. The marginal CMB solar dipole amplitude is 3362.7 ± 1.4 μK, where the error bar was derived directly from the posterior distribution without the need of any ad hoc instrumental corrections. We are currently not aware of any significant unmodeled systematic effects remaining in the Planck LFI data, and, for the first time, the 44 GHz channel is fully exploited in the current analysis. We argue that this framework can play a central role in the analysis of many current and future high-sensitivity CMB experiments, including LiteBIRD, and it will serve as the computational foundation of the emerging community-wide COSMOGLOBE effort, which aims to combine state-of-the-art radio, microwave, and submillimeter data sets into one global astrophysical model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Linking DSC/TGA to Cell Levels: Energetics, Evolved Gases, and Thermal Safety of NMC811‐Graphite Micro‐Cell

Thermochemical characterization of battery materials links intrinsic material properties to decomposition pathways, heat generation, and gas evolution that govern performance and safety. Despite extensive work on NMC811-Graphite, variability across partial configurations and the limited adoption of micro-cell architectures (cathode+anode+electrolyte+separator) hinder robust cell-scale interpretation. Accordingly, this work establishes a bottom-up, component-resolved methodology integrating DSC/TGA, evolved gas analysis (EGA), and in situ XRD to link decomposition pathways and energy release across partial and micro-cell configurations, providing a transferable assessment of safety and stability in emerging chemistries. In separator-free configurations, the gas–solid reaction between cathode-evolved O 2 and anode-leached Li dominates the net heat release (1139 J g −1 ). In contrast, in the micro-cell configuration, the separator hinders O 2 transport and alters the timing and pathways of other reactions, and reduces the net energy release to 618 J g −1 . Energy release was organized into defined temperature windows that provide a framework for a thermodynamic model combining quantified gas evolution with selected decomposition pathways and effective reaction enthalpies to estimate net specific energy release, with agreement between DSC and cell-level tests. Ex situ XPS of heat-treated samples extends post-mortem analysis to thermal-abuse regimes, supporting key pathway elements.

25 ENERGY STORAGE↗
DOE OSTI2021

Daily grass reference evapotranspiration with Meteosat Second Generation shortwave radiation and reference ET products

Here we assess the accuracy of estimating daily grass reference evapotranspiration (PM-ET o ) using daily shortwave radiation (R s ) and reference evapotranspiration (ET REF ) products provided by the Meteosat Second Generation (MSG) geostationary satellite delivered by the Satellite Applications Facility on Land Surface Analysis (LSA-SAF) framework. The accuracy of using reanalysis ERA5 shortwave radiation data (R s ERA5 ) provided by the European Center for Medium-Range Weather Forecasts (ECMWF) is also evaluated. The assessments were performed using observed weather variables at 37 weather stations distributed across continental Portugal, where climate conditions range from semi-arid to humid, and 12 weather stations located in Azores islands, characterized by humid, windy and often cloudy conditions. This study&#x27;s use of data from a variety of climate conditions contributed to a unique and innovative assessment of the usability of LSA-SAF and ERA5 products for ET o estimation. The first assessment focused on comparing LSA-SAF estimates of R s (R s LSA-SAF ) against ground stations (R s ground ). The results showed a good matching between the two R s data sets for continental Portugal but a tendency for R s LSA-SAF to under-estimate R s ground in the cloudy islands of Azores. ET o values computed using R s LSA-SAF data and observed temperature, humidity and wind speed (ET o LSA-SAF ) were then compared with PM-ET o estimates with ground-based data, which were used as benchmark; input data of temperature and humidity needed for PM-ET o were quality checked for surface aridity effects. It was observed that ET o LSA-SAF is strongly correlated with PM-ET o (R 2 > 0.97) for most locations in continental Portugal, with regression coefficient of a linear regression forced to the origin ranging between 0.95 and 1.05, mean root mean square error (RMSE) of 0.13 mm d –1 , and Nash and Sutcliff efficiency of modeling (EF) above 0.95. For most Azores locations, ET o LSA-SAF over-estimated PM-ET o . This is likely a consequence of the high spatio-temporal heterogeneity of weather conditions that occur in these oceanic islands together with the different footprints of satellite (averaged over the pixel) and station observations. Reanalysis ERA5 shortwave radiation data presented similar behavior to the LSA-SAF products, however with slightly lower accuracy. The daily LSA-SAF ET REF product (ET REF LSA-SAF ) was assessed and results have shown a good accuracy of this product, with acceptable RMSE and high EF values, for continental Portugal but a low accuracy for the Azores islands. A simplified bias correction approach was shown to improve both ET o derived from the LSA-SAF products, namely for Azores stations, which seem to be representative of smaller areas. The use of the FAO-PM temperature approach (PMT) was also assessed using the R s LSA-SAF and R s ERA5 data, which showed a superiority of the LSA-SAF product for ET o estimations (ET o PMT LSA-SAF ). No significant differences (p < 0.05) were observed in terms of the median value of the RMSE when adopting ET o PMT and ET REF LSA-SAF . Differently, results showed that using the R s LSA-SAF in the PMT approach (ET o PMT LSA-SAF ) produces significantly better RMSE results than ET o PMT and ET REF LSA-SAF . Overall, the performed assessment allows concluding that the use of R s LSA-SAF, and to a lesser extent the use of the R s ERA5 , highly improves the accuracy of computation of ET o when R s observations are not available, including when only temperature data are accessible. The use of the ET REF LSA-SAF product is a good alternative when observed weather data are not available.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Detailed biomass fast pyrolysis kinetics integrated to computational fluid dynamic (CFD) and discrete element modeling framework: Predicting product yields at the bench-scale

Fast pyrolysis is an intricate process due to the variability and anisotropy of lignocellulosic biomass and the complicated chemistry and physics during conversion in a bubbling fluidized bed reactor (BFBR). The complexity of biomass fast pyrolysis lends itself well to computational fluid dynamics (CFD) and discrete element (DEM) analysis, which promises to reduce experimental time and its associated cost. This work investigated switchgrass fast pyrolysis simulated by computational fluid dynamics coupled with a discrete element method to track individual reacting biomass particles throughout a bench-scale BFBR reactor. We accounted for the fast pyrolysis chemistry through a comprehensive reaction scheme with secondary cracking reactions. We performed a three-step reduction for secondary cracking reactions to convert the full cracking scheme into a reduced scheme easily incorporated into our model. We assessed the impact of operational conditions on the steady-state yields of liquid bio-oil, non-condensable gases (NCG), at 550 °C over a range of fluidization numbers (2 – 6 Umf), reported as a ratio to the minimum fluidization velocity (Umf). At steady-state, the volatile bio-oil yield had a range of 49.3–50.4 wt%. Levoglucosan was the primary volatile component present with 21 wt% of the bio-oil while water was the second largest with 20 wt%. The reduction of the secondary reaction schemes did not appreciably affect the overall yields of switchgrass pyrolysis compared to the full secondary scheme.

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

Subsurface Nitrogen Dissociation Kinetics in Lithium Metal from Metadynamics

The dissociation of molecular nitrogen in lithium is of interest for several promising technologies, such as the catalytic synthesis of ammonia in ambient or mild conditions. In this work we simulate nitrogen dissociation in the lithium BCC (110) surface at ambient and elevated temperatures using density functional theory (DFT) metadynamics simulations. The rate constants at temperatures of 300, 400, and 500 K are calculated by statistical analysis of the reaction time distributions from the accelerated simulations. This approach finds and estimates rate constants for transition pathways out of the initial state; the required input is the stable initial state and a reasonable choice of collective variable. A single collective variable is used in this case: the N–N distance. The results are robust to changes in metadynamics parameters, and the reaction time distributions follow the expected exponential distribution. We show that the metadynamics-derived rate constants are in agreement with results from the conventional harmonic approximation approach using a climbing image nudged elastic band (NEB) transition state search. The reaction barriers from metadynamics and the NEB/harmonic approximation agree to within 0.02–0.04 eV at all temperatures studied. This paper demonstrates that the harmonic approximation provides an accurate description of the rate constants for nitrogen dissociation in lithium metal, even at temperatures near or above the melting point of lithium, lending credence to previous and future theoretical studies using this approximation. Moreover, this work demonstrates a step toward the automated exploration and discovery of reaction mechanisms and associated rate constants for elementary surface-catalyzed reactions using DFT-based metadynamics.

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

Dynamics of Fungal and Bacterial Biomass Carbon in Natural Ecosystems: Site-level Applications of the CLM-Microbe Model

Explicitly representing microbial processes has been recognized as a key improvement to Earth system models for the realistic projections of soil carbon (C) and climate dynamics. The CLM-Microbe model builds upon the CLM4.5 and explicitly represents two major soil microbial groups, fungi and bacteria. Based on the compiled time-series data of fungal (FBC) and bacterial (BBC) biomass C from nine biomes, we parameterized and validated the CLM-Microbe model, and further conducted sensitivity analysis and uncertainty analysis for simulating C cycling. The model performance was evaluated with mean absolute error (MAE), root mean square error (RMSE), and coefficient of determination (R2) for relative change in FBC and BBC. The CLM-Microbe model is able to reasonably capture the seasonal dynamics of FBC and BBC across biomes, particularly for tropical/subtropical forest, temperate broadleaf forest, and grassland, with MAE < 0.49 for FBC and <0.36 for BBC and RMSE <0.52 FBC and <0.39 for BBC, while R2 values are relatively smaller in some biomes (e.g., shrub) due to small sample sizes. We found good consistencies between simulated and observed FBC (R 2 =0.70, P<0.001) and BBC (R 2 =0.26, P<0.05) on average across biomes, but the model is not able to fully capture the large variation in observed FBC and BBC. Sensitivity analysis shows the most critical parameters are turnover rate, carbon-to-nitrogen ratio of fungi and bacteria, and microbial assimilation efficiency. This study confirms that the explicit representation of soil microbial mechanisms enhances model performance in simulating C variables such as heterotrophic respiration and soil organic C density. The further application of the CLM-Microbe model would deepen our understanding of microbial contributions to the global C cycle.

54 ENVIRONMENTAL SCIENCES↗