Search NASA⌕ Search

SEARCH · Search NASA

Results for “reliability estimation”

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

The role of edge plasma parameters in H-mode density limit on the JET-ILW

A study of a dataset of JET H-mode plasma with the Be/W ITER-like wall (JET-ILW) shows that reaching the edge MHD ballooning limit leads to confinement degradation. However, unlike JET plasmas with a carbon wall (JET-C), the JET-ILW plasmas stay in a marginal dithering phase for a relatively long period, associated with a higher (≈20%) H-mode density limit (HDL) than JET-C equivalents. This suggests that ITER could be operated in H-mode with higher density than the scaling based on carbon wall devices, but likely with a dithering phase plasma with lower confinement. A new, reliable estimator for JET E r, min has been derived by combining HRTS measurements of pedestal gradient and edge-SOL decay lengths. JET radial E r ETB wells are observed in the range of -15 to -60 kV m -1 in high performance H-modes, consistent with previous CXRS results in ASDEX Upgrade. The results imply that a higher positive E × B shear in the near SOL plays a role in sustaining a marginal phase in JET-ILW which leads to a higher HDL than that in JET-C. The results of the JET-ILW dataset show agreement with the Goldston finite collisionality HD model for SOL broadening at high collisionality. A hypothesis for the dithering H-mode phase is proposed: as n e,SOL increases, ν *,SOL increases, SOL broadens, E r shear decreases, triggers L-mode; n e drops, ν *,SOL decreases, SOL becomes narrower, and E r shear increases, triggering H-mode, resulting in a cycle of H–L–H- oscillations. For burning plasma devices, such as ITER, operating just below the MHD limit for the dithering phase could be a promising regime for maximising core density, and fusion performance while minimising plasma-material interaction. The oscillatory signal during the dithering phase could be used as a precursor of undesirable plasma performance for control purposes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Deep neural network uncertainty quantification for LArTPC reconstruction

We evaluate uncertainty quantification (UQ) methods for deep learning applied to liquid argon time projection chamber (LArTPC) physics analysis tasks. As deep learning applications enter widespread usage among physics data analysis, neural networks with reliable estimates of prediction uncertainty and robust performance against overconfidence and out-of-distribution (OOD) samples are critical for their full deployment in analyzing experimental data. While numerous UQ methods have been tested on simple datasets, performance evaluations for more complex tasks and datasets are scarce. Here we assess the application of selected deep learning UQ methods on the task of particle classification using the PiLArNet monte carlo 3D LArTPC point cloud dataset. We observe that UQ methods not only allow for better rejection of prediction mistakes and OOD detection, but also generally achieve higher overall accuracy across different task settings. We assess the precision of uncertainty quantification using different evaluation metrics, such as distributional separation of prediction entropy across correctly and incorrectly identified samples, receiver operating characteristic curves (ROCs), and expected calibration error from observed empirical accuracy. We conclude that ensembling methods can obtain well calibrated classification probabilities and generally perform better than other existing methods in deep learning UQ literature.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

QSO photometric redshifts using machine learning and neural networks

ABSTRACT The scientific value of the next generation of large continuum surveys would be greatly increased if the redshifts of the newly detected sources could be rapidly and reliably estimated. Given the observational expense of obtaining spectroscopic redshifts for the large number of new detections expected, there has been substantial recent work on using machine learning techniques to obtain photometric redshifts. Here, we compare the accuracy of the predicted photometric redshifts obtained from deep learning (DL) with the k-nearest neighbour (kNN) and the decision tree regression (DTR) algorithms. We find using a combination of near-infrared, visible, and ultraviolet magnitudes, trained upon a sample of Sloan Digital Sky Survey quasi-stellar objects, that the kNN and DL algorithms produce the best self-validation result with a standard deviation of σΔz = 0.24 (σΔz(norm) = 0.11). Testing on various subsamples, we find that the DL algorithm generally has lower values of σΔz, in addition to exhibiting a better performance in other measures. Our DL method, which uses an easy to implement off-the-shelf algorithm with neither filtering nor removal of outliers, performs similarly to other, more complex, algorithms, resulting in an accuracy of Δz < 0.1 up to z ∼ 2.5. Applying the DL algorithm trained on our 70 000 strong sample to other independent (radio-selected) data sets, we find σΔz ≤ 0.36 (σΔz(norm) ≤ 0.17) over a wide range of radio flux densities. This indicates much potential in using this method to determine photometric redshifts of quasars detected with the Square Kilometre Array.

Curran, S. J.↗
DOE OSTI2025

Learning nuclear cross sections across the chart of nuclides with graph neural networks

We explore the use of deep learning techniques to learn how nuclear cross sections change as we add or remove protons and neutrons. As a proof of principle, we focus on the neutron-induced reactions in the fast energy regime. Our approach follows a two-stage learning framework. First, we apply representation learning to encode cross section data into a latent space using either variational autoencoders (VAEs) or implicit neural representations (INRs). Then, we train graph neural networks (GNNs) on the resulting embeddings to predict missing values across the nuclear chart by leveraging the topological structure of neighboring isotopes. We demonstrate accurate cross section predictions within a 9 × 9 block of missing nuclei. We also find that the optimal GNN training strategy depends on the type of latent representation used, with VAE embeddings performing best under end-to-end optimization in the original space, while INR embeddings achieve better results when the GNN is trained only in the latent space. Furthermore, using clustering algorithms, we map groups of latent vectors into regions of the nuclear chart and show that VAEs and INRs can discover some of the neutron magic numbers. These findings suggest that deep-learning models based on the representation encoding of cross sections combined with graph neural networks hold significant potential in augmenting nuclear theory models, e.g., by providing reliable estimates of covariances of cross sections, including cross-material covariances.

Machine learning↗
DOE OSTI2020

Dominant effect of anharmonicity on the equation of state and thermal conductivity of MgO under extreme conditions

Understanding heat transfer in the earth mantle region is of great scientific interest but highly challenging due to dissimilar conditions from our daily lives. Thermal conductivity of highly pressured and thermally activated materials can be completely different from that at ambient condition. In this study we calculate the equation of state and thermal conductivity of cubic crystalline MgO, a major component of the mantle, in such an environment. It is shown that the material properties are not accurately captured unless rigorous treatment of anharmonicity of phonon vibrational motion is explicitly applied. To predict renormalized phonon dispersion relations, interatomic force constants of high-order expansion are required, especially at elevated temperature. The anharmonicity plays an increasingly substantial role for regulating volume expansion as temperature rises, because the Gibbs free energy of MgO considerably deviates from what was predicted by the quasiharmonic approximation. We find that the lattice thermal conductivity of MgO and the heat transfer of the earth's mantle are reliably estimated only when both the temperature-dependent renormalization of phonon frequencies and four-phonon scattering effects are taken into consideration. We map the characteristics at Geotherm into contour diagrams in which equilibrium volumes and the thermal conductivities are described as functions of temperature and pressure. Furthermore, the thermal conductivity values of MgO in this study are smaller than those in previous theoretical studies due to the additional incorporation of four-phonon scatterings in the calculations. Specifically, deviation between the values in previous studies and in this study becomes larger as thermodynamic variables reach those near the core-mantle boundary. For instance, thermal conductivities of 52.15 and 23.67 Wm -1 K -1 are estimated at ambient and extreme conditions (T = 4000K and P = 135 GPa), respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Comparative study of He 6 β -decay based on different similarity-renormalization-group evolved chiral interactions

We report on a study of the Gamow-Teller matrix element contributing to ${}^6$He $\beta$-decay with similarity renormalization group (SRG) versions of momentum- and configuration-space two-nucleon interactions. These interactions are derived from two different formulations of chiral effective field theory ($\chi$EFT) -- without and with the explicit inclusion of $\Delta$-isobars. We consider evolution parameters $\Lambda_{\rm SRG}$ in the range between 1.2 and 2.0 fm$^{-1}$ and, for the $\Delta$-less case, also the unevolved (bare) interaction. The axial current contains one- and two-body terms, consistently derived at tree level (no loops) in the two distinct $\chi$EFT formulations we have adopted here. The ${}^6$He and ${}^6$Li ground-state wave functions are obtained from hyperspherical-harmonics (HH) solutions of the nuclear many-body problem. In $A\,$=$\,6$ systems, the HH method is limited at present to treat only two-body interactions and non-SRG evolved currents. Our results exhibit a significant dependence on $\Lambda_{\text{SRG}}$ of the contributions associated with two-body currents, suggesting that a consistent SRG-evolution of these is needed in order to obtain reliable estimates. We also show that the contributions from one-pion-exchange currents depend strongly on the model (chiral) interactions and on the momentum- or configuration-space cutoffs used to regularize them. These results might prove helpful in clarifying the origin of the sign difference recently found in No-Core-Shell-Model and Quantum Monte Carlo calculations of the ${}^6$He Gamow-Teller matrix element.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Elastic proton scattering off nonzero spin nuclei

In recent years, we constructed a microscopic optical potential (OP) for elastic nucleon-nucleus ( NA ) scattering using modern approaches based on chiral theories for the nucleon-nucleon ( NN ) interaction. The OP was derived at first order of the spectator expansion in Watson multiple scattering theory and its final expression was a folding integral between the NN t matrix and the nuclear density of the target. Two- and three-body forces are consistently included both in the target and in the projectile description. The purpose of this work is to apply our microscopic OP to nuclei characterized by a ground state of spin-parity quantum numbers J π ≠ 0 + . We extended our formalism to include the spin of the target nucleus. The full amplitudes of the NN reaction matrix are retained in the calculations starting from two- and three-body chiral forces. The microscopic OP can be applied in the energy range 100 ≤ E ≤ 350 MeV. Finally, we show a remarkable agreement with experimental data for the available observables and, simultaneously, provide reliable estimates for the theoretical uncertainties. This work paves the way toward a full microscopic approach to inelastic NA scattering, showing that the derivation of optical potentials between states with J π ≠ 0 + is completely under control.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

16 O ⁡(𝑒,𝑒′⁢𝛼)⁢ 12 C measurements and the 12 C ⁡(𝛼,𝛾)⁢ 16 O astrophysical reaction rate

The 12 C ⁡(𝛼,𝛾)⁢ 16 O reaction, an important component of stellar helium burning, has a key role in nuclear astrophysics. It has significant impact on the evolution and final state of heavy to low mass stars, maximum mass of stellar formed black holes and also shapes the elemental abundances resulting from nucleosynthesis in such stars. Providing a reliable estimate for the energy dependence of this reaction at stellar helium burning temperatures has been a longstanding and important goal. In this work, we study the role of potential new 𝐸⁢1 and 𝐸⁢2 measurements of the 16 O ⁡(𝑒,𝑒′⁢𝛼)⁢ 12 C reaction in reducing the overall uncertainty in the astrophysical 𝑆 factor for the 12 C ⁡(𝛼,𝛾)⁢ 16 O 𝐸⁢1 and 𝐸⁢2 ground state capture extrapolated to a stellar energy of 300 keV. A multilevel 𝑅-matrix analysis is used to make extrapolations of the 𝑆 𝐸⁢1 (300 keV) and 𝑆 𝐸⁢2 (300 keV) factors for the 12 C ⁡(𝛼,𝛾)⁢ 16 O reaction from existing ground state capture data. Bayesian analysis is used to quantify the uncertainties in the extrapolations for both the existing data alone and also when possible new experimental data are included. In particular, we consider a new experiment that would make use of a high-intensity low-energy electron beam that impinges on a windowless oxygen gas target as a means to determine the total 𝐸⁢1 and 𝐸⁢2 ground state cross sections for this reaction. Here we find that the new data could significantly reduce the 𝑆⁡(300) uncertainties. Splitting the new data into high and low energy regions shows that both low and high energy data are effective in reducing the uncertainty.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Model-Free Probabilistic Forecasting of Nodal Voltages in Distribution Systems

As the penetration of distributed energy resources (DERs) into distribution systems increases, so does the interest in forecasting relevant system variables to help mitigate the associated challenges. One such challenge is the more frequent occurrence of excessive voltages in distribution systems with higher shares of DERs. Accurate and reliable estimates together with forecasts of system states (i.e., nodal voltages) will therefore play a key role in improving the utilization of these variable and uncertain sources while mitigating potential operational risks. Whilst recent literature has explored machine learning (ML) methods for voltage estimation and their extrapolation for a short-time period into the future, few have taken uncertainty quantification into account, and these methods have not yet been translated into operations. This paper discusses the advantages offered by probabilistic voltage forecasts and proposes a non-parametric Bayesian method suitable for forecasting nodal voltages at short-term time horizons while accounting for uncertainties in load and distributed photovoltaic (PV) generation. We demonstrate the value of the proposed Gaussian process (GP) model for a case study using historical forecasts and observation data.

distribution system↗
DOE OSTI2021

Tissue scale agent-based simulation of premalignant progressions in Barrett’s esophagus

Barrett’s esophagus (BE) is a benign condition of the distal esophagus that initiates a multistage pathway to esophageal adenocarcinoma (EAC). Short of frequent intrusive (and costly) surveillance, effective screening for neoplasia in BE populations is yet to be established since progressors are rare and virtually undetectable without routine biopsies, which often sample only a small portion of the BE tissue. As a result, reliable estimation of the true prevalence of dysplasia in a BE population and evidence-based optimization of screening for at-risk individuals is challenging. Data-driven microsimulations, i.e., model-generated instances of disease history in a predefined virtual population, have found utility in the EAC screening literature as low-overhead alternatives to real-world hypothesis testing of optimal interventions for dysplasia. Despite the successes, computational limitations, paucity of knowledge and data on Barrett’s dysplasia, and the complexities of disease progression as a multiscale multiphysics process have hindered the treatment of disease progression in BE as a spatial process. Agent-based modeling of nucleation and proliferation processes in dysplasia warrants exploration in this context as an approximation that operates at a trade-off between computational tractability and precise representation of the composition and physics of the substrate (tissue). In this study, we describe spatially resolved simulations of premalignant progression toward EAC in a coarse-grained model of Barrett’s tissue that resolves the metaplastic tissue at a length scale of 0.42 mm (~3300 crypts/mm 2 ). Finally, the model is calibrated to reproduce historical high-grade dysplasia prevalence when model-generated patients are screened using the Seattle protocol.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

In-vivo Raman microspectroscopy reveals differential nitrate concentration in different developmental zones in Arabidopsis roots

Abstract Background Nitrate (NO 3 − ) is one of the two major forms of inorganic nitrogen absorbed by plant roots, and the tissue nitrate concentration in roots is considered important for optimizing developmental programs. Technologies to quantify the expression levels of nitrate transporters and assimilating enzymes at the cellular level have improved drastically in the past decade. However, a technological gap remains for detecting nitrate at a high spatial resolution. Using extraction-based methods, it is challenging to reliably estimate nitrate concentration from a small volume of cells (i.e., with high spatial resolution), since targeting a small or specific group of cells is physically difficult. Alternatively, nitrate detection with microelectrodes offers subcellular resolution with high cell specificity, but this method has some limitations on cell accessibility and detection speed. Finally, optical nitrate biosensors have very good ( in-vivo ) sensitivity (below 1 mM) and cellular-level spatial resolution, but require plant transformation, limiting their applicability. In this work, we apply Raman microspectroscopy for high-dynamic range in-vivo mapping of nitrate in different developmental zones of Arabidopsis thaliana roots in-situ . Results As a proof of concept, we have used Raman microspectroscopy for in-vivo mapping of nitrate content in roots of Arabidopsis seedlings grown on agar media with different nitrate concentrations. Our results revealed that the root nitrate concentration increases gradually from the meristematic zone (~ 250 µm from the root cap) to the maturation zone (~ 3 mm from the root cap) in roots grown under typical growth conditions used for Arabidopsis, a trend that has not been previously reported. This trend was observed for plants grown in agar media with different nitrate concentrations (0.5–10 mM). These results were validated through destructive measurement of nitrate concentration. Conclusions We present a methodology based on Raman microspectroscopy for in-vivo label-free mapping of nitrate within small root tissue volumes in Arabidopsis. Measurements are done in-situ without additional sample preparation. Our measurements revealed nitrate concentration changes from lower to higher concentration from tip to mature root tissue. Accumulation of nitrate in the maturation zone tissue shows a saturation behavior. The presented Raman-based approach allows for in-situ non-destructive measurements of Raman-active compounds.

Fernández González, Alma↗
DOE OSTI2024

A three-dimensional model of terrain-induced updrafts for movement ecology studies

Abstract Background Spatially explicit simulation models of animal movements through the atmosphere necessarily require a representation of the spatial and temporal variation of atmospheric conditions. In particular, for movements of soaring birds that rely extensively on vertical updrafts to avoid flapping flight, accurate and reliable estimation of the vertical component of wind is critical. The interaction between wind and complex terrain shapes both the horizontal and vertical wind fields, highlighting the need to model the coupling between local terrain features and atmospheric conditions at scales relevant to animal movement. Methods In this work, we propose a new empirical model for estimating the orographic updraft field. The model is developed using computational fluid dynamics simulations of canonical atmospheric conditions over moderately complex terrain. To isolate buoyancy and thermal effects, and focus on terrain-induced effects, we use only simulations of a neutrally stratified atmosphere to develop the model. The model, which we name Engineering Vertical Velocity Estimator (EVVE), is simple to implement and is a function of the underlying terrain elevation map, the desired height above ground level (AGL), and wind conditions at a reference height (80 m). We validate the model with data from the Alaiz mountain (Spain) field campaign. Results Compared to observations, the proposed improved model estimates the updrafts at 120 m AGL with a mean error of 0.11 m/s ( $$\sigma =0.28$$ σ = 0.28 m/s), compared to 0.85 m/s ( $$\sigma =0.58$$ σ = 0.58 m/s) for its baseline. For typical land-based wind turbine hub heights of 80 m AGL, the proposed model has a mean error of 0.04 m/s ( $$\sigma =0.25$$ σ = 0.25 m/s), compared to baseline 0.54 m/s ( $$\sigma =0.45$$ σ = 0.45 m/s) estimations. We illustrate an application of the model in movement ecology by comparing simulated tracks and presence maps of golden eagles ( Aquila chrysaetos ) moving across two distinct landscapes. The tracks and presence maps are obtained using a simple heuristic-based movement model, with the updraft field given by the proposed model and a wind vector-based estimation approach that is currently in wide use in movement ecology studies of raptors and other soaring birds. Conclusions We highlight that movement model results can be sensitive to the underlying orographic updraft model, especially in studies of fine-scale movements in regions of complex topography. We suggest adopting the proposed model rather than the wind vector estimation method for studies of soaring bird movements.

17 WIND ENERGY↗
DOE OSTI2021

Forecasting Marine Sediment Properties with Geospatial Machine Learning

Using a combination of geospatial machine learning prediction and sediment thermodynamic/physical modeling, we have developed a novel software workflow to create probabilistic maps of geoacoustic and geomechanical sediment properties of the global seabed. This new technique for producing reliable estimates of seafloor properties can better support Naval operations relying on sonar performance and seabed strength, can constrain models of shallow tomographic structure important for nuclear treaty compliance monitoring/detection, and can provide constraints on the distribution and inventory of shallow methane gas and gas hydrate accumulations on the continental shelves.

58 GEOSCIENCES↗
DOE OSTI2022

Convective Parameters Derived from Radiosonde Data (SONDEPARAM) Value-Added Product Report

Radiosondes provide fundamental observations of the vertical profile of atmospheric state (pressure, temperature, humidity, and winds), with important implications for subsequent studies on environmental controls on cloud conditions. Within convective cloud environments, there is an increasing demand for additional value-added products (VAPs) to facilitate the use of U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility radiosonde data sets. Such VAPs should provide quick and reliable estimates for several standard radiosonde parameters or quantities of interest using common assumptions, as well as open, flexible code for visualization and user interaction. The Convective Parameters Derived from Radiosonde Data (SONDEPARAM) VAP will apply several robust algorithms used in Wang et al. (2020) for the calculation of useful radiosonde convective cloud parameters, including the convective available potential energy (CAPE), convective inhibition (CIN), and other convective parameters, for several different assumptions regarding the initial parcel characteristics (i.e., surface-based, most unstable, mixed layer). These ARM VAP codes are developed in open, flexible Python formats, with the intention that these parameters/calculations will be incorporated into traditional ARM quick-look radiosonde plotting, yet associated with user-available codes for ease in user reproduction and assumption modification.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Marine carbon dioxide removal (mCDR) test site needs: A joint PNNL-PMEL workshop report

This workshop report summarizes testing and capability needs for marine Carbon Dioxide Removal (mCDR) research identified by researchers from DOE’s Pacific Northwest National Laboratory (PNNL) and NOAA’s Pacific Marine Environmental Laboratory (PMEL) as essential to advancing the state of the science and technology. A one-day workshop was held for 36 hybrid participants on June 2, 2023 at PNNL’s Marine and Coastal Sciences Laboratory in Sequim, WA. The workshop began with DOE and NOAA leadership offering perspectives on the mission alignment and long-term implementation of mCDR research within NOAA’s Office of Ocean and Atmospheric Research (OAR), DOE’s Water Power Technologies Office (WPTO), and DOE’s Office of Fossil Energy and Carbon Management (FECM). As identified by the participants and by program leadership at both agencies, many mCDR strategies are in early stages of development making it difficult to reliably estimate their impact. The essential next step in proving the efficacy, safety, and sustainability of these solutions and identifying market potential requires studies under field-relevant conditions.

13 HYDRO ENERGY↗
DOE OSTI2024

Williston Basin CORE-CM Initiative Final Report

The University of North Dakota Energy & Environmental Research Center (EERC) is leading the Williston Basin Carbon Ore, Rare Earth, and Critical Minerals (CORE-CM) Initiative to drive the expansion and transformation of coal and coal-based resource usage within the Williston Basin to produce rare-earth elements (REEs), CMs, and nonfuel carbon-based products (CBPs). This project is the first phase in a long-term program and set the stage for future work by assessing resource, market, technology, and infrastructure knowledge; identifying knowledge gaps; developing a series of plans to be carried out in future work; and initiating stakeholder engagement. Composed of several tasks, the project sought to identify, characterize, and assess several necessary aspects vital to make this future work a reality. The project’s fundamental task was to characterize the Williston Basin CORE-CM resources. Over 2500 samples from multiple sources were utilized to begin the assessment. Several locations were identified in western North Dakota where sample analysis identified the total REE (TREE) concentration as being over 500 parts per million (ppm), which is at a concentration level that would be suitable to consider for mining and extraction. Current operating coal mines have sufficient concentrations of TREEs for consideration. However, the current data across the basin are still not adequate to fully characterize REE and CM content nor give reliable estimates of the total resource potential. Waste stream reuse was also considered, and several streams were identified which ranged from potential energy sources to chemicals to material wastes. This includes streams that result from oil and gas production. These streams are not fully characterized, and further data are needed before they can be accurately assessed. Infrastructure within the Williston Basin is suitable for expansion of a new industry to mine, extract, and concentrate REEs and CMs. The development of this industry will not only preserve many existing jobs in the coal-mining industry but produce many new jobs. The supply chain for REEs and CMs is currently controlled outside of the United States in nations such as China, but the potential to develop the supply chain within the basin is considered possible. Processing of the mined materials for REEs and CMs needs further research. The technology and knowhow exist outside of the United States, and within the country much of the knowledge has been lost and must be regained. To develop the supply chain and regain lost processing technology, the creation of technology innovation centers (TICs) is crucial. The Williston Basin contains several similar centers and entrepreneurial assistance for other industries that can be applied in the development of REE and CM innovation centers. Education to develop the new skill sets required is also needed. Outreach is important for the development of the REE and CM industry within the basin. Understanding throughout federal and state governments, state agencies, industry, and resource end users is vital for the industry to form and grow. Through this project these groups have been contacted through bulletins, presentations, webinars, and annual symposiums. The report is a summary of the work conducted and throughout refers to a series of appendixes which contain more thorough and specific information about each section.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2026

Using PyBioNetFit to leverage qualitative and quantitative data in biological model parameterization and uncertainty quantification

Data generated in studies of cellular regulatory systems are often qualitative. For example, measurements of signaling readouts in the presence and absence of mutations may reveal a rank ordering of responses across conditions but not the precise extents of mutation-induced differences. Qualitative data are often ignored by mathematical modelers or are considered in an ad hoc manner, as in the study of Kocieniewski and Lipniacki (2013) [Phys Biol 10: 035006], which was focused on the roles of MEK isoforms in ERK activation. In this earlier study, model parameter values were tuned manually to obtain consistency with a combination of qualitative and quantitative data. This approach is not reproducible, nor does it provide insights into parametric or prediction uncertainties. Here, starting from the same data and the same ordinary differential equation (ODE) model structure, we generate formalized statements of qualitative observations, making these observations more reusable, and we improve the model parameterization procedure by applying a systematic and automated approach enabled by the software package PyBioNetFit. We also demonstrate uncertainty quantification (UQ), which was absent in the original study. Our results show that PyBioNetFit enables qualitative data to be leveraged, together with quantitative data, in parameterization of systems biology models and facilitates UQ. These capabilities are important for reliable estimation of model parameters and model analyses in studies of cellular regulatory systems and reproducibility.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Between Harmonic Crystal and Glass: Solids with Dimpled Potential-Energy Surfaces Having Multiple Local Energy Minima

Solids with dimpled potential-energy surfaces are ubiquitous in nature and, typically, exhibit structural (elastic or phonon) instabilities. Dimpled potentials are not harmonic; thus, the conventional quasiharmonic approximation at finite temperatures fails to describe anharmonic vibrations in such solids. At sufficiently high temperatures, their crystal structure is stabilized by entropy; in this phase, a diffraction pattern of a periodic crystal is combined with vibrational properties of a phonon glass. As temperature is lowered, the solid undergoes a symmetry-breaking transition and transforms into a lower-symmetry phase with lower lattice entropy. Here, we identify specific features in the potential-energy surface that lead to such polymorphic behavior; we establish reliable estimates for the relative energies and temperatures associated with the anharmonic vibrations and the solid–solid symmetry-breaking phase transitions. We show that computational phonon methods can be applied to address anharmonic vibrations in a polymorphic solid at fixed temperature. To illustrate the ubiquity of this class of materials, we present a range of examples (elemental metals, a shape-memory alloy, and a layered charge-density-wave system); we show that our theoretical predictions compare well with known experimental data.

36 MATERIALS SCIENCE↗