Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 1,585 records · Page 88

Statistical properties of filaments in the cosmic web

ABSTRACT In the context of the cosmological and constrained Exploring the Local Universe with the reConstructed Initial Density field (ELUCID) simulation, this study explores the statistical characteristics of filaments within the cosmic web, focussing on aspects such as the distribution of filament lengths and their radial density profiles. Using the classification of the cosmic web environment through the Hessian matrix of the density field, our primary focus is on how cosmic structures react to the two variables $R_{\rm s}$ and $\lambda _{\rm th}$. The findings show that the volume fractions of knots, filaments, sheets, and voids are highly influenced by the threshold parameter $\lambda _{\rm th}$, with only a slight influence from the smoothing length $R_{\rm s}$. The central axis of the cylindrical filament is pinpointed using the medial-axis thinning algorithm of the COsmic Web Skeleton (COWS) method. It is observed that median filament lengths tend to increase as the smoothing lengths increase. Analysis of filament length functions at different values of $R_{\rm s}$ indicates a reduction in shorter filaments and an increase in longer filaments as $R_{\rm s}$ increases, peaking around $2.5R_{\rm s}$. The study also shows that the radial density profiles of filaments are markedly affected by the parameters $R_{\rm s}$ and $\lambda _{\rm th}$, showing a valley at approximately $2R_{\rm s}$, with increases in the threshold leading to higher amplitudes of the density profile. Moreover, shorter filaments tend to have denser profiles than their longer counterparts.

Zhang, Youcai (ORCID:0000000319674091)↗

Intrinsic alignments of galaxies in the MassiveBlack-II simulation: Analysis of two-point statistics

The intrinsic alignment of galaxies with the large-scale density field in an important astrophysical contaminant in upcoming weak lensing surveys. We present detailed measurements of the galaxy intrinsic alignments and associated ellipticity-direction (ED) and projected shape (w g ₊) correlation functions for galaxies in the cosmological hydrodynamic MassiveBlack-II (MB-II) simulation. We carefully assess the effects on galaxy shapes, misalignment of the stellar component with the dark matter shape and two-point statistics of iterative weighted (by mass and luminosity) definitions of the (reduced and unreduced) inertia tensor. We find that iterative procedures must be adopted for a reliable measurement of the reduced tensor but that luminosity versus mass weighting has only negligible effects. Both ED and w g ₊ correlations increase in amplitude with subhalo mass (in the range of 10¹⁰ – 6.0 X 10¹⁴h⁻¹ M ⊙ ), with a weak redshift dependence (from z = 1 to z = 0.06) at fixed mass. At z ~ 0.3, we predict a w g ₊ that is in reasonable agreement with SDSS LRG measurements and that decreases in amplitude by a factor of ~ 5–18 for galaxies in the LSST survey. We also compared the intrinsic alignment of centrals and satellites, with clear detection of satellite radial alignments within the host halos. Finally, we show that w g ₊ (using subhalos as tracers of density and w δ (using dark matter density) predictions from the simulations agree with that of non-linear alignment models (NLA) at scales where the 2-halo term dominates in the correlations (and tabulate associated NLA fitting parameters). The 1-halo term induces a scale dependent bias at small scales which is not modeled in the NLA model.

79 ASTRONOMY AND ASTROPHYSICS↗

Reducing statistical uncertainties in simulated organ doses of phantoms immersed in water

In this study, methods are addressed to reduce the computational time to compute organ-dose rate coefficients using Monte Carlo techniques. Several variance reduction techniques are compared including the reciprocity method, importance sampling, weight windows and the use of the ADVANTG software package. For low-energy photons, the runtime was reduced by a factor of 10 5 when using the reciprocity method for kerma computation for immersion of a phantom in contaminated water. This is particularly significant since impractically long simulation times are required to achieve reasonable statistical uncertainties in organ dose for low-energy photons in this source medium and geometry. Although the MCNP Monte Carlo code is used in this paper, the reciprocity technique can be used equally well with other Monte Carlo codes.

61 RADIATION PROTECTION AND DOSIMETRY↗

Stokes-dependent droplet collection efficiency on a NACA 0012 airfoil from droplet-informed simulations with statistical overloading

Accurate modelling of ice accretion on aircraft wings requires analysing droplet impingement on the surface to optimize the design of ice-protection systems. We perform Euler–Lagrange simulations of a droplet-laden flow impinging on a NACA 0012 airfoil. Our study includes water droplets with eight discrete sizes ranging from 1 to 160 microns. We vary the free-stream velocity of the incoming airflow in the range 60 ≤ U ≤ 240 m s −1 and the chord length of the airfoil in the range 0.5 ≤ c ≤ 2 m. Due to the dilute nature of supercooled clouds, one-way coupling is used in the simulations. The effects of droplet breakup and collision are also neglected. To reduce the computational cost, we employ statistical overloading of droplets, allowing us to simulate millions of impinging droplets in a time span on the order of milliseconds. Our results show that the droplet collection efficiency, which measures the likelihood of droplet impingement on the airfoil surface, increases with droplet size and free-stream velocity but decreases with airfoil size. We demonstrate that collection efficiency, impingement velocity and impingement angle are primarily dictated by a single non-dimensional parameter, the droplet Stokes number. We also identify a critical stagnation-streamline Stokes number below which impingements do not occur and use it to estimate the minimum droplet size for impingement. In addition, we observe droplet behaviour to become Stokes number independent at large values of the Stokes number. This article is part of the theme issue ‘Heat and mass transfer in frost and ice’.

Science & Technology - Other Topics↗

Exploring quantum statistics for massive Dirac and Majorana neutrinos using spinor-helicity techniques

Recently, there has been interest in the applicability of quantum statistics to distinguish Dirac from Majorana neutrinos in multineutrino final states. In particular, debate has arisen over the validity of the Dirac-Majorana confusion theorem in these processes, i.e., that any distinction between the Dirac and Majorana processes goes to zero as the neutrino mass goes to zero. Here we approach this problem equipped with spinor-helicity methods generalized for massive Dirac and Majorana fermions. We explicitly calculate all helicity amplitudes, and their squares, for the decay of a light scalar particle to two neutrinos and two oppositely charged leptons. This allows us to pinpoint the crucial steps which could lead to claims of a violation of the confusion theorem. We show that, if the correct antisymmetrization of Dirac to Majorana amplitudes is used, identification of which is clear in this framework, and all relevant contributions are appropriately summed, a scalar decay into two charged leptons and two neutrinos satisfies the Dirac-Majorana confusion theorem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Braiding for the win: Harnessing braiding statistics in topological states to play quantum games

Nonlocal quantum games provide proof of principle that quantum resources can confer an advantage at certain tasks. They also provide a compelling way to explore the computational utility of phases of matter on quantum hardware. In a recent paper [O. Hart et al., Phys. Rev. Lett. 134, 130602 (2025)], we demonstrated that a toric code resource state conferred advantage at a certain nonlocal game, which remained robust to small deformations of the resource state. In this paper we demonstrate that this robust advantage is a generic property of resource states drawn from topological or fracton ordered phases of quantum matter. To this end, we illustrate how several other states from paradigmatic topological and fracton ordered phases can function as resources for suitably defined nonlocal games, notably the three-dimensional toric-code phase, the X-cube fracton phase, and the double-semion phase. The key in every case is to design a nonlocal game that harnesses the characteristic braiding processes of a quantum phase as a source of contextuality. We unify the strategies that take advantage of mutual statistics by relating the operators to be measured to order and disorder parameters of an underlying generalized symmetry-breaking phase transition. Additionally, by connecting the win probability to twist products, we show that success at the game serves as a many-body entanglement witness. Namely, if the players implement a perfect quantum strategy on large length scales, the quantum state they share cannot be connected to a trivial product state via a constant-depth local unitary circuit. Lastly, we massively generalize the family of games that admit perfect strategies when codewords of homological quantum error-correcting codes are used as resources.

Fractons↗

Kaon gluon parton distribution and momentum fraction from 2+1+1 lattice QCD with high statistics

We present a high-statistics lattice-QCD determination of the kaon gluon parton distribution function and gluon momentum fraction. We use clover valence fermion action to take 1,296,640 kaon-correlator measurements on a highly improved staggered quark ensemble with 𝑎 ≈ 0.12 fm and 310-MeV pion mass generated by the MILC Collaboration. A detailed investigation into the impact of gauge-link smearing on the gluonic matrix elements indicates that five steps of hypercubic smearing offer an effective balance between signal quality and preservation of long-distance physics. We report a nonperturbatively renormalized kaon gluon momentum fraction of Math output error at 𝜇 = 2 GeV in the Math output error scheme. Using reduced pseudo-Ioffe-time-distribution matrix elements and pseudo-parton-distribution-function (PDF) matching, we extract the kaon gluon PDF and compare with the prediction from the Dyson-Schwinger equation and with the pion PDF obtained from the same ensemble.

First-principles calculations↗

Relaxation mechanisms in a disordered system with Poisson-level statistics

Here we discuss the interplay between many-body localization and spin symmetry. To this end, we study the time evolution of several observables in the anisotropic t- j model. Like the Hubbard chain, the studied model contains charge and spin degrees of freedom, yet it has smaller Hilbert space and thus allows for numerical studies of larger systems. We compare the field disorder that breaks the $\mathbb{Z}_2$ spin symmetry and a potential disorder that preserves the latter symmetry. In the former case, sufficiently strong disorder leads to localization of all studied observables, at least for the studied system sizes. However, in the case of symmetry-preserving disorder, we observe that odd operators under the $\mathbb{Z}_2$ spin transformation relax towards the equilibrium value at relatively short timescales that grow only polynomially with the disorder strength. On the other hand, the dynamics of even operators and the level statistics within each symmetry sector are consistent with localization. Our results indicate that localization exists within each symmetry sector for symmetry-preserving disorder. Odd operators' apparent relaxation is due to their time evolution between distinct symmetry sectors.

36 MATERIALS SCIENCE↗

Statistical physics of multicomponent alloys using KKR-CPA

We apply variational principles from statistical physics and the Landau theory of phase transitions to multicomponent alloys using the multiple-scattering theory of Korringa-Kohn-Rostoker (KKR) and the coherent potential approximation (CPA). This theory is a multicomponent generalization of the S( 2 ) theory of binary alloys developed by G. M. Stocks, J. B. Staunton, D. D. Johnson and others. It is highly relevant to the chemical phase stability of high-entropy alloys as it predicts the kind and size of finite-temperature chemical fluctuations. In doing so it includes effects of rearranging charge and other electronics due to changing site occupancies. When chemical fluctuations grow without bound an absolute instability occurs and a second-order order-disorder phase transition may be inferred. The S( 2 ) theory is predicated on the fluctuation-dissipation theorem; thus we derive the linear response of the CPA medium to perturbations in site-dependent chemical potentials in great detail. The theory lends itself to a natural interpretation in terms of competing effects: entropy driving disorder and favorable pair interactions driving atomic ordering. Moreover, to further clarify interpretation we present results for representative ternary alloys CuAgAu, NiPdPt, RhPdAg, and CoNiCu within a frozen charge (or band-only) approximation. These results include the so-called Onsager mean field correction that extends the temperature range for which the theory is valid.

36 MATERIALS SCIENCE↗

Numerically exact full counting statistics of the nonequilibrium Anderson impurity model

The time-dependent full counting statistics of charge transport through an interacting quantum junction is evaluated from its generating function, controllably computed with the inchworm Monte Carlo method. Exact noninteracting results are reproduced; then, we continue to explore the effect of electron-electron interactions on the time-dependent charge cumulants, first-passage time distributions, and n-electron transfer distributions. We observe a crossover in the noise from Coulomb blockade to Kondo-dominated physics as the temperature is decreased. In addition, we uncover long-tailed spin distributions in the Kondo regime and analyze queuing behavior caused by correlations between single-electron transfer events

36 MATERIALS SCIENCE↗

Propagation of statistical uncertainties of Skyrme mass models to simulations of r -process nucleosynthesis

Uncertainties in nuclear models have a major impact on simulations that aim at understanding the origin of heavy elements in the universe through the rapid neutron capture process ( r process) of nucleosynthesis. Within the framework of the nuclear density functional theory, we use results of Bayesian statistical analysis to propagate uncertainties in the parameters of energy density functionals to the predicted r -process abundance pattern, by way not only of the nuclear masses but also through the influence of the masses on β -decay and neutron capture rates. We point out the importance of the nonequilibrium end stage of the r process in determining the width of the resulting abundance pattern uncertainty bands. We additionally make the first identifications of specific parameters of Skyrme-like energy density functionals which show tentative correlations with particular aspects of the r -process abundance pattern. While previous studies have explored the reduction in the abundance pattern uncertainties due to anticipated new measurements of neutron-rich nuclei, here we point out that an even larger reduction will occur when these new measurements are used to reduce the uncertainty of model predictions of masses, which are then propagated through to the abundance pattern. We make a quantitative prediction for how large this reduction will be.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Predicting the neutrinoless double-β-decay matrix element of 136 Xe using a statistical approach

Calculation of the nuclear matrix elements (NMEs) for double-β decay is of paramount importance for guiding experiments and for analyzing and interpreting the experimental data, especially for the search of the neutrinoless double β decay mode (0v⁢ββ). However, there are currently still large differences between the NME values calculated by different methods, hence a quantification of their uncertainties is very much required. Here, in this paper, we propose a statistical analysis of 0⁢v⁢ββ NME for the 136 Xe isotope, based on the interacting shell model, but using three independent effective Hamiltonians, emphasizing the range of the NMEs' most probable values and its correlations with observables that can be obtained from the existing nuclear data. Consequently, we propose a common probability distribution function for the 0⁢vββ NME, which has a range of (1.55–2.65) at 90% confidence level, with a mean value of 1.99 and a standard deviation of 0.37.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Solving the one-dimensional penetration problem for the fission channel in the statistical Hauser-Feshbach theory

Here we solve the Schrödinger equation for an arbitrary one-dimensional potential energy to calculate the transmission coefficient in the fission channel of compound nucleus reactions. We incorporate the calculated transmission coefficients into the statistical Hauser-Feshbach model calculation for neutron-induced reactions on 235,238 U and 239 Pu. The one-dimensional model reproduces the evaluated fission cross section data reasonably well considering the limited number of model parameters involved. A resonance-like structure appears in the transmission coefficient for a double-humped fission barrier shape that includes an intermediate well, which is understood to be a quantum mechanical effect in the fission channel. The calculated fission cross sections for the neutron-induced reactions on 235,238 U and 239 Pu all exhibit a similar structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Statistical error propagation in ab initio no-core full configuration calculations of light nuclei

We propagate the statistical uncertainty of experimental N N scattering data into the binding energy of 3 H and 4 He. Here, we also study the sensitivity of the magnetic moment and proton radius of the 3 H to changes in the N N interaction. The calculations are made with the no-core full configuration method in a sufficiently large harmonic oscillator basis. For those light nuclei we obtain Δ E stat (3H) = 0.015 MeV and Δ E stat ( 4 He) = 0.055 MeV .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Statistical Hauser-Feshbach theory with width-fluctuation correction including direct reaction channels for neutron-induced reactions at low energies

Here, a model to calculate particle-induced reaction cross sections with statistical Hauser-Feshbach theory including direct reactions is given. The energy average of the scattering matrix from the coupled-channels optical model is diagonalized by the transformation proposed by Engelbrecht and Weidenmüller [C. A. Engelbrecht and H. A. Weidenmüller, Phys. Rev. C 8, 859 (1973)]. The ensemble average of S -matrix elements in the diagonalized channel space is approximated by a model of Moldauer [P. A. Moldauer, Phys. Rev. C 12, 744 (1975)] using the newly parametrized channel degree-of-freedom ν a to better describe the Gaussian orthogonal ensemble (GOE) reference calculations. The Moldauer approximation is confirmed by a Monte Carlo study using a randomly generated S matrix, as well as the GOE threefold integration formula. Lastly, the method proposed is applied to the 238U(n,n') cross-section calculation in the fast-energy range, showing an enhancement in the inelastic scattering cross sections.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Estimation of M 1 scissors mode strength for deformed nuclei in the medium- to heavy-mass region by statistical Hauser-Feshbach model calculations

Radiative neutron capture is an important nuclear reaction whose accurate description is needed for many applications ranging from nuclear technology to nuclear astrophysics. The description of such a process relies on the Hauser-Feshbach theory which requires the nuclear optical potential, level density, and γ-strength function as model inputs. It has recently been suggested that the M 1 scissors mode may explain discrepancies between theoretical calculations and evaluated data. We explore statistical model calculations with the strength of the M 1 scissors mode estimated to be dependent on the nuclear deformation of the compound system. We show that the form of the M 1 scissors mode improves the theoretical description of evaluated data and the match to experiment in both the fission product and actinide regions. Since the scissors mode occurs in the range of a few keV to a few MeV, it may also impact the neutron capture cross sections of neutron-rich nuclei that participate in the rapid neutron capture process of nucleosynthesis. As a result, we comment on the possible impact to nucleosynthesis by evaluating neutron capture rates for neutron-rich nuclei with the M 1 scissors mode active.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Proton capture reaction cross section measurements on Er 162 as a probe of statistical model calculations

It is crucial to measure reaction cross sections relevant to the astrophysical γ process so that theoretical reaction rates can be tested and validated with experimental data. Here, the total cross sections for the 162 Er(p,γ) 163 Tm and the 162 Er(p,n) 162 Tm reactions have been measured by the activation method in center-of-mass energies from 3.973 to 8.944 MeV and from 5.962 to 8.944 MeV, respectively. The nucleus 162 Er is the heaviest p nuclide to be measured by the activation method using γ -ray spectroscopy, so far. It is important to note that the energy range for the (p,γ) reaction measurement covers a large fraction of the astrophysically relevant energy region between 2.71 and 5.34 MeV. The targets were prepared by evaporating 28.2 % isotopically enriched 162 Er 2 O 3 powder onto carbon backing foils, and bombarded with proton beams provided by the FN Tandem Accelerator at the University of Notre Dame. The reaction yields have been determined by the observed activity of produced radioactive isotopes, which was detected offline by a high-purity germanium detector. The results are presented and compared with calculations from two statistical model codes: NON-SMOKER and TALYS.

150 ≤ A ≤ 189↗

Statistical physics of black holes as quantum-mechanical systems

Some basic features of black-hole statistical mechanics are investigated, assuming that black holes respect the principles of quantum mechanics. Care is needed in defining an entropy S bh corresponding to the number of microstates of a black hole, given that the black hole interacts with its surroundings. An open question is then the relationship between this entropy and the Bekenstein-Hawking entropy S BH . For a wide class of models with interactions needed to ensure unitary quantum evolution, these interactions produce extra energy flux beyond that predicted by Hawking. Arguments are then presented that this results in an entropy S bh that is smaller than S BH . Correspondingly, in such scenarios equilibrium properties of black holes are modified. In this work, we examine questions of consistency of such an inequality; if it is not consistent, that provides significant constraints on models for quantum-mechanical black hole evolution.

79 ASTRONOMY AND ASTROPHYSICS↗