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

Hydrodynamic equations for electrons in graphene obtained from the maximum entropy principle

The maximum entropy principle is applied to the formal derivation of isothermal, Euler-like equations for semiclassical fermions (electrons and holes) in graphene. After proving general mathematical properties of the equations so obtained, their asymptotic form corresponding to significant physical regimes is investigated. In particular, the diffusive regime, the Maxwell-Boltzmann regime (high temperature), the collimation regime and the degenerate gas limit (vanishing temperature) are considered.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Multi-temperature model derived from state-to-state kinetics for hypersonic entry in Jupiter atmosphere

A state-to-state model of H{sub 2}/He plasmas coupling the master equations for internal distributions of heavy species with the transport equation for the free electrons has been used as a basis for implementing a multi-temperature kinetic model. In the multi-temperature model internal distributions of heavy particles are Boltzmann, the electron energy distribution function is Maxwell, and the rate coefficients of the elementary processes become a function of local temperatures associated to the relevant equilibrium distributions. The state-to-state and multi-temperature models have been compared in the case of a homogenous recombining plasma, reproducing the conditions met during supersonic expansion though converging-diverging nozzles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Shock response of He bubbles in single crystal Cu

With large-scale molecular dynamics simulations, we investigate shock response of He nanobubbles in single crystal Cu. For sufficient bubble size or internal pressure, a prismatic dislocation loop may form around a bubble in unshocked Cu. The internal He pressure helps to stabilize the bubble against plastic deformation. However, the prismatic dislocation loops may partially heal but facilitate nucleation of new shear and prismatic dislocation loops. For strong shocks, the internal pressure also impedes internal jetting, while a bubble assists local melting; a high speed jet breaks a He bubble into pieces dispersed among Cu. Near-surface He bubbles may burst and form high velocity ejecta containing atoms and small fragments, while the ejecta velocities do not follow the three-dimensional Maxwell-Boltzmann distributions expected for thermal equilibrium. The biggest fragment size deceases with increasing shock strength. With a decrease in ligament thickness or an increase in He bubble size, the critical shock strength required for bubble bursting decreases, while the velocity range, space extension and average velocity component along the shock direction, increase. Small bubbles are more efficient in mass ejecting. Compared to voids and perfect single crystal Cu, He bubbles have pronounced effects on shock response including bubble/void collapse, Hugoniot elastic limit (HEL), deformation mechanisms, and surface jetting. HEL is the highest for perfect single crystal Cu with the same orientations, followed by He bubbles without pre-existing prismatic dislocation loops, and then voids. Complete void collapse and shear dislocations occur for embedded voids, as opposed to partial collapse, and shear and possibly prismatic dislocations for He bubbles. He bubbles lower the threshhold shock strength for ejecta formation, and increase ejecta velocity and ejected mass.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Boltzmann-conserving classical dynamics in quantum time-correlation functions: “Matsubara dynamics”

We show that a single change in the derivation of the linearized semiclassical-initial value representation (LSC-IVR or “classical Wigner approximation”) results in a classical dynamics which conserves the quantum Boltzmann distribution. We rederive the (standard) LSC-IVR approach by writing the (exact) quantum time-correlation function in terms of the normal modes of a free ring-polymer (i.e., a discrete imaginary-time Feynman path), taking the limit that the number of polymer beads N → ∞, such that the lowest normal-mode frequencies take their “Matsubara” values. The change we propose is to truncate the quantum Liouvillian, not explicitly in powers of ħ{sup 2} at ħ{sup 0} (which gives back the standard LSC-IVR approximation), but in the normal-mode derivatives corresponding to the lowest Matsubara frequencies. The resulting “Matsubara” dynamics is inherently classical (since all terms O(ħ{sup 2}) disappear from the Matsubara Liouvillian in the limit N → ∞) and conserves the quantum Boltzmann distribution because the Matsubara Hamiltonian is symmetric with respect to imaginary-time translation. Numerical tests show that the Matsubara approximation to the quantum time-correlation function converges with respect to the number of modes and gives better agreement than LSC-IVR with the exact quantum result. Matsubara dynamics is too computationally expensive to be applied to complex systems, but its further approximation may lead to practical methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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↗