Search NASA⌕ Search

SEARCH · Search NASA

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Kinetic modeling of E-to-H mode transition in inductively coupled hydrogen plasmas

Radio Frequency (RF) Inductively Coupled Plasmas (ICPs) are widely known for their two discharge modes, i.e., H-mode and E-mode, where the dynamics of the plasmas are completely different from each other. We have performed a kinetic simulation of a hydrogen plasma discharge in order to clarify the discharge mechanism and the E-to-H transition of the RF ICPs. The numerical simulation results, such as the time variations of spatial distribution of electron density and the power dissipated in the plasma, show the characteristic changes of the plasma dynamics due to E-to-H mode transition. Especially, the drastic change during the mode transition has been observed in the time evolution of the electron energy distribution function (EEDF). The EEDF deviates from a Maxwellian distribution before/after the transition and the deviation is more significant in the E-mode phase. These results indicate the importance of kinetic modeling for the physical understanding of E-to-H transition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

The charge imbalance in ultracold plasmas

Ultracold plasmas are regarded as quasineutral but not strictly neutral. The results of charge imbalance in the expansion of ultracold plasmas are reported. The calculations are performed by a full molecular-dynamics simulation. The details of the electron velocity distributions are calculated without the assumption of electron global thermal equilibrium and Boltzmann distribution. Spontaneous evolutions of the charge imbalance from the initial states with perfect neutrality are given in the simulations. The expansion of outer plasma slows down with the charge imbalance. The influences of plasma size and parameters on the charge imbalance are discussed. The radial profiles of electron temperature are given for the first time, and the self-similar expansion can still occur even if there is no global thermal equilibrium. The electron disorder induced heating is also found in the simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

A comment on the position dependent diffusion coefficient representation of structural heterogeneity

Experimental studies of the variation of the mean square displacement (MSD) of a particle in a confined colloid suspension that exhibits density variations on the scale length of the particle diameter are not in agreement with the prediction that the spatial variation in MSD should mimic the spatial variation in density. The predicted behavior is derived from the expectation that the MSD of a particle depends on the system density and the assumption that the force acting on a particle is a point function of position. The experimental data are obtained from studies of the MSDs of particles in narrow ribbon channels and between narrowly spaced parallel plates and from new data, reported herein, of the radial and azimuthal MSDs of a colloid particle in a dense colloid suspension confined to a small circular cavity. In each of these geometries, a dense colloid suspension exhibits pronounced density oscillations with spacing of a particle diameter. We remove the discrepancy between prediction and experiment using the Fisher-Methfessel interpretation of how local equilibrium in an inhomogeneous system is maintained to argue that the force acting on a particle is delocalized over a volume with radius equal to a particle diameter. Our interpretation has relevance to the relationship between the scale of inhomogeneity and the utility of translation of the particle MSD into a position dependent diffusion coefficient and to the use of a spatially dependent diffusion coefficient to describe mass transport in a heterogeneous system.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2019

Voltage effects on the stability of Pd ensembles in Pd–Au/Au(111) surface alloys

The catalytic performance of multimetallic electrodes is often attributed to a beneficial combination of ligand, strain, and ensemble effects. Understanding the influence of the electrochemical environment on the stability of the alloy surface structure is thus a crucial component to the design of highly active and durable electrocatalysts. In this work, we study the effects of an applied voltage to electrocatalytic Pd–Au/Au(111) surface alloys in contact with a model continuum electrolyte. Using planewave density functional theory, two-dimensional cluster expansions are parameterized and used to simulate dilute Pd–Au surface alloys under electrochemical conditions via Metropolis Monte Carlo within an extended canonical ensemble. While Pd monomers are stable at all potentials considered, different extents of surface electrification are observed to promote the formation of Pd dimers and trimers, as well as clusters of Pd monomers. We find that the relative proportion of monomer, dimer, and trimer surface fractions is in good agreement with in situ scanning tunneling microscopy measurements. The further development and refinement of the approaches described herein may serve as a useful aid in the development of next-generation electrocatalysts.

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

Theory of the effect of external stress on the activated dynamics and transport of dilute penetrants in supercooled liquids and glasses

We generalize the self-consistent cooperative hopping theory for a dilute spherical penetrant or tracer activated dynamics in dense metastable hard sphere fluids and glasses to address the effect of external stress, the consequences of which are systematically established as a function of matrix packing fraction and penetrant-to-matrix size ratio. All relaxation processes speed up under stress, but the difference between the penetrant and matrix hopping (alpha relaxation) times decreases significantly with stress corresponding to less time scale decoupling. A dynamic crossover occurs at a critical “slaving onset” stress beyond which the matrix activated hopping relaxation time controls the penetrant hopping time. This characteristic stress increases (decreases) exponentially with packing fraction (size ratio) and can be well below the absolute yield stress of the matrix. Below the slaving onset, the penetrant hopping time is predicted to vary exponentially with stress, differing from the power law dependence of the pure matrix alpha time due to system-specificity of the stress-induced changes in the penetrant local cage and elastic barriers. An exponential growth of the penetrant alpha relaxation time with size ratio under stress is predicted, and at a fixed matrix packing fraction, the exponential relation between penetrant hopping time and stress for different size ratios can be collapsed onto a master curve. Direct connections between the short- and long-time activated penetrant dynamics and between the penetrant (or matrix) alpha relaxation time and matrix thermodynamic dimensionless compressibility are also predicted. The presented results should be testable in future experiments and simulations.

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

Water under hydrophobic confinement: entropy and diffusion

The properties of liquid water are known to change drastically in confined geometries. A most interesting and intriguing phenomenon is that the diffusion of water is found to be strongly enhanced by the proximity of a hydrophobic confining wall relative to the bulk diffusion. We report a molecular dynamics simulation using a classical water model investigating the water diffusion near a non-interacting smooth confining wall, which is assumed to imitate a hydrophobic surface, revealing a pronounced diffusion enhancement within several water layers adjacent to the wall. We present evidence that the observed diffusion enhancement can be accounted for, with a quantitative accuracy, using the universal scaling law for liquid diffusion that relates the diffusion rate to the excess entropy. These results show that the scaling law, which has so far only been used for the description of the diffusion in simple liquids, can successfully describe the diffusion in water. It is shown that the law can be used for the analysis of water dynamics under nanoscale hydrophobic confinement, which is currently a subject of intense research activity.

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

Collisional and collisionless shocks

Shock waves are one of the most common plasma phenomena. They play a significant role in astrophysical and magnetospheric environments, as well as in laboratory plasmas. The paper focuses on two issues that underline important similarities and dissimilarities between collisional and collisionless shocks. The two types of shocks are similar in that they can be formed by the overtaking mechanism. We demonstrate that geometrical structures appearing early in the process of shock formation from the generic smooth initial state are quite universal and similar in both cases. Characterization of these structures (dubbed here the “mussel-shell”) is presented. Alongside with these similarities, there exist also significant differences between the two types of shocks. The classical collisional shocks usually connect two well-defined equilibrium states (those before and after the shock transition). Each of these equilibrium states is characterized by the thermodynamic parameters of density, temperature, and pressure whose upstream and downstream values are related by the continuity of mass, momentum, and energy flux. In the collisionless plasmas, however, this description does not work: the final state can be any of the much broader class of states constrained only by the requirement of being stable with respect to collisionless plasma instabilities. To find a final state (which does not have to be Maxwellian) one now has to follow the evolution of the system through the whole transition, and a lot of universality is lost. In some situations even the separation of scales between the global flow and shock transition may be lost. Taken together, these phenomena reveal a fascinating interplay of hydrodynamics, statistical mechanics, and plasma physics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Photo-z outlier self-calibration in weak lensing surveys

Calibrating photometric redshift errors in weak lensing surveys with external data is extremely challenging. In this work, we show that both Gaussian and outlier photo-z parameters can be self-calibrated from the data alone. This comes at no cost for the neutrino masses, curvature and dark energy equation of state w0, but with a 65% degradation when both w0 and wa are varied. We perform a realistic forecast for the Vera Rubin Observatory (VRO) Legacy Survey of Space and Time (LSST) 3× 2 analysis, combining cosmic shear, projected galaxy clustering and galaxy - galaxy lensing. We confirm the importance of marginalizing over photo-z outliers. We examine a subset of internal cross-correlations, dubbed "null correlations", which are usually ignored in 3× 2 analyses. Despite contributing only ~ 10% of the total signal-to-noise, these null correlations improve the constraints on photo-z parameters by up to an order of magnitude. Using the same galaxy sample as sources and lenses dramatically improves the photo-z uncertainties too. Together, these methods add robustness to any claim of detected new Physics, and reduce the statistical errors on cosmology by 15% and 10% respectively. Finally, including CMB lensing from an experiment like Simons Observatory or CMB-S4 improves the cosmological and photo-z posterior constraints by about 10%, and further improves the robustness to systematics. To give intuition on the Fisher forecasts, we examine in detail several toy models that explain the origin of the photo-z self-calibration. Our Fisher code LaSSI (Large-Scale Structure Information), which includes the effect of Gaussian and outlier photo-z, shear multiplicative bias, linear galaxy bias, and extensions to LCDM, is publicly available at \href{https://github.com/EmmanuelSchaan/LaSSI}{https://github.com/EmmanuelSchaan/LaSSI}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Modelling Z→TT processes in ATLAS with T-embedded Z →μμ data

We describe the concept, technical realisation and validation of a largely data-driven method to model events with Z→ττ decays. In Z→μμ events selected from proton-proton collision data recorded at √s=8 TeV with the ATLAS experiment at the LHC in 2012, the Z decay muons are replaced by τ leptons from simulated Z →ττ decays at the level of reconstructed tracks and calorimeter cells. The τ lepton kinematics are derived from the kinematics of the original muons. Thus, only the well-understood decays of the Z boson and τ leptons as well as the detector response to the τ decay products are obtained from simulation. All other aspects of the event, such as the Z boson and jet kinematics as well as effects from multiple interactions, are given by the actual data. This so-called τ-embedding method is particularly relevant for Higgs boson searches and analyses in ττ final states, where Z→ττ decays constitute a large irreducible background that cannot be obtained directly from data control samples. In this paper, we discuss the relevant concepts based on the implementation used in the ATLAS Standard Model H→ττ analysis of the full datataset recorded during 2011 and 2012.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Energy calibration of germanium detectors for the Majorana Demonstrator

The Majorana Demonstrator was a search for neutrinoless double-beta decay (0νββ) in the 76 Ge isotope. It was staged at the 4850-foot level of the Sanford Underground Research Facility (SURF) in Lead, SD. The experiment consisted of 58 germanium detectors housed in a low background shield and was calibrated once per week by deploying a 228 Th line source for 1 to 2 hours. The energy scale calibration determination for the detector array was automated using custom analysis tools. We describe the offline procedure for calibration of the Demonstrator germanium detectors, including the simultaneous fitting of multiple spectral peaks, estimation of energy scale uncertainties, and the automation of the calibration procedure.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Critical assessment of the recent report on the gigaparsec-scale correlation of the orientations of large quasar groups

ABSTRACT Recently, it was reported that large quasar groups (LQGs) identified from the Sloan Digital Sky Survey (SDSS) data release seven catalogue are not randomly oriented but preferentially aligned or orthogonal over scales 1–2 Gpc. To confirm this claim, I reproduced the same LQG sample and performed Sobolev tests of uniformity on the LQG orientation axes in the redshift space. Contrary to the original report based on the bimodal distribution of the LQG position angles in the sky, I found no departure from uniformity in the distribution of the LQG orientation axes. I also examined whether the LQGs are physical structures using a statistically more reliable data set constructed from the SDSS data release 16 (DR16) large-scale structure (LSS) quasar catalogue. Considering the Gaussian primordial density fluctuations and shot noise, I estimated the mass density contrasts of the LQGs from the number counts of the DR16 LSS quasars and found that most of the LQGs do not trace statistically significant high-density regions. I conclude that the LQG sample is a collection of unphysical chance associations and should not be used for any cosmological studies.

Fujii, Hirokazu (ORCID:0000000161473512)↗
DOE OSTI2022

The local wavenumber model for computation of turbulent mixing

We present an overview of the current status in the development of a two-point spectral closure model for turbulent flows, known as the local wavenumber (LWN) model. The model is envisioned as a practical option for applications requiring multi-physics simulations in which statistical hydrodynamics quantities such as Reynolds stresses, turbulent kinetic energy, and measures of mixing such as density-correlations and mix-width evolution, need to be captured with relatively high fidelity. In this review, we present the capabilities of the LWN model since it was first formulated in the early 1990s, for computations of increasing levels of complexity ranging from homogeneous isotropic turbulence, inhomogeneous and anisotropic single-fluid turbulence, to two-species mixing driven by buoyancy forces. The review concludes with a discussion of some of the more theoretical considerations that remain in the development of this model.

42 ENGINEERING↗
DOE OSTI2023

Theorem on extensive spectral degeneracy for systems with rigid higher symmetries in general dimensions

We establish, in the spirit of the Lieb-Schultz-Mattis theorem, lower bounds on the spectral degeneracy of quantum systems with higher (gaugelike) symmetries with rather generic physical boundary conditions in an arbitrary number of spatial dimensions. Contrary to applying twists or equivalent adiabatic operations, we exploit the effects of modified boundary conditions. We find that when a general choice of boundary geometry is immaterial in determining spectral degeneracies while approaching the thermodynamic limit, systems that exhibit rigid noncommuting gaugelike symmetries, such as the orbital compass model, must have an exponential (in the size of the boundary) degeneracy of each of their spectral levels.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Entropy and Seebeck signals meet on the edges

Here, we explore the electronic entropy per particle s and Seebeck coefficient S in zigzag graphene ribbons. Pristine and edge-doped ribbons are considered using tight-binding models to inspect the role of edge states in the observed thermal transport properties. As a band gap opens when the ribbons are doped at one or both edges, due to asymmetric edge potentials, we find that s and S signals are closely related to each other: both develop sharp dip-peak line shapes as the chemical potential lies in the gap, while the ratio s/S exhibits a near-constant value equal to the elementary charge e at low temperatures. This constant ratio suggests that S can be seen as the transport differential entropy per charge, as suggested by some authors. Our calculations also indicate that measurement of s and S may be useful as a spectroscopic probe of different electronic energy scales involved in such quantities in gapped materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Bipartite mutual information in classical many-body dynamics

Information theoretic measures have helped to sharpen our understanding of many-body quantum states. As perhaps the most well-known example, the entanglement entropy (or more generally, the bipartite mutual information) has become a powerful tool for characterizing the dynamical growth of quantum correlations. By contrast, although computable, the bipartite mutual information (MI) is almost never explored in classical many particle systems; this owes in part to the fact that computing the MI requires keeping track of the evolution of the full probability distribution, a feat which is rarely done (or thought to be needed) in classical many-body simulations. Here, we utilize the MI to analyze the spreading of information in 1D elementary cellular automata (CA). Broadly speaking, we find that the behavior of the MI in these dynamical systems exhibits a few different types of scaling that roughly correspond to known CA universality classes. Of particular note is that we observe a set of automata for which the MI converges parametrically slowly to its thermodynamic value. We develop a microscopic understanding of this behavior by analyzing a two-species model of annihilating particles moving in opposite directions. Furthermore, our work suggests the possibility that information theoretic tools such as the MI might enable a more fine-grained characterization of classical many-body states and dynamics.

Cellular automata↗
DOE OSTI2019

Reexamining the proton-radius problem using constrained Gaussian processes

The “proton radius puzzle” refers to an 8-year-old problem that highlights major inconsistencies in the extraction of the charge radius of the proton from muonic Lamb-shift experiments as compared against experiments using elastic electron scattering. For the latter approach, the determination of the charge radius involves an extrapolation of the experimental form factor to zero momentum transfer. To estimate the proton radius, a novel and powerful nonparametric method based on a constrained Gaussian process is introduced. Furthermore, the constrained Gaussian process models the electric form factor as a function of the momentum transfer. Within a Bayesian paradigm, we develop a model flexible enough to fit the data without any parametric assumptions on the form factor. The Bayesian estimation is guided by imposing only two physical constraints on the form factor: (a) its value at zero momentum transfer (normalization) and (b) its overall shape, assumed to be a monotonically decreasing function of the momentum transfer. Variants of these assumptions are explored to assess their impact.By adopting both constraints and incorporating the whole range of experimental data available we extracted a charge radius of r p = 0.845 ± 0.001 fm, consistent with the muonic experiment. Nevertheless, we show that within our model the extracted radius depends on both the assumed constraints and the range of experimental data used to fit the Gaussian process. For example, if only low-momentum-transfer data are used, relaxing the normalization constraint provides a value compatible with the larger electronic value. We have presented a novel technique to estimate the proton radius from electron-scattering data based on a constrained Gaussian process. We demonstrated that the impact of imposing sensible physical constraints on the form factor is substantial. Also critical is the range of the experimental data used in the extrapolation. We are hopeful that as the technique gets refined, together with the anticipated new results from the PRad experiment, we will get closer to a resolution of the puzzle.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Finite-temperature many-body perturbation theory in the canonical ensemble

Benchmark data are presented for the zeroth- through third-order many-body perturbation corrections to the electronic Helmholtz energy, internal energy, and entropy in the canonical ensemble in a wide range of temperature. They are determined as numerical λ-derivatives of the respective quantities computed by thermal full configuration interaction with a perturbation-scaled Hamiltonian, $\hat{H} = \hat{H}_0 + λ\hat{V}$. Sum-over-states analytical formulas for up to the third-order corrections to these properties are also derived as analytical λ-derivatives. These formulas, which are verified by exact numerical agreement with the benchmark data, are given in terms of the Hirschfelder–Certain degenerate perturbation energies and should be valid for both degenerate and nondegenerate reference states at any temperature down to zero. Further, the results in the canonical ensemble are compared with the same in the grand canonical ensemble.

74 ATOMIC AND MOLECULAR PHYSICS↗