Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical physics”

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,063 records · Page 59

FlameNEST: explicit profile likelihoods with the Noble Element Simulation Technique

We present FlameNEST, a framework providing explicit likelihood evaluations in noble element particle detectors using data-driven models from the Noble Element Simulation Technique. FlameNEST provides a way to perform statistical analyses on real data with no dependence on large, computationally expensive Monte Carlo simulations by evaluating the likelihood on an event-by-event basis using analytic probability elements convolved together in a single TensorFlow multiplication. Furthermore, this robust framework creates opportunities for simple inter-collaboration analyses which will be fundamental for the future of experimental dark matter physics.

47 OTHER INSTRUMENTATION↗

Using Dense Array Waveform Correlations to Build a Velocity Model with Stochastic Variability

In an ongoing effort to improve 3D seismic-wave propagation modeling for frequencies up to 10 Hz, we used cross correlations between vertical-component waveforms from an underground chemical explosion to estimate the statistical properties of small-scale velocity heterogeneities. The waveforms were recorded by a dense 2D seismic array deployed during the Source Physics Experiments for event number 5 (SPE-5) in a series of six underground chemical explosions, conducted at the Nevada National Security Site. The array consisted of 996 geophones with a 50–100 m grid spacing, deployed at the SPE site at the north end of the Yucca Flat basin. The SPE were conducted to investigate the generation and propagation of seismic and acoustic waves from underground explosions. Comparisons of decay rates of waveform cross correlations as function of interstation distance, computed for observed and synthetic seismograms from the SPE-5 chemical explosion, were used to constrain statistical properties of correlated stochastic velocity perturbations representing small-scale heterogeneities added to a geology-based velocity model of the Yucca Flat basin. Using comparisons between recorded and simulated waveform cross correlations, we were able to recover sets of statistical properties of small-scale velocity perturbations in the velocity model that produce the best-fit between the recorded and simulated ground motion. The stochastic velocity fluctuations in the velocity model that produced the smallest misfits have a horizontal correlation distance of between 400 and 800 m, a vertical correlation distance between 100 and 200 m, and a standard deviation of 10% from the nominal model velocity in the alluvium basin layers. They also have a horizontal correlation distance of 1000 m, a vertical correlation distance of 250 m, and a standard deviation of 6% in the underlying and consolidated sedimentary layers, up to a depth of 4 km. Comparisons between observed and simulated wavefields were used to assess the proposed small-scale heterogeneity enhancements to the Yucca Flat basin model. We found that adding a depth-resolved stochastic variability to the geology-based velocity model improves the overall performance of ground-motion simulations of an SPE-5 explosion in the modeled frequency range up to 10 Hz. The results may be applicable to other similar basins.

58 GEOSCIENCES↗

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↗

An approximately 300 Mpc void of rich clusters of galaxies

Kirshner et al. (1981) have reported an apparent absence of relatively bright galaxies over an enormous volume in the direction of Bootes. This result, if statistically significant, has important implications for theories of galaxy clustering and formation. In connection with the study of the correlation function of rich clusters of galaxies, Bahcall and Soneira (1982) have noticed a large area in the sky, extending approximately 100 deg across, which is completely void of nearby, but not distant, rich Abell clusters. This apparent absence of rich clusters may represent a physical void. If this interpretation is correct, the void extends approximately 150 Mpc/h along the line of sight. The sample studied is Abell's (1958) statistical sample of clusters of galaxies. The distribution of rich clusters is discussed, taking into account observed versus expected characteristics. Attention is also given to aspects of absorption and galaxy counts in the considered region.

Bahcall, N. A.↗

Earth mover’s distance as a measure of CP violation

We introduce a new unbinned two sample test statistic sensitive to CP violation utilizing the optimal transport plan associated with the Wasserstein (earth mover’s) distance. The efficacy of the test statistic is shown via two examples of CP asymmetric distributions with varying sample sizes: the Dalitz distributions of B 0 → K + π – π 0 and of D 0 → π + π – π 0 decays. The windowed version of the Wasserstein distance test statistic is shown to have comparable sensitivity to CP violation as the commonly used energy test statistic, but also retains information about the localized distributions of CP asymmetry over the Dalitz plot. For large statistic datasets we introduce two modified Wasserstein distance based test statistics — the binned and the sliced Wasserstein distance statistics, which show comparable sensitivity to CP violation, but improved computing time and memory scalings. Finally, general extensions and applications of the introduced statistics are discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Novel approaches in hadron spectroscopy

The last two decades have witnessed the discovery of a myriad of new and unexpected hadrons. The future holds more surprises for us, thanks to new-generation experiments. Understanding the signals and determining the properties of the states requires a parallel theoretical effort. To make full use of available and forthcoming data, a careful amplitude modeling is required, together with a sound treatment of the statistical uncertainties, and a systematic survey of the model dependencies. We review the contributions made by the Joint Physics Analysis Center to the field of hadron spectroscopy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Machine learning approaches for structural and thermodynamic properties of a Lennard-Jones fluid

Predicting the functional properties of many molecular systems relies on understanding how atomistic interactions give rise to macroscale observables. However, current attempts to develop predictive models for the structural and thermodynamic properties of condensed-phase systems often rely on extensive parameter fitting to empirically selected functional forms whose effectiveness is limited to a narrow range of physical conditions. Here, we illustrate how these traditional fitting paradigms can be superseded using machine learning. Specifically, we use the results of molecular dynamics simulations to train machine learning protocols that are able to produce the radial distribution function, pressure, and internal energy of a Lennard-Jones fluid with increased accuracy in comparison to previous theoretical methods. The radial distribution function is determined using a variant of the segmented linear regression with the multivariate function decomposition approach developed by Craven et al. [J. Phys. Chem. Lett. 11, 4372 (2020)]. The pressure and internal energy are determined using expressions containing the learned radial distribution function and also a kernel ridge regression process that is trained directly on thermodynamic properties measured in simulation. The presented results suggest that the structural and thermodynamic properties of fluids may be determined more accurately through machine learning than through human-guided functional forms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reversible to irreversible transitions for cyclically driven particles on periodic obstacle arrays

Here, we examine the collective dynamics of disks moving through a square array of obstacles under cyclic square wave driving. Below a critical density, we find that the system organizes into a reversible state in which the disks return to the same positions at the end of every drive cycle. Above this density, the dynamics are irreversible and the disks do not return to the same positions after each cycle. The critical density depends strongly on the angle θ between the driving direction and a symmetry axis of the obstacle array, with the highest critical densities appearing at commensurate angles such as θ = 0° and θ = 45° and the lowest critical densities falling at θ = arctan (0.618), the inverse of the golden ratio, where the flow is the most degenerate. As the density increases, the number of cycles required to reach a reversible state grows as a power law with an exponent near ν = 1.36, similar to what is found in periodically driven colloidal and superconducting vortex systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Calculating adsorption isotherms using the two-phase thermodynamic method and molecular dynamics simulations

We describe the calculation of adsorption isotherms from molecular dynamics simulations based on the two-phase thermodynamic (2PT) model. The 2PT model developed for bulk fluid phases treats the gas-like components as hard spheres (HSs), which correctly recovers the limiting behaviors of unconfined fluids. We showed that this treatment, however, does not always lead to the correct zero-loading behavior in strongly confining systems. For methane adsorption into zeolite MFI, the HS reference state underestimates entropy by up to 20% at low loadings and leads to an order-of-magnitude increase in the adsorption onset pressure. To fix these issues, we propose the use of ideal adsorbed gas (IAG) as the gas reference model, the properties of which can be computed using the Widom insertion method on an empty adsorbent. We further describe three routes to compute adsorption isotherms from the Helmholtz free energy at different loadings. Comparing against established Monte Carlo (MC) methods, we found that the adsorption isotherms obtained using the IAG reference state agrees to within 40%, which corresponds to deviations of <5% in adsorption free energy. The isotherms calculated using the HS reference state underestimate the adsorption uptake at low to medium loadings in strongly confining systems, but its accuracy improves at higher loadings and as the pore size increases relative to the sorbate diameter. The methods described here provide an alternative approach for computing adsorption isotherms when MC simulations in an open ensemble are undesirable and enable a direct comparison of computed adsorption thermodynamics with experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Supervised machine learning-based multivariate regression of parallel closures for a high-collisionality deuterium-carbon plasma

Many plasmas of interest in laboratory experiments and space consist of multiple ion species. In tokamak edge plasmas, for instance, ionized impurities expelled from the vessel wall influence plasma transport. When describing multi-species plasmas using fluid equations, we need accurate closure relations to close the set of fluid equations. In this study, we introduce the development of fitting formulas for parallel closures using supervised machine learning, in conjunction with the recent closure theory, considering multi-ion collisions and arbitrary ion temperatures. We apply this approach to a high-collisionality deuterium-carbon plasma and demonstrate its effectiveness. As a result, the machine learning-based method for developing practical and accurate closures can be extended to a wider range of plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Thermal mean-field theories

Several closely related ab initio thermal mean-field theories for fermions, both well-established and new ones, are compared with one another at the formalism level and numerically. The theories considered are Fermi–Dirac theory; thermal Hartree–Fock (HF) theory; two modifications of the thermal single-determinant and the first-order finite-temperature many-body perturbation theory based on a zero-temperature or thermal HF reference. Furthermore, thermal full-configuration-interaction theory is used as the benchmark.

74 ATOMIC AND MOLECULAR PHYSICS↗

Analysis of streaked images of x-ray self-emission in laser-driven spherical implosions

Imaging of x-ray self-emission provides a powerful in situ measurement of the spatial and temporal evolution of high-energy-density plasmas. However, interpretation of these measurements requires detailed understanding of the data-generating process. This work presents a case study in the interpretation of x-ray self-emission data for the specific application of streaked one-dimensional slit imaging of spherical laser-driven implosions. A comprehensive generative model of the streaked slit-imaging diagnostic is developed including detailed treatments of the radiation transfer, photometrics, and photostatistics associated with the measurement. The model is used to generate realistic synthetic streaked images and to analyze experimental streaked images to extract important physical quantities of interest. An example analysis of streaked images from implosion experiments on the OMEGA laser is presented, where the model developed in this work is used to constrain the trajectory and peak velocity of the implosion using Bayesian inference.

Bayesian inference↗

Scattering theory of frequency-entangled biphoton states facilitated by cavity polaritons

The use of quantum light to probe exciton properties in semiconductor and molecular nanostructures typically occurs in the low-intensity regime. A substantial enhancement of exciton-photon coupling can be achieved with photonic cavities, where excitons hybridize with cavity modes to form polariton states. Here, to provide a theoretical framework for interpreting emerging experimental efforts in this direction, we develop a scattering theory describing the interaction of frequency-entangled photon pairs with cavity polariton and bipolariton states under various coupling regimes. Employing the Tavis-Cummings model in combination with our scattering approach, we present a quantitative analysis of how the interaction of the entangled photon pair with the polariton or bipolariton modifies its joint spectral amplitude (JSA). Specifically, we examine the effects of the cavity-mode steady-state population, exciton-cavity coupling strength, and different forms of the input photon JSA. Our results show that the entanglement entropy of the scattered photons is highly sensitive to the interplay between the input JSA and the spectral line shapes of the polariton resonances, emphasizing the cavity filtering effects. We suggest that biphoton-scattering quantum light spectroscopy best serves as a sensitive probe of polariton and bipolariton states in the photon-vacuum cavity state. Our approach is not only robust to various regimes of cavity-exciton coupling, but also amenable to extensions beyond the Tavis-Cummings model, enabling the representation of a broad class of molecular systems and solid state quantum materials.

36 MATERIALS SCIENCE↗

Strangeness-correlations on the pseudocritical line in ( 2 + 1 )-flavor QCD

We present some lattice QCD results on first ( χ 1 i ) and second ( χ 2 i ) cumulants of and correlations ( χ 11 i j ) among net baryon-number ( B ), strangeness ( S ) and electric charge ( Q ) along the pseudocritical line [ T p c ( μ B ) ] in the temperature ( T )–baryon chemical potential ( μ B ) phase diagram of ( 2 + 1 )-flavor QCD. We point out that violations of sum rules among second order cumulants, which hold in the isospin symmetric limit of vanishing electric charge chemical potential, are small along the T p c ( μ B ) for the entire range of μ B covered in the RHIC beam energy scan. For the strangeness neutral matter produced in heavy-ion collisions this leads to a close relation between χ 11 B S and χ 11 Q S . We compare lattice QCD results for χ 11 B S / χ 2 S along the T p c ( μ B ) line with preliminary experimental measurements of χ 11 B S / χ 2 S for collision energies 7.7 GeV ≤ s N N ≤ 62.4 GeV . While we find good agreements for s N N ≥ 39 GeV , differences are sizeable at smaller values of s N N . Moreover, we compare lattice QCD results for the ratio of the strangeness ( μ S ) to baryon ( μ B ) chemical potentials, which define a strangeness neutral system with fixed electric charge to baryon number density, with experimental results obtained by the STAR collaboration for μ S / μ B using strange baryon yields on the freeze-out line. Finally, we determine the baryon chemical potential at the freeze-out ( μ B f ) by comparing χ 1 B / χ 2 B along the T p c ( μ B ) with the experimentally measured net-proton cumulants χ 1 p / χ 2 p . We find that { μ B f , T p c ( μ B f ) } are consistent with the freeze-out parameters of the statistical-model fits to experimentally measured hadron yields for s N N ≥ 11.5 GeV . Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Quantum Jarzynski Equality in Open Quantum Systems from the One-Time Measurement Scheme

In open quantum systems, a clear distinction between work and heat is often challenging, and extending the quantum Jarzynski equality to systems evolving under general quantum channels beyond unitality remains an open problem in quantum thermodynamics. In this Letter, we introduce well-defined notions of guessed quantum heat and guessed quantum work, by exploiting the one-time measurement scheme, which only requires an initial energy measurement on the system alone. We derive a modified quantum Jarzynski equality and the principle of maximum work with respect to the guessed quantum work, which requires the knowledge of the system only. We further show the significance of guessed quantum heat and work by linking them to the problem of quantum hypothesis testing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗