Bayesian inverse uncertainty quantification of the physical model parameters for the spallation neutron source first target station
Not Available
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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.
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↗
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↗
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↗
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↗
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)↗
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↗
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↗
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↗
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↗
Critical point particle number fluctuations from molecular dynamics
Not Available
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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↗
Weak scaling of the contact distance between two fluctuating interfaces with system size
A pair of flat parallel surfaces, each freely diffusing along the direction of their separation, will eventually come into contact. If the shapes of these surfaces also fluctuate, then contact will occur when their centers-of-mass remain separated by a nonzero distance ℓ. An example of such a situation is the motion of interfaces between two phases at conditions of thermodynamic coexistence, and in particular the annihilation of domain wall pairs under periodic boundary conditions. Here we present a general approach to calculate the probability distribution of the contact distance ℓ and determine how its most likely value ℓ* depends on the surfaces' lateral size L. Using the Edward-Wilkinson equation as a model for interfaces, we demonstrate that ℓ* scales weakly with system size, i.e., the dependence of ℓ* on L for both (1+1)- and (2+1)-dimensional interfaces is such that lim L→∞ (ℓ*/L) = 0. In particular, for (2+1)-dimensional interfaces ℓ* is an algebraic function of logL, a result that is confirmed by computer simulations of slab-shaped domains formed under periodic boundary conditions. Overall, this weak scaling implies that such domains remain topologically intact until ℓ becomes very small compared to the lateral size of the interface, contradicting expectations from equilibrium thermodynamics.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Supervised and unsupervised machine learning of structural phases of polymers adsorbed to nanowires
Here, we identify configurational phases and structural transitions in a polymer nanotube composite by means of machine learning. We employ various unsupervised dimensionality reduction methods, conventional neural networks, as well as the confusion method, an unsupervised neural-network-based approach. We find neural networks are able to reliably recognize all configurational phases that have been found previously in experiment and simulation. Furthermore, we locate the boundaries between configurational phases in a way that removes human intuition or bias. This could be done before only by relying on preconceived, ad hoc order parameters.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Free-energy estimates from nonequilibrium trajectories under varying-temperature protocols
The Jarzynski equality allows the calculation of free-energy differences using values of work measured from nonequilibrium trajectories. The number of trajectories required to accurately estimate free-energy differences in this way grows sharply with the size of work fluctuations, motivating the search for protocols that perform desired transformations with minimum work. However, protocols of this nature can involve varying temperature, to which the Jarzynski equality does not apply. Here, we derive a variant of the Jarzynski equality that applies to varying-temperature protocols, and show that it can have better convergence properties than the standard version of the equality. We derive this modified equality and the associated fluctuation relation within the framework of Markovian stochastic dynamics, complementing related derivations done within the framework of Hamiltonian dynamics.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Universal energy-speed-accuracy trade-offs in driven nonequilibrium systems
The connection between measure theoretic optimal transport and dissipative nonequilibrium dynamics provides a language for quantifying nonequilibrium control costs, leading to a collection of thermodynamic speed limits, which rely on the assumption that the target probability distribution is perfectly realized. This is almost never the case in experiments or numerical simulations, so here we address the situation in which the external controller is imperfect. We obtain a lower bound for the dissipated work in generic nonequilibrium control problems that (1) is asymptotically tight and (2) matches the thermodynamic speed limit in the case of optimal driving. Along with analytically solvable examples, we refine this imperfect driving notion to systems in which the controlled degrees of freedom are slow relative to the nonequilibrium relaxation rate, and identify independent energy contributions from fast and slow degrees of freedom. Furthermore, we develop a strategy for optimizing minimally dissipative protocols based on optimal transport flow matching, a generative machine learning technique. Furthermore, this latter approach ensures the scalability of both the theoretical and computational framework we put forth. Crucially, we demonstrate that we can compute the terms in our bound numerically using efficient algorithms from the computational optimal transport literature and that the protocols we learn saturate the bound.
59 BASIC BIOLOGICAL SCIENCES↗
Unifying Quantum and Classical Speed Limits on Observables
Not Available
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Shortcuts in Stochastic Systems and Control of Biophysical Processes
Not Available
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗