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 1,117 records · Page 62

Low-frequency Carbon Radio Recombination Lines. I. Calculations of Departure Coefficients

In the first paper of this series, we study the level population problem of recombining carbon ions. We focus our study on high quantum numbers, anticipating observations of carbon radio recombination lines to be carried out by the Low Frequency Array. We solve the level population equation including angular momentum levels with updated collision rates up to high principal quantum numbers. We derive departure coefficients by solving the level population equation in the hydrogenic approximation and including low-temperature dielectronic capture effects. Our results in the hydrogenic approximation agree well with those of previous works. When comparing our results including dielectronic capture, we find differences that we ascribe to updates in the atomic physics (e.g., collision rates) and to the approximate solution method of the statistical equilibrium equations adopted in previous studies. A comparison with observations is discussed in an accompanying article, as radiative transfer effects need to be considered.

79 ASTRONOMY AND ASTROPHYSICS↗

Simulating Z 2 lattice gauge theory on a quantum computer

The utility of quantum computers for simulating lattice gauge theories is currently limited by the noisiness of the physical hardware. Various quantum error mitigation strategies exist to reduce the statistical and systematic uncertainties in quantum simulations via improved algorithms and analysis strategies. We perform quantum simulations of Z 2 gauge theory with matter to study the efficacy and interplay of different error mitigation methods: readout error mitigation, randomized compiling, rescaling, and dynamical decoupling. We compute Minkowski correlation functions in this confining gauge theory and extract the mass of the lightest spin-1 state from fits to their time dependence. Quantum error mitigation extends the range of times over which our correlation function calculations are accurate by a factor of 6 and is therefore essential for obtaining reliable masses.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The Statistical Mechanics of Ideal MHD Turbulence

Turbulence is a universal, nonlinear phenomenon found in all energetic fluid and plasma motion. In particular. understanding magneto hydrodynamic (MHD) turbulence and incorporating its effects in the computation and prediction of the flow of ionized gases in space, for example, are great challenges that must be met if such computations and predictions are to be meaningful. Although a general solution to the "problem of turbulence" does not exist in closed form, numerical integrations allow us to explore the phase space of solutions for both ideal and dissipative flows. For homogeneous, incompressible turbulence, Fourier methods are appropriate, and phase space is defined by the Fourier coefficients of the physical fields. In the case of ideal MHD flows, a fairly robust statistical mechanics has been developed, in which the symmetry and ergodic properties of phase space is understood. A discussion of these properties will illuminate our principal discovery: Coherent structure and randomness co-exist in ideal MHD turbulence. For dissipative flows, as opposed to ideal flows, progress beyond the dimensional analysis of Kolmogorov has been difficult. Here, some possible future directions that draw on the ideal results will also be discussed. Our conclusion will be that while ideal turbulence is now well understood, real turbulence still presents great challenges.

Shebalin, John V.↗

Modeling Protein–Protein and Protein–Ligand Interactions by the ClusPro Team in CASP16

ABSTRACT In the CASP16 experiment, our team employed hybrid computational strategies to predict both protein–protein and protein–ligand complex structures. For protein–protein docking, we combined physics‐based sampling—using ClusPro FFT docking and molecular dynamics—with AlphaFold (AF)‐based sampling, followed by AF‐based refinement. Our method produced numerous high‐accuracy complex models, including cases where AF alone failed, underscoring the critical role of physics‐based sampling alongside deep learning‐based refinement. For protein–ligand docking, we integrated the ClusPro LigTBM template‐based approach with a machine learning‐based confidence model for rescoring. The method preserves conserved interaction fragments derived from homologous complexes, followed by local resampling using physics‐based sampling and a diffusion model. Our template‐based strategy achieved a mean lDDT‐PLI of 0.69 across 233 targets, which was highly competitive. These results demonstrate that combining physics‐based modeling with AI‐driven refinement can significantly enhance the accuracy of both protein–protein and protein–ligand structure predictions.

Ashizawa, Ryota [Department of Applied Mathematics↗

Theoretical study of various nonlinear phenomena in plasma systems and scaling of magneto-inertial-fusion targets

Plasma physics is an exciting field of study with a wide variety of nonlinear processes that come into play. Examples of such processes include the interaction of small-scale turbulence with large-scale plasma structures and the nonlinear saturation of plasma instabilities, for example those of magneto-hydrodynamical nature. During this Truman LDRD project, I studied a collection of nonlinear problems that are of interest to the field of plasma physics. This LDRD report summarizes four main research accomplishments. First, a new statistical model for describing inhomogeneous drift-wave turbulence inter- acting with zonal flows was developed. This new model includes the effects of nonlinear wave-wave collisions, which are expected to change the spectrum of the underlying DW turbulence and therefore the generation of zonal flows. Second, a new mathematical formalism was proposed to systematically apply the non- linear WKB approximation to general field theories, including those often used in fluid dynamics. This formalism represents an interesting tool for studying physical systems that show an explicit scale separation. Third, a weakly nonlinear model was developed to describe the magneto-Rayleigh-Taylor instability. This instability is of paramount importance to understand as it can reduce the performance of magnetic-inertial-fusion (MIF) platforms. The developed models captures the effects of harmonic generation and saturation of the linear growth of the instability. Finally, a framework was proposed for scaling magneto-inertial fusion (MIF) targets to larger pulsed-power drivers. From this framework, a set of scaling rules were derived that conserve the physical regimes of MIF systems when scaling up in peak current. By doing so, deleterious nonlinear processes that affect MIF performance may be kept at bay.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Confronting Models with Data: The GEWEX Cloud Systems Study

The GEWEX Cloud System Study (GCSS; GEWEX is the Global Energy and Water Cycle Experiment) was organized to promote development of improved parameterizations of cloud systems for use in climate and numerical weather prediction models, with an emphasis on the climate applications. The strategy of GCSS is to use two distinct kinds of models to analyze and understand observations of the behavior of several different types of clouds systems. Cloud-system-resolving models (CSRMs) have high enough spatial and temporal resolutions to represent individual cloud elements, but cover a wide enough range of space and time scales to permit statistical analysis of simulated cloud systems. Results from CSRMs are compared with detailed observations, representing specific cases based on field experiments, and also with statistical composites obtained from satellite and meteorological analyses. Single-column models (SCMs) are the surgically extracted column physics of atmospheric general circulation models. SCMs are used to test cloud parameterizations in an un-coupled mode, by comparison with field data and statistical composites. In the original GCSS strategy, data is collected in various field programs and provided to the CSRM Community, which uses the data to "certify" the CSRMs as reliable tools for the simulation of particular cloud regimes, and then uses the CSRMs to develop parameterizations, which are provided to the GCM Community. We report here the results of a re-thinking of the scientific strategy of GCSS, which takes into account the practical issues that arise in confronting models with data. The main elements of the proposed new strategy are a more active role for the large-scale modeling community, and an explicit recognition of the importance of data integration.

Randall, David↗

Coupled-oscillator theory of dispersion and Casimir-Polder interactions

We address the question of the applicability of the argument theorem (of complex variable theory) to the calculation of two distinct energies: (i) the first-order dispersion interaction energy of two separated oscillators, when one of the oscillators is excited initially and (ii) the Casimir-Polder interaction of a ground-state quantum oscillator near a perfectly conducting plane. We show that the argument theorem can be used to obtain the generally accepted equation for the first-order dispersion interaction energy, which is oscillatory and varies as the inverse power of the separation r of the oscillators for separations much greater than an optical wavelength. However, for such separations, the interaction energy cannot be transformed into an integral over the positive imaginary axis. If the argument theorem is used incorrectly to relate the interaction energy to an integral over the positive imaginary axis, the interaction energy is non-oscillatory and varies as r{sup −4}, a result found by several authors. Rather remarkably, this incorrect expression for the dispersion energy actually corresponds to the nonperturbative Casimir-Polder energy for a ground-state quantum oscillator near a perfectly conducting wall, as we show using the so-called “remarkable formula” for the free energy of an oscillator coupled to a heat bath [G. W. Ford, J. T. Lewis, and R. F. O’Connell, Phys. Rev. Lett. 55, 2273 (1985)]. A derivation of that formula from basic results of statistical mechanics and the independent oscillator model of a heat bath is presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Extraction of $\sin^2\theta^{\rm lept}_{\rm eff}$ and indirect measurement of $m_w$ from the 9 fb$^{-1}$ full run ii sample of $\mu^+\mu^-$ events at cdf

We report on the extraction of $\sin^2\theta^{\rm lept}_{\rm eff}$ and indirect measurement of the mass of the W boson from the forward-backward asymmetry of $\mu^+\mu^-$ events in the $Z$ boson mass region. The data sample collected by the CDF detector corresponds to the full 9 fb$^{-1}$ run II sample. We measure $\sin^2 \theta^{\rm lept}_{\rm eff} = 0.2315 \pm 0.0010$,$ \sin^2 \theta_W = 0.2233 \pm 0.0009$ and $M_W ({\rm indirect}) = 80.365 \pm 0.047 \;{\rm GeV}/c^2$, where each uncertainty includes both statistical and systematic contributions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The origin of quark color

Some 50 years ago, just for fun, I began playing around with different types of particle statistics. Those investigations led to a surprising application.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The size of the Proton: Recent Results from the PRad Experiment

Two new extremely high precision measurements of the proton rms charge radius performed in 2010-2013 with muonic hydrogen atom demonstrated up to six standard deviations smaller values than the accepted average from all previous experiments performed with different methods (scattering and atomic spectroscopy) on ordinary hydrogen. This discrepancy triggered the well-known ?proton radius puzzle? in hadronic physics. To address this puzzle, the PRad collaboration in spring of 2016 performed a novel magnetic-spectrometer-free ep-scattering experiment in Hall B at Jefferson Lab accumulating high statistics and a rich experimental data set. The specifics of the PRad experiment and the final physics results, including the extracted proton radius, as well as, the current status of the puzzle will be presented and discussed in this talk.

Gasparian, Ashot↗

Objective determination of image end-members in spectral mixture analysis

Spectral mixture analysis was shown to be a powerful, multifaceted tool for analysis of multi- and hyper-spectral data. The essence of the first phase of the approach is to determine a set of image end-members that best account for the spectral variance in an image cube within a constrained, linear least squares mixing model. The selection of the image end-member is usually achieved using a priori knowledge and successive trial and error solutions to refine the total number and physical location of the end-members. However, in many situations a more objective method of determining these essential components is desired. The problem of image end-member determination was approached objectively by using the inherent variance of the data. Unlike purely statistical methods such as factor analysis, this approach derives solutions that conform to a physically realistic model.

Tompkins, Stefanie↗

Objective determination of image end-members in spectral mixture analysis of AVIRIS data

Spectral mixture analysis has been shown to be a powerful, multifaceted tool for analysis of multi- and hyper-spectral data. Applications of AVIRIS data have ranged from mapping soils and bedrock to ecosystem studies. During the first phase of the approach, a set of end-members are selected from an image cube (image end-members) that best account for its spectral variance within a constrained, linear least squares mixing model. These image end-members are usually selected using a priori knowledge and successive trial and error solutions to refine the total number and physical location of the end-members. However, in many situations a more objective method of determining these essential components is desired. We approach the problem of image end-member determination objectively by using the inherent variance of the data. Unlike purely statistical methods such as factor analysis, this approach derives solutions that conform to a physically realistic model.

Tompkins, Stefanie↗

Cosmological Probes of Fundamental Physics (Final Report)

This report summarizes the research performed by Hirata’s group under award # DE-SC0006624, which contains both theoretical and observational parts. The theoretical part includes studies of the physics of the reionization era, the period in cosmic history when the first stars and galaxies emitted ultraviolet radiation and turned the intergalactic medium from neutral to ionized. We investigated the ways in which future observations of this era using the 21 cm line of hydrogen could probe inflation (through the effects of gravitational waves on the statistics of this line) and primordial magnetic fields. We also investigated how the physics of the reionization era could impact measurements of the Lyman-α forest using current (BOSS/eBOSS) and near-future (DESI) instruments. It also included a detailed study of sterile neutrinos as a dark matter candidate, including the early Universe modeling required to predict the distribution of sterile neutrinos on small scales. The observational arm of the project included preparations for ground- and space-based dark energy surveys (including PFS and WFIRST), as well as measurements of the growth of structure from the Sloan Digital Sky Survey.

79 ASTRONOMY AND ASTROPHYSICS↗

Kinetic theory of particle-in-cell simulation plasma and the ensemble averaging technique

Abstract We derive the kinetic theory of fluctuations in physically and numerically stable particle-in-cell (PIC) simulations of electrostatic plasmas. The starting point is the single-time correlation at the start of the simulation between the statistical fluctuations of the weighted densities of macroparticle centers in the plasma particle phase-space. The fluctuations are associated with different initial conditions, typically due to the random initial conditions (in velocity space) of the macroparticles/simulation plasma, assigned according to their initial distribution of probability. The single-time correlations at all time steps and in each spatial grid cell are then determined from the Laplace–Fourier transforms of the discretized Klimontovich-like equation for the macroparticles and Maxwell’s equations for the fields, as computed by modern PIC codes. We recover the expressions for the electrostatic field and the plasma particle density fluctuation autocorrelation spectra as well as the kinetic equations describing the average evolution of PIC-simulated plasma particles, first derived by Langdon (1970b Proc. 4th Conf. Numerical Simulation of Plasmas ) using a test macroparticle approach perturbing a discretized Vlasovian plasma and then averaging the obtained physical quantity over the initial macroparticle velocity distribution. We generalize and extend these results to the modern algorithms in PIC codes using arbitrary macroparticle weights. Analytical estimates of statistical fluctuation amplitudes are derived as a function of the plasma simulation parameters, using the central limit theorem in the limit of a large number of macroparticles per cell. The theory is then used to analyze the ensemble averaging technique of PIC simulations where statistical averages are performed over ensembles of PIC simulations, modeling the same plasma physics problem but using different statistical realizations of the initial distribution functions of the macroparticles. This method is illustrated by linear Landau damping uncovering (from noise, which is usually considered numerical) the physical fluctuations driven by a single small amplitude electrostatic wave perturbing a PIC simulation plasma in equilibrium.

fluctuations correlations↗

Real-time calibration of a feedback trap

Feedback traps use closed-loop control to trap or manipulate small particles and molecules in solution. They have been applied to the measurement of physical and chemical properties of particles and to explore fundamental questions in the non-equilibrium statistical mechanics of small systems. These applications have been hampered by drifts in the electric forces used to manipulate the particles. Although the drifts are small for measurements on the order of seconds, they dominate on time scales of minutes or slower. Here, we show that a recursive maximum likelihood (RML) algorithm can allow real-time measurement and control of electric and stochastic forces over time scales of hours. Simulations show that the RML algorithm recovers known parameters accurately. Experimental estimates of diffusion coefficients are also consistent with expected physical properties.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Open quantum systems and random matrix theory

A simple model for open quantum systems is analyzed with RMT. The system is coupled to the continuum in a minimal way. In this paper we see the effect of opening the system on the level statistics, in particular the level spacing, width distribution and Δ{sub 3}(L) statistic are examined as a function of the strength of this coupling. The usual super-radiant state is observed, and it is seen that as it is formed, the level spacing and Δ{sub 3}(L) statistic exhibit the signatures of missed levels.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Flow field topology of transient mixing driven by buoyancy

Transient mixing driven by buoyancy occurs through the birth of a symmetric Rayleigh-Taylor morphology (RTM) structure for large length scales. Beyond its critical bifurcation the RTM structure exhibits self-similarity and occurs on smaller and smaller length scales. The dynamics of the RTM structure, its nonlinear growth and internal collision, show that its genesis occurs from an explosive bifurcation which leads to the overlap of resonance regions in phase space. This event shows the coexistence of regular and chaotic regions in phase space which is corroborated with the existence of horseshoe maps. A measure of local chaos given by the topological entropy indicates that as the system evolves there is growth of uncertainty. Breakdown of the dissipative RTM structure occurs during the transition from explosive to catastrophic bifurcation; this event gives rise to annihilation of the separatrices which drives overlap of resonance regions. The global bifurcation of explosive and catastrophic events in phase space for the large length scale of the RTM structure serves as a template for which mixing occurs on smaller and smaller length scales. Copyright 2004 American Institute of Physics.

Nonlinear Dynamics↗