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 73 records · Page 4

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↗

Observation of a prethermal discrete time crystal

Extending the framework of statistical physics to the nonequilibrium setting has led to the discovery of previously unidentified phases of matter, often catalyzed by periodic driving. However, preventing the runaway heating that is associated with driving a strongly interacting quantum system remains a challenge in the investigation of these newly discovered phases. In this work, we utilize a trapped-ion quantum simulator to observe the signatures of a nonequilibrium driven phase without disorder—the prethermal discrete time crystal. Here, the heating problem is circumvented not by disorder-induced many-body localization, but rather by high-frequency driving, which leads to an expansive time window where nonequilibrium phases can emerge. Floquet prethermalization is thus presented as a general strategy for creating, stabilizing, and studying intrinsically out-of-equilibrium phases of matter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Improving the Parameterization of Cloud and Rain Microphysics in E3SM using Novel Observationally-Constrained Bayesian Approach (Final Technical Report)

In this project, we sought to develop new cloud and rain microphysics frameworks within the Energy Exascale Earth System Model (E3SM). This work encompassed two primary avenues of research: 1) Further development of a Bayesian-based scheme called BOSS (Bayesian Observationally-constrained Statistical-physical Scheme) to represent cloud and rain microphysics, testing it in realistic high-resolution cloud models, and implementing it in E3SM; 2) Development of a methodology utilizing machine learning to enable computationally tractable use of tractable use of Markov chain Monte Carlo sampling for Bayesian parameter estimation in Earth system and cloud models. In this project, we adapted the BOSS microphysics scheme, originally formulated for rain-only, to include all liquid-phase microphysical processes for cloud and rain, in particular the processes that mediate between these two categories, for example the conversion from cloud to rain through collision and coalescence of drops. We constrained the scheme via comparison and testing against a detailed model that explicitly represents the evolution of cloud and rain particles, called a bin microphysics scheme.

54 ENVIRONMENTAL SCIENCES↗

Comment on “Review of experimental studies of secondary ice production” by Korolev and Leisner (2020)

Abstract. This is a comment on the review by Korolev and Leisner (2020, hereafter KL2020). The only two laboratory/field studies ever to measure the breakup in ice–ice collisions for in-cloud conditions were negatively criticised by KL2020, as were our subsequent theoretical and modelling studies informed by both studies. First, hypothetically, even without any further laboratory experiments, such theoretical and modelling studies would continue to be possible, based on classical mechanics and statistical physics. They are not sensitive to the accuracy of lab data for typical situations, partly because the nonlinear explosive growth of ice concentrations continues until some maximum concentration is reached. To a degree, the same final concentration is expected regardless of the fragment number per collision. Second, there is no evidence that both lab/field observational studies characterising fragmentation in ice–ice collisions are either mutually conflicting or erroneous such that they cannot be used to represent this breakup in numerical models, contrary to the review. The fact that the ice spheres of one experiment were hail sized (2 cm) is not a problem if a universal theoretical formulation, such as ours, with fundamental dependencies, is informed by it. Although both lab/field studies involved head-on collisions, rotational kinetic energy for all collisions generally is only a small fraction of the initial collision kinetic energy (CKE) anyway. Although both lab/field experiments involved fixed targets, that is not a problem since the fixing of the target is represented via CKE in any energy-based formulation such as ours. Finally, scaling analysis suggests that the breakup of ice during sublimation can make a significant contribution to ice enhancement in clouds, again contrary to the impression given by the review.

54 ENVIRONMENTAL SCIENCES↗

Stratospheric warmings in the Southern Hemisphere deduced from satellite radiation data, 1969-73

The investigation is largely based on radiation data from the vertical temperature profile radiometers on the satellite NOAA 2. The physical-statistical interpretation of the radiance data is considered. Aspects regarding the nonuniqueness of the minor warming in August 1973 are discussed. It is pointed out that conditions leading to major warmings in the Northern Hemisphere were not satisfied in the Southern Hemisphere during the period from 1969 to 1973.

Quiroz, R. S.↗

A remark about pointed bubbles

The polymer expansion is a formal algebraic identity between a partition function and logarithm in statistical physics problems. The expansion gives a systematic method to control the free energy or to establish exponential tree-graph decay of connected correlations. Here, the convergence properties of the polymer expansion are analyzed in connection with three practical examples, including: intersecting bonds in chemical polymer chains; a connected closed hypersurface built from the (d-1)-faces of the d-dimensional unit cubes; and the set of Feynamn diagrams in the perturbation series of the Euclidean field theory partition function Z. The example of connected polymer chains is generalized to apply to other lattice models, including n-state Ising models at high temperature; short range lattice gases at high temperature; and weak coupling lattice field and gauge theories.

Garabedian, P. R.↗

Photo-thermo-dissociation. I - A general mechanism for destroying molecules

Photothermodissociation (PTD) of an isolated molecule is the name of the process when a photon is absorbed by an electronic transition, the energy is transferred to the vibrations (hot molecule) and then an atom or a fragment is ejected. A statistical physics treatment is presented that goes beyond the thermal approximation and fully takes into account the isolation of the system. A model is given that permits an explicit calculation of the ejection rates of hydrogen and carbon fragments from aromatic molecules as a function of their size and internal energy. Then, the lifetime of species against PTD in different astronomical radiation fields will result when a precise dynamical evolution of their internal energy can be established.

Leger, A.↗

Confidence bounds on structural reliability

Different approaches for quantifying physical, statistical, and model uncertainties associated with the distribution parameters which are aimed at determining structural reliability are described. Confidence intervals on the distribution parameters of the input random variables are estimated using four algorithms to evaluate uncertainty of the response. Design intervals are evaluated using either Monte Carlo simulation or an iterative approach. A first order approach can be used to compute a first approximation of the design interval, but its accuracy is not satisfactory. The regression approach which combines the iterative approach with Monte Carlo simulation is capable of providing good results if the performance function can be accurately represented using regression analysis. It is concluded that the design interval-based approach seems to be quite general and takes into account distribution and model uncertainties.

Mehta, S. R.↗

Dynamics of Granular Materials

Granular materials exhibit a rich variety of dynamical behavior, much of which is poorly understood. Fractal-like stress chains, convection, a variety of wave dynamics, including waves which resemble capillary waves, l/f noise, and fractional Brownian motion provide examples. Work beginning at Duke will focus on gravity driven convection, mixing and gravitational collapse. Although granular materials consist of collections of interacting particles, there are important differences between the dynamics of a collections of grains and the dynamics of a collections of molecules. In particular, the ergodic hypothesis is generally invalid for granular materials, so that ordinary statistical physics does not apply. In the absence of a steady energy input, granular materials undergo a rapid collapse which is strongly influenced by the presence of gravity. Fluctuations on laboratory scales in such quantities as the stress can be very large-as much as an order of magnitude greater than the mean.

Behringer, Robert P.↗

Contrail-Cirrus Studies at FARS

The three primary goals have been to: (1) increase our contrail remote sensing database collected from the Facility for Atmospheric Remote Sensing (FARS) by expanding our Project FIRE Extended Time Observations (ETO) program to include contrail persistence studies: (2) continue our retrospective analyses of the statistical, physical, and radiative properties of contrails derived from our 10-year ETO cirrus cloud dataset to examine their potential direct impact for regional climate change: and, (3) prepare our mobile remote sensing systems and participate in the SUCCESS field campaign from the DOE Southern Great Plains ARM CART site.

Sassen, Kenneth↗

Multi-Satellite Estimates of Land-Surface Properties for Determination of Energy and Water Budgets

Using the WETNET database, existing methods for the estimation of surface wetness from SSM/I data have been assessed and further developed. A physical-statistical method for optimal estimation of daily surface heat flux and Bowen ratio on the mesoscale has been developed and tested. This method is based on observations of daytime planetary boundary layer (PBL) growth from operational ravansonde and daytime land-surface temperature amplitude from Geostationary Operational Environmental (GOES) satellites. The mesoscale patterns of these heat fluxes have been compared with an AVHRR-based vegetation index and surface wetness (separately estimated from SSM/I and in situ observations). Cases of the 1988 Midwest drought and a surface/atmosphere moisture gradient (dry-line) in the southern Plains were studied. The analyses revealed significant variations in sensible heat flux (S(sub 0), and Bowen ratio, B(sub 0)) associated with vegetation cover and antecedent precipitation. Relationships for surface heat flux (and Bowen ratio) from antecedent precipitation and vegetation index have been developed and compared to other findings. Results from this project are reported in the following reviewed literature.

Menzel, W. Paul↗

Engineering of Novel Biocolloid Suspensions

Colloidal suspensions are materials with a variety of uses from cleaners and lubricants to food, cosmetics, and coatings. In addition, they can be used as a tool for testing the fundamental tenets of statistical physics. Colloidal suspensions can be synthesized from a wide variety of materials, and in the form of monodisperse particles, which can self-assemble into highly ordered colloidal crystal structures. As such they can also be used as templates for the construction of highly ordered materials. Materials design of colloids has, to date, relied on entropic self-assembly, where crystals form as result of lower free energy due to a transition to order. Here, our goal is to develop a completely new method for materials fabrication using colloidal precursors, in which the self-assembly of the ordered colloidal structures is driven by a highly controllable, attractive interaction. This will greatly increase the range of potential structures that can be fabricated with colloidal particles. In this work, we demonstrate that colloidal suspensions can be crosslinked through highly specific biological crosslinking reactions. In particular, the molecules we use are protein-carbohydrate interactions derived from the immune system. This different driving force for self-assembly will yield different and novel suspensions structures. Because the biological interactions are heterotypic (A binding to B), this chemical system can be used to make binary alloys in which the two colloid subpopulations vary in some property - size, density, volume fraction, magnetic susceptibility, etc. An additional feature of these molecules which is unique - even within the realm of biological recognition - is that the molecules bind reversibly on reasonable time-scales, which will enable the suspension to sample different configurations, and allow us to manipulate and measure the size of the suspension dynamically. Because of the wide variety of structures that can be made from these novel colloids, and because the suspension structure can be altered dynamically, we believe this biocolloid system will yield a novel set of materials with many technological applications, including sensors (both biological and non-biological), optical filters and separation media.

Hammer, D. A.↗

Ising model of cardiac thin filament activation with nearest-neighbor cooperative interactions

We have developed a model of cardiac thin filament activation using an Ising model approach from equilibrium statistical physics. This model explicitly represents nearest-neighbor interactions between 26 troponin/tropomyosin units along a one-dimensional array that represents the cardiac thin filament. With transition rates chosen to match experimental data, the results show that the resulting force-pCa (F-pCa) relations are similar to Hill functions with asymmetries, as seen in experimental data. Specifically, Hill plots showing (log(F/(1-F)) vs. log [Ca]) reveal a steeper slope below the half activation point (Ca(50)) compared with above. Parameter variation studies show interplay of parameters that affect the apparent cooperativity and asymmetry in the F-pCa relations. The model also predicts that Ca binding is uncooperative for low [Ca], becomes steeper near Ca(50), and becomes uncooperative again at higher [Ca]. The steepness near Ca(50) mirrors the steep F-pCa as a result of thermodynamic considerations. The model also predicts that the correlation between troponin/tropomyosin units along the one-dimensional array quickly decays at high and low [Ca], but near Ca(50), high correlation occurs across the whole array. This work provides a simple model that can account for the steepness and shape of F-pCa relations that other models fail to reproduce.

Evaluation Studies↗

Fractal mechanisms and heart rate dynamics. Long-range correlations and their breakdown with disease

Under healthy conditions, the normal cardiac (sinus) interbeat interval fluctuates in a complex manner. Quantitative analysis using techniques adapted from statistical physics reveals the presence of long-range power-law correlations extending over thousands of heartbeats. This scale-invariant (fractal) behavior suggests that the regulatory system generating these fluctuations is operating far from equilibrium. In contrast, it is found that for subjects at high risk of sudden death (e.g., congestive heart failure patients), these long-range correlations break down. Application of fractal scaling analysis and related techniques provides new approaches to assessing cardiac risk and forecasting sudden cardiac death, as well as motivating development of novel physiologic models of systems that appear to be heterodynamic rather than homeostatic.

NASA Discipline Cardiopulmonary↗

A Comparative Study of Probability Collectives Based Multi-agent Systems and Genetic Algorithms

We compare Genetic Algorithms (GA's) with Probability Collectives (PC), a new framework for distributed optimization and control. In contrast to GA's, PC-based methods do not update populations of solutions. Instead they update an explicitly parameterized probability distribution p over the space of solutions. That updating of p arises as the optimization of a functional of p. The functional is chosen so that any p that optimizes it should be p peaked about good solutions. The PC approach works in both continuous and discrete problems. It does not suffer from the resolution limitation of the finite bit length encoding of parameters into GA alleles. It also has deep connections with both game theory and statistical physics. We review the PC approach using its motivation as the information theoretic formulation of bounded rationality for multi-agent systems. It is then compared with GA's on a diverse set of problems. To handle high dimensional surfaces, in the PC method investigated here p is restricted to a product distribution. Each distribution in that product is controlled by a separate agent. The test functions were selected for their difficulty using either traditional gradient descent or genetic algorithms. On those functions the PC-based approach significantly outperforms traditional GA's in both rate of descent, trapping in false minima, and long term optimization.

Huang, Chien-Feng↗

The X-43A Flight Research Program: Lessons Learned on the Road to Mach 10

During an aerospace engineer's undergraduate studies, he or she will attend classes in aerodynamics, thermodynamics, structures, stability and control, dynamics, design, propulsion, and computer science, along with the related courses in mathematics, physics, statistics, and chemistry required to understand the material. Upon graduation, the new engineer will have acquired a basic knowledge of how to build an aerospace vehicle. What only comes through experience, however, is the understanding of the inevitable imperfect process through which an aerospace vehicle is built. This is the adventure of turning a basic concept into functional hardware. Engineers working on a project must often deal with ambiguous situations. They are routinely asked by management to provide risk assessments of a project, yet even after careful analysis uncertainties remain. The project must be accomplished within finite limits of time and money. The question an engineer answers is whether the solution to potential problem is worth the cost and schedule delay, or if the solution might actually be worse than the problem it is meant to solve. Review protocols are established to ensure that an unknown has not been overlooked. But these cannot protect against an unknown unknown. Examples of these situations can be found in the history of the X-43A Hyper-X (Hypersonic Experiment) program. In this NASA project, a supersonic combustion ramjet (scramjet) engine was flight tested on a subscale vehicle. The X-43A Hyper-X Research Vehicle (HXRV) was launched from a B-52B mothership, then boosted to the test speed by a modified Pegasus rocket first stage, called the Hyper-X Launch Vehicle (HXLV). Once at the proper speed and altitude, the X-43A separated from the booster, stabilized itself, and then the engine test began. Although wind-tunnel scramjet engine tests had begun in the late 1950s, before the Hyper-X program there had never been an actual in-flight test of such an engine integrated with an appropriate airframe. Thus, while the scramjet had successfully operated in the artificial airflow of wind tunnels, the concept had yet to be proven in real air. These conditions meant changes in density and temperature, as well as changes in angle of attack and sideslip of a free-flying vehicle. A wind tunnel is limited in its ability to simulate these subtle factures, which have a major impact on almost any vehicle, but especially that of a scramjet's performance. The Hyper-X project was to provide a real-world benchmark of the ground test data. The full scale X-43A engine would be operated in the wind tunnel, and then flown, and the data from its operation would then be compared with projections. If these matched, the wind tunnel data would be considered a reliable design tool for future scramjet. If there were significant differences, the reasons for these would have to be identified. Until such information was available, scramjets would lack the technological maturity to be considered for future space launch or high-speed atmospheric flight vehicles.

Peebles, Curtis↗