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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Applications of a Rayleigh-Taylor model to direct-drive laser fusion

Here, this paper presents a simple physics-based model for the interpretation of key metrics in laser direct drive. The only input parameters required are target scale, in-flight aspect ratio, and beam-to-target radius, and the importance of each has been quantified with a tailored set of cryogenic implosion experiments. These analyses lead to compact and accurate predictions of the fusion yield and areal density as a function of hydrodynamic stability, and suggest new ways to take advantage of direct drive. To provide examples, we will discuss how the inferred mix width behaves relative to theory, then show how it could be exploited to perform a direct drive implosion with a Lawson metric or Χ no α of 0.24 ± 0.02—using a novel parameter space at high velocities and beam radii on the OMEGA laser—that projects to ignition at a laser energy ≤2.0 MJ.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Fredkin Staircase: An Integrable System with a Finite-Frequency Drude Peak

We introduce and explore an interacting integrable cellular automaton, the Fredkin staircase, that lies outside the existing classification of such automata, and has a structure that seems to lie beyond that of any existing Bethe-solvable model. The Fredkin staircase has two families of ballistically propagating quasiparticles, each with infinitely many species. Despite the presence of ballistic quasiparticles, charge transport is diffusive in the d.c. limit, albeit with a highly non- gaussian dynamic structure factor. Remarkably, this model exhibits persistent temporal oscillations of the current, leading to a delta-function singularity (Drude peak) in the a.c. conductivity at nonzero frequency. Here, we analytically construct an extensive set of operators that anticommute with the time-evolution operator; the existence of these operators both demonstrates the integrability of the model and allows us to lower-bound the weight of this finite-frequency singularity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Noisy intermediate-scale quantum algorithms

A universal fault-tolerant quantum computer that can efficiently solve problems such as integer factorization and unstructured database search requires millions of qubits with low error rates and long coherence times. While the experimental advancement toward realizing such devices will potentially take decades of research, noisy intermediate-scale quantum (NISQ) computers already exist. These computers are composed of hundreds of noisy qubits, i.e., qubits that are not error corrected, and therefore perform imperfect operations within a limited coherence time. In the search for achieving quantum advantage with these devices, algorithms have been proposed for applications in various disciplines spanning physics, machine learning, quantum chemistry, and combinatorial optimization. The overarching goal of such algorithms is to leverage the limited available resources to perform classically challenging tasks. In this review, a thorough summary of NISQ computational paradigms and algorithms is provided. The key structure of these algorithms and their limitations and advantages are discussed. Finally, a comprehensive overview of various benchmarking and software tools useful for programming and testing NISQ devices is additionally provided.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Dataset for manuscript "Equipartition and the temperature of maximum density of TIP4P/2005 water"

We simulate TIP4P/2005 water in the temperature range of 257 K to 318 K with time-steps 0.25, 0.50, 1.00, 2.00, and 4.00 fs. The density-temperature behavior obtained using 0.25 or 0.50 fs are in excellent agreement with each other but differ from those obtained using time-steps that have been shown earlier to lead to a breakdown of equipartition. The temperature of maximum density (TMD) is 277.15 K with time-step 0.25 or 0.50 fs, but is shifted to progressively lower values for longer time-steps, a trend that holds for different thermostat/barostat combinations. Enhancing the water-water dispersion interaction, as has been recommended for simulating disordered proteins in TIP4P/2005, degrades the description of the liquid-vapor phase envelope. We present a simple physically transparent reasoning to highlight the separation of the time-scales between translational and rotational motion. We also develop a metric, Chi, that we term the equipartition anomaly, to detect equipartition violations in simulations that include molecules that are treated as rigid objects. Calculating Chi is shown to be straightforward and sensitive to equipartition violations. A key takeaway from this study is that using sufficiently short time-steps (less than or equal to 0.5 fs) to preserve equipartition is essential for obtaining meaningful liquid water properties and for producing reliable simulation data, as correct-ensemble sampling is fundamental to ensure reproducibility across codes and simulation alogrithms. The included dataset provides the raw data used in the preparation of the graphs noted in the manuscript.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Machine Learning for Neutron Resonance Evaluations [Slides]

The performance of nuclear reactors and other nuclear systems depends on a precise understanding of the neutron interaction cross sections for materials used in these systems. These cross sections exhibit a resonance structure whose shape is determined in part by the angular momentum quantum numbers of the resonances. The correct assignment of the quantum numbers of neutron resonances is therefore of paramount importance. In this presentation, we describe the application of machine learning to automate the quantum number assignments. Scikit-learn classifiers were trained on simulated resonance data whose statistical properties were chosen to mimic real data. We explored the use of several physics (and random matrix theory)-motivated features for training the classifiers, including the nearest neighbor spacing distribution, cumulative level distribution, and channel width distributions. Initial results demonstrated that we can determine resonance spin groups somewhat reliably. We are now investigating the application of our approach to 52 Cr resonance data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Glassy interphases reinforce elastomeric nanocomposites by enhancing percolation-driven volume expansion under strain

For nearly a century, introduction of nanoparticles to elastomers has yielded extraordinarily tough nanocomposites that are critical to technologies from actuators to tires. The mechanisms by which this reinforcement occurs have nevertheless remained a central open question in material science. One widely debated hypothesis posits that strong interactions between polymer and particles induce "glassy bridges" that cement particles into a cohesive percolating network that resists elongation. Here, molecular dynamics simulations show that glassy particle shells do not primarily provide elongational cohesion. Instead, they amplify an underlying mechanism wherein competition between filler and elastomer networks causes the elastomer's volume to increase on deformation. This induces contributions from the elastomer's bulk modulus, which is of order 1000 times larger than its Young's modulus. These findings establish a unified understanding of low-strain reinforcement in filled elastomers as emanating from volumetric competition between coexisting particulate and elastomeric networks. This reframes and unifies our understanding of low-strain reinforcement, provides a clear-cut diagnostic for the presence of glassy bridging, and offers a new design principle for tough elastomeric nanocomposites.

Computational Physics (physics.comp-ph)↗
DOE OSTI2020

Characterizing Dislocation Defects in Dark Field X-ray Microscopy Images of Bulk Aluminum

Material defects play a large role in material response under shock loading, yet the process of how these defects initiate, propagate, and annihilate is not well understood. Using a newly developed diagnostic, dark field X-ray microscopy (DFXM), we can now visualize the behavior of dislocation defects in materials at the mesoscale under varying conditions. Using data from the European Synchrotron Radiation Facility, we apply a variety of image processing techniques to capture relevant features to locate and characterize size and orientation of dislocation defects in bulk aluminum DFXM images. Beyond simple visualization, this analysis drives statistical characterization of the defects and their dynamic response to further improve the relevant physics models.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Improving Solution Accuracy and Convergence for Stochastic Physics Parameterizations with Colored Noise

Stochastic parameterizations are used in numerical weather prediction and climate modeling to help capture the uncertainty in the simulations and improve their statistical properties. Convergence issues can arise when time integration methods originally developed for deterministic differential equations are applied naively to stochastic problems. In previous studies, it has been demonstrated that a correction term, known in stochastic analysis as the Itô correction, can help improve solution accuracy for various deterministic numerical schemes and ensure convergence to the physically relevant solution without substantial computational overhead. The usual formulation of the Itô correction is valid only when the stochasticity is represented by white noise. In this study, a generalized formulation of the Itô correction is derived for noises of any color. The formulation is applied to a test problem described by an advection–diffusion equation forced with a spectrum of fast processes. We present numerical results for cases with both constant and spatially varying advection velocities to show that, for the same time step sizes, the introduction of the generalized Itô correction helps to substantially reduce time integration error and significantly improve the convergence rate of the numerical solutions when the forcing term in the governing equation is rough (fast varying); alternatively, for the same target accuracy, the generalized Itô correction allows for the use of significantly longer time steps and, hence, helps to reduce the computational cost of the numerical simulation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Quantum Liouville theorem based on Haar measure

Liouville theorem (L theorem) reveals robust incompressibility of the distribution function in phase space, given arbitrary potentials. However, its quantum generalization, Wigner flow, is compressible, i.e., L theorem is only conditionally true (e.g., for perfect Harmonic potential). Here, we develop quantum L theorem (rigorous incompressibility) for arbitrary potentials (interacting or not) in Hamiltonians. Haar measure, instead of symplectic measure dp$\bigwedge$dq used in Wigner’s scheme, plays a central role. The argument is based on general measure theory, independent of specific spaces or coordinates. Comparison of classical and quantum is made: for instance, here we address why Haar measure and metric preservation do not work in the classical case. Applications of the theorems in statistics, topological phase transition, ergodic theory, etc., are discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Modified Brink-Axel hypothesis for astrophysical Gamow-Teller transitions

Weak interaction charged current transition strengths from highly excited nuclear states are fundamental ingredients for accurate modeling of compact object composition and dynamics, but are difficult to obtain either from experiment or theory. For lack of alternatives, calculations have often fallen back upon a generalized Brink-Axel hypothesis, that is, assuming the strength function (transition probability) is independent of the initial nuclear state but depends only upon the transition energy and the weak interaction properties of the parent nucleus ground state. Here we present numerical evidence for a modified “local” Brink-Axel hypothesis for Gamow-Teller transitions for pf-shell nuclei relevant to astrophysical applications. Specifically, while the original Brink-Axel hypothesis does not hold globally, strength functions from initial states nearby in energy are similar within statistical fluctuations. This agrees with previous work on strength function moments. Using this modified hypothesis, we can tackle strength functions at previously intractable initial energies, using semiconverged initial states at arbitrary excitation energy. Our work provides a well-founded method for computing accurate thermal weak transition rates for medium-mass nuclei at temperatures occurring in stellar cores near collapse. We finish by comparing results to previous calculations of astrophysical rates.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Distinct universality classes of diffusive transport from full counting statistics

The hydrodynamic transport of local conserved densities furnishes an effective coarse-grained description of the dynamics of a many-body quantum system. However, the full quantum dynamics contains much more structure beyond the simplified hydrodynamic description. Here we show that systems with the same hydrodynamics can nevertheless belong to distinct dynamical universality classes, as revealed by new classes of experimental observables accessible in synthetic quantum systems, which can, for instance, measure simultaneous site-resolved snapshots of all of the particles in a system. Specifically, we study the full counting statistics of spin transport, whose first moment is related to linear-response transport, but the higher moments go beyond. We present an analytic theory of the full counting statistics of spin transport in various integrable and nonintegrable anisotropic one-dimensional spin models, including the XXZ spin chain. We find that spin transport, while diffusive on average, is governed by a distinct non-Gaussian dynamical universality class in the models considered. We consider a setup in which the left and right half of the chain are initially created at different magnetization densities, and consider the probability distribution of the magnetization transferred between the two half-chains. We derive a closed-form expression for the probability distribution of the magnetization transfer, in terms of random walks on the half-line. We show that this distribution strongly violates the large-deviation form expected for diffusive chaotic systems, and explain the physical origin of this violation. Here, we discuss the crossovers that occur as the initial state is brought closer to global equilibrium. Our predictions can directly be tested in experiments using quantum gas microscopes or superconducting qubit arrays.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Lattice QCD Calculation of the Subtraction Function in Forward Compton Amplitude

The subtraction function plays a pivotal role in calculations involving the forward Compton amplitude, which is crucial for predicting the Lamb shift in muonic atoms, as well as the proton-neutron mass difference. In this Letter, we present a lattice QCD calculation of the subtraction function using two domain wall fermion gauge ensembles near the physical pion mass. We utilize a recently proposed subtraction point, demonstrating its advantage in mitigating statistical and systematic uncertainties by eliminating the need for ground-state subtraction. Our results reveal significant contributions from N⁢π intermediate states to the subtraction function. Incorporating these contributions, we compute the proton, neutron, and nucleon isovector subtraction functions at photon momentum transfer Q 2 ϵ[0,2] GeV 2 . For the proton subtraction function, we compare our lattice results with chiral perturbation theory prediction at low Q 2 and with the results from the perturbative operator-product expansion at high Q 2 . Finally, using these subtraction functions as input, we determine their contribution to two-photon exchange effects in the Lamb shift and isovector nucleon electromagnetic self-energy.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Wormholes and black hole microstates in AdS/CFT

It has long been known that the coarse-grained approximation to the black hole density of states can be computed using classical Euclidean gravity. In this work we argue for another entry in the dictionary between Euclidean gravity and black hole physics, namely that Euclidean wormholes describe a coarse-grained approximation to the energy level statistics of black hole microstates. To do so we use the method of constrained instantons to obtain an integral representation of wormhole amplitudes in Einstein gravity and in full-fledged AdS/CFT. These amplitudes are non-perturbative corrections to the two-boundary problem in AdS quantum gravity. The full amplitude is likely UV sensitive, dominated by small wormholes, but we show it admits an integral transformation with a macroscopic, weakly curved saddle-point approximation. The saddle is the “double cone” geometry of Saad, Shenker, and Stanford, with fixed moduli. In the boundary description this saddle appears to dominate a smeared version of the connected two-point function of the black hole density of states, and suggests level repulsion in the microstate spectrum. Using these methods we further study Euclidean wormholes in pure Einstein gravity and in IIB supergravity on Euclidean AdS 5 × S 5 . We address the perturbative stability of these backgrounds and study brane nucleation instabilities in 10d supergravity. In particular, brane nucleation instabilities of the Euclidean wormholes are lifted by the analytic continuation required to obtain the Lorentzian spectral form factor from gravity. Our results indicate a factorization paradox in AdS/CFT.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Self-calibrating the look-elsewhere effect: fast evaluation of the statistical significance using peak heights

In experiments where one searches a large parameter space for an anomaly, one often finds many spurious noise-induced peaks in the likelihood. This is known as the look-elsewhere effect, and must be corrected for when performing statistical analysis. Here, this paper introduces a method to calibrate the false alarm probability (FAP), or p-value, for a given dataset by considering the heights of the highest peaks in the likelihood. Specifically, we derive an equation relating the global p-value to the rank and height of local maxima. In the simplest form of self-calibration, the look-elsewhere-corrected $\chi^2$ of a physical peak is approximated by the $\chi^2$ of the peak minus the $\chi^2$ of the highest noise-induced peak, with accuracy improved by considering lower peaks. In contrast to alternative methods, this approach has negligible computational cost as peaks in the likelihood are a byproduct of every peak-search analysis. We apply to examples from astronomy, including planet detection, periodograms, and cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

The Completed SDSS-IV Extended Baryon Oscillation Spectroscopic Survey: Baryon Acoustic Oscillations with Ly$α$ Forests

Here, we present a measurement of baryonic acoustic oscillations (BAO) from Lyman-α (Lyα) absorption and quasars at an effective redshift $z = 2.33$ using the complete extended Baryonic Oscillation Spectroscopic Survey (eBOSS). The sixteenth and final eBOSS data release (SDSS DR16) contains all data from eBOSS and its predecessor, the Baryonic Oscillation Spectroscopic Survey (BOSS), providing 210,005 quasars with $z_q > 2.10$ that are used to measure Lyα absorption. We measure the BAO scale both in the auto-correlation of Ly$α$ absorption and in its cross correlation with 341,468 quasars with redshift $z_q > 1.77$. Apart from the statistical gain from new quasars and deeper observations, the main improvements over previous work come from more accurate modeling of physical and instrumental correlations and the use of new sets of mock data. Combining the BAO measurement from the auto- and cross-correlation yields the constraints of the two ratios $D_H(z = 2.33)/r_d = 8.99 ± 0.19$ and $D_M(z = 2.33)/r_d = 37.5 ± 1.1$, where the error bars are statistical. These results are within $1.5σ$ of the prediction of the flat-ΛCDM cosmology of Planck (2016). The analysis code, picca, the catalog of the flux-transmission field measurements, and the $Δ\chi^2$ surfaces are publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Artificial spin ice: Paths forward

Artificial spin ice systems are lithographically prepared arrays of interacting nanoscale magnetic moments with collective behavior resulting from the chosen array geometry. These many-body systems are unusual in that their simple constituent elements can be configured to design their interactions with exquisite control. They can then be probed experimentally over a remarkably large range of time scales and length scales, including imaging of the individual moments. The study of artificial spin ice has broadened well beyond the original focus on systems with connection to spin ice and ice models, and now the term is used in reference to many other structures that exhibit a much wider range of physical phenomena. Here in this perspective, we review progress in the field of artificial spin ice since its inception more than a decade ago. We then discuss prospects for future directions – continuing the exploration of nanomagnetism, statistical spin models, and the potential for technological devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗