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 109 records · Page 6

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↗

Path integral Monte Carlo and density functional molecular dynamics simulations of warm dense MgSiO 3

In order to provide a comprehensive theoretical description of MgSiO 3 at extreme conditions, we combine results from path integral Monte Carlo and density functional molecular dynamics simulations and generate a consistent equation of state for this material. We consider a wide range of temperature and density conditions from 10 4 to 10 8 K and from 0.321 to 64.2 g/cm -3 (0.1- to 20-fold the ambient density). We study how the L and K shell electrons are ionized with increasing temperature and pressure. We derive the shock Hugoniot curve and compare with experimental results. Our Hugoniot curve is in good agreement with the experiments, and we predict a broad compression maximum that is dominated by the K shell ionization of all three nuclei while the peak compression ratio of 4.70 is obtained when the Si and Mg nuclei are ionized. To conclude, we analyze the heat capacity and structural properties of the liquid.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

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↗

Applying the Scientific Method of Cybersecurity Research

The cyber environment has rapidly evolved from a curiosity to an essential component of the contemporary world. As the cyber environment has expanded and become more complex, so have the nature of adversaries and styles of attacks. Today, cyber incidents are an expected part of life. As a result, cybersecurity research emerged to address adversarial attacks interfering with or preventing normal cyber activities. Historical response to cybersecurity attacks is heavily skewed to tactical responses with an emphasis on rapid recovery. While threat mitigation is important and can be time critical, a knowledge gap exists with respect to developing the science of cybersecurity. Such a science will enable the development and testing of theories that lead to understanding the broad sweep of cyber threats and the ability to assess trade-offs in sustaining network missions while mitigating attacks. The Asymmetric Resilient Cybersecurity Initiative at Pacific Northwest National Laboratory is a multi-year, multi-million dollar investment to develop approaches for shifting the advantage to the defender and sustaining the operability of systems under attack. The initiative established a Science Council to focus attention on the research process for cybersecurity. The Council shares science practices, critiques research plans, and aids in documenting and reporting reproducible research results. The Council members represent ecology, economics, statistics, physics, computational chemistry, microbiology and genetics, and geochemistry. This paper reports the initial work of the Science Council to implement the scientific method in cybersecurity research. The second section describes the scientific method. The third section in this paper discusses scientific practices for cybersecurity research. Section four describes initial impacts of applying the science practices to cybersecurity research.

Tardiff, Mark F.↗

Optimal structure and parameter learning of Ising models

Reconstruction of the structure and parameters of an Ising model from binary samples is a problem of practical importance in a variety of disciplines, ranging from statistical physics and computational biology to image processing and machine learning. The focus of the research community shifted toward developing universal reconstruction algorithms that are both computationally efficient and require the minimal amount of expensive data. Here, we introduce a new method, interaction screening, which accurately estimates model parameters using local optimization problems. The algorithm provably achieves perfect graph structure recovery with an information-theoretically optimal number of samples, notably in the low-temperature regime, which is known to be the hardest for learning. Here, the efficacy of interaction screening is assessed through extensive numerical tests on synthetic Ising models of various topologies with different types of interactions, as well as on real data produced by a D-Wave quantum computer. Finally, this study shows that the interaction screening method is an exact, tractable, and optimal technique that universally solves the inverse Ising problem.

97 MATHEMATICS AND COMPUTING↗

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↗

Emergent inequality and self-organized social classes in a network of power and frustration

We propose a simple agent-based model on a network to conceptualize the allocation of limited wealth among more abundant expectations at the interplay of power, frustration, and initiative. Concepts imported from the statistical physics of frustrated systems in and out of equilibrium allow us to compare subjective measures of frustration and satisfaction to collective measures of fairness in wealth distribution, such as the Lorenz curve and the Gini index. We find that a completely libertarian, law-of-the-jungle setting, where every agent can acquire wealth from or lose wealth to anybody else invariably leads to a complete polarization of the distribution of wealth vs. opportunity. This picture is however dramatically ameliorated when hard constraints are imposed over agents in the form of a limiting network of transactions. There, an out of equilibrium dynamics of the networks, based on a competition between power and frustration in the decision-making of agents, leads to network coevolution. The ratio of power and frustration controls different dynamical regimes separated by kinetic transitions and characterized by drastically different values of equality. It also leads, for proper values of social initiative, to the emergence of three self-organized social classes, lower, middle, and upper class. Their dynamics, which appears mostly controlled by the middle class, drives a cyclical regime of dramatic social changes.

97 MATHEMATICS AND COMPUTING↗

Equation of State and Damage in Polyethylene

The dynamic response of polymers differs significantly from those of metals, upon which many of the National Laboratories' deformation, damage, and failure models are based. Their moduli, yield strength, and damage characteristics are highly strain rate-, temperature-, and phase-dependent, requiring models that encompass a wide range of phenomena including some not in equilibrium. Recently, Los Alamos developed the Glassy Amorphous Polymer (GAP) 1 model [1] to address limitations in existing models of polymer deformation. GAP captures both volumetric (equation of state) and deviatoric (shear) response, including a non-equilibrium component to the former (a feature determined to be crucial in capturing the low-pressure, viscoelastic response to impact loading). GAP has already been applied to polymers such as PMMA, PTFE, epoxy, and Kel-F 800, but with an emphasis on impact response as opposed to damage or failure. The current effort was launched to address this gap in predictive capability. For reasons that will be made clear, semi-crystalline polyethylene (PE) was chosen to serve as a model system for parameterization and validation. PE (-C 2 H 4- )n is one of the most widely used polymers in industrial and engineering contexts, chiey due to the versatility of its mechanical response. This response can be tuned through network and chain structure, degree of crystallinity, and molecular weight. PE is found in several forms including low density (LDPE), high density (HDPE), and ultra-high molecular weight (UHMWPE). The focus here was on HDPE and UHMWPE, of pedigree described in the following section. Materials were well-characterized prior to study and are representative of semi-crystalline polymers of interest to DOE and DoD. Semi-crystalline PE undergoes a glass transition at low temperature (-35°C) and melts across a range of moderate temperatures (~80-180°C), depending on its structure. It is typically inert chemically, has low strength and high ductility, and the high strength and anisotropy of UHMWPE ber, in particular, have driven its use in engineering, impact, and armor applications. Surprisingly little is known, however, about the influence of PE's crystalline structure and associated phase transitions (including melt) on its response to dynamic compression. A broad suite of experiments was used to calibrate the GAP model for HDPE and UHMWPE. Section IV examines the effects of tensile strain on the structure and integrity of PE crystalline domains. These data were used to inform the preliminary damage model described in Section XII, whose roots lie in statistical physics and network theory. The viscoelastic and plastic components of GAP rely heavily on the stress-strain data of Section VI, which also include dynamic extrusion and Taylor anvil experiments used to validate the damage model. The thermal data of Section X provide crucial inputs to the equilibrium EOS in GAP, as well as the much broader range SESAME EOS whose construction is outlined in Section XI. Section VII details plate impact experiments characterizing the low-pressure shock locus and failure (spall) using in situ electromagnetic gauges. A previously reported \cusp" in the principal Hugoniot near 0.5 kbar was confirmed, and a multi-wave structure was observed over a limited input stress range above the cusp. This cusp is believed due to solid-solid phase transitions associated with the crystalline domains of the polymer.

36 MATERIALS SCIENCE↗

Data-driven Alignment of the Tracker for the Mu2e Experiment

An alignment of the Mu2e Tracker is investigated using the Mu2e Offline Software and Millepede-II. Simulated misalignments of Tracker planes were inserted between Monte-Carlo event generation and before reconstruction using randomly generated $(x,y,z)$-translations and rotations of the rigid-body planes in the detector model. Cosmic ray muons, with no magnetic field present in the Detector Solenoid, will follow a straight-line trajectory. A maximum-likelihood fit is used to reconstruct a track and determine its parameters. Tracks traversing at least four planes, with ten registered straw hits, and a maximum time residual of $20$ ns are selected to be used in an alignment sample. The cuts were chosen to discard tracks with low information content, such as those with high levels of multiple scattering or those with low numbers of measurements contributing to a single track fit. Millepede-II produces the resulting alignment constants according to a simultaneous least-squares fit of the alignment parameters and the track parameters for all input tracks. The improved alignment constants are used to repeat the process and recover an increased number of suitable tracks, further improving the estimate of the alignment constants in subsequent iterations. A performance study of the process was conducted on a misaligned Tracker in the case of translation-only, and translation plus rotational misalignments. In both cases, the correct alignment constants were recovered to within statistical error after three global fit iterations. Additionally, a continuous integration system was developed for the \texttt{Mu2e/Offline} GitHub repository, introducing a convenient automated testing process for all proposed changes (``pull requests"). The tests developed included a full build of the software, followed by tests of the validation routines, geometry overlap checking, and multi-threading mode. Additionally, code quality checks and low-statistics physics comparison checks a re also available.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

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.↗