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 55 records · Page 3

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↗

Two Mesoporous Domains Are Better Than One for Catalytic Deconstruction of Polyolefins

Catalytic hydrogenolysis of polyolefins into valuable liquid, oil, or wax-like hydrocarbon chains for second-life applications is typically accompanied by the hydrogen-wasting co-formation of low value volatiles, notably methane, that increase greenhouse gas emissions. Catalytic sites confined at the bottom of mesoporous wells, under conditions in which the pore exerts the greatest influence over the mechanism, are capable of producing less gases than unconfined sites. A new architecture was designed to emphasize this pore effect, with the active platinum nanoparticles embedded between linear, hexagonal mesoporous silica and gyroidal cubic MCM-48 silica (mSiO 2 /Pt/MCM-48). This catalyst deconstructs polyolefins selectively into ~C 20 –C 40 paraffins and cleaves C–C bonds at a rate (TOF = 4.2 ± 0.3 s –1 ) exceeding that of materials lacking these combined features while generating negligible volatile side products including methane. The time-independent product distribution is consistent with a processive mechanism for polymer deconstruction. In contrast to time- and polymer length-dependent products obtained from non-porous catalysts, mSiO 2 /Pt/MCM-48 yields a C 28 -centered Gaussian distribution of waxy hydrocarbons from polyolefins of varying molecular weight, composition, and physical properties, including low-density polyethylene, isotactic polypropylene, ultrahigh-molecular-weight polyethylene, and mixtures of multiple, post-industrial polyolefins. Here, coarse-grained simulation reveals that the porous-core architecture enables the paraffins to diffuse away from the active platinum site, preventing secondary reactions that produce gases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Shear-induced lift force on spheres in a viscous linear shear flow at finite volume fractions

Several studies have shown a significant increase in drag on a distribution of solid spherical particles within a fluid with increasing particle volume fraction. As a result, many empirical drag laws accounting for the dependence on the Reynolds number and volume fraction can be found in the literature. This study investigates the possibility of a similar effect of the particle volume fraction on the mean hydrodynamic lift force on randomly distributed spherical particles in a linear shear flow. Particle-resolved direct numerical simulations are performed to evaluate the mean lift force, and the results are compared with the case of an isolated particle in a linear shear flow for the same Reynolds number and shear rate. The mean lift force acting on the particles appears to remain nearly the same as that on an isolated particle. However, due to the influence of neighboring particles, there is a substantial force variation in transverse directions on each individual particle, whose magnitude is comparable to the mean drag force. The distribution of drag force in a linear shear flow is shown to be nearly the same as in a uniform flow at the same volume fraction and Reynolds number. A simple stochastic model based on a Gaussian distribution is presented for the lift force variation, and its performance is compared to the prediction of the deterministic pairwise interaction extended point-particle model.

42 ENGINEERING↗

Normalizing flows for domain adaptation when identifying Λ hyperon events

Here this study focuses on the application of a normalizing flow as a method of domain adaptation when classifying physics data. Normalizing flows offer a way to transform data points between two different distributions. The present study investigates a novel method of transforming latent representations of physics data to a normal distribution and then to a physics distribution again. The final distribution models a simulated distribution. After being transformed, the data can be classified by a neural network trained on labeled simulation data. The present study succeeds in training two normalizing flows that can transform between data (or simulation) and a Gaussian distribution.

47 OTHER INSTRUMENTATION↗

Forecasting cosmological constraints from the weak lensing magnification of type Ia supernovae measured by the Nancy Grace Roman Space Telescope

We report that the weak lensing magnification of type Ia supernovae (SNe Ia) is sensitive to the clustering of matter and provides an independent cosmological probe complementary to SN Ia distance measurements. The Nancy Grace Roman Space Telescope is uniquely sensitive to this measurement as it can discover high redshift SNe Ia and measure them with high precision. We present a methodology for reconstructing the probability distribution of the weak lensing magnification μ of SNe Ia, p(μ), from observational data, and using it to constrain cosmological parameters. We find that the reconstructed p(μ) can be fitted accurately by a stretched Gaussian distribution and used to measure the variance of μ, ξ μ , which can be compared to theoretical predictions in a likelihood analysis. Applying our methodology to a set of realistically simulated SNe Ia expected from the Roman Space Telescope, we find that using the weak lensing magnification of the SNe Ia constrains a combination of matter density Ω m and matter clustering amplitude σ 8 . SN Ia distances alone lead to a better than 1% measurement of Ωm. The combination of SN Ia weak lensing magnification and distance measurements result in a ~10 % measurement on σ 8 . The SNe Ia from Roman will be powerful in constraining the cosmological model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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↗

Disorder-induced local strain distribution in Y-substituted TmVO 4

We report an investigation of the effect of substitution of Y for Tm in Tm 1-x ⁢Y x VO 4 via low-temperature heat capacity measurements, with the yttrium content x varying from 0 to 0.997. Because the Tm ions support a local quadrupolar (nematic) moment, they act as reporters of the local strain state in the material, with the splitting of the ion's non-Kramers crystal field ground state proportional to the quadrature sum of the in-plane tetragonal symmetry-breaking transverse and longitudinal strains experienced by each ion individually. Analysis of the heat capacity, therefore, provides detailed insights into the distribution of local strains that arise as a consequence of the chemical substitution. These local strains suppress long-range quadrupole order for x > 0.22, and result in a broad Schottky-like feature for higher concentrations. Heat capacity data are compared to expectations for a distribution of uncorrelated (random) strains. For dilute Tm concentrations, the heat capacity cannot be accounted for by randomly distributed strains, demonstrating the presence of significant strain correlations between sites. For intermediate Tm concentrations, these correlations must still exist, but the data cannot be distinguished from that which would be obtained from a two-dimensional Gaussian distribution. The crossover between these limits is discussed in terms of the interplay of key lengthscales in the substituted material. Furthermore, the central result of this work, namely that local strains arising from chemical substitution are not uncorrelated, has implications for the range of validity of theoretical models based on random effective fields that are used to describe such chemically substituted materials, particularly when electronic nematic correlations are present.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

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↗

How “hot” are hotspots: Statistically localizing the high-activity areas on soil and rhizosphere images

The topic of microbial hotspots in soil requires not only visualizing their spatial distribution and biochemical analyses, but also statistical approaches to identify these hotspots and separate them from the surrounding activities (background). We hypothesized that each hotspot type (e.g. enzyme activities in the rhizosphere, root exudation, localization of herbicide accumulation) is a result of local process driven by biotic and/or abiotic factors, and the process rates in the hotspots are much faster than those in the soil background. We further hypothesized that the background and hotspot activities in soil belong to different statistical distributions. Consequently, hotspot determination should be based on statistical separation of activities significantly higher than the background. We analyzed for the statistical distributions of grey values on three groups of published images: 1) 14 C images of carbon input by roots into the rhizosphere, 2) 14 C glyphosate accumulation in the plant, and 3) zymogram of leucine aminopeptidase activity in rooted soil. The two Gaussian distributions were fit (the first representing the background, the second the hotspots) to the distribution of grey values in the images, the parameters (means and standard deviations, SD) of the fitted distributions were calculated, and the background was removed. Thus, we identified hotspots as areas outside of the Mean+2SD image intensity (corresponding to the upper ~ 2.5% of activity, being over 97.5% of background values) and finally, visualized images of solely hotspot locations. Finally, these results were compared with previously used decisions on hotspot intensity thresholding (i.e. Top-25% and 17 standard thresholding approaches in ImageJ) and discussed the advantages of the Mean+2SD as well as Mean+3SD approaches. These advantages include: i) simple unification of the thresholding approach for several imaging methods with various principles of activity distribution, ii) identification of hotspots with various activity levels, iii) analysis of “time-specific” hotspots in temporal sequences of images. Compared with 17 standard thresholding methods, we concluded that objectively elucidating and separating the hotspots should be based on statistical distribution analysis, e.g. using the Mean+2SD or Mean+3SD approaches. Furthermore, this simple Mean+2SD approach delivered suitable results for three groups of images and so, helps to understand the processes responsible for the highest activities and elucidate hotspots.

59 BASIC BIOLOGICAL SCIENCES↗

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↗

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↗

Pseudorapidity distributions of charged particles in pp($\bar p$), p(d)A and AA collisions using Tsallis thermodynamics

The pseudorapidity distributions of charged particles measured in p + p($\overline{\text{p}}$) collisions for energies ranging from $\sqrt{{s}_{\mathrm{N}\mathrm{N}}}=23.6$ GeV to 13 TeV and A + A collisions at RHIC and LHC are investigated in the fireball model with Tsallis thermodynamics. We assume that the rapidity axis is populated with fireballs following q-Gaussian distribution and the charged particles follow the Tsallis distribution in the fireball. We also extend the fireball model to asymmetric collision systems, i.e. d + Au collisions at $\sqrt{{s}_{\mathrm{N}\mathrm{N}}}=200$ GeV and p + Pb collisions at $\sqrt{{s}_{\mathrm{N}\mathrm{N}}}=5.02$ TeV, by taking into account the asymmetric geometry configuration. The model can fit the experimental data well for all the collision systems and centralities investigated. The collision energy and centrality dependence of the model parameters for the symmetric (asymmetric) collision system, i.e. the central position y 0 (y 0a , y 0A ) and its width σ (σ a , σ A ) of the fireball distribution, are also investigated and discussed. Furthermore, our results suggest that the fireball model with Tsallis thermodynamics can be used as a universal framework for the pseudorapidity distributions of charged particles in high energy collisions at RHIC and LHC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Evaluating Space Object Conjunction Probabilities Using Characteristic Function Inversion

This report discusses an approach to computing the probability of a conjunction between two space objects in the short-term encounter scenario. A conjunction is defined here as an event where the miss distance between the objects is less than some specified value. The scenario assumptions are that the motion of the objects is linear, their positions are Gaussian distributed, and their velocities are known and constant. Under these assumptions, the squared-miss distance is shown to have the generalized chi-square distribution. An established statistical technique called characteristic function inversion is employed to evaluate the distribution and obtain conjunction probabilities. The method is closely related to a recent approach based on moment generating function inversion, and a qualitative comparison of the approaches is provided. Last, the method is tested on several benchmark test cases where it agrees with numerical integration on the cases with conjunction probabilities above 10 –12 . However, the exact probability in these cases is usually not needed and this probability can be bounded above using an independent Gaussian approximation. Overall, the report shows how to compute conjunction probabilities using a standard statistical method, though numerical integration seems to perform equally well.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Novel Data-based Stochastic Distribution Control for Non-Gaussian Stochastic Systems

In this note, we present a novel data-based approach to investigate the non-Gaussian stochastic distribution control problem. As the motivation of this note, the existing methods have been summarised regarding to the drawbacks, for example, neural network weights training for unknown stochastic distribution and so on. To overcome these disadvantages, a new transformation for dynamic probability density function is given by kernel density estimation using interpolation. Based upon this transformation, a representative model has been developed while the stochastic distribution control problem has been transformed into an optimisation problem. Then, data-based direct optimisation and identification-based indirect optimisation have been proposed. In addition, the convergences of the presented algorithms are analysed and the effectiveness of these algorithms has been evaluated by numerical examples. In summary, the contributions of this note are as follows: 1) a new data-based probability density function transformation is given; 2) the optimisation algorithms are given based on the presented model; and 3) a new research framework is demonstrated as the potential extensions to the existing stochastic distribution control.

42 ENGINEERING↗

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),↗

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↗

On the Stochastic Stability of Deep Markov Models

Deep Markov models (DMM) are generative models which are scalable and expressive generalization of Markov models for representation, learning, and inference problems. DMMs using deep neural networks to parametrize the transition of Markov probability distributions have recently been shown to provide more expressiveness in modeling sequential data and dynamical system responses. However, the fundamental stochastic stability guarantees of such models have not been thoroughly investigated. In this paper, we present a rigorous analytical method to prove the necessary and sufficient conditions of DMM's stochastic stability. This task is achieved by spectral analysis of the efficiently computed Jacobians of probabilistic maps modeled by deep neural networks. We make theoretical connections between the eigenvalues of neural network's weights and the different activation function types used on the stability and overall dynamic behavior of DMMs with Gaussian distributions. We empirically substantiate our theoretical results on stochastic stability and eigenvalue spectra via several numerical experiments. Formal stability guarantees of DMMs can substantially improve their robustness and trustworthiness, necessary for reliable use in safety-critical real-world applications.

Drgona, Jan↗

Assessment of Local Observation of Atomic Ordering in Alloys via the Radial Distribution Function: A Computational and Experimental Approach

As a powerful analytical technique, atom probe tomography (APT) has the capacity to acquire the spatial distribution of millions of atoms from a complex sample. However, extracting information at the Ångstrom-scale on atomic ordering remains a challenge due to the limits of the APT experiment and data analysis algorithms. The development of new computational tools enable visualization of the data and aid understanding of the physical phenomena such as disorder of complex crystalline structures. Here, we report progress towards this goal using two steps. We describe a computational approach to evaluate atomic ordering in the crystal structure by generating radial distribution functions (RDF). Atomic ordering is rendered as the Fractional Cumulative Radial Distribution Function (FCRDF) which allows for greater visibility of local compositions at short range in the structure. Further, we accommodate in the analysis additional parameters such as uncertainty in the atomic coordinates and the atomic abundance to ascertain short-range ordering in APT data sets. We applied the FCRDF analysis to synthetic and experimental APT data sets for Ni 3 Al. The ability to observe a signal of atomic ordering consistent with the known L1 2 crystal structure is heavily dependent on spatial uncertainty, irrespective of abundance. Detection of atomic ordering is subject to an upper limit of spatial uncertainty of atoms described with Gaussian distributions with a standard deviation of 1.3 Å. The FCRDF analysis was also applied to the APT data set for a six-component alloy, Al 1.3 CoCrCuFeNi. In this case, we are currently able to visualize elemental segregation at the nanoscale, though unambiguous identification of atomic ordering at the Ångstrom (nearest-neighbor) scale remains a goal.

36 MATERIALS SCIENCE↗