SEARCH · Engineering Papers

Results for “parameterization”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 937 records · Page 52

All sources NASA NTRS DOE OSTI

DOE OSTI2018

A bulk adjusted linear combination of atomic orbitals (BA-LCAO) approach for nanoparticles

We describe a bulk adjusted linear combination of atomic orbitals (BA-LCAO) approach for nanoparticles. In this method, we apply a many-body scaling function (in similar manner as in the environment-modified total energy based tight-binding method) to the DFT-derived diatomic AO interaction potentials (like in the conventional orbital-based density-functional tight binding approach) strictly according to atomic valences acquired naturally in a bulk structure. This modification, (a) facilitates all atom orbital-based electronic structure calculations of charge carrier dynamics in nanoscale structures with a molecular acceptor, and (b) allows to closely match high-level density functional calculation data (previously adjusted to the available experimental findings) for bulk structures. To advance practical application of the BA-LCAO approach we parameterize the Hamiltonian of wurtzite CdSe by fitting its band structure to a high-level DFT reference, corrected for experimentally measured band edges. Here, unlike in conventional DFTB approach, we: (1) use hydrogen-like AOs for the basis as exact atomic eigenfunctions, while orbital energies of which are taken from experimentally measured ionization potentials, and (2) parameterize the many-body scaling functions rather than the atomic wavefunctions. Development of this approach and parameters is guided by our aim to devise a method capable of simultaneously treating the problems of (i) interfacial electron/hole transfer between finite, variable size nanoparticles and electron scavenging molecules, and (ii) high-energy electronic transitions (Auger transitions) that mediate multi-exciton decay in quantum dots. Electronic structure results are described for CdSe quantum dots of various sizes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2013

Implications of Low Volatility SOA and Gas-Phase Fragmentation Reactions on SOA Loadings and their Spatial and Temporal Evolution in the Atmosphere

Recent laboratory and field measurements by a number of groups show that secondary organic aerosol (SOA) evaporates orders of magnitude slower than traditional models assume. In addition, chemical transport models using volatility basis set (VBS) SOA schemes neglect gas-phase fragmentation reactions, which are known to be extremely important. In this work, we present modeling studies to investigate the implications of non-evaporating SOA and gas-phase fragmentation reactions. Using the 3-D chemical transport model, WRF-Chem, we show that previous parameterizations, which neglect fragmentation during multi-generational gas-phase chemistry of semi-volatile/inter-mediate volatility organics ("aging SIVOC"), significantly over-predict SOA as compared to aircraft measurements downwind of Mexico City. In sharp contrast, the revised models, which include gas-phase fragmentation, show much better agreement with measurements downwind of Mexico City. We also demonstrate complex differences in spatial SOA distributions when we transform SOA to non-volatile secondary organic aerosol (NVSOA) to account for experimental observations. Using a simple box model, we show that for same amount of SOA precursors, earlier models that do not employ multi-generation gas-phase chemistry of precursors ("non-aging SIVOC"), produce orders of magnitude lower SOA than "aging SIVOC" parameterizations both with and without fragmentation. In addition, traditional absorptive partitioning models predict almost complete SOA evaporation at farther downwind locations for both "non-aging SIVOC" and "aging SIVOC" with fragmentation. In contrast, in our revised approach, SOA transformed to NVSOA implies significantly higher background concentrations as it remains in particle phase even under highly dilute conditions. This work has significant implications on understanding the role of multi-generational chemistry and NVSOA formation on SOA evolution in the atmosphere.

Shrivastava, ManishKumar B.↗

DOE OSTI2022

Effects of random forest modeling decisions on biogeochemical time series predictions

Abstract Random forests (RF) are an increasingly popular machine learning approach used to model biogeochemical processes in the Earth system. While RF models are robust to many assumptions that complicate deterministic models, there are several important parameterization decisions for appropriate use and optimal model fit. We explored the role that parameter decisions, including training/testing data splitting strategies, variable selection, and hyperparameters play on RF goodness‐of‐fit by constructing models using 1296 unique parameter combinations to predict concentrations of nitrate, a key nutrient for biogeochemical cycling in aquatic ecosystems. Models were built on long‐term, publicly available water quality and meteorology time series collected by the National Estuarine Research Reserve monitoring network for two contrasting ecosystems representing freshwater and brackish estuaries. We found that accounting for temporal dependence when splitting data into training and testing subsets was key for avoiding over‐estimation of model predictive power. In addition, variable selection, the ratio of training to testing data, and to a lesser degree, variables per split and number of trees, were significant parameters for optimizing RF goodness‐of‐fit. We also explored how model parameter decisions influenced interpretation of the relative importance of predictors to the model, and model predictor‐dependent variable relationships, with results suggesting that both data structure and model parameterization influence these factors. Because much of the current RF literature is written for the computational and statistical science communities, the primary goal of this study is to provide guidelines for aquatic scientists new to machine learning to apply RF techniques appropriately to aquatic biogeochemical datasets.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2017

Quantitative impact of changes in marrow cellularity, skeletal size, and bone mineral density on active marrow dosimetry based upon a reference model

Purpose The hematopoietically active tissues of skeletal bone marrow are a prime target for computational dosimetry given potential risks of leukemia and, at higher dose levels, acute marrow toxicity. The complex three‐dimensional geometry of trabecular spongiosa, however, complicates schema for dose assessment in such a way that only a few reference skeletal models have been developed to date, and which are based upon microimaging of a limited number of cadaveric bone spongiosa cores. The question then arises as to what degree of accuracy is achievable from reference skeletal dose models when applied to individual patients or specific exposed populations? Methods Patient variability in marrow dosimetry were quantified for three skeletal sites — the ribs, lumbar vertebrae, and cranium — for the beta‐emitters 45 Ca, 153 Sm, and 90 Y, and the alpha‐particle emitters 223 Ra, 219 Rn, and 215 Po, the latter two being the immediate progeny of the former. For each radionuclide and bone site, three patient parameters were altered from their values in the reference model: (1) bone size as a surrogate for patient stature, (2) marrow cellularity as a surrogate for age‐ or disease‐related changes in marrow adiposity, and (3) the trabecular bone volume fraction as a surrogate for bone mineral density. Marrow dose variability is expressed as percent differences in the radionuclide S value given by the reference model and the patient‐parameterized model. The impact of radionuclide biokinetics on marrow dosimetry was not considered. Results Variations in overall bone size play a very minor role in active marrow dose variability. Marrow cellularity is a significant factor in dose variability for active marrow self‐irradiation, but it plays no role for radionuclides localized to the trabecular bone matrix. Variations in trabecular bone volume fractions impact the active marrow dose variability for short‐range particle emitters 45 Ca, 223 Ra, 219 Rn, and 215 Po in the vertebrae and ribs, skeletal sites with small spongiosa proportions of trabecular bone. In the cranium, with its relative high proportion of trabecular bone, significant differences in marrow dosimetry from the reference model were noted for all radionuclides. Conclusions Skeletal models of active marrow dosimetry should be more fully parameterized to permit closer matching to patient bone density and marrow cellularity, particularly when considering short‐range particle emitters localized to either the bone trabeculae or active marrow, respectively.

Geyer, Amy M.↗

DOE OSTI2023

Diurnal cycle of precipitation over the tropics and central United States: intercomparison of general circulation models

Diurnal precipitation is a fundamental mode of variability that climate models have difficulty in accurately simulating. Here, in this work, the diurnal cycle of precipitation (DCP) in participating climate models from the Global Energy and Water Exchanges' DCP project is evaluated over the tropics and central United States. Common model biases such as excessive precipitation over the tropics, too frequent light-to-moderate rain, and the failure to capture propagating convection in the central United States still exist. Over the central United States, the issues of too weak rainfall intensity in climate runs is well improved in their hindcast runs with initial conditions from numerical weather prediction analyses. But the improvement is minimal over the central Amazon. Incorporating the role of the large-scale environment in convective triggering processes helps resolve the phase-locking issue in many models where precipitation often incorrectly peaks near noon due to maximum insolation over land. Allowing air parcels to be lifted above the boundary layer improves the simulation of nocturnal precipitation which is often associated with the propagation of mesoscale systems. Including convective memory in cumulus parameterizations acts to suppress light-to-moderate rain and promote intense rainfall; however, it also weakens the diurnal variability. Simply increasing model resolution (with cumulus parameterizations still used) cannot fully resolve the biases of low-resolution climate models in DCP. The hierarchy modeling framework from this study is useful for identifying the missing physics in models and testing new development of model convective processes over different convective regimes.

58 GEOSCIENCES↗

DOE OSTI2017

Twofold symmetries of the pure gravity action

Here, we recast the action of pure gravity into a form that is invariant under a twofold Lorentz symmetry. To derive this representation, we construct a general parameterization of all theories equivalent to the Einstein-Hilbert action up to a local field redefinition and gauge fixing. We then exploit this freedom to eliminate all interactions except those exhibiting two sets of independently contracted Lorentz indices. The resulting action is local, remarkably simple, and naturally expressed in a field basis analogous to the exponential parameterization of the nonlinear sigma model. The space of twofold Lorentz invariant field redefinitions then generates an infinite class of equivalent representations. By construction, all off-shell Feynman diagrams are twofold Lorentz invariant while all on-shell tree amplitudes are automatically twofold gauge invariant. We extend our results to curved spacetime and calculate the analogue of the Einstein equations. Finally, while these twofold invariances are hidden in the canonical approach of graviton perturbation theory, they are naturally expected given the double copy relations for scattering amplitudes in gauge theory and gravity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

SHAPER: can you hear the shape of a jet?

The identification of interesting substructures within jets is an important tool for searching for new physics and probing the Standard Model at colliders. Many of these substructure tools have previously been shown to take the form of optimal transport problems, in particular the Energy Mover’s Distance (EMD). In this work, we show that the EMD is in fact the natural structure for comparing collider events, which accounts for its recent success in understanding event and jet substructure. We then present a Shape Hunting Algorithm using Parameterized Energy Reconstruction (SHAPER), which is a general framework for defining and computing shape-based observables. SHAPER generalizes N-jettiness from point clusters to any extended, parametrizable shape. This is accomplished by efficiently minimizing the EMD between events and parameterized manifolds of energy flows representing idealized shapes, implemented using the dual-potential Sinkhorn approximation of the Wasserstein metric. We show how the geometric language of observables as manifolds can be used to define novel observables with built-in infrared-and-collinear safety. We demonstrate the efficacy of the SHAPER framework by performing empirical jet substructure studies using several examples of new shape-based observables.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2017

Studies of Zγ production in association with a high-mass dijet system in pp collisions at $\sqrt{s}=8$ TeV with the ATLAS detector

The production of a Z boson and a photon in association with a high-mass dijet system is studied using 20.2 fb -1 of proton-proton collision data at a centre-of-mass energy of $\sqrt{s}$ = 8 TeV recorded with the ATLAS detector in 2012 at the Large Hadron Collider. Final states with a photon and a Z boson decaying into a pair of either electrons, muons, or neutrinos are analysed. Electroweak and total pp → Zγjj cross-sections are extracted in two fiducial regions with different sensitivities to electroweak production processes. Quartic couplings of vector bosons are studied in regions of phase space with an enhanced contribution from pure electroweak production, sensitive to vector-boson scattering processes VV → Zγ. Finally, no deviations from Standard Model predictions are observed and constraints are placed on anomalous couplings parameterized by higher-dimensional operators using effective field theory.

43 PARTICLE ACCELERATORS↗

DOE OSTI2021

Conformal symmetries for extremal black holes with general asymptotic scalars in STU supergravity

We present a construction of the most general BPS black holes of STU supergravity (N = 2 supersymmetric D = 4 supergravity coupled to three vector super-multiplets) with arbitrary asymptotic values of the scalar fields. These solutions are obtained by acting with a subset of the global symmetry generators on STU BPS black holes with zero values of the asymptotic scalars, both in the U-duality and the heterotic frame. The solutions are parameterized by fourteen parameters: four electric and four magnetic charges, and the asymptotic values of the six scalar fields. We also present BPS black hole solutions of a consistently truncated STU supergravity, which are parameterized by two electric and two magnetic charges and two scalar fields. These latter solutions are significantly simplified, and are very suitable for further explicit studies. We also explore a conformal inversion symmetry of the Couch-Torrence type, which maps any member of the fourteen-parameter family of BPS black holes to another member of the family. Furthermore, these solutions are expected to be valuable in the studies of various swampland conjectures in the moduli space of string compactifications.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2014

Infrared consistency and the weak gravity conjecture

The weak gravity conjecture (WGC) asserts that an Abelian gauge theory coupled to gravity is inconsistent unless it contains a particle of charge q and mass m such that q ≥ m/m Pl . This criterion is obeyed by all known ultraviolet completions and is needed to evade pathologies from stable black hole remnants. In this paper, we explore the WGC from the perspective of low-energy effective field theory. Below the charged particle threshold, the effective action describes a photon and graviton interacting via higher-dimension operators. We derive infrared consistency conditions on the parameters of the effective action using i ) analyticity of light-by-light scattering, ii ) unitarity of the dynamics of an arbitrary ultraviolet completion, and iii ) absence of superluminality and causality violation in certain non-trivial backgrounds. For convenience, we begin our analysis in three spacetime dimensions, where gravity is non-dynamical but has a physical effect on photon-photon interactions. We then consider four dimensions, where propagating gravity substantially complicates all of our arguments, but bounds can still be derived. Operators in the effective action arise from two types of diagrams: those that involve electromagnetic interactions (parameterized by a charge-to-mass ratio q/m) and those that do not (parameterized by a coefficient γ). In conclusion, infrared consistency implies that q/m is bounded from below for small γ.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

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 OSTI2017

The atmospheric hydrologic cycle in the ACME v0.3 model

We examine the global water cycle characteristics in the Accelerated Climate Modeling for Energy v0.3 model (a close relative to version 5.3 of the Community Atmosphere Model) in atmosphere-only simulations spanning the years 1980–2005. We evaluate the simulations using a broad range of observational and reanalysis datasets, examine how the simulations change when the horizontal resolution is increased from 1° to 0.25, and compare the simulations against models participating in the the Atmosphere Model Intercomparison Project of the 5th Coupled Model Intercomparison Project (CMIP5). Particular effort has been made to evaluate the model using the best available observational estimates and verifying model biases with additional datasets when differences are known to exist among the observations. Regardless of resolution, the model exhibits several biases: global-mean precipitation, evaporation, and precipitable water are too high, light precipitation occurs too frequently, and the atmospheric residence time of water is too short. Many of these biases are shared by the multi-model mean climate of models participating in CMIP5. The reasons behind regional biases in precipitation are discussed by examining how different fields, such as local evaporation and transport of water vapor, contribute to the bias. Although increasing the horizontal resolution does not drastically change the water cycle, it does lead to a few differences: an increase in global mean precipitation rate, an increase in the fraction of total precipitation that falls over land, more frequent heavy precipitation (>30 mm day -1 ), and a decrease in precipitable water. One of the most notable changes is the shift of precipitation produced by the convective parameterization to that produced by the large-scale microphysics parameterization. We analyze how changes in moisture and circulation with resolution contribute to this shift in the precipitation partitioning. Because changing horizontal resolution requires some re-tuning, the effect of that tuning was evaluated by performing an additional simulation at 1 but using the tunings from the 0.25 simulation. In conclusion, the evaluation shows that the more frequent heavy precipitation, the decrease in precipitable water, and the shift from convective to large-scale precipitation are predominantly due to resolution changes, while tuning changes have a major influence on the global mean precipitation and the land/ ocean partitioning of precipitation.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Region and cloud regime dependence of parametric sensitivity in E3SM atmosphere model

Abstract The Department of Energy (DOE)’s Energy Exascale Earth System Model (E3SM), including its atmosphere model (EAM), has many relatively new features. In a previous study we conducted a systematic parametric sensitivity analysis for EAM based on short, perturbed parameter ensemble (PPE) simulations, mainly focusing on global mean climate features and metrics. While parameter values in global climate models are generally invariant in space and time, model response to parameters perturbation may vary by regions and climate regimes, which motivates the need to better understand the EAM model behaviors and physics at regional scale and process level. In this study, using the same set of PPE simulations and a similar sensitivity analysis framework, we identify parameters that cause largest sensitivities over different regions and compare model responses in fast atmospheric processes to the parameters across different cloud regimes for several important cloud-related fidelity metrics. We find that cloud forcing has opposite response to some parameters over mid-latitude vs. tropical land. We also analyze how the parametric sensitivity varies as stratocumulus transitions to shallow convection and to deep convection over ocean. Low cloud forcing and shortwave cloud forcing in the subtropical eastern Pacific are most sensitive to the parameters controlling the width of the probability density function (PDF) of the subgrid vertical velocity ( w’ ) ( gamma ) and the damping of the w’ skewness ( c8 ) near the coast but become more sensitive to the parameter affecting the damping of the w’ variance ( c1 ) further offshore. Detailed interpretation of the spatial dependence of parametric sensitivity is provided. We also investigate how the parametric sensitivity evolves with prediction duration. This study improves our process-level understanding of cloud physics and parameterization and provides insights for developing more advanced regime-aware parameterization schemes in global climate model.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2019

Non-stationary Geostatistical Modeling: A Case Study Comparing LVA Estimation Frameworks

Incorporating locally varying anisotropy (LVA) in geostatistical modeling improves estimates for structurally complex domains where a single set of anisotropic parameters modeled globally do not account for all geological features. In this work, the properties of two LVA-geostatistical modeling frameworks are explored through application to a complexly folded gold deposit in Ghana. The inference of necessary parameters is a significant requirement of geostatistical modeling with LVA; this work focuses on the case where LVA orientations, derived from expert geological interpretation, are used to improve the grade estimates. The different methodologies for inferring the required parameters in this context are explored. The results of considering different estimation frameworks and alternate methods of parameterization are evaluated with a cross-validation study, as well as visual inspection of grade continuity along select cross sections. Results show that stationary methodologies are outperformed by all LVA techniques, even when the LVA framework has minimal guidance on parameterization. Findings also show that additional improvements are gained by considering parameter inference where the LVA orientations and point data are used to infer the local range of anisotropy. Considering LVA for geostatistical modeling of the deposit considered in this work results in better reproduction of curvilinear geological features.

58 GEOSCIENCES↗

DOE OSTI2016

A Bayesian Machine Learning Model for Estimating Building Occupancy from Open Source Data

Understanding building occupancy is critical to a wide array of applications including natural hazards loss analysis, green building technologies, and population distribution modeling. Due to the expense of directly monitoring buildings, scientists rely in addition on a wide and disparate array of ancillary and open source information including subject matter expertise, survey data, and remote sensing information. These data are fused using data harmonization methods which refer to a loose collection of formal and informal techniques for fusing data together to create viable content for building occupancy estimation. In this paper, we add to the current state of the art by introducing the Population Data Tables (PDT), a Bayesian based informatics system for systematically arranging data and harmonization techniques into a consistent, transparent, knowledge learning framework that retains in the final estimation uncertainty emerging from data, expert judgment, and model parameterization. PDT probabilistically estimates ambient occupancy in units of people/1000ft2 for over 50 building types at the national and sub-national level with the goal of providing global coverage. The challenge of global coverage led to the development of an interdisciplinary geospatial informatics system tool that provides the framework for capturing, storing, and managing open source data, handling subject matter expertise, carrying out Bayesian analytics as well as visualizing and exporting occupancy estimation results. We present the PDT project, situate the work within the larger community, and report on the progress of this multi-year project.Understanding building occupancy is critical to a wide array of applications including natural hazards loss analysis, green building technologies, and population distribution modeling. Due to the expense of directly monitoring buildings, scientists rely in addition on a wide and disparate array of ancillary and open source information including subject matter expertise, survey data, and remote sensing information. These data are fused using data harmonization methods which refer to a loose collection of formal and informal techniques for fusing data together to create viable content for building occupancy estimation. In this paper, we add to the current state of the art by introducing the Population Data Tables (PDT), a Bayesian model and informatics system for systematically arranging data and harmonization techniques into a consistent, transparent, knowledge learning framework that retains in the final estimation uncertainty emerging from data, expert judgment, and model parameterization. PDT probabilistically estimates ambient occupancy in units of people/1000 ft 2 for over 50 building types at the national and sub-national level with the goal of providing global coverage. The challenge of global coverage led to the development of an interdisciplinary geospatial informatics system tool that provides the framework for capturing, storing, and managing open source data, handling subject matter expertise, carrying out Bayesian analytics as well as visualizing and exporting occupancy estimation results. We present the PDT project, situate the work within the larger community, and report on the progress of this multi-year project.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2018

Predicting the Electric Field Effect on the Lateral Interactions Between Adsorbates: O/Fe(100) from First Principles

A density functional theory parameterized lattice gas cluster expansion of oxygen on an Fe(100) surface is developed in the presence and absence of an applied positive and negative external electric field, characterizing the heterogeneity in oxygen's surface distribution and the effect of an external electric field on the lateral interactions between the adsorbates. We show that the presence of a negative electric field tends to weaken both attractive and repulsive interactions while the positive electric field tends to strengthen these interactions, altering surface distributions and ground state configurations. Since lateral interactions have been shown to play a critical role in defining catalytic behavior, the application of an applied electric field has the potential to be a useful tool in adjusting critical chemical properties of Fe-based hydrodeoxygenation catalysts, by decreasing the repulsive interactions between the adspecies in the presence of a negative field and thereby mitigating the formation of an oxide.

Bray, Jacob E.↗

DOE OSTI2023

Neural Network-Based Variational Methods for Solving Quadratic Porous Medium Equations in High Dimensions

Here, in this paper, we propose and study neural network-based methods for solutions of high-dimensional quadratic porous medium equation (QPME). Three variational formulations of this nonlinear PDE are presented: a strong formulation and two weak formulations. For the strong formulation, the solution is directly parameterized with a neural network and optimized by minimizing the PDE residual. It can be proved that the convergence of the optimization problem guarantees the convergence of the approximate solution in the $L^1$ sense. The weak formulations are derived following (Brenier in Examples of hidden convexity in nonlinear PDEs, 2020) which characterizes the very weak solutions of QPME. Specifically speaking, the solutions are represented with intermediate functions who are parameterized with neural networks and are trained to optimize the weak formulations. Extensive numerical tests are further carried out to investigate the pros and cons of each formulation in low and high dimensions. This is an initial exploration made along the line of solving high-dimensional nonlinear PDEs with neural network-based methods, which we hope can provide some useful experience for future investigations.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Advancing quantitative models of soil microbiology, ecology, and biochemistry

Rapid advances in soil biology are increasingly reflected in mathematical models. From molecular to global scales, models contribute to fundamental understanding and prediction of critical soil processes, such as carbon sequestration. Recent models explicitly incorporate microbial control over soil carbon and nutrient cycling, an approach that offers greater biological realism but raises new challenges in model parameterization and stability. Although they have proliferated at community and ecosystem scales, microbial-explicit models have yet to be incorporated into most Earth system models. Applying soil models across scales requires integration with data to support parameterization and validation using approaches such as Bayesian data assimilation. In conclusion, to catalyze future progress, we recommend breaking down barriers between modeling and empirical disciplines while broadening access to the computational infrastructure that supports soil biological modeling.

Xie, Wally↗