Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Random fields”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

The direct and indirect drivers shaping RNA viral communities in grassland soil

Recent studies have revealed diverse RNA viral communities in soils. Yet, how environmental factors influence soil RNA viruses remains largely unknown. Here, we recovered RNA viral communities from 24 metatranscriptomes sequenced from grassland soils managed under a range of environmental conditions including 1) water content: 100% and 25% water holding capacity, 2) plant presence: planted with tall wheatgrass (Thinopyrum ponticum) and bare soil, 3) cultivar type: Alkar and Jose, and 4) soil depth: 0-5 cm and 15-25 cm. The recovered RNA viral communities were novel with nearly one-third of the RNA viral contigs uniquely detected in the studied grassland. The classified RNA viral contigs are mostly known as eukaryotic RNA viruses (74.7%) belonging to Phyla Duplornaviricota, Kitrinoviricota, Lenarviricota, and Pisuviricota. Eukaryotic RNA viruses of Family Mitoviridae as well as their natural hosts, Fungi, are one of the most dominant taxa. Consistent with the results of nonmetric multidimensional scaling analysis, the four environmental conditions (water content, plant presence, cultivar, and soil depth) significantly influence the assemblages of soil RNA viral communities as suggested by the correlation analysis and the random forest algorithm. The modularity analysis of the factor network and the structural equation modeling further support the hierarchical associations among the four environmental factors and the community factors representing the co-existing eukaryotic, prokaryotic, and RNA viral communities. The soil water content, plant presence, and type of cultivar demonstrate a significant positive impact on eukaryotic RNA viral richness directly as well as indirectly on eukaryotic RNA viral abundance via influencing the co-existing eukaryotic members in this soil. Our data also provide statistical support for the negative influence of soil depth on soil eukaryotic richness and abundances resulting in its indirect impact on soil eukaryotic RNA viral communities. This study provides field-relevant information on how environmental and community factors collectively shape soil RNA communities and contribute to ecological understanding of RNA viral survival under various environmental conditions and virus-host interactions in soil.

Wu, Ruonan↗
DOE OSTI2023

Operator Relaxation and the Optimal Depth of Classical Shadows

Classical shadows are a powerful method for learning many properties of quantum states in a sample-efficient manner, by making use of randomized measurements. Here we study the sample complexity of learning the expectation value of Pauli operators via “shallow shadows,” a recently proposed version of classical shadows in which the randomization step is effected by a local unitary circuit of variable depth t. Here we show that the shadow norm (the quantity controlling the sample complexity) is expressed in terms of properties of the Heisenberg time evolution of operators under the randomizing (“twirling”) circuit—namely the evolution of the weight distribution characterizing the number of sites on which an operator acts nontrivially. For spatially contiguous Pauli operators of weight k, this entails a competition between two processes: operator spreading (whereby the support of an operator grows over time, increasing its weight) and operator relaxation (whereby the bulk of the operator develops an equilibrium density of identity operators, decreasing its weight). From this simple picture we derive (i) an upper bound on the shadow norm which, for depth t~log⁡(k), guarantees an exponential gain in sample complexity over the t=0 protocol in any spatial dimension, and (ii) quantitative results in one dimension within a mean-field approximation, including a universal subleading correction to the optimal depth, found to be in excellent agreement with infinite matrix product state numerical simulations. Our Letter connects fundamental ideas in quantum many-body dynamics to applications in quantum information science, and paves the way to highly optimized protocols for learning different properties of quantum states.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Quantifying nitrogen loss hotspots and mitigation potential for individual fields in the US Corn Belt with a metamodeling approach

The high productivity in the US Corn Belt is largely enabled by the consumption of millions of tons of manufactured fertilizer. Excessive application of nitrogen (N) fertilizer has been pervasive in this region, and the unrecovered N eventually escaped from croplands in forms of nitrous oxide (N 2 O) emission and N leaching. Mitigating these negative impacts is hindered by a lack of practical information on where to focus and how much mitigation potential to expect. At a large scale, process-based crop models are the primary tools for predicting variables required by decision making, but their applications are prohibited by expensive computational and data storage costs. To overcome these challenges, we built a series of metamodels to learn the key mechanisms regarding the carbon (C) and N cycle from a well-validated process-based biogeochemical model, ecosys. The trained metamodel captures over 98% of the variability of the ecosys simulated outputs for 99 randomly selected counties in Iowa, Illinois, and Indiana. To identify hotspots with high mitigation potential, we introduce net societal benefit (NSB) as an indicator for synthesizing the loss in yield and social benefits through emissions and pollutants avoided. Our results show that reducing N fertilizer by 10% leads to 9.8% less N 2 O emissions and 9.6% less N leaching at the cost of 4.9% more SOC depletion and 0.6% yield reduction over the study region. The estimated total annual NSB is $\$395$ M (uncertainty ranges from $\$114$ M to $\$1271$ M), including $\$334$ from social benefits (uncertainty ranges from $\$46$ M to $\$1076$ M), $\$100$ M from saving fertilizer (uncertainty ranges from $\$13$ M to $\$455$ M), and –$\$40$ M due to yield changes (uncertainty ranges from –$\$261$ M to $\$69$ M). For the median scenario, we noted that 20% of the study area accounts for nearly 50% of the NSB, and thus represent hotspot locations for targeted mitigation. Although the uncertainty range suggests that developing such a high-resolution framework is not yet settled and the scenario based estimations are not appropriate to inform the management practices for individual farmers, our efforts shed light on the new generation of analytical tools for life cycle assessment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Daejeon16 interaction with contact-term corrections for heavy nuclear systems

The Daejeon16 two-nucleon interaction is employed in many-body approaches based on the mean-field approximation. The perturbative character of Daejeon16 is verified by comparing results for 16 O from the Hartree–Fock (HF) approximation and from the no-core shell model and by examining the magnitude of perturbative corrections to the HF energy in light and heavy nuclei. Here, in order to approximately describe energies and radii across the nuclear chart, a phenomenological correction in the form of a two-plus-three-nucleon contact interaction is introduced. With fitted parameters we achieve a very good description of medium-mass nuclei in terms of energy and size and also in terms of the centroid energy of the giant monopole resonance and the dipole polarizability calculated within the random-phase approximation. Our results provide further justification for the use of Daejeon16 augmented with phenomenological corrections as an effective interaction of perturbative character in a variety of applications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Phonon-induced collective modes in spin-orbit coupled polar metals

We study the interplay between collective electronic and lattice modes in polar metals in an applied magnetic field aligned with the polar axis. Static spin-orbit coupling leads to the appearance of a particle-hole spin-flip continuum that is gapped at low energies in a finite field. We find that a weak spin-orbit assisted coupling between electrons and polar phonons induces the emergence of electronic collective modes. The strength of the applied magnetic field tunes the number of modes and their energies, which can lie both above and below the particle-hole continuum. For a range of field values, we identify Fano-type interference between the electronic continuum and phonons. Here, we show that signatures of these collective modes can be observed in optical spectroscopy experiments, and we provide the corresponding theoretical predictions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

ESR Dipole Power Supply Current Ripple and Noise Specifications

This note presents key findings for the ESR main magnet dipole power supplies (PS), where we find the current ripple specification to be close to or beyond the state-of-the-art. These specifications originate from beam-beam considerations, with the requirement to limit the ripple-induced hadron emittance growth to below 10%/hour. Beam dynamics that drive this PS ripple specification arise from the beam motions at the Interaction Point (IP). The frequency of the motions can be separated into "low", compared to the betatron frequency, and "high", i.e. around the betatron frequency and harmonics. In terms of the driving frequency, "low" implies ƒ<<ƒ 0 v x,y and "fast" means ƒ≈{ƒ 0 v x,y , ƒ 0 (1- v x,y , etc.}, where ƒ 0 =1/ T 0 =78.2 kHz is the revolution frequency, and v x,y are the fractional parts of the betatron tunes. Frequencies higher than ƒ 0 /2 will be folded back due to the particles sampling the field once per turn. To provide flexibility for future lattice adjustments and working point variations, we do not consider the tunes as fixed. Instead, we assume a certain margin and allow them to potentially fall within the range of 0.1< v x,y <0.5. In other words, the high-frequency region spans approximately from 8kHz to 40 kHz. Consequently,, we definte the dipole PS ripple in two distinct frequency ranges: the low-frequency range of [1-8000] Hz and the high-frequency range of [8-40] kHz. For the physics effects we analyzed in this note, there is no distinction between the ripple (which can be approximately reproduced in the frequency domain) and the random noise if both have some power within the frequency bandwidth of interest. Therefore, while we will use the term "ripple" for short, it should always be understood that we are referring to "ripple plus noise." Except for the lower end of the low-frequency range, the impact of the rippling PS current on the beam will be considerably reduced due to the eddy currents induced in the walls of the vacuum chamber. We will account for this effect in the PS ripple specifications to follow. The remaining sections of this note are structure as follows: Section II outlines the beam-beam physics requirements for the positional stability of the beam at the IP. Section III describes the anticipated shielding effect of the eddy currents induced in the vacuum chamber. Section IV derives the ripple requirement for the low-frequency range by propagating the closed orbit ripple resulting from the rippling dipoles to the IP (relevant lattice simulation results are summarized in the Appendix). In Section V, we present the analytical criterion for the ripple in the high-frequency range by considering resonant oscillations of the electron beam around a stable closed orbit near the betatron frequency. Finally, Section VI provides a summary of our findings and discusses related work.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

SDSS-IV MaNGA: calibration of astrophysical line-widths in the Hα region using HexPak observations

ABSTRACT We have re-observed $\rm \sim$40 low-inclination, star-forming galaxies from the MaNGA survey (σ ∼ 65 km s−1) at ∼6.5 times higher spectral resolution (σ ∼ 10 km s−1) using the HexPak integral field unit on the WIYN 3.5-m telescope. The aim of these observations is to calibrate MaNGA’s instrumental resolution and to characterize turbulence in the warm interstellar medium and ionized galactic outflows. Here we report the results for the Hα region observations as they pertain to the calibration of MaNGA’s spectral resolution. Remarkably, we find that the previously reported MaNGA line-spread-function (LSF) Gaussian width is systematically underestimated by only 1 per cent. The LSF increase modestly reduces the characteristic dispersion of H ii regions-dominated spectra sampled at 1–2 kpc spatial scales from 23 to 20 km s−1 in our sample, or a 25 per cent decrease in the random-motion kinetic energy. This commensurately lowers the dispersion zeropoint in the relation between line-width and star-formation rate surface-density in galaxies sampled on the same spatial scale. This modest zero-point shift does not appear to alter the power-law slope in the relation between line-width and star-formation rate surface-density. We also show that adopting a scheme whereby corrected line-widths are computed as the square root of the median of the difference in the squared measured line width and the squared LSF Gaussian avoids biases and allows for lower signal-to-noise data to be used reliably.

Chattopadhyay, Sabyasachi↗
DOE OSTI2021

A Spatiotemporal Prediction Model for Regional Scheduling of Shared Bicycles Based on the INLA Method

Dock-less bicycle-sharing programs have been widely accepted as an efficient mode to benefit health and reduce congestions. And modeling and prediction has always been a core proposition in the field of transportation. Most of the existing demand prediction models for shared bikes take regions as research objects; therefore, a POI-based method can be a beneficial complement to existing research, including zone-level, OD-level, and station-level techniques. Point of interest (POI) is the location description of spatial entities, which can reflect the cycling route characteristics for both commuting and noncommuting trips to a certain extent, and is also the main generating point and attraction point of shared-bike travel flow. In this study, we make an effort to model a POI-level cycling demand with a Bayesian hierarchical method. The proposed model combines the integrated nested Laplace approximation (INLA) and random partial differential equation (SPDE) to cope with the huge computation in the modeling process. In particular, we have adopted the dock-less bicycle-sharing rental records of Mobike as a case study to validate our method; the study area was one of the fastest growing urban districts in Shanghai in August 2016. The operation results show that the method can help better understand, measure, and characterize spatiotemporal patterns of bike-share ridership at the POI level and quantify the impact of the spatiotemporal effect on bicycle-sharing use.

Yu, Zhuoran↗
DOE OSTI2024

NoRA: A Tensor Network Ansatz for Volume-Law Entangled Equilibrium States of Highly Connected Hamiltonians

Motivated by the ground state structure of quantum models with all-to-all interactions such as mean-field quantum spin glass models and the Sachdev-Ye-Kitaev (SYK) model, we propose a tensor network architecture which can accomodate volume law entanglement and a large ground state degeneracy. We call this architecture the non-local renormalization ansatz (NoRA) because it can be viewed as a generalization of MERA, DMERA, and branching MERA networks with the constraints of spatial locality removed. We argue that the architecture is potentially expressive enough to capture the entanglement and complexity of the ground space of the SYK model, thus making it a suitable variational ansatz, but we leave a detailed study of SYK to future work. We further explore the architecture in the special case in which the tensors are random Clifford gates. Here the architecture can be viewed as the encoding map of a random stabilizer code. We introduce a family of codes inspired by the SYK model which can be chosen to have constant rate and linear distance at the cost of some high weight stabilizers. We also comment on potential similarities between this code family and the approximate code formed from the SYK ground space.

Physics↗
DOE OSTI2023

Observations and Machine-Learned Models of Near-Surface Permafrost along the Koyukuk River, Alaska, USA

This dataset contains GeoTIFs (raster) and GeoPackages (vector) that map observations of near-surface permafrost and not-permafrost from a field campaign conducted near the village of Huslia, AK along the Koyukuk River and its floodplain in July 2018. These data were collected as part of a campaign to understand if and how permafrost impacts riverbank erosion. This problem cannot be assessed without knowing where permafrost exists. Permafrost was observed via frost probing (to a maximum depth of one meter), coring (to a maximum depth of two meters) and bank/bar excavations. An additional boat survey was performed wherein expert (Joel Rowland) judgment assessed the presence or absence of distinctive permafrost features (e.g., overhanging tundra mats, thermoerosional niching, ice wedges, active drainage of ice melt from soils). This dataset also contains the input features and results of two machine learning models (random forest and convolutional neural network) that extrapolate the observations to the full floodplain that may be useful for building, testing, or validating other machine-learned permafrost models. Permafrost data are provided as georasters of the same shape and geovectors (polylines/polygons) and are all projected into EPSG:32605. All data can be visualized with a GIS (QGIS, ArcGIS, etc.).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Complexity: from Turbulence to Accelerators to Poetry

Though it may seem intuitive at times, the term 'complexity' is multifaceted, often encompassing various interacting aspects and components that follow local rules. This interaction leads to non-linearity, randomness, collective dynamics, hierarchy, and emergence. The phenomenological complexity found in large accelerators has recently been addressed by the Snowmass’21 Future Colliders Implementation Task Force (ITF). In this presentation, we will provide a concise overview of the mathematical constructs employed to assess hierarchical complexity, summarize the key findings from the ITF analysis, and offer a few examples demonstrating how this approach can be applied to analyze certain literaryclassics.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

High-fidelity retrieval from instantaneous line-of-sight returns of nacelle-mounted lidar including supervised machine learning

Abstract. Wind turbine applications that leverage nacelle-mounted Doppler lidar are hampered by several sources of uncertainty in the lidar measurement, affecting both bias and random errors. Two problems encountered especially for nacelle-mounted lidar are solid interference due to intersection of the line of sight with solid objects behind, within, or in front of the measurement volume and spectral noise due primarily to limited photon capture. These two uncertainties, especially that due to solid interference, can be reduced with high-fidelity retrieval techniques (i.e., including both quality assurance/quality control and subsequent parameter estimation). Our work compares three such techniques, including conventional thresholding, advanced filtering, and a novel application of supervised machine learning with ensemble neural networks, based on their ability to reduce uncertainty introduced by the two observed nonideal spectral features while keeping data availability high. The approach leverages data from a field experiment involving a continuous-wave (CW) SpinnerLidar from the Technical University of Denmark (DTU) that provided scans of a wide range of flows both unwaked and waked by a field turbine. Independent measurements from an adjacent meteorological tower within the sampling volume permit experimental validation of the instantaneous velocity uncertainty remaining after retrieval that stems from solid interference and strong spectral noise, which is a validation that has not been performed previously. All three methods perform similarly for non-interfered returns, but the advanced filtering and machine learning techniques perform better when solid interference is present, which allows them to produce overall standard deviations of error between 0.2 and 0.3 m s−1, or a 1 %–22 % improvement versus the conventional thresholding technique, over the rotor height for the unwaked cases. Between the two improved techniques, the advanced filtering produces 3.5 % higher overall data availability, while the machine learning offers a faster runtime (i.e., ∼ 1 s to evaluate) that is therefore more commensurate with the requirements of real-time turbine control. The retrieval techniques are described in terms of application to CW lidar, though they are also relevant to pulsed lidar. Previous work by the authors (Brown and Herges, 2020) explored a novel attempt to quantify uncertainty in the output of a high-fidelity lidar retrieval technique using simulated lidar returns; this article provides true uncertainty quantification versus independent measurement and does so for three techniques rather than one.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Understanding vortex dynamics in CaK(Fe,Ni) 4 As 4 and Ba(Fe,Co) 2 As 2 single crystals under the influence of random point disorder

We report on the influence of doping on vortex dynamics in 3 MeV proton-irradiated single crystals of CaK(Fe 1–x Ni x ) 4 As 4 (1144, x = 0.015, 0.025, and 0.03) and Ba(Fe 1–x Co x ) 2 As 2 (x = 0.04, 0.062, 0.066 and 0.074). Non-irradiated crystals of the 1144 system display superconducting critical temperatures ranging from 31 K for x = 0.015–20.5 K, as doping increases to 0.03. On the other hand, pristine crystals of the 122 system show T c values between 14.6 and 23.6 K, with the maximum T c occurring at intermediate doping levels. The fluence was set at 3 × 10 16 p cm –2 , resulting in a decrease in the T c by around 1.5 K for all samples and significantly affecting the vortex dynamics by reducing the flux creep relaxation compared to previously reported values for unirradiated crystals. Parameters such as vortex pinning energy U 0 and the glassy exponent μ dependencies on doping and magnetic field strength are identified. For the 1144 system, U 0 reaches values approaching 500 K for small fields in samples with T c = 29.3 K (x = 0.015), systematically decreasing to around 200 K as T c falls below 20 K. Furthermore, U 0 decreases as the field increases to 3 T for the same sample, varying from approximately 250 K to 100 K as T c decreases. These changes are typically accompanied by modifications in μ, gradually increasing from values around 1 towards 1.5, corresponding to small bundle relaxation in the collective creep theory. Despite differences in the substitutional disorder and magnetic phase diagram with respect to the 1144 system, the results for 122 single crystals follow a similar tendency in which U 0 usually reduces and μ increase rise as the applied magnetic field is increased. Due to moderate U 0 in these systems (few hundreds of kelvins), the resulting decay of persistent current at liquid helium temperatures is primarily determined by a balance between U 0 and bundle size contribution. Furthermore, these findings provide valuable insights for potential applications of these systems, particularly in the context of intrinsic superconducting parameters and the resulting pinning landscape.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Quantum optimal control of ten-level nuclear spin qudits in Sr 87

We study the ability to implement unitary maps on states of the I = 9 / 2 nuclear spin in Sr 87 , a d = 10 dimensional (qudecimal) Hilbert space, using quantum optimal control. Through a combination of nuclear spin resonance and a tensor ac Stark shift, by solely modulating the phase of a radio-frequency magnetic field, the system is quantum controllable. Alkaline-earth-metal atoms, such as Sr 87 , have a very favorable figure of merit for such control due to narrow intercombination lines and the large hyperfine splitting in the excited states. We numerically study the quantum speed limit, optimal parameters, and the fidelity of arbitrary state preparation and full SU(10) maps, including the presence of decoherence due to optical pumping induced by the light-shifting laser. We also study the use of robust control to mitigate some dephasing due to inhomogeneities in the light shift. We find that with an rf Rabi frequency of Ω rf and 0.5% inhomogeneity in the the light shift we can prepare an arbitrary Haar-random state in a time T = 4.5 π / Ω rf with average fidelity ( F ψ ) = 0.9992 , and an arbitrary Haar-random SU(10) map in a time T = 24 π / Ω rf with average fidelity ( F U ) = 0.9923 .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Understanding Growth Dynamics and Yield Prediction of Sorghum Using High Temporal Resolution UAV Imagery Time Series and Machine Learning

Unmanned aerial vehicles (UAV) carrying multispectral cameras are increasingly being used for high-throughput phenotyping (HTP) of above-ground traits of crops to study genetic diversity, resource use efficiency and responses to abiotic or biotic stresses. There is significant unexplored potential for repeated data collection through a field season to reveal information on the rates of growth and provide predictions of the final yield. Generating such information early in the season would create opportunities for more efficient in-depth phenotyping and germplasm selection. This study tested the use of high-resolution time-series imagery (5 or 10 sampling dates) to understand the relationships between growth dynamics, temporal resolution and end-of-season above-ground biomass (AGB) in 869 diverse accessions of highly productive (mean AGB = 23.4 Mg/Ha), photoperiod sensitive sorghum. Canopy surface height (CSM), ground cover (GC), and five common spectral indices were considered as features of the crop phenotype. Spline curve fitting was used to integrate data from single flights into continuous time courses. Random Forest was used to predict end-of-season AGB from aerial imagery, and to identify the most informative variables driving predictions. Improved prediction of end-of-season AGB (RMSE reduction of 0.24 Mg/Ha) was achieved earlier in the growing season (10 to 20 days) by leveraging early- and mid-season measurement of the rate of change of geometric and spectral features. Early in the season, dynamic traits describing the rates of change of CSM and GC predicted end-of-season AGB best. Late in the season, CSM on a given date was the most influential predictor of end-of-season AGB. The power to predict end-of-season AGB was greatest at 50 days after planting, accounting for 63% of variance across this very diverse germplasm collection with modest error (RMSE 1.8 Mg/ha). End-of-season AGB could be predicted equally well when spline fitting was performed on data collected from five flights versus 10 flights over the growing season. This demonstrates a more valuable and efficient approach to using UAVs for HTP, while also proposing strategies to add further value.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Framework of compressive sensing and data compression for 4D-STEM

Four-dimensional Scanning Transmission Electron Microscopy (4D-STEM) is a powerful technique for high-resolution and high-precision materials characterization at multiple length scales, including the characterization of beam-sensitive materials. However, the field of view of 4D-STEM is relatively small, which in absence of live processing is limited by the data size required for storage. Furthermore, the rectilinear scan approach currently employed in 4D-STEM places a resolution- and signal-dependent dose limit for the study of beam sensitive materials. Improving 4D-STEM data and dose efficiency, by keeping the data size manageable while limiting the amount of electron dose, is thus critical for broader applications. Here we introduce a general method for reconstructing 4D-STEM data with subsampling in both real and reciprocal spaces at high fidelity. The approach is first tested on the subsampled datasets created from a full 4D-STEM dataset, and then demonstrated experimentally using random scan in real-space. The same reconstruction algorithm can also be used for compression of 4D-STEM datasets, leading to a large reduction (100 times or more) in data size, while retaining the fine features of 4D-STEM imaging, for crystalline samples.

4D-STEM↗
DOE OSTI2022

Model-Form Epistemic Uncertainty Quantification for Modeling with Differential Equations: Application to Epidemiology

Modeling real-world phenomena to any degree of accuracy is a challenge that the scientific research community has navigated since its foundation. Lack of information and limited computational and observational resources necessitate modeling assumptions which, when invalid, lead to model-form error (MFE). The work reported herein explored a novel method to represent model-form uncertainty (MFU) that combines Bayesian statistics with the emerging field of universal differential equations (UDEs). The fundamental principle behind UDEs is simple: use known equational forms that govern a dynamical system when you have them; then incorporate data-driven approaches – in this case neural networks (NNs) – embedded within the governing equations to learn the interacting terms that were underrepresented. Utilizing epidemiology as our motivating exemplar, this report will highlight the challenges of modeling novel infectious diseases while introducing ways to incorporate NN approximations to MFE. Prior to embarking on a Bayesian calibration, we first explored methods to augment the standard (non-Bayesian) UDE training procedure to account for uncertainty and increase robustness of training. In addition, it is often the case that uncertainty in observations is significant; this may be due to randomness or lack of precision in the measurement process. This uncertainty typically manifests as “noisy” observations which deviate from a true underlying signal. To account for such variability, the NN approximation to MFE is endowed with a probabilistic representation and is updated using available observational data in a Bayesian framework. By representing the MFU explicitly and deploying an embedded, data-driven model, this approach enables an agile, expressive, and interpretable method for representing MFU. In this report we will provide evidence that Bayesian UDEs show promise as a novel framework for any science-based, data-driven MFU representation; while emphasizing that significant advances must be made in the calibration of Bayesian NNs to ensure a robust calibration procedure.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Machine learning study of magnetism in uranium-based compounds

Actinide and lanthanide-based materials display exotic properties that originate from the presence of itinerant or localized f electrons and include unconventional superconductivity and magnetism, hidden order, and heavy-fermion behavior. Due to the strongly correlated nature of the 5f electrons, magnetic properties of these compounds depend sensitively on applied magnetic field and pressure, as well as on chemical doping. However, precise connection between the structure and magnetism in actinide-based materials is currently unclear. In this investigation, we established such structure-property links by assembling and mining two datasets that aggregate, respectively, the results of high-throughput density functional theory simulations and experimental measurements for the families of uranium- and neptunium-based binary compounds. Various regression algorithms were utilized to identify correlations among accessible attributes (features or descriptors) of the material systems and predict their cation magnetic moments and general forms of magnetic ordering. Descriptors representing compound structural parameters and cation f-subshell occupation numbers were identified as most important for accurate predictions. The best machine learning model developed employs the random forest regression algorithm. It can predict both spin and orbit moment size with root-mean-square error of 0.17 μ B and 0.19 μ B , respectively. Lastly, the random forest classification algorithm is used to predict the ordering (paramagnetic, ferromagnetic, and antiferromagnetic) of such systems with 76% accuracy.

36 MATERIALS SCIENCE↗