Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “functional modeling”

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

Tradeoffs and Synergies in Tropical Forest Root Traits and Dynamics for Nutrient and Water Acquisition: Field and Modeling Advances

Vegetation processes are fundamentally limited by nutrient and water availability, the uptake of which is mediated by plant roots in terrestrial ecosystems. While tropical forests play a central role in global water, carbon, and nutrient cycling, we know very little about tradeoffs and synergies in root traits that respond to resource scarcity. Tropical trees face a unique set of resource limitations, with rock-derived nutrients and moisture seasonality governing many ecosystem functions, and nutrient versus water availability often separated spatially and temporally. Root traits that characterize biomass, depth distributions, production and phenology, morphology, physiology, chemistry, and symbiotic relationships can be predictive of plants’ capacities to access and acquire nutrients and water, with links to aboveground processes like transpiration, wood productivity, and leaf phenology. In this review, we identify an emerging trend in the literature that tropical fine root biomass and production in surface soils are greatest in infertile or sufficiently moist soils. We also identify interesting paradoxes in tropical forest root responses to changing resources that merit further exploration. For example, specific root length, which typically increases under resource scarcity to expand the volume of soil explored, instead can increase with greater base cation availability, both across natural tropical forest gradients and in fertilization experiments. Also, nutrient additions, rather than reducing mycorrhizal colonization of fine roots as might be expected, increased colonization rates under scenarios of water scarcity in some forests. Efforts to include fine root traits and functions in vegetation models have grown more sophisticated over time, yet there is a disconnect between the emphasis in models characterizing nutrient and water uptake rates and carbon costs versus the emphasis in field experiments on measuring root biomass, production, and morphology in response to changes in resource availability. Closer integration of field and modeling efforts could connect mechanistic investigation of fine-root dynamics to ecosystem-scale understanding of nutrient and water cycling, allowing us to better predict tropical forest-climate feedbacks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Linear magnetoconductivity as a probe of time-reversal symmetry breaking

Several optical experiments have shown that in magnetic materials, the principal axes of response tensors can rotate as an odd function of an applied magnetic field. Here we offer a microscopic explanation of this effect, and we propose a closely related dc transport phenomenon—an off-diagonal symmetric conductivity, linear and odd in a magnetic field, which we refer to as linear magnetoconductivity (LMC). Although LMC has the same functional dependence on a magnetic field as the Hall effect, its origin is fundamentally different: LMC requires time-reversal symmetry to be broken even before a magnetic field is applied, and is therefore a sensitive probe of magnetism. We demonstrate LMC in three different ways: via a tight-binding toy model, a density functional theory calculation on MnPSe3, and a semiclassical treatment. The third approach identifies two distinct mechanisms yielding LMC: momentum-dependent band magnetization and Berry curvature. Finally, we propose an experimental geometry suitable for detecting LMC, and we demonstrate its applicability using Landauer-Büttiker simulations. Our results emphasize the importance of measuring the full conductivity tensor in magnetic materials, and they introduce LMC as a new transport probe of symmetry.

Sunko, V↗
DOE OSTI2024

PSInet: a new global water potential network

Abstract Given the pressing challenges posed by climate change, it is crucial to develop a deeper understanding of the impacts of escalating drought and heat stress on terrestrial ecosystems and the vital services they offer. Soil and plant water potential play a pivotal role in governing the dynamics of water within ecosystems and exert direct control over plant function and mortality risk during periods of ecological stress. However, existing observations of water potential suffer from significant limitations, including their sporadic and discontinuous nature, inconsistent representation of relevant spatio-temporal scales and numerous methodological challenges. These limitations hinder the comprehensive and synthetic research needed to enhance our conceptual understanding and predictive models of plant function and survival under limited moisture availability. In this article, we present PSInet (PSI—for the Greek letter Ψ used to denote water potential), a novel collaborative network of researchers and data, designed to bridge the current critical information gap in water potential data. The primary objectives of PSInet are as follows. (i) Establishing the first openly accessible global database for time series of plant and soil water potential measurements, while providing important linkages with other relevant observation networks. (ii) Fostering an inclusive and diverse collaborative environment for all scientists studying water potential in various stages of their careers. (iii) Standardizing methodologies, processing and interpretation of water potential data through the engagement of a global community of scientists, facilitated by the dissemination of standardized protocols, best practices and early career training opportunities. (iv) Facilitating the use of the PSInet database for synthesizing knowledge and addressing prominent gaps in our understanding of plants’ physiological responses to various environmental stressors. The PSInet initiative is integral to meeting the fundamental research challenge of discerning which plant species will thrive and which will be vulnerable in a world undergoing rapid warming and increasing aridification.

Forestry↗
DOE OSTI2022

Modification at the C2'-O-Position with 2-Methylbenzothiophene Induces Unique Structural Changes and Thermal Transitions on Duplexes of RNA and DNA

Oligonucleotides can be chemically modified for a variety of applications that include their use as biomaterials, in therapeutics, or as tools to understand biochemical processes, among others. This work focuses on the functionalization of oligonucleotides of RNA and DNA (12- or 14-nucleotides long) with methylbenzothiophene (BT), at the C2'-O-position, which led to unique structural features. Circular dichroism (CD) analyses showed that positioning the BT units on one strand led to significant thermal destabilization, while duplexes where each strand contained 4-BT rings formed a distinct arrangement with cooperativity/interactions among the modifications (evidenced from the appearance of a band with positive ellipticity at 235 nm). Interestingly, the structural arrays displayed increased duplex stabilization (>10 °C higher than the canonical analogue) as a function of [Na + ] with an unexpected structural rearrangement at temperatures above 50 °C. Density functional theory–polarizable continuum model (DFT-PCM) calculations were carried out, and the analyses were in agreement with induced structural changes as a function of salt content. A model was proposed where the hydrophobic surface allows for an internal nucleobase rearrangement into a more thermodynamically stable structure, before undergoing full denaturation, with increased heat. While this behavior is not common, B- to Z-form duplex transitions can occur and are dependent on parameters that were probed in this work, i.e., temperature, nature of modification, or ionic content. To take advantage of this phenomenon, we probed the ability of the modified duplexes to be recognized by Zα (an RNA binding protein that targets Z-form RNA) via electrophoretic analysis and CD. Interestingly, the protein did not bind to canonical duplexes of DNA or RNA; however, it recognized the modified duplexes, in a [monovalent/divalent salt] dependent manner. Overall, the findings describe methodology to attain unique structural motifs of modified duplexes of DNA or RNA, and control their behavior as a function of salt concentration. While their affinity to RNA binding proteins, and the corresponding mechanism of action, requires further exploration, the tunable properties can be of potential use to study this, and other, types of modifications. The novel arrays that formed, under the conditions described herein, provide a useful way to explore the structure and behavior of modified oligonucleotides, in general.

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

Climate Drives Modeled Forest Carbon Cycling Resistance and Resilience in the Upper Great Lakes Region, USA

Forests dominate the global terrestrial carbon budget, but their ability to continue doing so in the face of a changing climate is uncertain. A key uncertainty is how forests will respond (resistance) and recover from (resilience) to rising levels of disturbance of varying intensities. This knowledge gap can optimally be addressed by integrating manipulative field experiments with ecophysiological modeling. We used the ED-2.2 model to project carbon fluxes for a northern temperate deciduous forest subjected to a real-world disturbance severity manipulation experiment. ED-2.2 was run for 150 years, starting from near bare ground in 1900, and subjected to three disturbance treatments under an ensemble of climate conditions. Both disturbance severity and climate strongly affected carbon fluxes, and interacted with one another. We then calculated metrics of functional ecosystem stability. Modeled gross primary production exhibited a mean resistance of -0.13, -0.23, and -0.42 (unitless) in response to increasing disturbance severity of 45%, 65%, and 85% mortality, respectively; and resilience of 0.002, 0.004, 0.008 (yr-1) respectively. This pattern held true for other carbon fluxes, meaning that in general simulated GPP declined more initially but recovered more quickly from higher mortality. Notably, however, heterotrophic respiration responded more slowly to disturbance, and its highly variable response was affected by different drivers. This work provides insight into how future conditions might affect the functional stability of mature forests in this region under ongoing climate change and changing disturbance regimes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

2‐Aminobenzenethiol‐Functionalized Silver‐Decorated Nanoporous Silicon Photoelectrodes for Selective CO 2 Reduction

Abstract A molecularly thin layer of 2‐aminobenzenethiol (2‐ABT) was adsorbed onto nanoporous p‐type silicon (b‐Si) photocathodes decorated with Ag nanoparticles (Ag NPs). The addition of 2‐ABT alters the balance of the CO 2 reduction and hydrogen evolution reactions, resulting in more selective and efficient reduction of CO 2 to CO. The 2‐ABT adsorbate layer was characterized by Fourier transform infrared (FTIR) spectroscopy and modeled by density functional theory calculations. Ex situ X‐ray photoelectron spectroscopy (XPS) of the 2‐ABT modified electrodes suggests that surface Ag atoms are in the +1 oxidation state and coordinated to 2‐ABT via Ag−S bonds. Under visible light illumination, the onset potential for CO 2 reduction was −50 mV vs. RHE, an anodic shift of about 150 mV relative to a sample without 2‐ABT. The adsorption of 2‐ABT lowers the overpotentials for both CO 2 reduction and hydrogen evolution. A comparison of electrodes functionalized with different aromatic thiols and amines suggests that the primary role of the thiol group in 2‐ABT is to anchor the NH 2 group near the Ag surface, where it serves to bind CO 2 and also to assist in proton transfer.

Kan, Miao↗
DOE OSTI2024

Addressing bias in bagging and boosting regression models

As artificial intelligence (AI) becomes widespread, there is increasing attention on investigating bias in machine learning (ML) models. Previous research concentrated on classification problems, with little emphasis on regression models. This paper presents an easy-to-apply and effective methodology for mitigating bias in bagging and boosting regression models, that is also applicable to any model trained through minimizing a differentiable loss function. Our methodology measures bias rigorously and extends the ML model's loss function with a regularization term to penalize high correlations between model errors and protected attributes. We applied our approach to three popular tree-based ensemble models: a random forest model (RF), a gradient-boosted model (GBT), and an extreme gradient boosting model (XGBoost). We implemented our methodology on a case study for predicting road-level traffic volume, where RF, GBT, and XGBoost models were shown to have high accuracy. Despite high accuracy, the ML models were shown to perform poorly on roads in minority-populated areas. Our bias mitigation approach reduced minority-related bias by over 50%.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Capacitor Voltage Ripple Estimation and Optimal Sizing of Modular Multi-Level Converters for Variable-Speed Drives

This article deals with the sizing of submodule capacitors for modular multilevel converters (MMCs) used for variable speed drives for medium voltage and high voltage electric machines. Variable speed drives are required to run over a wide frequency range from zero to a few hundred Hz or even over 1 kHz. Under low frequency operation (<30 Hz), a large current ripple flowing through the MMC submodule capacitors will cause significant capacitor voltage fluctuation. This article proposes an optimal submodule capacitor selection method based on the common-mode voltage injection. The proposed sizing model includes two functions: 1) accurately estimating the capacitor voltage ripple based on the drive system parameters; and 2) identifying the minimum capacitance of the submodule capacitors based on the allowed maximum capacitor voltage fluctuations. Simulation and experimental results show that the error between the actual and estimated values of the voltage fluctuation is within 3% using the proposed model.

42 ENGINEERING↗
DOE OSTI2023

Quantifying the effect of 3D models on moment tensor results using synthetic data

Moment tensors provide vital information on seismic source properties for seismic events. Moment tensors require seismic wavespeed models to compute Green’s functions, which measure the impulse response between a given source and receiver. Traditionally, researchers have used one-dimensional velocity models to calculate Green’s functions since 1D Green’s functions are computationally cheap to compute. Local 1D velocity models can also accurately model waveforms at short distances (< 500 km). However, 1D velocity models do not account for lateral heterogeneity, which can cause significant misfit in tectonically complex regions such as the Middle East (Covellone and Savage, 2012). Green’s functions calculated using 3D seismic wavespeed models have been shown to perform better in tectonically complex regions (Covellone and Savage, 2012; Kintner and Modrak, 2022), so we are interested in quantifying the effect of considering 3D structure on moment tensor inversion results. The Middle East is an ideal study area for a synthetic moment tensor test for two reasons. Firstly, the Middle East is a tectonically complex region that has been heavily studied. Secondly, the Middle East has significant seismic activity throughout the region, but imperfect station coverage due to limited open data through large swaths of the domain. The tectonic complexity and uneven station coverage will test real-world performance even in a synthetic experiment.

58 GEOSCIENCES↗
DOE OSTI2021

Computational Investigation of the Role of Active Site Heterogeneity for aSupported Organovanadium(III) Hydrogenation Catalyst

A crucial consideration for supported heterogeneous catalysts is the nonuniformity of the active sites, particularly for supported organometallic catalysts. Standard spectroscopic techniques, such as X-ray absorption spectroscopy, reflect the nature of the most populated sites, which are often intrinsically structurally distinct from the most active catalytic sites. Additionally, with computational models, often, only a few representative structures are used to depict catalytic active sites on a surface, even though there are numerous observable factors of surface heterogeneity that contribute to the kinetically favorable active species. A previously reported study on the mechanism of a surface organovanadium-(III) catalyst [(SiO 2 )V III (Mes)(THF)] for styrene hydrogenation yielded two possible mechanisms: heterolytic cleavage and redox cycling. These two mechanistic scenarios are challenging to differentiate experimentally since the kinetic readouts of the catalyst are identical. To showcase the importance of modeling surface heterogeneity and its effect on catalytic activity, density functional theory (DFT) computational models of a series of potential active sites of [(SiO 2 )V III (Mes)(THF)] for the reaction pathways are applied in combination with kinetic Monte Carlo (kMC) simulations. Computed results were then compared to the previously reported experimental kinetic study: (1) DFT free-energy reaction pathways indicated the likely active site and pathway for styrene hydrogenation, a heterolytic cleavage pathway requiring a bare tripodal vanadium site. (2) From the kMC simulations, a mixture of different bond lengths from the support oxygen to the metal center was required to qualitatively describe the experimentally observed kinetic aspects of a supported organovanadium(III) catalyst for olefin hydrogenation. This work underscores the importance of modeling surface heterogeneity in computational catalysis.

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

Leveraging prior mean models for faster Bayesian optimization of particle accelerators

Tuning particle accelerators is a challenging and time-consuming task that can be automated and carried out efficiently using suitable optimization algorithms, such as model-based Bayesian optimization techniques. One of the major advantages of Bayesian algorithms is the ability to incorporate prior information about beam physics and historical behavior into the model used to make control decisions. In this work, we examine incorporating prior accelerator physics information into Bayesian optimization algorithms by utilizing fast executing, neural network models trained on simulated or historical datasets as prior mean functions in Gaussian process models. We show that in ideal cases, this technique substantially increases convergence speed to optimal solutions in high-dimensional tuning parameter spaces. Additionally, we demonstrate that even in non-ideal cases, where prior models of beam dynamics do not exactly match experimental conditions, the use of this technique can still enhance convergence speed. Finally, we demonstrate how these methods can be used to improve optimization in practical applications, such as transferring information gained from beam dynamics simulations to online control of the LCLS injector, and transferring knowledge gained from experimental measurements across different operating modes, such as accelerating different ion species at the ATLAS heavy ion accelerator.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Reaction Mechanism of Vapor-Phase Formic Acid Decomposition over Platinum Catalysts: DFT, Reaction Kinetics Experiments, and Microkinetic Modeling

A combination of periodic density functional theory (DFT, PW91-GGA) calculations, reaction kinetics experiments, and mean-field microkinetic modeling is used to derive insights on the reaction mechanism and determine the nature of the active site under reaction conditions for the vapor-phase decomposition of formic acid (FA, HCOOH) over Pt/C catalysts. Microkinetic models formulated using DFT energetics derived on the clean Pt(100) and Pt(111) required large parameter adjustments to reproduce the experimentally measured apparent activation energies and reaction orders. Further, these models predicted high surface coverage of adsorbed carbon monoxide (CO*), inconsistent with the environment of the active site in the DFT calculations on the clean surfaces. Consequently, we reperformed DFT calculations for the entire reaction network on partially CO*-covered (4/9 monolayer, ML) Pt(111) and Pt(100). The resultant microkinetic models, with thermochemistry and kinetics explicitly dependent on CO* coverage, were able to reproduce the experimentally determined activation energies and reaction orders, in addition to being self-consistent in CO* coverage. Our results suggest that Pt(100) is likely poisoned by CO* under typical reaction conditions and does not contribute significantly to the experimentally observed reactivity. Instead, we find that Pt(111) better represents the active site for FA decomposition reaction on Pt/C catalysts. The optimized model on 4/9 ML CO*-covered Pt(111) suggests that the reaction occurs via the carboxyl (COOH*) intermediate and that the spectator CO*-assisted pathways play a significant role under reaction conditions. Lastly, this study underscores the importance of spectator species on the energetics and the mechanism of a catalytic reaction and their key role in developing a model that better addresses the nature of the active site under realistic catalytic reaction conditions.

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

Modeling quasielastic interactions of monoenergetic kaon decay-at-rest neutrinos

Monoenergetic muon neutrinos at 236 MeV are readily produced in intense medium-energy proton facilities (≳2–3GeV) when a positive kaon decays at rest (KDAR; K + →μ + ν μ ). These neutrinos provide a unique opportunity to both study the neutrino interaction and probe the nucleus with a monoenergetic weak-interaction-only tool. We present cross-section calculations for quasielastic scattering of these 236-MeV neutrinos off 12 C and 40 Ar, paying special attention to low-energy aspects of the scattering process. Our model takes the description of the nucleus in a mean-field approach as the starting point, where we solve Hartree-Fock equations using a Skyrme type nucleon-nucleon interaction. Thereby, we introduce long-range nuclear correlations by means of a continuum random-phase approximation (CRPA) framework where we solve the CRPA equations using a Green's function method. The model successfully describes (e,e') data on 12 C and 40 Ca in the kinematic region that overlaps with the KDAR ν μ phase space. In addition to these results, we present future prospects for precision KDAR cross-section measurements and applications of our calculations in current and future experiments that will utilize these neutrinos.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Electronic Structure Theory and Novel Materials

This grant supported research on electronic structure and materials theory, with focus on three main issues: (i) novel techniques to deal with correlation in the electronic ground-state, (ii) topological materials, (iii) the phase diagram of lattice spin models. Regarding (i), we applied to the homogeneous electron liquid an approach that we previously developed in the context of molecular systems. In this scheme the electronic occupation probabilities and the natural spin orbitals are used to construct an approximate two-body density matrix for the electronic ground-state. Regarding (ii) we used standard electronic structure methods based on density functional theory to model topological materials and interpret experimental observations. Finally, regarding (iii) we further developed a numerical approach to compute the renormalized couplings within real space renormalization group theory in the context of lattice spin models. The main findings were the following. (i) We found that with our approximate two-body density matrix, which works well for small molecules, is not sufficiently accurate for condensed phase systems. Missing a systematic way of improving on the adopted approximations, we decided not to pursue this approach. (ii) We performed two studies. In one, we investigated the influence of Te defects on the topological properties of a WTe2 monolayer, finding that while Te vacancies, even in modest concentration, destroy the topological character, Te adatoms do not, consistent with a recent experiment. In another study, we predicted Weyl semimetal character and strong anomalous Hall effect in the Heusler compensated ferrimagnet Ti2MnAl. (iii) We developed a new Monte Carlo method to do real space renormalization group calculations for lattice spin models. We subsequently extended the scheme to deal with lattice spin models in presence of quenched disorder, finding that the approach can distinguish systems with finite and strong disorder. In the finite disorder case, the method allows one to find with good approximation the critical coupling distribution and the critical exponents.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Supervised learning and the finite-temperature string method for computing committor functions and reaction rates

A central object in the computational studies of rare events is the committor function. Though costly to compute, the committor function encodes complete mechanistic information of the processes involving rare events, including reaction rates and transition-state ensembles. Under the framework of transition path theory, Rotskoff et al. [Proceedings of the 2nd Mathematical and Scientific Machine Learning Conference, Proceedings of Machine Learning Research (PLMR, 2022), Vol. 145, pp. 757–780] proposes an algorithm where a feedback loop couples a neural network that models the committor function with importance sampling, mainly umbrella sampling, which collects data needed for adaptive training. Here, in this work, we show additional modifications are needed to improve the accuracy of the algorithm. The first modification adds elements of supervised learning, which allows the neural network to improve its prediction by fitting to sample-mean estimates of committor values obtained from short molecular dynamics trajectories. The second modification replaces the committor-based umbrella sampling with the finite-temperature string (FTS) method, which enables homogeneous sampling in regions where transition pathways are located. We test our modifications on low-dimensional systems with non-convex potential energy where reference solutions can be found via analytical or finite element methods, and show how combining supervised learning and the FTS method yields accurate computation of committor functions and reaction rates. We also provide an error analysis for algorithms that use the FTS method, using which reaction rates can be accurately estimated during training with a small number of samples. The methods are then applied to a molecular system in which no reference solution is known, where accurate computations of committor functions and reaction rates can still be obtained.

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

Investigating late-stage particle production in pp collisions with balance functions

Balance functions have been regarded in the past as a method of investigating the late-stage hadronization found in the presence of a strongly-coupled medium. They are also used to constrain mechanisms of particle production in large and small collision systems. Charge balance functions for inclusive and identified particle pairs are reported as a function of charged particle multiplicity in proton–proton collisions simulated with the PYTHIA8 and the EPOS4 models. The charge balance functions of inclusive, pion, kaon, and proton pairs exhibit amplitudes and shapes that depend on particle species and differ significantly in the two models due to the different particle production mechanisms implemented in PYTHIA and EPOS. The shapes and amplitudes also evolve with multiplicity in both models. In addition, the evolution of the longitudinal rms width and that of balance functions integrals with multiplicity (and average transverse momentum) feature significant differences in the two models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Haplotype‐resolved genome assembly of Populus tremula × P. alba reveals aspen‐specific megabase satellite DNA

SUMMARY Populus species play a foundational role in diverse ecosystems and are important renewable feedstocks for bioenergy and bioproducts. Hybrid aspen Populus tremula × P. alba INRA 717‐1B4 is a widely used transformation model in tree functional genomics and biotechnology research. As an outcrossing interspecific hybrid, its genome is riddled with sequence polymorphisms which present a challenge for sequence‐sensitive analyses. Here we report a telomere‐to‐telomere genome for this hybrid aspen with two chromosome‐scale, haplotype‐resolved assemblies. We performed a comprehensive analysis of the repetitive landscape and identified both tandem repeat array‐based and array‐less centromeres. Unexpectedly, the most abundant satellite repeats in both haplotypes lie outside of the centromeres, consist of a 147 bp monomer PtaM147, frequently span >1 megabases, and form heterochromatic knobs. PtaM147 repeats are detected exclusively in aspens (section Populus ) but PtaM147‐like sequences occur in LTR‐retrotransposons of closely related species, suggesting their origin from the retrotransposons. The genomic resource generated for this transformation model genotype has greatly improved the design and analysis of genome editing experiments that are highly sensitive to sequence polymorphisms. The work should motivate future hypothesis‐driven research to probe into the function of the abundant and aspen‐specific PtaM147 satellite DNA.

Zhou, Ran↗