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 55 records · Page 3

Time-domain theory of transient heat conduction in the local limit

Ultrafast and nanoscale heat conduction demands a unified theoretical framework that rigorously bridges macroscopic transport equations with microscopic material properties derived from statistical physics. Existing empirical generalizations of Fourier's law often lack a solid microscopic foundation, failing to connect observed non-Fourier behavior with underlying atomic-scale mechanisms. In this work, we present a time-domain theory of transient heat conduction rooted in Zwanzig's statistical theory of irreversible processes. Central to this framework is the time-domain transport function $\overleftrightarrow{𝑍}$⁡(𝑡) defined through equilibrium time-correlation functions of heat fluxes. This function generalizes the conventional concept of steady-state thermal conductivity, governing the transition of conduction dynamics from onset second sound type wave propagation at finite speeds to diffusion-dominated behavior across broad temporal and spatial scales. Unlike phonon hydrodynamic models that rely on mesoscopic constructs such as phonon drift velocity, our approach provides a quantitative and microscopic description of intrinsic memory effects in transient heat fluxes and applies universally to bulk materials at any temperature or length scale. By integrating atomistic-scale first-principles calculations with continuum-level macroscopic equations, this framework offers a robust foundation for numerical simulations of transient temperature fields. Furthermore, it facilitates the interpretation and design of transient thermal grating experiments using nanometer-scale heat sources and ultrafast laser systems in the extreme ultraviolet and x-ray wavelength ranges, advancing our understanding of heat dissipation dynamics.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Water vapor profile retrieval simulation studies for the HIRS/MSU and AMTS/MSU sounders

A combined physical statistical method for the retrieval of atmospheric water vapor column densities has been developed at the Goddard Laboratory for Atmospheric Sciences. In this method differences between observed and calculated clear column brightness temperatures are used in a regression formalism to obtain estimates of integrated column densities in given atmospheric layers. The physics of the radiative transfer problem is contained in the calculated clear column brightness temperatures and in accounting for the effects of clouds on the observed radiances. The regression matrix is equivalent to a set of empirical sensitivity factors.

Reuter, D.↗

ORION: Operational FoRecastIng Of INduced Seismicity

The Operational Forecasting of Induced Seismicity toolkit “ORION” (ORION) is an open-source, observation-based ensemble forecasting toolkit which is geared towards helping operators understand the seismic hazard (i.e., probabilistic assessment of the magnitude and frequency of induced seismic events) at a site. ORION analyzes how the seismic hazard evolves during injection and suggests possible mitigation strategies to employ if an earthquake that exceeds certain threshold is observed. Through its ensemble modeling approach, ORION leverages the benefits of statistical-, physics-, and machine learning-based forecasting methodologies, while reducing the impact of each model’s respective limitations. The ORION toolkit consists of an easy-to-use GUI interface that affords a user as much or as little interaction as desired. Advanced capabilities allow the user to upload local, high-precision earthquake catalogs, projected injection profiles and/or spatiotemporal estimates of pressure/stress, and to tune various model parameters. ORION will then provide a spatial and temporal ensemble forecast of seismicity defined as the probability of exceedance of a given earthquake magnitude over a forecast period. Additionally, ORION will provide probability distribution of the statistically derived maximum possible earthquake magnitude that may be expected. Finally, ORION will provide suggested operational management strategies (e.g. reduce injection volumes at specific wells) based on the level of hazard.

NRAP,NRAP Tools,ORION,Operational Forcasting of In↗

Nonlinearities in Magnetic Confinement, Ionospheric Physics, and Population Explosion Leading to Profile Resilience Нелінійності в магнетному утриманні, фізиці іоносфери та процесі демографічного вибуху, які приводять до стійкості профілю

Nonlinearities play an important role in many fields. In the field of thermonuclear fusion, they are involved in questions such as profile resilience and fluid closure. A nonlinear phenomenon common to both fusion and astrophysical planets is the generation of zonal flows. These flows play a significant role in determining the level of turbulence and fluid closure in fusion. The effects of resonance broadening and nonlinearities are investigated, specifically focusing on the case of nonlinear instability that has appeared in drift waves. Similarities and differences between our systems are discussed, with population explosion and the dynamics of nonlinear systems for drift waves by different states in profile resilience described with great precision. The aim of our study is to put our fluid model for drift waves in tokamaks within the wider framework of statistical physics principles. This reinforces our belief in the broad application of our drift wave model, which encompasses current tokamaks, ITER, and the fusion pilot plant.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quantum machine learning for chemistry and physics

Machine learning (ML) has emerged as a formidable force for identifying hidden but pertinent patterns within a given data set with the objective of subsequent generation of automated predictive behavior. In recent years, it is safe to conclude that ML and its close cousin, deep learning (DL), have ushered in unprecedented developments in all areas of physical sciences, especially chemistry. Not only classical variants of ML, even those trainable on near-term quantum hardwares have been developed with promising outcomes. Such algorithms have revolutionized materials design and performance of photovoltaics, electronic structure calculations of ground and excited states of correlated matter, computation of force-fields and potential energy surfaces informing chemical reaction dynamics, reactivity inspired rational strategies of drug designing and even classification of phases of matter with accurate identification of emergent criticality. In this review we shall explicate a subset of such topics and delineate the contributions made by both classical and quantum computing enhanced machine learning algorithms over the past few years. We shall not only present a brief overview of the well-known techniques but also highlight their learning strategies using statistical physical insight. The objective of the review is not only to foster exposition of the aforesaid techniques but also to empower and promote cross-pollination among future research in all areas of chemistry which can benefit from ML and in turn can potentially accelerate the growth of such algorithms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamics of nonequilibrium magnons in gapped Heisenberg antiferromagnets

Nonequilibrium dynamics in spin systems is a topic currently under intense investigation as it provides fundamental insights into thermalization, universality, and exotic transport phenomena. While most of the studies have been focused on ideal closed quantum many-body systems such as ultracold atomic quantum gases and one-dimensional spin chains, driven-dissipative Bose gases in steady states away from equilibrium in classical systems also lead to intriguing nonequilibrium physics. Here, in this work, we theoretically investigate out-of-equilibrium dynamics of magnons in a gapped Heisenberg quantum antiferromagnet based on Boltzmann transport theory. We show that, by treating scattering terms beyond the relaxation-time approximation in the Boltzmann transport equation, energy and particle number conservation mandate that nonequilibrium magnons cannot relax to equilibrium, but decay to other nonequilibrium stationary states. The only decay channel for these stationary states back to equilibrium is through the nonconserving interactions (i.e., changing particle number and/or energy within the magnon system) such as boundary or magnon-phonon scattering. At low temperatures, these nonconserving interactions are much slower processes than intrinsic magnon-magnon interaction in a gapped spin system. Using magnon-phonon interaction as a quintessential type of nonconserving interaction, we then propose that nonequilibrium steady states of magnons can be maintained and tailored using periodic driving at frequencies faster than relaxation due to phonon interactions. These findings reveal a class of classical material systems that are suitable platforms to study nonequilibrium statistical physics and macroscopic phenomena such as classical Bose-Einstein condensation of quasiparticles and magnon supercurrents that are relevant for spintronic applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Ultrafast radiographic imaging and tracking: An overview of instruments, methods, data, and applications

Ultrafast radiographic imaging and tracking (U-RadIT) use state-of-the-art ionizing particle and light sources to experimentally study sub-nanosecond transients or dynamic processes in physics, chemistry, biology, geology, materials science and other fields. These processes are fundamental to modern technologies and applications, such as nuclear fusion energy, advanced manufacturing, communication, and green transportation, which often involve one mole or more atoms and elementary particles, and thus are challenging to compute by using the first principles of quantum physics or other forward models. One of the central problems in U-RadIT is to optimize information yield through, e.g. high-luminosity X-ray and particle sources, efficient imaging and tracking detectors, novel methods to collect data, and large-bandwidth online and offline data processing, regulated by the underlying physics, statistics, and computing power. We review and highlight recent progress in: (a.) Detectors such as high-speed complementary metal-oxide semiconductor (CMOS) cameras, hybrid pixelated array detectors integrated with Timepix4 and other application-specific integrated circuits (ASICs), and digital photon detectors; (b.) U-RadIT modalities such as dynamic phase contrast imaging, dynamic diffractive imaging, and four-dimensional (4D) particle tracking; (c.) U-RadIT data and algorithms such as neural networks and machine learning, and (d.) Applications in ultrafast dynamic material science using XFELs, synchrotrons and laser-driven sources. Hardware-centric approaches to U-RadIT optimization are constrained by detector material properties, low signal-to-noise ratio, high cost and long development cycles of critical hardware components such as ASICs. Interpretation of experimental data, including comparisons with forward models, is frequently hindered by sparse measurements, model and measurement uncertainties, and noise. Alternatively, U-RadIT make increasing use of data science and machine learning algorithms, including experimental implementations of compressed sensing. Machine learning and artificial intelligence approaches, refined by physics and materials information, may also contribute significantly to data interpretation, uncertainty quantification and U-RadIT optimization.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Sequence2Self: Self-supervised image sequence denoising of pixel-level spray breakup morphology

Optical imaging of fast and transient phenomena such as the turbulent breakup of liquid sprays exhibit low signal-to-noise ratios due to the limited illumination intensity relative to the short exposure time. Image denoising is required to facilitate physical studies over these data but is challenging due to the absence of clean ground-truths and the stringency of the denoising task (e.g., strong and complex noise, limited resolution, preserving physical fidelity), preventing supervised and existing un-/self-supervised deep learning methods. To this end, Sequence2Self (Seq2S) is proposed, an extension of Self2Self (S2S) to image sequences that leverages both the signal’s spatial and temporal correlation. Seq2S is demonstrated on time-resolved x-ray phase contrast imaging of liquid jet fuel sprays in a gas turbine combustor, which possesses all of challenges detailed above. Experiments are conducted across four fuels with different breakup morphology using various state-of-the-art methods. Overall, many of the methods failed and Seq2S was most successful: (1) Accurate spray structures were reconstructed with consistent evolution across frames void of artifacts. (2) The performance was robust, invariant to the hyperparameter choice. (3) Computational time is short and can be made eligible for real-time denoising. In particular, the images denoised by Seq2S showed spray droplet diameter distributions with near-zero Kullback–Leibler divergence (0.01 ± 0.01) to a cleaner reference, whereas the second best method yielded 0.06 ± 0.03. In conclusion, this suggests that Seq2S can be reliably used prior to subsequent quantitative spray analyses as it retains (if not, improves) the statistical physical properties of the data.

97 MATHEMATICS AND COMPUTING↗

Distributions in CFT. Part II. Minkowski space

CFTs in Euclidean signature satisfy well-accepted rules, such as the convergent Euclidean OPE. It is nowadays common to assume that CFT correlators exist and have various properties also in Lorentzian signature. Some of these properties may represent extra assumptions, and it is an open question if they hold for familiar statistical-physics CFTs such as the critical 3d Ising model. Here we consider Wightman 4-point functions of scalar primaries in Lorentzian signature. We derive a minimal set of their properties solely from the Euclidean unitary CFT axioms, without using extra assumptions. We establish all Wightman axioms (temperedness, spectral property, local commutativity, clustering), Lorentzian conformal invariance, and distributional convergence of the s-channel Lorentzian OPE. This is done constructively, by analytically continuing the 4-point functions using the s-channel OPE expansion in the radial cross-ratios ρ, ρ¯. We prove a key fact that |ρ|, |ρ¯| < 1 inside the forward tube, and set bounds on how fast |ρ|, |ρ¯| may tend to 1 when approaching the Minkowski space.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The endpoint of partial deconfinement

We study the matrix quantum mechanics of two free hermitian N × N matrices subject to a singlet constraint in the microcanonical ensemble. This is the simplest example of a theory that at large N has a confinement/deconfinement transition. In the microcanonical ensemble, it also exhibits partial deconfinement with a Hagedorn density of states. We argue that the entropy of these configurations, based on a combinatorial counting of Young diagrams, are dominated by Young diagrams that have the VKLS shape. When the shape gets to the maximal depth allowed for a Young diagram of SU(N), namely N, we argue that the system stops exhibiting the Hagedorn behavior. The number of boxes (energy) at the transition is N 2 /4, independent of the charge of the state.

1/N Expansion↗

Entanglement Cost for Infinite-Dimensional Physical Systems

We prove that the entanglement cost equals the regularized entanglement of formation for any infinite-dimensional quantum state ρ ΑΒ with finite quantum entropy on at least one of the subsystems A or B. This generalizes a foundational result in quantum information theory that was previously formulated only for operations and states on finite-dimensional systems. The extension to infinite-dimensional systems is nontrivial because the conventional tools for establishing both the direct and converse bounds, i.e., strong typicality, monotonicity, and asymptotic continuity, are no longer directly applicable. To address this problem, we construct a new entanglement dilution protocol for infinite-dimensional states implementable by local operations and a finite amount of one-way classical communication (one-way LOCC), using weak and strong typicality multiple times. We also prove the optimality of this protocol among all protocols, even under infinite-dimensional separable operations, by developing an argument based on alternative forms of monotonicity and asymptotic continuity of the entanglement of formation for infinite-dimensional states. Along the way, we derive a new integral representation for the quantum entropy of infinite-dimensional states, which we believe to be of independent interest. Our results allow us to fully characterize an important operational entanglement measure—the entanglement cost—for all infinite-dimensional physical systems.

Complexity↗

Record acceleration of the two-dimensional Ising model using a high-performance wafer-scale engine

The versatility and wide-ranging applicability of the Ising model, originally introduced to study phase transitions in magnetic materials, have made it a cornerstone in statistical physics and a valuable tool for evaluating the performance of emerging computer hardware. Here, we present a novel implementation of the two-dimensional Ising model on Cerebras Wafer-Scale Engine (WSE) – a revolutionary processor that is opening new frontiers in computing. In our deployment of the checkerboard algorithm, we optimized the Ising model to take advantage of the unique WSE architecture. Specifically, we employed a compressed bit representation storing 16 spins on each int16 word, and efficiently distributed the spins over the processing units enabling seamless weak scaling and limiting communications to only immediate neighboring units. Our implementation can handle up to 754 simulations in parallel, achieving an aggregate of over 61.8 trillion flip attempts per second for Ising models with up to 200 million spins. This represents a gain of up to 148 times over previously reported single-devices with a highly optimized implementation on NVIDIA V100 and up to 88 times in productivity compared to NVIDIA H100. Our findings highlight the significant potential of the WSE in scientific computing, particularly in the field of materials modeling.

Ising model↗

Electrons matter in ferromagnetic criticality

When elemental iron is cooled from high temperatures, magnetic moments carried by its d-electrons begin a continuous transition to a state of parallel aligned moments below its Curie temperature TC of about 1040 K. As established in the mid-1900s, the statistical physics of this second order phase transition is universal, independent of the ordering temperature, material structure and chemical composition as well as the character of moment-bearing electrons. Finally, the thermodynamics of the phase transition is determined solely by symmetries of the transition’s order parameter – magnetization in the case of a ferromagnet.

36 MATERIALS SCIENCE↗

Feasibility of a multigroup Boltzmann–Fokker–Planck solution for electron beam dose calculations

Legacy nuclear-reactor Boltzmann solvers start clinical deployment as an alternative to Monte Carlo (MC) codes and Fermi–Eyges semiemprical models in radiation oncology treatment planning. Today’s certified clinical solvers are limited to photon beams. In this paper, ELECTR, a state-of-the-art multigroup electron cross sections generation module in NJOY is presented and validated against Lockwood’s calorimetric measurements, EGS-nrc and GEANT-4 for 1–20 MeV unidirectional electron beams. The nuclear-reactor DRAGON-5 solver is upgraded to access the library and solve the Boltzmann–Fokker–Planck (BFP) equation. A variety of heterogeneous radiotherapy and radiosurgery phantom configurations were used for validation purpose. Case studies include a thorax benchmark, that of a typical breast Intra-Operative Radiotherapy and a high-heterogeneity patient-like benchmark. For all beams, 100% of the water voxels satisfied the American Association of Physicists in Medicine accuracy criterion for a BFP-MC dose error below 2%. At least, 97.0% of adipose, muscle, bone, lung, tumor and breast voxels satisfied the 2% criterion. The average BFP-MC relative error was about 0.56% for all voxels, beams and materials combined. By irradiating homogeneous slabs from Z = 1 (hydrogen) to Z = 99 (einsteinium), we reported performance and defects of the CEPXS mode [US. Sandia National Lab., SAND-89-1685] in ELECTR for the entire periodic table. For all Lockwood’s benchmarks, NJOY-DRAGON dose predictions are within the experimental data precision for 98% of voxels.

42 ENGINEERING↗

Robust finite-temperature many-body scarring on a quantum computer

Mechanisms for suppressing thermalization in disorder-free many-body systems, such as Hilbert space fragmentation and quantum many-body scars, have recently attracted much interest in foundations of quantum statistical physics and potential quantum information processing applications. However, their sensitivity to realistic effects such as finite temperature remains largely unexplored. Here, we have utilized IBM's Kolkata quantum processor to demonstrate an unexpected robustness of quantum many-body scars at finite temperatures when the system is prepared in a thermal Gibbs ensemble. We identify such robustness in the PXP model, which describes quantum many-body scars in experimental systems of Rydberg atom arrays and ultracold atoms in tilted Bose-Hubbard optical lattices. By contrast, other theoretical models which host exact quantum many-body scars are found to lack such robustness and their scarring properties quickly decay with temperature. Our study sheds light on the important differences between scarred models in terms of their algebraic structures, which impacts their resilience to finite temperature. Published by the American Physical Society 2024

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Quantum many-body simulations of the two-dimensional Fermi-Hubbard model in ultracold optical lattices

Understanding quantum many-body states of correlated electrons is one main theme in modern condensed-matter physics. Given that the Fermi-Hubbard model, the prototype of correlated electrons, was recently realized in ultracold optical lattices, it is highly desirable to have controlled numerical methodology to provide precise finite-temperature results upon doping to directly compare with experiments. Here, we demonstrate the exponential tensor renormalization group (XTRG) algorithm [Chen et al., Phys. Rev. X 8, 031082 (2018)], complemented by independent determinant quantum Monte Carlo, offers a powerful combination of tools for this purpose. XTRG provides full and accurate access to the density matrix and thus various spin and charge correlations, down to an unprecedented low temperature of a few percent of the tunneling energy. Finally, we observe excellent agreement with ultracold fermion measurements at both half filling and finite doping, including the sign-reversal behavior in spin correlations due to formation of magnetic polarons, and the attractive hole-doublon and repulsive hole-hole pairs that are responsible for the peculiar bunching and antibunching behaviors of the antimoments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Numerical fluid dynamics for FRG flow equations: Zero-dimensional QFTs as numerical test cases. I. The O ( N ) model

The functional renormalization group (FRG) approach is a powerful tool for studies of a large variety of systems, ranging from statistical physics over the theory of the strong interaction to gravity. The practical application of this approach relies on the derivation of so-called flow equations, which describe the change of the quantum effective action under the variation of a coarse-graining parameter. In the present work, we discuss in detail a novel approach to solve such flow equations. This approach relies on the fact that RG equations can be rewritten such that they exhibit similarities with the conservation laws of fluid dynamics. This observation can be exploited in different ways. First of all, we show that this allows to employ powerful numerical techniques developed in the context of fluid dynamics to solve RG equations. In particular, it allows us to reliably treat the emergence of nonanalytic behavior in the RG flow of the effective action as it is expected to occur in studies of, e.g., spontaneous symmetry breaking. Second, the analogy between RG equations and fluid dynamics offers the opportunity to gain novel insights into RG flows and their interpretation in general, including the irreversibility of RG flows. Further, we work out this connection in practice by applying it to zero-dimensional quantum-field theoretical models. The generalization to higher-dimensional models is also discussed. Our findings are expected to help improving future FRG studies of quantum field theories in higher dimensions both on a qualitative and quantitative level.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Observing Quantum Measurement Collapse as a Learnability Phase Transition

During a quantum measurement, superpositions of states with different observable properties probabilistically collapse into one with a sharp value of the measured observable. In macroscopic quantum systems, this collapse arises via a continuous measurement-induced phase transition (MIPT) at a critical value of the strength of interaction with the measurement apparatus. MIPTs lie outside established paradigms for equilibrium or nonequilibrium critical phenomena and delineate distinct, stable dynamical and computational phases of matter. Quantum computers enable programmable simulation of the interaction of a measurement apparatus with a dynamical quantum system, to explore MIPT phenomena over a range of system sizes while retaining quantum coherence. Yet, existing experimental protocols rely on fundamentally nonscalable postselection techniques or direct classical simulation of quantum circuits. Here, we report the scalable observation of finite-size scaling evidence for an observable-sharpening MIPT in monitored quantum circuits in a chain of Yb + 171 ions in Quantinuum’s H1-1 trapped-ion quantum processor. By leveraging an equivalent description as a statistical physics problem, we implement scalable classical algorithms to infer the value of the measured observable from a single experimental shot. This technique enables a truly scalable protocol to observe observable-sharpening MIPTs in generic classes of circuits that cannot be directly classically simulated and also provides enhanced means to detect and suppress errors in the quantum simulation. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗