Search NASA⌕ Search

SEARCH · Search NASA

Results for “Particle 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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Zero-pairing and zero-temperature limits of finite-temperature Hartree-Fock-Bogoliubov theory

Recently, variational Hartree-Fock-Bogoliubov (HFB) mean-field equations were shown to possess a mathematically well-defined zero-pairing limit, independently of the closed- or open-shell character of the system under consideration. This limit is non-trivial for open-shell systems such that HFB theory does {\it not} reduce to the Hartree-Fock (HF) formalism in all cases. The present work extends this analysis to finite-temperature HFB (FTHFB) theory by investigating the behavior of this more general formalism in the combined zero-temperature and zero-pairing limits. The zero-pairing and zero-temperature limits of the FTHFB statistical density operator constrained to carry an arbitrary (integer) number of particles A on average is worked out analytically and realized numerically using a two-nucleon interaction. While the FTHFB density operator reduces to the projector corresponding to a pure HF Slater determinant for closed-shell nuclei, the FTHFB formalism does not reduce to the HF theory in all cases in the zero-temperature and zero-pairing limits, i.e. for open-shell nuclei. However, the fact that a nucleus can be of open-shell character in these joint limits is necessarily the result of some symmetry restrictions. Whenever it is the case, the non-trivial description obtained for open-shell systems is shown to depend on the order with which both limits are taken, i.e. the two limits do not commute for these systems. When the zero-temperature limit is performed first, the FTHFB density operator is demoted to a projector corresponding to a pure state made out of a linear combination of a finite number of Slater determinants with different (even) numbers of particles. When the zero-pairing limit is performed first, the FTHFB density operator remains a statistical mixture of a finite number of Slater determinants with both even and odd particle numbers. While the entropy (pairing density) is zero in the first (second) case, it does not vanish in the second (first) case in spite of the temperature (pairing) tending towards zero. The difference between both limits can have striking consequences for the (thermal) expectation values of observables. For instance, the particle-number variance does not vanish in either case and has limiting values that differ by a factor of two in both cases. In conclusion, while in the textbook situation associated with closed-shell nuclei Hartree-Fock-Bogoliubov (finite-temperature Hartree-Fock) theory reduces to Hartree-Fock theory in the zero-pairing (zero-temperature) limit, the present analysis demonstrates that a non trivial and unexpected limit is obtained for this formalism in open-shell systems. This result sheds a new light on certain aspects of this otherwise very well-studied many-body formalism.

20 ≤ A ≤ 38↗
DOE OSTI2022

Origin of nontopological soliton dark matter: solitosynthesis or phase transition

This work demonstrates that nontopological solitons with large global charges and masses, even above the Planck scale, can form in the early universe and dominate the dark matter abundance. In solitosynthesis, solitons prefer to grow as large as possible under equilibrium dynamics when an initial global charge asymmetry is present. Their abundance is set by when soliton formation via particle fusion freezes out, and their charges are set by the time it takes to accumulate free particles. This work improves the estimation of both quantities, and in particular shows that much larger-charged solitons form than previously thought. The results are estimated analytically and validated numerically by solving the coupled Boltzmann equations. Without solitosynthesis, phase transitions can still form solitons from particles left inside false-vacuum pockets and determine their present-day abundance and properties. Even with zero charge asymmetry, solitons formed in this way can have very large charges on account of statistical fluctuations in the numbers of (anti)particles inside each pocket.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Characterizing the expected behavior of non-Poissonian template fitting

We have performed a systematic study of the statistical behavior of non-Poissonian template fitting (NPTF), a method designed to analyze and characterize unresolved point sources in general counts datasets. Here, in this paper, we focus on the properties and characteristics of the Fermi-LAT gamma-ray data set. In particular, we have simulated and analyzed gamma-ray sky maps under varying conditions of exposure, angular resolution, pixel size, energy window, event selection, and source brightness. We describe how these conditions affect the sensitivity of NPTF to the presence of point sources, for inner-galaxy studies of point sources within the Galactic Center excess, and for the simplified case of isotropic emission. We do not find opportunities for major gains in sensitivity from varying these choices, within the range available with current Fermi-LAT data. We provide an analytic estimate of the NPTF sensitivity to point sources for the case of isotropic emission and perfect angular resolution, and find good agreement with our numerical results for that case.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Introducing the Video In Situ Snowfall Sensor (VISSS)

The open-source Video In Situ Snowfall Sensor (VISSS) is introduced as a novel instrument for the characterization of particle shape and size in snowfall. The VISSS consists of two cameras with LED backlights and telecentric lenses that allow accurate sizing and combine a large observation volume with relatively high pixel resolution and a design that limits wind disturbance. VISSS data products include various particle properties such as maximum extent, cross-sectional area, perimeter, complexity, and sedimentation velocity. Initial analysis shows that the VISSS provides robust statistics based on up to 10 000 unique particle observations per minute. Comparison of the VISSS with the collocated PIP (Precipitation Imaging Package) and Parsivel instruments at Hyytiälä, Finland, shows excellent agreement with the Parsivel but reveals some differences for the PIP that are likely related to PIP data processing and limitations of the PIP with respect to observing smaller particles. The open-source nature of the VISSS hardware plans, data acquisition software, and data processing libraries invites the community to contribute to the development of the instrument, which has many potential applications in atmospheric science and beyond.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Data From Experiments on Bubbling Fluidization of Zeolite in a Rectangular Bubbling Fluidized Bed

Fluidization experiments were conducted in a lab-scale rectangular bubbling fluidized bed with the objective of generating a high-quality dataset for model validation and artificial intelligence/machine learning (AI/ML) training. Zeolite was chosen as the bed material, and the fluidizing medium was air as supplied by a compressor. Three different flow rates at the inlet were chosen such that the particles were fluidized but not elutriated from the system. The test matrix involved randomization and replicates to provide uncertainty estimates as well as four different batches of zeolite as the bed material. The quantities of interest obtained from this study were statistics of differential pressures, interface heights, and particle velocities. Considering all the components of the elaborate test plan, the results obtained were consistent and reproducible. Characterization tests were performed to estimate particle properties including size, density, coefficient of friction, coefficient of restitution, and minimum fluidization velocity. In addition, the angle of repose from granular discharge experiments has been reported to account for rolling friction, though its effect on the overall process is expected to be negligible.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Study of η′ → π+π−l+l− decays at BESIII

Abstract With a sample of (10087±44)×10 6 J/ψevents accumulated with the BESIII detector, we analyze the decaysη′→ π + π − l + l − (l=e, μ) via the processJ/ψ → γη′. The branching fractions are measured to be$$ \mathcal{B} $$ B (η′→ π + π − e + e − ) = (2.45±0.02(stat.)±0.08(syst.))×10 −3 and$$ \mathcal{B} $$ B (η′→ π + π − μ + μ − ) = (2.16±0.12(stat.)±0.06(syst.))×10 −5 , and the ratio is$$ \frac{\mathcal{B}\left({\eta}^{\prime}\to {\pi}^{+}{\pi}^{-}{e}^{+}{e}^{-}\right)}{\mathcal{B}\left({\eta}^{\prime}\to {\pi}^{+}{\pi}^{-}{\mu}^{+}{\mu}^{-}\right)}=113.4\pm 0.9\left(\textrm{stat}.\right)\pm 3.7\left(\textrm{syst}.\right) $$ B η ′ → π + π − e + e − B η ′ → π + π − μ + μ − = 113.4 ± 0.9 stat . ± 3.7 syst . . In addition, by combining theη′ →π + π − e + e − andη′ →π + π − μ + μ − decays, the slope parameter of the electromagnetic transition form factor is measured to beb η′ = 1.30 ± 0.19 (GeV/c 2 ) −2 , which is consistent with previous measurements from BESIII and theoretical predictions from the VMD model. The asymmetry in the angle between theπ + π − andl + l − decay planes, which has the potential to reveal theCP-violation originating from an unconventional electric dipole transition, is also investigated. The asymmetry parameters are determined to be$$ {\mathcal{A}}_{CP}\left({\eta}^{\prime}\to {\pi}^{+}{\pi}^{-}{e}^{+}{e}^{-}\right)=\left(-0.21\pm 0.73\left(\textrm{stat}.\right)\pm 0.01\left(\textrm{syst}.\right)\right)\% $$ A CP η ′ → π + π − e + e − = − 0.21 ± 0.73 stat . ± 0.01 syst . % and$$ {\mathcal{A}}_{CP}\left({\eta}^{\prime}\to {\pi}^{+}{\pi}^{-}{\mu}^{+}{\mu}^{-}\right)=\left(0.62\pm 4.71\left(\textrm{stat}.\right)\pm 0.08\left(\textrm{syst}.\right)\right)\% $$ A CP η ′ → π + π − μ + μ − = 0.62 ± 4.71 stat . ± 0.08 syst . % , implying that no evidence ofCP-violation is observed at the present statistics. Finally, an axion-like particle is searched for via the decayη′ →π + π − a, a→e + e − , and upper limits of the branching fractions are presented for the mass assumptions of the axion-like particle in the range of 0−500 MeV/c 2 .

Physics↗
DOE OSTI2024

Dynamic data-driven multiscale modeling for predicting the degradation of a 316L stainless steel nuclear cladding material

Here, we have developed a long short-term memory stacked ensemble (LSTM-SE) surrogate modeling approach that can provide rapid predictions of microstructural evolution and the resultant mechanical properties of American Iron and Steel Institute (AISI) 316L series stainless steel (316LSS) fuel cladding under conditions of varying temperature and radiation dose rate. To acquire training data, we developed and implemented a kinetic Monte Carlo (KMC) model to simulate precipitation kinetics of M 23 C 6 , γ', and G phases within SS316L cladding. Experimentally reported precipitation kinetics of SS316L in literature were linked to the kinetic parameters of the simulated precipitation in our KMC model. The model was then used to simulate microstructure evolution under synthetically generated treatments of varying temperature and radiation dose rate, for periods of up to 3000 hours. Changes in volume fraction, number density, and particle size of precipitates were recorded, and particle area fractions were correlated using statistical methods to develop the surrogate model. Simultaneously, the mechanical properties of the simulated microstructures were evaluated using microstructure-based finite element method (FEM) analysis to determine the elastic modulus, yield stress, ultimate tensile strength, and elongation to failure of the aged microstructures. Using this approach, our surrogate model can predict precipitation behavior within 0.25% volume fraction and mechanical properties within 6% relative error from the values predicted by the KMC and FEM models using 50 training simulations as input. The trained recurrent neural network-based model can return estimations of precipitation kinetics and mechanical properties ~1000 times faster than the physics-based codes. This work demonstrates, as a proof of concept, that reactor material service lifetimes under variable service conditions can be predicted for a statistics-based model from a practicably obtainable dataset.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Orbit-averaging and deposition accuracy for runaway electron beams in hybrid kinetic-MHD simulations of the runaway plateau

We develop a new procedure that combines the kinetic orbit runaway electron code (KORC) and the NIMROD extended-magnetohydrodynamic code to simulate runaway electrons (REs) in the post-disruption plateau. KORC integrates guiding-center orbits, with a barycentric-based binary search strategy providing initial guesses for the Newton–Raphson logical-to-physical coordinate inversion, ensuring reliable particle-to-mesh mapping in NIMROD, whose fields remain static for the present study. Samples are drawn in accord with experimental parallel current profiles of RE beams during the plateau phase. Deposition in NIMROD is verified through comparison with a Python-based finite-element code that ensures periodicity in the poloidal direction and continuity at the magnetic axis. Accurate representation of near-axis fields requires finer mesh resolution to prevent under- and overshoots in current density from orbit inaccuracies. Yet, at a fixed particle count, increasing mesh resolution amplifies statistical noise in the deposited fields. An orbit-averaging method accumulates partial current deposits over multiple kinetic steps and reduces the statistical noise with little added computational cost. By coupling kinetic routines from KORC directly into the NIMROD codebase, these developments lay essential groundwork for future self-consistent KORC–NIMROD coupling.

Algorithms and data structure↗
DOE OSTI2022

Phase diagram of superconductivity in the integer quantum Hall regime

An interplay between pairing and topological orders has been predicted to give rise to superconducting states supporting exotic emergent particles, such as Majorana particles obeying non-Abelian braid statistics. We consider a system of spin polarized electrons on a Hofstadter lattice with nearest-neighbor attractive interaction and solve the mean-field Bogoliubov–de Gennes equations in a self-consistent fashion, leading to gauge-invariant observables and a rich phase diagram as a function of the chemical potential, the magnetic field, and the interaction. As the strength of the attractive interaction is increased, the system first makes a transition from a quantum Hall phase to a skyrmion lattice phase that is fully gapped in the bulk but has topological chiral edge current, characterizing a topologically nontrivial state. This is followed by a vortex phase in which the vortices carrying Majorana modes form a lattice; the spectrum contains a low-energy Majorana band arising from the coupling between neighboring vortex-core Majorana modes but does not have chiral edge currents. For some parameters, a dimer vortex lattice occurs with no Majorana band. The experimental feasibility and the observable consequences of skyrmions as well as Majorana modes are indicated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Statistical mechanical model for crack growth

Analytic relations that describe crack growth are vital for modeling experiments and building a theoretical understanding of fracture. Upon constructing an idealized model system for the crack and applying the principles of statistical thermodynamics, it is possible to formulate the rate of thermally activated crack growth as a function of load, but the result is analytically intractable. In this report an asymptotically correct theory is used to obtain analytic approximations of the crack growth rate from the fundamental theoretical formulation. These crack growth rate relations are compared to those that exist in the literature and are validated with respect to Monte Carlo calculations and experiments. The success of this approach is encouraging for future modeling endeavors that might consider more complicated fracture mechanisms, such as inhomogeneity or a reactive environment.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dynamic Data-Driven Multiscale Modeling for Predicting Microstructure and Mechanical Properties of 316L Stainless Steel

We have developed a long short-term memory stacked ensemble (LSTM-SE) surrogate modeling approach that can provide rapid predictions of microstructural evolution and the resultant mechanical properties of American Iron and Steel Institute (AISI) 316L series stainless steel (316LSS) fuel cladding under conditions of varying temperature and radiation dose rate. To acquire training data, we developed and implemented a kinetic Monte Carlo (KMC) model to simulate precipitation kinetics of M 23 C 6 , γ', and G phases within SS316L cladding. Experimentally reported precipitation kinetics of SS316L in literature were linked to the kinetic parameters of the simulated precipitation in our KMC model. The model was then used to simulate microstructure evolution under synthetically generated treatments of varying temperature and radiation dose rate, for periods of up to 3000 hours. Changes in volume fraction, number density, and particle size of precipitates were recorded, and particle area fractions were correlated using statistical methods to develop the surrogate model. Simultaneously, the mechanical properties of the simulated microstructures were evaluated using microstructure-based finite element method (FEM) analysis to determine the elastic modulus, yield stress, ultimate tensile strength, and elongation to failure of the aged microstructures. Using this approach, our surrogate model can predict precipitation behavior within 0.25% volume fraction and mechanical properties within 6% relative error from the values predicted by the KMC and FEM models using 50 training simulations as input. The trained recurrent neural network-based model can return estimations of precipitation kinetics and mechanical properties ~1000 times faster than the physics-based codes. This work demonstrates, as a proof of concept, that reactor material service lifetimes under variable service conditions can be predicted for a statistics-based model from a practicably obtainable dataset.

36 MATERIALS SCIENCE↗
DOE OSTI2021

General Quantum Field Theory of Flavor Mixing and Oscillations

We review the canonical transformation in quantum physics known as the Bogoliubov transformation and present its application to the general theory of quantum field mixing and oscillations with an arbitrary number of mixed particles with either boson or fermion statistics. The mixing relations for quantum states are derived directly from the definition of mixing for quantum fields and the unitary inequivalence of the Fock space of energy and flavor eigenstates is shown by a straightforward algebraic method. The time dynamics of the interacting fields is then explicitly solved and the flavor oscillation formulas are derived in a unified general formulation with emphasis on antiparticle content and effect introduced by nontrivial flavor vacuum.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Optimization Studies of Radiation Shielding for the PIP-II Project at Fermilab

The PIP-II project at Fermilab, which includes an 800-MeV superconducting LINAC, demands rigorous radiation shielding optimization to meet safety requirements. We updated the MARS geometry model to reflect new magnet and collimator designs and introduced high-resolution detector planes to better capture radiation field distributions. To overcome the significant computational demands, we implemented a well-known branching technique that drastically reduced simulation runtimes while maintaining statistical integrity. This was achieved through particle splitting and the application of Russian Roulette techniques. Additionally, new graphical tools were created to streamline data visualization and MARS code usability.

Makovec, Alajos [Fermilab] (ORCID:0000000286157492↗
DOE OSTI2021

Understanding protein-complex assembly through grand canonical maximum entropy modeling

Inside a cell, heterotypic proteins assemble in inhomogeneous, crowded systems where the abundance of these proteins vary with cell types. While some protein complexes form putative structures that can be visualized with imaging, there are far more protein complexes that are yet to be solved because of their dynamic associations with one another. Nevertheless, it is possible to infer these protein complexes through a physical model. However, it is often not clear to physicists what kind of data from biology is necessary for such a modeling endeavor. Here, we aim to model these clusters of coarse-grained protein assemblies from multiple subunits through the constraints of interactions among the subunits and the chemical potential of each subunit. We obtained the constraints on the interactions among subunits from the known protein structures. We inferred the chemical potential that dictates the particle number distribution of each protein subunit from the knowledge of protein abundance from experimental data. Guided by the maximum entropy principle, we formulated an inverse statistical mechanical method to infer the distribution of particle numbers from the data of protein abundance as chemical potentials for a grand canonical multicomponent mixture. Using grand canonical Monte Carlo simulations, we captured a distribution of high-order clusters in a protein complex of succinate dehydrogenase with four known subunits. The complexity of hierarchical clusters varies with the relative protein abundance of each subunit in distinctive cell types such as lung, heart, and brain. When the crowding content increases, we observed that crowding stabilizes emergent clusters that do not exist in dilute conditions. We, therefore, proposed a testable hypothesis that the hierarchical complexity of protein clusters on a molecular scale is a plausible biomarker of predicting the phenotypes of a cell.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Towards an understanding of particle-scale flaws and microstructure evolution in cold spray via accumulation of single particle impacts

Cold spray coatings are the sum of countless individual bonding events between single particles impacting on top of one another at high velocities. Thus, the collective behavior of microparticles must be considered to elucidate the origins of coating flaws at the scale of the particles and larger, or the dynamic evolution of the overall coating microstructure. Laser-induced particle impact testing (LIPIT) has been extensively used to study single-particle impacts, and in this work is adapted to study the accumulation of numerous particles with knowledge of each individual particle's impact parameters (particle size, velocity). The method reproducibly deposits stacks of gold particles (>20 particles) with different characteristic spectra of impact velocity. The observation of impact-induced erosion lets us define a critical velocity for material-build-up that is higher than that for single-particle bonding. The quantitative single-particle data are analyzed in a correlative manner to the structure and flaws in the resulting stacks, providing some first statistical connections between, e.g., strain and recrystallization, or aberrant particle characteristics and defects. The results highlight opportunities for the study of many-particle phenomena in microparticle impact—from interaction of particles in cold spray to multi-step erosion processes—with a quantitative view of the behavior of single particles.

Erosion↗
DOE OSTI2025

A radar view of ice microphysics and turbulence in Arctic cloud systems

Ice microphysical processes are inherently complex because of their sensitivity to temperature and humidity, the diversity of ice crystal habits, and their interaction with supercooled liquid water (SCL) and turbulence. Long-term surface-based radar observations have been systematically used to unravel the different processes that affect ice particle growth. In this study, we present a statistical analysis of 6.5 years of Ka-band radar observations in Arctic cloud systems, combined with thermodynamic profiles derived from radiosonde measurements. For the first time, ice particle growth and sublimation – diagnosed from vertical gradients of radar reflectivity and mean Doppler velocity – are systematically mapped across a broad range of temperature and moisture conditions. These vertical gradients correspond closely to saturation levels relative to ice and exhibit a strong temperature dependence in supersaturated regions. Notably, distinct signatures near −15 °C are indicative of dendritic growth. Turbulence, quantified via the eddy dissipation rate (EDR), is most frequently observed in regions containing SCL. The co-occurrence of SCL and elevated turbulence results in significantly enhanced ice particle growth compared to conditions in which either is present alone. This work provides new observational constraints that are critical for improving the representation of ice microphysics in atmospheric models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Collective mode Brownian dynamics: A method for fast relaxation of statistical ensembles

Sampling equilibrium configurations of correlated systems of particles with long relaxation times (e.g., polymeric solutions) using conventional molecular dynamics and Monte Carlo methods can be challenging. This is especially true for systems with complicated, extended bond network topologies and other interactions that make the use and design of specialized relaxation protocols infeasible. We introduce a method based on Brownian dynamics simulations that can reduce the computational time it takes to reach equilibrium and draw decorrelated samples. Importantly, the method is completely agnostic to the particle configuration and the specifics of interparticle forces. In particular, we develop a mobility matrix that excites non-local, collective motion of N particles and can be computed efficiently in $\mathcal{O}(N)$ time. Particle motion in this scheme is computed by integrating the overdamped Langevin equation with an Euler–Maruyama scheme, in which Brownian displacements are drawn efficiently using a low-rank representation of the mobility matrix in position and wave space. We demonstrate the efficacy of the method with various examples from the realm of soft condensed matter and release a massively parallel implementation of the code as a plugin for the open-source package HOOMD-blue [J. A. Anderson et al., J. Comput. Phys. 227, 5342 (2008) and J. Glaser et al., Comput. Phys. Commun. 192, 97 (2015)] which runs on graphics processing units

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Analysis of thermal radiation of a gas-particulate-cloud in a pressurized oxy-coal combustor

Pressurized oxy-combustion (POC) aims to reduce carbon emissions from coal-fired power plants. In the POC process, pulverized coal is burned at elevated pressure in an O2-CO2 environment. The elevated pressure and high CO2 and H2O concentrations strongly impact thermal radiation, thereby significantly distinguishing POC system from conventional atmospheric-pressure air-fired and oxy-fuel combustion. Consequently, a thorough understanding of thermal radiation in POC is required to promote the development of novel combustors. The present computational work analyzes thermal radiation of the gaseous and particle phases in the pilot-scale POC combustor in Washington University in St. Louis. Particularly, the radiation characteristics of the gaseous and particulate cloud were estimated by employing the statistical narrow band model and the large-particle model. It is found that thermal radiation of a particulate cloud dominates in the combustor under the conditions of a furnace temperature of 1500 K and no substantial particle loss. Another important result is that the gas-and-particulate cloud can be approximately treated as a graybody under the same conditions. To be more specific, the results of spectral radiation intensity of a gas comprised of 40% (vol) H2O and 60% CO2 show that when the pressure is increased to 15 bar, and the radiation pathlength to 100 cm, the spectral radiation profile of the gas phase approaches that of a blackbody at the respective temperature. In addition, the emissivity of the particulate cloud has been evaluated as a function of the particle concentration and diameter by means of the large particle model. It is shown that the emissivity grows with the particle concentration but decreases with the particle size for the same mass of a particulate cloud. The emissivity of a particulate cloud in oxy- combustion exceeds that in air-fired combustion. This work could be used to validate the assumption of the gray gas model adopted in numerical simulations.

Li, Lei↗