Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “GAUSSIAN DISTRIBUTION”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Impact of super-Gaussian electron distributions on plasma K-shell emission

Electron distributions in laser-produced plasmas will be driven toward a super-Gaussian distribution due to inverse bremsstrahlung absorption [Langdon, Phys. Rev. Lett. 44, 575 (1980)]. Both theoretical and experimental evidence suggest that fundamental plasma properties are altered by the super-Gaussian distribution. Here, this paper examines how the super-Gaussian distribution affects the ionization balance and K-shell emission of atomic plasmas, utilizing approximate formulas and detailed collisional-radiative simulations. While the impact on plasma ionization is small, K-shell spectra can be significantly modified. Based on these findings, we demonstrate that K-shell spectroscopy can be used to infer super-Gaussian or other similar nonequilibrium electron distributions.

Atomic spectra↗
DOE OSTI2026

Algorithm to extract direction in 2D discrete distributions and a continuous Frobenius norm

In this study, we present a novel algorithm for determining directionality in 2D distributions of discrete data. We compare a reference dataset with a known direction to a measured dataset with an unknown direction by the Frobenius norm of the difference (FND) to find the unknown direction. To generalize this concept, we develop a continuous Frobenius norm of the difference (CFND) as a continuous analog of the FND and derive its analytical expression. By relating fitted and normalized 2D Gaussian distributions, we show that the CFND approximates the FND, and we validate this relationship with computer simulations. We find that a first-order approximation of the CFND between two similar Gaussian distributions takes the form of an absolute sine function, offering a simple analytical form with potential applications in specialized areas such as segmented inverse beta decay neutrino detectors, astronomy, machine learning, and more. Our methodology consists of modeling a 2D Gaussian distribution, binning the data into a histogram, and encoding it as a square matrix. Rotating this matrix around its geometric center and comparing it to a measured dataset using the FND gives us rotational data that we fit with an absolute sine function. The location of the minimum of this fit is the angle closest to the true angle of the direction in the measured dataset. We present the derivation and discuss initial applications of the CFND in our novel algorithm, demonstrating its success in approximating directionality in 2D distributions.

Physics↗
DOE OSTI2025

Algorithm to extract direction in 2D discrete distributions and a continuous Frobenius norm

In this study, we present a novel algorithm for determining directionality in 2D distributions of discrete data. We compare a reference dataset with a known direction to a measured dataset with an unknown direction by the Frobenius norm of the difference (FND) to find the unknown direction. To generalize this concept, we develop a continuous Frobenius norm of the difference (CFND) as a continuous analog of the FND and derive its analytical expression. By relating fitted and normalized 2D Gaussian distributions, we show that the CFND approximates the FND, and we validate this relationship with computer simulations. We find that a first-order approximation of the CFND between two similar Gaussian distributions takes the form of an absolute sine function, offering a simple analytical form with potential for specialized applications in segmented inverse beta decay (IBD) neutrino detectors, astronomy, machine learning, and more. Although this method may easily extend to 3D scalar fields, our focus here is on 2D real-valued fields as it directly applies to directionality. Our methodology consists of modeling a 2D Gaussian distribution, binning the data into a histogram, and encoding it as a square matrix. Rotating this matrix around its geometric center and comparing it to a measured dataset using the FND gives us rotational data that we fit with an absolute sine function. The location of the minimum of this fit is the angle closest to the true angle of the direction in the measured dataset. We present the derivation and discuss initial applications of the CFND in our novel algorithm, demonstrating its success in approximating directionality in 2D distributions.

Data Analysis, Statistics and Probability (physics↗
DOE OSTI2026

MAGIC: M arching Cubes Isosurface Uncertainty Visualization for G auss i an Uncertain Data With Spatial C orrelation

Here, in this paper, we study the propagation of data uncertainty through the marching cubes algorithm for isosurface visualization for correlated uncertain data. Consideration of correlation has been shown paramount for avoiding errors in uncertainty quantification and visualization in multiple prior studies. Although the problem of isosurface uncertainty with spatial data correlation has been previously addressed, there are two major limitations to prior treatments. First, there are no analytical formulations for uncertainty quantification of isosurfaces when the data uncertainty is characterized by a Gaussian distribution with spatial correlation. Second, as a consequence of the lack of analytical formulations,existing techniques resort to a Monte Carlo sampling approach, which is expensive and difficult to integrate into visualization tools. To address these limitations, we present a closed-form framework to efficiently derive uncertainty in marching cubes level-sets for Gaussian uncertain data with spatial correlation (MAGIC). To derive closed-form solutions, we leverage the Hinkley's derivation on the ratio of Gaussian distributions. With our analytical framework, we achieve a significant speed-up and enhanced accuracy of uncertainty quantification over classical Monte Carlo methods. We further accelerate our analytical solutions using many-core processors to achieve speed-ups up to 585× and integrability with production visualization tools for broader impact. We demonstrate the effectiveness of our correlation-aware uncertainty framework through experiments on meteorology, urban flow, and astrophysics simulation datasets.

Gaussian↗
DOE OSTI2025

Discrete generative diffusion models without stochastic differential equations: A tensor network approach

Diffusion models (DMs) are a class of generative machine learning methods that sample a target distribution by transforming samples of a trivial (often Gaussian) distribution using a learned stochastic differential equation. In standard DMs, this is done by learning a “score function” that reverses the effect of adding diffusive noise to the distribution of interest. Here we consider the generalisation of DMs to lattice systems with discrete degrees of freedom, and where noise is added via Markov chain jump dynamics. We show how to use tensor networks (TNs) to efficiently define and sample such “discrete diffusion models” (DDMs) without explicitly having to solve a stochastic differential equation. We show the following: (i) by parametrising the data and evolution operators as TNs, the denoising dynamics can be represented exactly; (ii) the auto-regressive nature of TNs allows to generate samples efficiently and without bias; (iii) for sampling Boltzmann-like distributions, TNs allow to construct an efficient learning scheme that integrates well with Monte Carlo. We illustrate this approach to study the equilibrium of two models with non-trivial thermodynamics, the d = 1 constrained Fredkin chain and the d = 2 Ising model. Published by the American Physical Society 2025

Causer, Luke (ORCID:0000000194243473)↗
DOE OSTI2025

Bulk synthesis of radiation resistant W – Ti – Cr – V compositionally complex alloys

Refractory compositionally complex alloys are candidate material systems for next generation advanced nuclear reactors. This work showcases the first successful bulk synthesis of low activation W–Ti based refractory compositionally complex alloys using arc melting and provides insights on using additive manufacturing for these compositions using directed energy deposition. Both techniques produce equiaxed grains composed of a tungsten matrix with Ti–V–Cr dendritic boundaries. The arc melted specimen possesses a multi-modal grain size distribution, while the directed energy deposition specimen possesses a more gaussian distribution of grain size. Both arc melted and directed energy deposition specimens demonstrate high thermal stability up to 900 °C, as well as promising radiation resistance with low loop formation and the presence of homogeneously distributed helium cavities maintaining small diameters at ≥10 dpa under simultaneous light (helium) and heavy (krypton) ion irradiation at 900 °C.

Arc melting (AM)↗
DOE OSTI2026

LANL-TAMUS Graduate Fellowship Report

The goal of this project is to conduct analysis of 147 Sm neutron transmission (n,tot) data at DICER alongside previously collected neutron capture (n,γ) data and report updated resonance parameters and spins. Measurement and analysis of neutron-capture data, without neutron transmission data, has been conducted before and showed that the gamma widths, Γ γ , and neutron widths, Γn, both deviate from the distributions predicted by well-established theory. Γ γ is expected to follow a narrow Gaussian distribution independent of energy but instead a broadening of the distribution around 350 eV incident neutron energy was observed (P. E. Koehler et al. Phys. Rev. C 108, 142502 (2012)). Γ n is predicted to follow a chi-squared distribution with one degree of freedom known as the Porter-Thomas distribution but this prediction also failed above 350 eV incident neutron energy (P. E. Koehler et al. Phys. Rev. C 76, 025804 (2007)). This project seeks to reassess this past measurement by including neutron transmission data and simultaneously analyzing them alongside the neutron capture data to extract the most accurate measure of Γ γ and Γ n ever taken for the resonances of 147 Sm and see if theory holds or fails under a more refined measurement with newer data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Photometry of Outer Solar System Objects from the Dark Energy Survey. II. A Joint Analysis of Trans-Neptunian Absolute Magnitudes, Colors, Light Curves and Dynamics

For the 696 trans-Neptunian objects (TNOs) with absolute magnitudes 5.5 < H r < 8.2 detected in the Dark Energy Survey, we characterize the relationships between their dynamical state and physical properties—namely H r , indicating size; colors, indicating surface composition; and flux variation semiamplitude A, indicating asphericity and surface inhomogeneity. We seek “birth” physical distributions that can recreate these parameters in every dynamical class. We show that the observed colors of these TNOs are consistent with two Gaussian distributions in griz space, “near-infrared bright” (NIRB) and “near-infrared faint” (NIRF), presumably an inner and outer birth population, respectively. We find a model in which both the NIRB and NIRF H r and A distributions are independent of current dynamical states, supporting their assignment as birth populations. All objects are consistent with a common rolling p(H r ), but NIRF objects are significantly more variable. Cold classicals (CCs) are purely NIRF, while hot classical (HC), scattered, and detached TNOs are consistent with ≈ 70% NIRB and the resonance NIRB fractions show significant variation. The NIRB components of the HCs and of some resonances have broader inclination distributions than the NIRFs, i.e. their current dynamics retains information about birth location. We find evidence for radial stratification within the birth NIRB population, in that HC NIRBs are on average redder than detached or scattered NIRBs; a similar effect distinguishes CCs from other NIRFs. We estimate total object counts and masses of each class within our H r range. These results will strongly constrain models of the outer solar system.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Solving high-dimensional inverse problems using amortized likelihood-free inference with noisy and incomplete data

Here, we present a likelihood-free probabilistic inversion method based on normalizing flows for high-dimensional inverse problems. The proposed method is composed of two complementary networks: a summary network for data compression and an inference network for parameter estimation. The summary network encodes raw observations into a fixed-size vector of summary features, while the inference network generates samples of the approximate posterior distribution of the model parameters based on these summary features. The posterior samples are produced in a deep generative fashion by sampling from a latent Gaussian distribution and passing these samples through an invertible transformation. We construct this invertible transformation by sequentially alternating conditional invertible neural network and conditional neural spline flow layers. The summary and inference networks are trained simultaneously. We apply the proposed method to an inversion problem in groundwater hydrology to estimate the posterior distribution of the log-conductivity field conditioned on spatially sparse time-series observations of the system’s hydraulic head responses. The conductivity field is represented with 706 degrees of freedom in the considered problem. Comparison with the likelihood-based iterative ensemble smoother PEST-IES method demonstrates that the proposed method accurately estimates the parameter posterior distribution and the observations’ predictive posterior distribution at a fraction of the inference time of PEST-IES.

conditional invertible neural network↗
DOE OSTI2025

Energy migration and scintillation kinetics in compositionally complex (Gd 1/4 Y 1/4 Tb 1/4 Lu 1/4 ) 3 Al 5 O 12 :Ce single crystal scintillator

It is well-established that compositional tuning through binary admixture can improve scintillation performance in several materials systems, including Ce-activated garnets. Although recent work on ternary or quaternary cation admixture shows promise, the impact of this increased compositional complexity on thermal stability and carrier-defect dynamics has not been addressed. Here, we investigate a compositionally complex garnet, (Gd 1/4 Y 1/4 Tb 1/4 Lu 1/4 ) 3 Al 5 O 12 :Ce (GYTLAG), grown by the Czochralski method using temperature-dependent photoluminescence (PL), PL decay, and thermoluminescence (TL). PL and PL decay measurements support a thermally activated Tb 3+ -Ce 3+ energy transfer, where Tb 3+ emission dominates below 60 K, but Ce 3+ emission increases from 20-300 K. Thermal quenching of Ce 3+ emission occurs around T 50 = 508 K, with an activation energy of 0.6 eV. TL and wavelength-resolved TL spectra from 20–500 K show that GYTLAG contains similar trap groups to LuAG but with a broader distribution of glow peaks below room temperature, possibly caused by quaternary cation mixing. A combination of dose dependence, partial cleaning and initial rise, and glow curve fitting to a first order continuous Gaussian distribution model are used to understand the contribution of electronic point defects to scintillation decay and afterglow at room temperature. Furthermore, these results inform how increased compositional complexity influences recombination dynamics in garnet scintillators.

Compositionally complex↗
DOE OSTI2026

Exploring sodium dynamics in the dilute oxygen regime of mixed oxy-sulfide NaPSiSO glasses

All-solid-state‑sodium batteries have great potential to lower the cost of grid-scale energy storage systems. However, they require a low cost, high conductivity solid electrolyte to become a reality. The electrochemical properties of mixed oxy-sulfide glasses NaPSiSO exceed those of either the pure sulfide or the pure oxide glass alone. In this work, we have tested models of sodium ion dynamics in the NaPSiSO series: zNa 2 S+(1-z)[(1-y)[(1-x)SiS 2 + xPS 5/2 ] + yNaPO 3 ] where 0.57 ≤ z ≤ 0.588, 0 ≤ x ≤ 0.35, and 0.15 ≤ y ≤ 0.35 using 23 Na nuclear magnetic resonance to directly probe the local fields and their fluctuations at the nuclear sites. We determined the distribution of activation energies in these materials arising from distributions of local environments and provided quantitative information about the connectivity of the ionic pathways through which the sodium ions move. We found that the sodium ion conduction dynamics were described as an ionic percolation through a Gaussian distribution of energy barriers and that the composition dependence of the ionic conductivity in the dilute oxygen regime is well-described by the Nernst-Einstein relation modified for percolation.

25 ENERGY STORAGE↗
DOE OSTI2026

Statistics of base polytopes in F-theory

We propose a new statistical ensemble of toric bases for elliptic Calabi-Yaus used in F-theory models, by focusing on only the convex hull of the base, i.e., the base polytope. This physically motivated coarse-graining greatly simplifies the combinatorial complexity of the part of the 4d F-theory landscape with toric bases. We develop a Monte Carlo approach that randomly samples the base polytopes within fixed boxes, with proper statistical weights. We first apply the algorithm to the set of 2d base polytopes, generating an enlarged set of toric 2d bases that include certain types of codimension-two (4,6) points, and we validate our approach against exact numbers. We then explore the set of 3d base polytopes which fit in a set of “maximal” 3d boxes, and estimate the total number of inequivalent 3d base polytopes to be 10 85 –10 90 . We provide statistical data such as the distribution of non-Higgsable gauge groups on these bases. Amusingly, a similar method can also be applied to generate reflexive polytopes in various dimensions. In both the reflexive and base polytope cases, the number of relevant polytopes obeys a Gaussian distribution as a function of the number of vertices, which can be understood in terms of other results on random polytopes in the math literature.

Differential and algebraic geometry↗
DOE OSTI2025

Search for electroweak-scale dijet resonances using trigger-level analysis with the ATLAS detector in 132 fb −1 of 𝑝⁢𝑝 collisions at $\sqrt{𝑠}$ = 13 TeV

This article reports on a search for dijet resonances using 132 fb −1 of 𝑝⁢𝑝 collision data recorded at $\sqrt{𝑠}$ = 13 TeV by the ATLAS detector at the Large Hadron Collider. The search is performed solely on jets reconstructed within the ATLAS trigger to overcome bandwidth limitations imposed on conventional single-jet triggers, which would otherwise reject data from decays of sub-TeV dijet resonances. Collision events with two jets satisfying transverse momentum thresholds of 𝑝 T ≥ 85 GeV and jet rapidity separation of |𝑦*| <0.6 are analysed for dijet resonances with invariant masses from 375 to 1800 GeV. A data-driven background estimate is used to model the dijet mass distribution from multijet processes. No significant excess above the expected background is observed. Upper limits are set at 95% confidence level on coupling values for a benchmark leptophobic axial-vector 𝑍′ model and on the production cross section for a new resonance contributing a Gaussian-distributed line-shape to the dijet mass distribution.

hadron colliders↗
DOE OSTI2025

Uncertainty quantification of a physics-informed model based on sparse identification of a Thermal Energy Distribution System

Integrated energy systems (IES)s are crucial for enhancing the economy and efficiency of power generation sources (e.g., nuclear energy) necessary to unleash American energy dominance. These systems can be integrated with thermal energy storage (TES) and intermittent renewable energies to optimize overall energy use, peak-load regulation, and demand-side responses. However, the stabilization of energy generation, transport, and utilization introduces operational complexities that exceed the challenges of managing each sub-component individually. Currently, though IESs rely on human operators for efficiency and stability, reducing human error risk and enhancing performance through automation is highly desirable. Recent advances at Idaho National Laboratory have demonstrated successful control of the Thermal Energy Distributed System (TEDS). However, the automatic control system depends on a deterministic Sparse Identification of Nonlinear Dynamics with Control (SINDyC) model, which are trained based on simulation data from physics-based simulations. Because of uncertainties in physics-based simulation, SINDyC model results in large discrepancies against experimental data and cannot be reliably used in automatic control. In this paper, we present an innovative approach to address these discrepancies by quantifying uncertainties and developing a more robust model. We first generated trajectories by using first-principles physics codes to encapsulate the experiment. Next, we trained thousands of models by randomly sampling these trajectories. We then collapsed all those models into one probabilistic SINDyC by fitting a multivariate Gaussian distribution onto the resulting coefficient’s distribution. Despite its simplicity, our approach successfully produced 95% confidence intervals that captured the experimental trajectories. It even did so with a higher probability and better U-pooling score across six of the seven relevant quantities of interest (QoIs), as compared to other classical approaches. In conclusion, ongoing research is focusing on generating new experimental trajectories to validate this approach, and on employing Bayesian calibration to refine parametric uncertainties and guide future model development efforts.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Machine learning-assisted profiling of a kinked ladder polymer structure using scattering

Ladder polymers consisting of fused rings in the backbone have very limited conformational freedom, which results in very different properties from traditional linear polymers. However, accurately determining their size and chain conformations from solution scattering remains a challenge. Their chain conformations of kinked ladder polymers are largely governed by the structures and relative orientations or configurations of the repeat units, unlike conventional polymer chains whose bending angles between repeat units follow a unimodal Gaussian distribution. Meanwhile, traditional scattering models for polymer chains do not account for these unique structural features. This work introduces a novel approach that integrates machine learning with Monte Carlo simulations to construct a model that can describe the geometry of a type of kinked CANAL ladder polymers. We first develop a Monte Carlo simulation model for sampling the configuration space of CANAL ladder polymers, where each repeat unit is modeled as a biaxial segment. Then, we establish a machine learning-assisted scattering analysis framework based on Gaussian Process Regression. Finally, we conduct small-angle neutron scattering experiments on a CANAL ladder polymer solution to apply our approach. Our method uncovers structural features of such ladder polymers that conventional methods fail to capture.

Ding, Lijie [Oak Ridge National Laboratory (ORNL),↗
DOE OSTI2025

Exploring the energy landscape of RBMs: reciprocal space insights into bosons, hierarchical learning and symmetry breaking

Deep generative models have become ubiquitous due to their ability to learn and sample from complex distributions. Despite the proliferation of various frameworks, the relationships among these models remain largely unexplored, a gap that hinders the development of a unified theory of AI learning. In this work, we address two central challenges: clarifying the connections between different deep generative models and deepening our understanding of their learning mechanisms. We focus on Restricted Boltzmann Machines (RBMs), a class of generative models known for their universal approximation capabilities for discrete distributions. By introducing a reciprocal space formulation for RBMs, we reveal a connection between these models, diffusion processes, and systems of coupled bosons. Our analysis shows that at initialization, the RBM operates at a saddle point, where the local curvature is determined by the singular values of the weight matrix, whose distribution follows the Marc̆enko-Pastur law and exhibits rotational symmetry. During training, this rotational symmetry is broken due to hierarchical learning, where different degrees of freedom progressively capture features at multiple levels of abstraction. This leads to a symmetry breaking in the energy landscape, reminiscent of Landau’s theory. This symmetry breaking in the energy landscape is characterized by the singular values and the weight matrix eigenvector matrix. We derive the corresponding free energy in a mean-field approximation. We show that in the limit of infinite size RBM, the reciprocal variables are Gaussian distributed. Our findings indicate that in this regime, there will be some modes for which the diffusion process will not converge to the Boltzmann distribution. To illustrate our results, we trained replicas of RBMs with different hidden layer sizes using the MNIST dataset. Our findings not only bridge the gap between disparate generative frameworks but also shed light on the fundamental processes underpinning learning in deep generative models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Feature learning and generalization in deep networks with orthogonal weights

Fully-connected deep neural networks with weights initialized from independent Gaussian distributions can be tuned to criticality, which prevents the exponential growth or decay of signals propagating through the network. However, such networks still exhibit fluctuations that grow linearly with the depth of the network, which may impair the training of networks with width comparable to depth. We show analytically that rectangular networks with tanh activations and weights initialized from the ensemble of orthogonal matrices have corresponding preactivation fluctuations which are independent of depth, to leading order in inverse width. Moreover, we demonstrate numerically that, at initialization, all correlators involving the neural tangent kernel (NTK) and its descendants at leading order in inverse width—which govern the evolution of observables during training—saturate at a depth of ~20, rather than growing without bound as in the case of Gaussian initializations. We speculate that this structure preserves finite-width feature learning while reducing overall noise, thus improving both generalization and training speed in deep networks with depth comparable to width. We provide some experimental justification by relating empirical measurements of the NTK to the superior performance of deep non-linear orthogonal networks trained under full-batch gradient descent on the MNIST and CIFAR-10 classification tasks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Statistical Analysis of Intertube Tunneling Contacts in the Macroscopic Electrical Conductivity of Carbon Nanotube Fibers

Here, this study investigates the influence of tunneling contact resistances between carbon nanotubes (CNTs) on electron transport and electrical conductivity of macroscopic carbon nanofibers (CNFs), which profoundly impacts the performance of CNT thin film electronics, CNF electron emitters and cathodes, and energy conversion and storage devices. Utilizing a self-consistent electrical contact model coupling a transmission line model with tunneling current, we calculate the contact resistances of a plethora of CNT-CNT contacts within a CNF fiber, which consists of aligned, densely packed CNTs. A statistical analysis is conducted, using Gaussian distributions to account for variations in contact lengths, tunneling gap distances, and single CNT aspect ratios, to calculate the CNT-CNT contact resistance and the overall resistance of CNT fiber. By scaling our model to a macroscopic level, our results are in good agreement with experimental measurements. Our calculation suggests that while increasing the contact overlap length diminishes individual CNT-CNT contact resistance, it could paradoxically increase macroscopic CNT fiber resistance for a given constant CNF mass density, which is due to that fact that a larger overlap length allows more CNTs to pack along an electrical conduction path per unit length, leading to more tunneling contact junctions connected in series and thus less number of parallel conduction paths within the fiber cross section. Increasing tunneling gap distance increases both individual contact and overall fiber resistance. This research provides a simple design tool for tailoring CNT fiber electrical properties to promote real-world applications using CNTs or similar low-dimensional materials.

36 MATERIALS SCIENCE↗