Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model Inference”

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 919 records · Page 51

Probing Sizes and Shapes of Nobelium Isotopes by Laser Spectroscopy

Until recently, ground-state nuclear moments of the heaviest nuclei could only be inferred from nuclear spectroscopy, where model assumptions are required. Here, laser spectroscopy in combination with modern atomic structure calculations is now able to probe these moments directly, in a comprehensive and nuclear- model-independent way. In this work, we reported on unique access to the differential mean-square charge radii of 252,253,254 No, and therefore to changes in nuclear size and shape. State-of-the-art nuclear density functional calculations describe well the changes in nuclear charge radii in the region of the heavy actinides, indicating an appreciable central depression in the deformed proton density distribution in 252,254 No isotopes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Magnetoreflection spectroscopy of monolayer transition-metal dichalcogenide semiconductors in pulsed magnetic fields

The authors describe recent experimental efforts to perform polarization-resolved optical spectroscopy of monolayer transition-metal dichalcogenide semiconductors in very large pulsed magnetic fields to 65 T. The experimental setup and technical challenges are discussed in detail, and temperature-dependent magnetoreflection spectra from atomically thin tungsten disulphide are presented. The data clearly reveal not only the valley Zeeman effect in these two-dimensional semiconductors but also the small quadratic exciton diamagnetic shift from which the very small exciton size can be directly inferred. Lastly, the authors present model calculations that demonstrate how the measured diamagnetic shifts can be used to constrain estimates of the exciton binding energy in this new family of monolayer semiconductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A Near-Complete Haplotype-Phased Genome of the Dikaryotic Wheat Stripe Rust Fungus Puccinia striiformis f. sp. tritici Reveals High Interhaplotype Diversity

A long-standing biological question is how evolution has shaped the genomic architecture of dikaryotic fungi. To answer this, high-quality genomic resources that enable haplotype comparisons are essential. Short-read genome assemblies for dikaryotic fungi are highly fragmented and lack haplotype-specific information due to the high heterozygosity and repeat content of these genomes. Here, we present a diploid-aware assembly of the wheat stripe rust fungus Puccinia striiformis f. sp. tritici based on long reads using the FALCON-Unzip assembler. Transcriptome sequencing data sets were used to infer high-quality gene models and identify virulence genes involved in plant infection referred to as effectors. This represents the most complete Puccinia striiformis f. sp. tritici genome assembly to date (83 Mb, 156 contigs, N 50 of 1.5 Mb) and provides phased haplotype information for over 92% of the genome. Comparisons of the phase blocks revealed high interhaplotype diversity of over 6%. More than 25% of all genes lack a clear allelic counterpart. When we investigated genome features that potentially promote the rapid evolution of virulence, we found that candidate effector genes are spatially associated with conserved genes commonly found in basidiomycetes. Yet, candidate effectors that lack an allelic counterpart are more distant from conserved genes than allelic candidate effectors and are less likely to be evolutionarily conserved within the P. striiformis species complex and Pucciniales. In summary, this haplotype-phased assembly enabled us to discover novel genome features of a dikaryotic plant-pathogenic fungus previously hidden in collapsed and fragmented genome assemblies.

54 ENVIRONMENTAL SCIENCES↗

A Bayesian Approach for Characterizing and Mitigating Gate and Measurement Errors

Various noise models have been developed in quantum computing study to describe the propagation and effect of the noise that is caused by imperfect implementation of hardware. Identifying parameters such as gate and readout error rates is critical to these models. We use a Bayesian inference approach to identify posterior distributions of these parameters such that they can be characterized more elaborately. By characterizing the device errors in this way, we can further improve the accuracy of quantum error mitigation. Experiments conducted on IBM’s quantum computing devices suggest that our approach provides better error mitigation performance than existing techniques used by the vendor. Also, our approach outperforms the standard Bayesian inference method in some scenarios.

97 MATHEMATICS AND COMPUTING↗

Fast Adaptive Neural Control of Resonant Extraction at Fermilab

We present progress on the development of a machine learning (ML) regulation system for third-order resonant extraction of the beam delivered to the Mu2e experiment at Fermilab. We consider classical and ML-based controllers optimized on semi-analytic simulations and provide performance comparisons for several models. Additionally, we discuss the efficiency of each model in training, which has implications for future work on adaptive control. We also discuss progress on developing optimized implementations of ML models for edge-based inference.

Whitbeck, A. [Fermilab] (ORCID:0000000342245164)↗

Improved Sampling Algorithms in the Risk-Informed Safety Margin Characterization Toolkit

The RISMC approach is developing advanced set of methodologies and algorithms in order to perform Probabilistic Risk Analyses (PRAs). In contrast to classical PRA methods, which are based on Event-Tree and Fault-Tree methods, the RISMC approach largely employs system simulator codes applied to stochastic analysis tools. The basic idea is to randomly perturb (by employing sampling algorithms) timing and sequencing of events and internal parameters of the system codes (i.e., uncertain parameters) in order to estimate stochastic parameters such as core damage probability. This approach applied to complex systems such as nuclear power plants requires to perform a series of computationally expensive simulation runs given a large set of uncertain parameters. These types of analysis are affected by two issues. Firstly, the space of the possible solutions (a.k.a., the issue space or the response surface) can be sampled only very sparsely, and this precludes the ability to fully analyze the impact of uncertainties on the system dynamics. Secondly, large amounts of data are generated and tools to generate knowledge from such data sets are not yet available. This report focuses on the first issue and in particular employs novel methods that optimize the information generated by the sampling process by sampling unexplored and risk-significant regions of the issue space: adaptive (smart) sampling algorithms. They infer system response from surrogate models constructed from existing samples and predict the most relevant location of the next sample. It is therefore possible to understand features of the issue space with a small number of carefully selected samples. In this report, we will present how it is possible to perform adaptive sampling using the RISMC toolkit and highlight the advantages compared to more classical sampling approaches such Monte-Carlo. We will employ RAVEN to perform such statistical analyses using both analytical cases but also another RISMC code: RELAP-7.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Characterizing Fracture Property Using Resistivity Measured at Different Frequencies

The objective was to develop geophysical approaches to detecting and evaluating the fractures created or existing in EGS and other geothermal reservoirs by measuring the resistivity at different frequencies. This project has been divided into two phases: Phase I (first year): Proof of Concept – develop the resistivity approach and verify the effect of frequency on the resistivity in rocks with artificial or natural fractures over a wide range of frequencies. Phase II: Prototyping Part 1 (second year): measure the resistivity in rocks with fractures of different apertures, different length, and different configurations at different frequencies. Part 2 (third year): develop mathematical models and the resistivity method; infer the fracture properties using the measured resistivity data.

15 GEOTHERMAL ENERGY↗

Exploiting Cross-sensitivity by Bayesian Decoding of Mixed Potential Sensor Arrays

LANL mixed-potential electrochemical sensor (MPES) device arrays were coupled with advanced Bayesian inference treatment of the physical model of relevant sensor-analyte interactions. We demonstrated that our approach could be used to uniquely discriminate the composition of ternary gas sensors with three discreet MPES sensors with an average error of less than 2%. We also observed that the MPES exhibited excellent stability over a year of operation at elevated temperatures in the presence of test gases.

47 OTHER INSTRUMENTATION↗

Hierarchical Linear Regression (Rev. 1)

We describe an algorithm for making inferences in a hierarchical linear model. The model allows different slopes and intercepts for different groups, and ties together the different slopes by representing them as samples from a higher level distribution. (The intercepts are treated as statistically independent.) Our model and code permit many options for the variance of this hierarchical distribution for the slopes, including gamma distributions on the precision, t-distributions on the standard error, pooled slopes (all equal), and independent slopes. The algorithm is described in the paper; an implementation of the algorithm in R is available on request.

97 MATHEMATICS AND COMPUTING↗

Analysis of the effect of spherical aberrations on a solenoid scan

Solenoid scans are often used to determine the beam emittance of an intense relativistic electron beam (IREB). In a solenoid scan, the beam size is measured at a location downstream of the solenoid as a function of the strength of the solenoid magnetic field. The emittance is then inferred using a beam optics model of the transport of the beam through the solenoid magnet to the imaging location. This requires precise knowledge of the field strength and shape as a function of the current supplied to the solenoid. The beam optics model for an IREB needs to include an accurate model of the magnetic field. Unfortunately, most beam optics models do not include the effect of spherical aberrations which are inherent in a solenoid magnet. Spherical aberrations are manifested as an essentially quadratic increase in the solenoid magnetic field as a function of radius. One consequence of this is that the focal length will be dependent on the beam radius. This note examines the sensitivity of a typical solenoid scan performed on DARHT axis 2 to spherical aberrations. The sensitivity to spherical aberrations is examined as a function of the beam envelope size and the drift distance from the solenoid magnet to the imaging station.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Geant4 Event Biasing and Fast Simulation

Geant4 offers advanced event biasing techniques to significantly accelerate simulations involving rare events. Various biasing methods, such as leading particle selection, cross-section biasing, radioactive decay enhancement, and bremsstrahlung splitting, enable efficient event sampling, though they require careful handling. Additionally, Geant4 provides a Fast Simulation Interface, allowing the replacement of standard processes in specific region and for selected particles, enabling faster execution or external code integration. Applications of fast simulation include electromagnetic shower modeling in calorimeters, machine learning inference, and offloading tasks to specialized hardware like GPUs, making Geant4 a powerful tool for computationally demanding simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Fast Adaptive Neural Control of Resonant Extraction at Fermilab

We present progress on the development of a machine learning (ML) regulation system for third-order resonant extraction of the beam delivered to the Mu2e experiment at Fermilab. We consider classical and ML-based controllers optimized on semi-analytic simulations and provide perfor- mance comparisons for several models. Additionally, we discuss the efficiency of each model in training, which has implications for future work on adaptive control. We also discuss progress on developing optimized implementations of ML models for edge-based inference.

Whitbeck, A. [Fermilab]↗

Clustering of DESI galaxies split by thermal Sunyaev-Zeldovich effect

The thermal Sunyaev-Zeldovich (tSZ) effect is associated with galaxy clusters - extremely large and dense structures tracing the dark matter with a higher bias than isolated galaxies. We propose to use the tSZ data to separate galaxies from redshift surveys into distinct subpopulations corresponding to different densities and biases independently of the redshift survey systematics. Leveraging the information from different environments, as in density-split and density-marked clustering, is known to tighten the constraints on cosmological parameters, like $\Omega_m$, $\sigma_8$ and neutrino mass. We use data from the Dark Energy Spectroscopic Instrument (DESI) and the Atacama Cosmology Telescope (ACT) in their region of overlap to demonstrate informative tSZ splitting of Luminous Red Galaxies (LRGs). We discover a significant increase in the large-scale clustering of DESI LRGs corresponding to detections starting from 1-2 sigma in the ACT DR6 + Planck tSZ Compton-$y$ map, below the cluster candidate threshold (4 sigma). We also find that such galaxies have higher line-of-sight coordinate (and velocity) dispersions and a higher number of close neighbors than both the full sample and near-zero tSZ regions. We produce simple simulations of tSZ maps that are intrinsically consistent with galaxy catalogs and do not include systematic effects, and find a similar pattern of large-scale clustering enhancement with tSZ effect significance. Moreover, we observe that this relative bias pattern remains largely unchanged with variations in the galaxy-halo connection model in our simulations. This is promising for future cosmological inference from tSZ-split clustering with semi-analytical models. Thus, we demonstrate that valuable cosmological information is present in the lower signal-to-noise regions of the thermal Sunyaev-Zeldovich map, extending far beyond the individual cluster candidates.

Astronomy data analysis↗

On the robustness of the pendulum model for large-amplitude longitudinal oscillations in prominences

Large-amplitude longitudinal oscillations (LALOs) in prominences are spectacular manifestations of solar activity. In such events nearby energetic disturbances induce periodic motions on filaments with displacements comparable to the size of the filaments themselves and with velocities larger than 20 km s{sup −1}. The pendulum model, in which the gravity projected along a rigid magnetic field is the restoring force, was proposed to explain these events. However, it can be objected that in a realistic situation where the magnetic field reacts to the mass motion of the heavy prominence, the simplified pendulum model could be no longer valid. We have performed nonlinear time-dependent numerical simulations of LALOs considering a dipped magnetic field line structure. In this work we demonstrate that for even relatively weak magnetic fields the pendulum model works very well. We therefore validate the pendulum model and show its robustness, with important implications for prominence seismology purposes. With this model it is possible to infer the geometry of the dipped field lines that support the prominence.

79 ASTRONOMY AND ASTROPHYSICS↗

XRISM Finds the Changing-look Active Galactic Nucleus NGC 1365 in an Extended Low State: A Dense, Highly Ionized Outflow Obscures the Central Source

We present the first XRISM/Resolve observations of the active galactic nucleus NGC 1365, obtained in 2024 February and July. NGC 1365 is known for rapid transitions between Compton-thick and Compton-thin states, along with strong absorption from a highly ionized wind. During our observations, the source was found in a persistent low-flux state, characterized by a decrease in hard-X-ray luminosity and significant line-of-sight obscuration. In this state, XRISM/Resolve reveals clear Fe XXV and Fe XXVI absorption lines together with, for the first time in this source, corresponding emission lines. These features may arise either from reemission from a photoionized wind (P Cygni profile) or from collisionally ionized gas associated with outflow-driven shocks in the interstellar medium. We estimate the wind launch radius to be approximately 10 16 cm (∼10 4 R g ), consistent with the location of the X-ray broadline region. We also resolve a broadened Fe Kα line by σ ∼ 1300 km s −1 , placing it at similar scales to the wind, consistent with radii inferred from disk-broadening models and the variability of the Fe Kα broad line. The similarity of the Fe Kα profile to the Hβ wing and broad Paα width indicates that the X-ray-emitting region is likely cospatial with the optical/infrared broadline region and originates from the same gas.

Astronomy and AstroPhysics↗

Probabilistic estimation of depth-resolved profiles of soil thermal diffusivity from temperature time series

Abstract. Improving the quantification of soil thermal and physical properties is key to achieving a better understanding and prediction of soil hydro-biogeochemical processes and their responses to changes in atmospheric forcing. Obtaining such information at numerous locations and/or over time with conventional soil sampling is challenging. The increasing availability of low-cost, vertically resolved temperature sensor arrays offers promise for improving the estimation of soil thermal properties from temperature time series, and the possible indirect estimation of physical properties. Still, the reliability and limitations of such an approach need to be assessed. In the present study, we develop a parameter estimation approach based on a combination of thermal modeling, sliding time windows, Bayesian inference, and Markov chain Monte Carlo simulation to estimate thermal diffusivity and its uncertainty over time, at numerous locations and at an unprecedented vertical spatial resolution (i.e., down to 5 to 10 cm vertical resolution) from soil temperature time series. We provide the necessary framework to assess under which environmental conditions (soil temperature gradient, fluctuations, and trend), temperature sensor characteristics (bias and level of noise), and deployment geometries (sensor number and position) soil thermal diffusivity can be reliably inferred. We validate the method with synthetic experiments and field studies. The synthetic experiments show that in the presence of median diurnal fluctuations ≥ 1.5 ∘C at 5 cm below the ground surface, temperature gradients > 2 ∘C m−1, and a sliding time window of at least 4 d the proposed method provides reliable depth-resolved thermal diffusivity estimates with percentage errors ≤ 10 % and posterior relative standard deviations ≤ 5 % up to 1 m depth. Reliable thermal diffusivity under such environmental conditions also requires temperature sensors to be spaced precisely (with accuracy to a few millimeters), with a level of noise ≤ 0.02 ∘C, and with a bias defined by a standard deviation ≤ 0.01 ∘C. Finally, the application of the developed approach to field data indicates significant repeatability in results and similarity with independent measurements, as well as promise in using a sliding time window to estimate temporal changes in soil thermal diffusivity, as needed to potentially capture changes in bulk density or water content.

54 ENVIRONMENTAL SCIENCES↗

Tensor-Dictionary Learning with Deep Kruskal-Factor Analysis

We introduce new dictionary learning methods for tensor-variate data of any order. We represent each data item as a sum of Kruskal decomposed dictionary atoms within the framework of beta-process factor analysis (BPFA). Our model is nonparametric and can infer the tensor-rank of each dictionary atom. This Kruskal-Factor Analysis (KFA) is a natural generalization of BPFA. We also extend KFA to a deep convolutional setting and develop online learning methods. We test our approach on image processing and classification tasks achieving state of the art results for 2D & 3D inpainting and Caltech 101. The experiments also show that atom-rank impacts both overcompleteness and sparsity.

Stevens, Andrew J.↗

Bayesian inference for plasmonic nanometrology

Here, we introduce a Bayesian method for the characterization of plasmonic nanoparticles, which is applicable to both near- and far-field problems. Designed to combine data generated from any photon-plasmon interaction experiment with physically motivated theoretical models, our approach leverages state-of-the-art Markov chain Monte Carlo sampling techniques and returns parameter estimates on nanometric scales. Simulated spectral data sets, describing resonant scattering of photons from ellipsoidal and toroidal nanoparticles, are explored as concrete examples of our approach, with the resulting Bayesian estimates showing excellent agreement with the ground truth, even under conditions of high statistical noise. By incorporating Bayes factors into the method as well, we reveal how model selection can determine which one of competing geometric shapes better explains the observed data. Our comprehensive nanometrology procedure can be tailored to a variety of light-particle interaction models, and its reliance on Bayesian inference furnishes automatic uncertainty quantification. In addition to applicability to a host of plasmonic configurations such as nanoparticle dimers, trimers, and array studies, it is proposed that the presented analysis can be extended to the quantum regime, where nonclassical photon statistics may provide additional insight for inference of scatterer properties.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗