Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ordered statistics”

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 199 records · Page 11

How Well Do We Know the Neutron-Matter Equation of State at the Densities Inside Neutron Stars? A Bayesian Approach with Correlated Uncertainties

Here, we introduce a new framework for quantifying correlated uncertainties of the infinite-matter equation of state derived from chiral effective field theory (𝜒⁢EFT ). Bayesian machine learning via Gaussian processes with physics-based hyperparameters allows us to efficiently quantify and propagate theoretical uncertainties of the equation of state, such as 𝜒⁢EFT truncation errors, to derived quantities. We apply this framework to state-of-the-art many-body perturbation theory calculations with nucleon-nucleon and three-nucleon interactions up to fourth order in the 𝜒⁢EFT expansion. This produces the first statistically robust uncertainty estimates for key quantities of neutron stars. We give results up to twice nuclear saturation density for the energy per particle, pressure, and speed of sound of neutron matter, as well as for the nuclear symmetry energy and its derivative. At nuclear saturation density, the predicted symmetry energy and its slope are consistent with experimental constraints.

79 ASTRONOMY AND ASTROPHYSICS↗

Misanthropic entropy and renormalization as a communication channel

A central physical question is the extent to which infrared (IR) observations are sufficient to reconstruct a candidate ultraviolet (UV) completion. We recast this question as a problem of communication, with messages encoded in field configurations of the UV being transmitted to the IR degrees of freedom via a noisy channel specified by renormalization group (RG) flow, with noise generated by coarse graining/decimation. We present an explicit formulation of these considerations in terms of lattice field theory, where we show that the “misanthropic entropy” — the mutual information obtained from decimating neighbors — encodes the extent to which information is lost in marginalizing over/tracing out UV degrees of freedom. Our considerations apply both to statistical field theories as well as density matrix renormalization of quantum systems, where in the quantum case, the statistical field theory analysis amounts to a leading-order approximation. As a concrete example, we focus on the case of the 2D Ising model, where we show that the misanthropic entropy detects the onset of the phase transition at the Ising model critical point.

Physics↗

First Measurement of 87 Rb( α , xn ) Cross Sections at Weak r -process Energies in Supernova ν -driven Ejecta to Investigate Elemental Abundances in Low-metallicity Stars

Observed abundances of Z ∼ 40 elements in metal-poor stars vary from star to star, indicating that the rapid and slow neutron capture processes may not contribute alone to the synthesis of elements beyond iron. The weak r-process was proposed to produce Z ∼ 40 elements in a subset of old stars. Thought to occur in the ν-driven ejecta of a core-collapse supernova, ( α, xn ) reactions would drive the nuclear flow toward heavier masses at T = 2−5 GK. However, current comparisons between modeled and observed yields do not bring satisfactory insights into the stellar environment, mainly due to the uncertainties of the nuclear physics inputs where the dispersion in a given reaction rate often exceeds 1 order of magnitude. Involved rates are calculated with the statistical model where the choice of an α -optical-model potential ( α OMP) leads to such a poor precision. The first experiment on 87 Rb( α, xn ) reactions at weak r -process energies is reported here. Total inclusive cross sections were assessed at E c.m. = 8.1−13 MeV (3.7−7.6 GK) with the active target MUlti-Sampling Ionization Chamber. With an N = 50 seed nucleus, the measured values agree with statistical model estimates using the α OMP Atomki-V2. A reevaluated reaction rate was incorporated into new nucleosynthesis calculations, focusing on ν-driven ejecta conditions known to be sensitive to this specific rate. These conditions were found to fail to reproduce the lighter heavy element abundances in metal-poor stars.

79 ASTRONOMY AND ASTROPHYSICS↗

Bayesian Inference of Hybrid Star Properties from Future High-precision Measurements of Their Radii

Future high-precision X-ray and gravitational-wave observations of neutron stars (NSs) are expected to constrain NS radii with uncertainties as small as σ ≃ 0.1 km. Such unprecedented precision offers a unique opportunity to extract new information about the nature and equation of state (EOS) of supradense matter in NS cores. Using mock radius data with uncertainties ranging from σ = 1.0 to 0.1 km, together with a flexible meta-model NS EOS that allows for a first-order hadron–quark phase transition, we perform a Bayesian statistical analysis to assess the impact of radius measurements on EOS constraints. We find that high-precision radius measurements, particularly for massive NSs, significantly tighten constraints on the hadron–quark transition density ρ t , the quark matter mass fraction in NS cores, and several parameters characterizing the EOS of supranuclear hadronic matter, although the degree of improvement depends on the assumed prior range of ρ t . In contrast, even with the highest precision considered, NS radii—including those of massive stars—remain largely insensitive to the stiffness of quark matter, independent of the measurement accuracy or the prior range adopted for ρ t .

Neutron stars↗

The QCD phase diagram and statistics friendly distributions

The preliminary STAR data for proton cumulants for central collisions at \sqrt{s}=7.7 GeV are consistent with a two-component proton multiplicity distribution. We show that this two-component distribution is statistics friendly in that factorial cumulants of surprisingly high orders may be extracted with a relatively small number of events. As a consequence the two-component model can be tested and verified right now with the presently available STAR data from the first phase of the RHIC beam energy scan.

Koch, Volker↗

Towards a Comparative Assessment of Data-Driven Process Models in Health Information Technology

Process mining for conformance analysis focuses on comparing a reference process model against a data-driven process model that is generated via log files from information technology systems. While this approach is helpful when there is an existing process model in an organization, it leaves the question of what to do in the absence of a complete reference process model unanswered. In this paper, we present a comparative assessment approach that combines process mining, process mapping for dimensionality reduction, and statistical analysis. Our goal is to find similarities and dissimilarities in data-driven process models among U.S. Veterans Health Administration (VHA) facilities to assess process conformance among different healthcare facilities, which can help assess the standardization of care. We illustrate our approach by applying it to two clinical radiology order process models generated by two similar facilities. Our results demonstrate statistical similarities in the standardization of care among those two facilities.

Klasky, Hilda↗

Superconductivity in a system of interacting spinful semions

Noninteracting particles obeying certain fractional statistics have been predicted to exhibit superconductivity. We discuss the issue in an attractively interacting system of spinful semions on a lattice by numerically investigating the presence of off-diagonal long-range order at zero temperature. For this purpose we construct a Hubbard model wherein two semions with opposite spin can virtually coincide while maintaining consistency with the fractional braiding statistics. Clear off-diagonal long-range order is seen in the strong coupling limit, consistent with the expectation that a pair of semions obeys Bose statistics. We find that the semion system behaves similarly to a system of fermions with the same attractive Hubbard U interaction for a wide range of U, suggesting that semions also undergo a BCS to BEC crossover as a function of U.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

The impact of detection rate changes and correlations on random-coincidence background measurements

Coincidence detection of multiple particles emitted during an experiment can yield a new depth of understanding of the underlying process under study. However, the probability of detecting particles that are generated from the same physical event within a given coincidence time window is generally much lower than that of detecting particles that appear in the same coincidence time window, but were not created from the same physical event, and are therefore detected randomly in coincidence with each other. Thus, accurate and precise methods of measuring this random-coincidence background are essential for a wide variety of fields of science. A method to determine this background directly using the data themselves without any additional experimental run time or fake signals introduced in the data was recently established (O’Donnell, 2016). This method yields a statistical uncertainty on the random-coincidence background that is orders of magnitude smaller than that of the true coincidence data, though the potential for systematic errors of backgrounds from this method was never explored. In this work, we discuss common varieties of correlated and uncorrelated changes in the detection rates of each particle detected in an experiment. Here we demonstrate here that a correlation between particle detection rates from, for example, an incident particle beam that initiates a physical process of interest, creates systematic errors in the random-coincidence background measurement. We also discuss the impact of a variety of other realistic scenarios for rate changes in experiments. Lastly, a method is introduced to correct for errors in the random-coincidence background from any source, yielding an optimization between statistical precision and eliminating potential lingering systematic errors.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Stochastic Vector Techniques in Ground-State Electronic Structure

Herein we review a suite of stochastic vector computational approaches for studying the electronic structure of extended condensed matter systems. These techniques help reduce algorithmic complexity, facilitate efficient parallelization, simplify computational tasks, accelerate calculations, and diminish memory requirements. While their scope is vast, we limit our study to ground-state and finite temperature density functional theory (DFT) and second-order many-body perturbation theory. More advanced topics, such as quasiparticle (charge) and optical (neutral) excitations and higher-order processes, are covered elsewhere. We start by explaining how to use stochastic vectors in computations, characterizing the associated statistical errors. Next, we show how to estimate the electron density in DFT and discuss effective techniques to reduce statistical errors. Finally, we review the use of stochastic vectors for calculating correlation energies within the second-order Møller-Plesset perturbation theory and its finite temperature variational form. Example calculation results are presented and used to demonstrate the efficacy of the methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Band-limited photodetection of temporal coherence

The quantum theory of optical coherence plays a ubiquitous role in identifying optical emitters. An unequivocal identification, however, presumes that the photon number statistics is resolved from timing uncertainties. We demonstrate from first principle that the observed nth-order temporal coherence is a n-fold convolution of the instrument responses and the expected coherence. The consequence is detrimental in which the photon number statistics is masked from the unresolved coherence signatures. The experimental investigations are thus far consistent with the theory developed. We envision the present theory will mitigate the false identification of optical emitters and enlarge the coherence deconvolution to an arbitrary order.

42 ENGINEERING↗

Iterative Reconstruction for Multimodal Neutron Tomography

Here, we describe a unified framework for model-based iterative 3-D reconstruction of multimodal neutron transmission, hydrogen-scatter, and induced-fission images from low resolution data recorded using 14.1-MeV neutrons and the associated-particle imaging (API) technique. The framework, which was developed to facilitate use in challenging field-deployment scenarios, is centered around physics-based system models and a total variation (TV) constrained implementation of the simultaneous iterative reconstruction technique (SIRT). Modified to solve a statistically weighted least squares (WLS) problem, the SIRT algorithm is accelerated using ordered subsets and Nesterov’s momentum for which we derive a near-optimal value of the governing Lipschitz constant. The approach enables the reconstruction of images that are high resolution compared to the acquired data and is robust to both limited statistics and a limited number of projection angles. Moreover, the framework is fast enough to be practical. Example images are provided that demonstrate both the ability to perform fast-neutron imaging of high-atomic-number materials with low radiation dose and the benefit of multimodal neutron imaging to identify key materials.

Hydrogen scatter↗

A Preliminary Investigation for the Use of Digital Gamma-Gamma Coincidence Spectrometry to Determine 239 Pu

Gamma-ray spectrometry is one of the most effective ways to determine the activity of 239 Pu, depending on activity levels. However, often high backgrounds in complex spectra with low amounts of 239 Pu can increase detection limits. The effectiveness of the use of gamma-gamma spectrometry in the characterization of 239 Pu was studied for the first time. Using the XIA Pixie-16 digital pulse processor, the gamma-gamma system was studied to establish unshielded background radiation, source to detector distances, Poisson statistics and gating of several gamma-rays from 239 Pu. The results clearly show that the distance between gamma ray spectrometers has a much bigger effort in counting statistics for the gamma-gamma coincidence detectors than a single detector. In addition, measurements demonstrated the background radiation to be virtually negligible, meaning gammagamma spectroscopy can still be very effective without requiring the usual lead shielding. Other measurements were taken to assure Poisson statistics were attained in the digital system. Preliminary measurements revealed 3 orders of magnitude background reduction for the measurement of 239 Pu using gated gamma-rays.

61 RADIATION PROTECTION AND DOSIMETRY↗

The effect of local chemical ordering on Frank-Read source activation in a refractory multi-principal element alloy

In this work, we investigate the operation of Frank-Read (FR) sources in a refractory multi-principal element alloy (MPEA). Simulations of discrete dislocation motion in MPEAs is enabled by the development of a phase field dislocation dynamics model that treats the atomic-scale fluctuations in lattice energies across the glide plane, present due to local ordering in the chemical composition within the nominally random MPEA atomic structure. We consider, through simulation, a range of length scales over which ordering occurs, varying from short-range lengths, a few times dislocation core width, to long-range lengths, an order of magnitude longer than the core. Characteristic of this body-centered cubic MPEA, the simulations also include screw/edge character dependence in glide resistance, as informed by atomic scale simulation. The critical stresses to activate the source for the same source size are found to be statistically distributed, as a direct consequence of the underlying variation in lattice energy. Furthermore, analysis of critical state for activating edge and screw FR sources in the MPEA reveals that FR source operation occurs via a two-step mechanism, involving athermal kink-pair formation, unlike the conventional FR source operation in a material with no composition fluctuations. This mechanism lowers the average critical stress required to activate the FR source and causes the statistical dispersion in critical stress to depend on the range of composition ordering. More importantly, it leads to a more severe dependence of source strength on FR source length than predicted by line tension alone.

36 MATERIALS SCIENCE↗

Theory of topological defects and textures in two-dimensional quantum orders with spontaneous symmetry breaking

In this article, we consider two-dimensional (2d) quantum many-body systems with long-range orders, where the only gapless excitations in the spectrum are Goldstone modes of spontaneously broken continuous symmetries. To understand the interplay between classical long-range order of local order parameters and quantum order of long-range entanglement in the ground states, we study the topological point defects and textures of order parameters in such systems. We show that the universal properties of point defects and textures are determined by the remnant symmetry enriched topological order in the symmetry-breaking ground states with a nonfluctuating order parameter, and provide a classification for their properties based on the inflation-restriction exact sequence. We highlight a few phenomena revealed by our theory framework. First, in the absence of intrinsic topological orders, we show a connection between the symmetry properties of point defects and textures to deconfined quantum criticality. Second, when the symmetry-breaking ground state has intrinsic topological orders, we show that the point defects can permute different anyons when braided around. They can also obey projective fusion rules in the sense that multiple vortices can fuse into an Abelian anyon, a phenomenon for which we coin “defect fractionalization.” Finally, we provide a formula to compute the fractional statistics and fractional quantum numbers carried by textures (skyrmions) in Abelian topological orders.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Bose–Einstein condensation of a two-magnon bound state in a spin-1 triangular lattice

In ordered magnets, the elementary excitations are spin waves (magnons), which obey Bose–Einstein statistics. Similarly to Cooper pairs in superconductors, magnons can be paired into bound states under attractive interactions. The Zeeman coupling to a magnetic field is able to tune the particle density through a quantum critical point, beyond which a ‘hidden order’ is predicted to exist. Here, in this work, we report direct observation of the Bose–Einstein condensation of the two-magnon bound state in Na 2 BaNi(PO 4 ) 2 . Comprehensive thermodynamic measurements confirmed the two-dimensional Bose–Einstein condensation quantum critical point at the saturation field. Inelastic neutron scattering experiments were performed to establish the microscopic model. An exact solution revealed stable two-magnon bound states that were further confirmed by electron spin resonance and nuclear magnetic resonance experiments, demonstrating that the quantum critical point is due to the pair condensation, and the phase below the saturation field is likely the long-sought-after spin nematic phase.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Statistical uncertainty of fission matrix eigenvalues using perturbation theory

Eigenvalue search of high dominance ratio systems may be slow to converge. The fission matrix element is defined by its element (FM){sub ij}, which are the probability for a neutron born in cell i to create a fission in cell j for a spatial mesh of n{sub i}*n{sub j} elements. Fission matrices are used in Monte Carlo criticality simulations to enhance computing speed, but also to find higher order eigenvalues. However, few studies have been made on the link between statistical uncertainties of fission matrix elements and eigenvalues uncertainties. Thus, dominance ratio statistical uncertainties remain unknown. This paper uses a new generalized perturbation theory (GPT) method to estimate sensitivities of eigenvalues to fission matrix elements and then to calculate dominance ratio uncertainties.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Accelerated Depth Computation for Surface Boxplots with Deep Learning

Functional depth is a well-known technique used to derive descriptive statistics (e.g., median, quartiles, and outliers) for 1D data. Surface boxplots extend this concept to ensembles of images, helping scientists and users identify representative and outlier images. However, the computational time for surface boxplots increases cubically with the number of ensemble members, making it impractical for integration into visualization tools. In this paper, we propose a deep-learning solution for efficient depth prediction and computation of surface boxplots for time-varying ensemble data. Our deep learning framework accurately predicts member depths in a surface boxplot, achieving average speedups of 6X on a CPU and 15X on a GPU for the 2D Red Sea dataset with 50 ensemble members compared to the traditional depth computation algorithm. Our approach achieves at least a 99% level of rank preservation, with order flipping occurring only at pairs with extremely similar depth values that pose no statistical differences. This local flipping does not significantly impact the overall depth order of the ensemble members.

Han, Mengjiao↗

A compute-bound formulation of Galerkin model reduction for linear time-invariant dynamical systems

This work aims to advance computational methods for projection-based reduced-order models (ROMs) of linear time-invariant (LTI) dynamical systems. For such systems, current practice relies on ROM formulations expressing the state as a rank-1 tensor (i.e., a vector), leading to computational kernels that are memory bandwidth bound and, therefore, ill-suited for scalable performance on modern architectures. This weakness can be particularly limiting when tackling many-query studies, where one needs to run a large number of simulations. This work introduces a reformulation, called rank-2 Galerkin, of the Galerkin ROM for LTI dynamical systems which converts the nature of the ROM problem from memory bandwidth to compute bound. We present the details of the formulation and its implementation, and demonstrate its utility through numerical experiments using, as a test case, the simulation of elastic seismic shear waves in an axisymmetric domain. We quantify and analyze performance and scaling results for varying numbers of threads and problem sizes. In conclusion, we present an end-to-end demonstration of using the rank-2 Galerkin ROM for a Monte Carlo sampling study. We show that the rank-2 Galerkin ROM is one order of magnitude more efficient than the rank-1 Galerkin ROM (the current practice) and about 970 times more efficient than the full-order model, while maintaining accuracy in both the mean and statistics of the field.

97 MATHEMATICS AND COMPUTING↗