Search NASA⌕ Search

SEARCH · Search NASA

Results for “statistical fluctuations”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Particle-number fluctuations near the critical point of nuclear matter

Equation of state with the quantum statistics corrections is used for particle-number fluctuations ω of the isotopically symmetric nuclear matter with interparticle van der Waals and Skyrme local density interactions. Here, the fluctuations, ω ∝ 1/K, are analytically derived through the isothermal incompressibility K at first order over a small quantum-statistics parameter. Our approximate analytical results appear to be in good agreement with the results of accurate numerical calculations. These results are also close to those obtained by using more accurate Tolman and Rowlinson expansions of the incompressibility K near the critical point. More general formula for fluctuations ω, improved at the critical point, was obtained for a finite particle-number average $\langle$N$\rangle$ by neglecting, for simplicity, small quantum statistics effects. It is shown that for a large dimensionless parameter, α ∝ K 2 $\langle$N$\rangle$/K", where K" is the second derivative of the incompressibility K as function of the average particle density n, far from the critical point (α >> 1), one finds the traditional asymptote, ω ∝ 1/K, for the fluctuations ω. For a small parameter, α << 1, near the critical point, where K = 0 and α = 0, one obtains another asymptote of ω. These fluctuations, having a maximum near the critical point as function of the average density n, for finite values of $\langle$N$\rangle$ are finite and relatively small, in contrast to the results of the traditional calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

coh3

CoH3 (CoH ver.3) is an optical model, exciton pre-equilibrium, and Hauser-Feshbach statistical model code, which calculates nuclear reaction cross sections for medium to heavy targets in the keV to MeV energy region. This program is written in standard C++, divided into approximately 200 source and header files. CoH solves the Schroedinger equation for optical potentials defined in the code, and calculates differential elastic scattering, reaction, and total cross sections, for neutron, proton, deuteron, triton, 3He, and alpha-particle. Deformed optical potentials are solved with the coupled-channels method, in which the ground state rotational band members, or vibrational phonon states are coupled. The optical model gives particle transmission coefficients that are fed into the statistical model calculations. CoH includes the pre-equilibrium model (exciton model), the direct/semidirect capture model, and the multi-stage Hauser-Feshbach statistical decay with width fluctuation correction based on the Gaussian orthogonal ensemble. For weakly coupled levels, the DWBA (distorted wave Born approximation) method is used to calculate the direct inelastic scattering process to the excited states.

Kawano, Toshihiko↗
DOE OSTI2025

Stochastic frequency fluctuation super-resolution imaging

The inherent non-linearity of intensity correlation functions can be used to spatially distinguish identical emitters beyond the diffraction limit, as achieved, for example, in super-resolution optical fluctuation imaging (SOFI). Here, we propose a complementary concept based on spectral correlation functions, termed spectral fluctuation super-resolution (SFSR) imaging. Through theoretical and computational analysis, we show that spatially resolving time-frequency correlation functions in the image plane can improve the imaging resolution by a factor of $\sqrt2$ in most cases and up to twofold for strictly two emitters. This improvement is achieved by quantifying the degree of correlation in spectral fluctuations across the spatial domain. Experimentally, SFSR can be implemented using a combination of interferometry and photon-correlation measurements. The method works for non-blinking emitters and stochastic spectral fluctuations with arbitrary temporal statistics. This suggests its utility in super-resolution microscopy of quantum emitters at low temperatures, where spectral diffusion is often more pronounced than emitter blinking.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Comment on “Brownian motion of droplets induced by thermal noise”

We simulate phase separated fluids using the Cahn-Hillard fluctuating hydrodynamic (CH-FHD) model and measure the statistical properties of capillary waves generated by thermal fluctuations. Our measurements are in good agreement with stochastic lubrication theory and molecular dynamics simulations but differ significantly from recent CH-FHD results by Zhang et al. [Phys. Rev. E 109, 024208 (2024)2470-004510.1103/PhysRevE.109.024208]. Specifically, we find that capillary wave statistics at thermodynamic equilibrium are independent of transport properties, namely viscosity and species diffusion.

Bell, J B↗
DOE OSTI2022

Reynolds number dependence of Lagrangian dispersion in direct numerical simulations of anisotropic magnetohydrodynamic turbulence

Large-scale magnetic fields thread through the electrically conducting matter of the interplanetary and interstellar medium, stellar interiors and other astrophysical plasmas, producing anisotropic flows with regions of high-Reynolds-number turbulence. It is common to encounter turbulent flows structured by a magnetic field with a strength approximately equal to the root-mean-square magnetic fluctuations. In this work, direct numerical simulations of anisotropic magnetohydrodynamic (MHD) turbulence influenced by such a magnetic field are conducted for a series of cases that have identical resolution, and increasing grid sizes up to $2048^3$ . The result is a series of closely comparable simulations at Reynolds numbers ranging from 1400 up to 21 000. We investigate the influence of the Reynolds number from the Lagrangian viewpoint by tracking fluid particles and calculating single-particle and two-particle statistics. The influence of Alfvénic fluctuations and the fundamental anisotropy on the MHD turbulence in these statistics is discussed. Single-particle diffusion curves exhibit mildly superdiffusive behaviours that differ in the direction aligned with the magnetic field and the direction perpendicular to it. Competing alignment processes affect the dispersion of particle pairs, in particular at the beginning of the inertial subrange of time scales. Scalings for relative dispersion, which become clearer in the inertial subrange for a larger Reynolds number, can be observed that are steeper than indicated by the Richardson prediction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

(U) A General-Purpose Code for Correlated Sampling Using Batch Statistics with MCNP6 for Fixed-Source Problems

Correlated sampling can be used to reduce the uncertainty of a difference of tallies by taking advantage of the negative covariance term in the sandwich formula. Booth first showed how correlated sampling can be applied with batch statistics using MCNP’s tally fluctuation chart (TFC) to reduce the uncertainty of a difference of tallies in fixed-source problems. Booth presented a problem in which a 1273% uncertainty in a difference was reduced to 8% by accounting for correlations. Researchers He and Su recently studied correlated sampling using the TFC in MCNP version 5. They determined that the code did not print enough digits in the TFC tally means for accurate batch statistics in some cases. After modifying the source code, they concluded that “correlated sampling can yield a standard deviation of about one magnitude smaller than that predicted by the direct, un-correlated simulation when the changes in system response are small (say about 1%), which is equivalent to saving in CPU time by a factor of 100. Such saving [sic] becomes less significant as the change in system response becomes larger.” He and Su provided the formulas needed to apply batch statistics to compute the correlated uncertainty of a difference of tallies. In this report, we follow up on their work by providing the formulas needed to apply batch statistics to compute the correlated uncertainty of a ratio of tallies and of a difference of two tallies divided by a third tally. We extend these formulas to differences and ratios of ratios. These formulas are applied to reduce the uncertainty associated with calculating a relative sensitivity. He and Su did not investigate the accuracy of their correlated sampling uncertainty estimates. We use their test problems and evaluate the accuracy of the uncertainty estimates by comparing with results obtained from random sampling, and, in simple cases, with theoretical values of the “exact” uncertainties. We find that the uncertainties obtained from batch statistics are accurate as long as at least 100 batches are used. We present a new computer code, COSUBS (COrrelated Sampling Using Batch Statistics), that reads MCNP6 TFCs and applies correlated sampling using batch statistics for the tally combinations that the user specifies. COSUBS is a very general tool that compares all TFCs for a base case and one or two perturbed cases. It computes uncertainties for ratios if given only a base case. This report is organized as follows. The equations to apply batch statistics to the difference of random tallies are reviewed in Sec. II. Section III presents the equations for applying batch statistics to a ratio of random tallies; this is useful for computing relative sensitivities using a one-sided finite difference and the relative sensitivity using the differential operator method. Section IV presents the equations for applying batch statistics to a difference of two random tallies divided by a third; this is useful for computing a relative sensitivities using a central difference. Section V presents the equations for applying batch statistics to a difference of two ratios with four random tallies. Section VI presents the equations for applying batch statistics to a one-sided finite difference estimate of the relative sensitivity of a ratio (this uses four random tallies). Section VII presents the equations for applying batch statistics to a central difference estimate of the relative sensitivity of a ratio (this uses six random tallies). Section VIII presents the equations for applying batch statistics to a sum of random tallies. Section IX discusses how to apply batch statistics using MCNP6. Section X presents COSUBS, describing its command-line options and logic. Sections XI through XVI present numerical results for various test problems. Section XVII is a summary and conclusions. Appendix A derives the theoretical Monte Carlo tally variance given certain assumptions; these variances are used to verify the batch statistics for some of the problems. Appendix B lists the MCNP6 input for the unperturbed example problem. Appendix C presents modifications made to MCNP6.3 to support this work.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Experimental Measurements of the Wake of a Sphere at Subcritical Reynolds Numbers

This work experimentally investigated the flow phenomena and vortex structures in the wake of a sphere located in a water loop at Reynolds numbers of Re = 850, 1,250, and 1,700. Velocity fields in the wake region were obtained by applying the time-resolved stereoscopic particle image velocimetry (TR-SPIV) technique. From the acquired TR-SPIV velocity vector fields, the statistical values of mean and fluctuating velocities were computed. Spectral analysis, two-point velocity–velocity cross-correlation, proper orthogonal decomposition (POD) and vortex identification analyses were also performed. The velocity fields show a recirculation region that decreases in length with an increase of Reynolds numbers. The power spectra from the spectral analysis had peaks corresponding to a Strouhal number of St = 0.2, which is a value commonly found in the literature studies of flow over a sphere. The two-point cross-correlation analysis revealed elliptical structures in the wake, with estimated integral length scales ranging between 12% and 63% of the sphere diameter. The POD analysis revealed the statistically dominant flow structures that captured the most flow kinetic energy. It is seen that the flow kinetic energy captured in the smaller scale flow structures increased as Reynolds number increased. The POD modes contained smaller structure as the Reynolds number increased and as mode order increased. In addition, spectral analysis performed on the POD temporal coefficients revealed peaks corresponding to St = 0.2, similar to the spectral analysis on the fluctuating velocity. The ability of POD to produce low-order reconstructions of the flow was also utilized to facilitate vortex identification analysis, which identified average vortex sizes of 0.41D for Re1, 0.33D for Re2, and 0.32D for Re3.

Engineering↗
DOE OSTI2021

Data-driven learning of Mori–Zwanzig operators for isotropic turbulence

Developing reduced-order models for turbulent flows, which contain dynamics over a wide range of scales, is an extremely challenging problem. In statistical mechanics, the Mori–Zwanzig (MZ) formalism provides a mathematically exact procedure for constructing reduced-order representations of high-dimensional dynamical systems, where the effects due to the unresolved dynamics are captured in the memory kernel and orthogonal dynamics. Turbulence models based on MZ formalism have been scarce due to the limited knowledge of the MZ operators, which originates from the difficulty in deriving MZ kernels for complex nonlinear dynamical systems. In this work, we apply a recently developed data-driven learning algorithm, which is based on Koopman's description of dynamical systems and Mori's linear projection operator, on a set of fully resolved isotropic turbulence datasets to extract the Mori–Zwanzig operators. With data augmentation using known turbulence symmetries, the extracted Markov term, memory kernel, and orthogonal dynamics are statistically converged and the generalized fluctuation–dissipation relation can be verified. The properties of the memory kernel and orthogonal dynamics, and their dependence on the choices of observables are investigated to address the modeling assumptions that are commonly used in MZ-based models. A series of numerical experiments are then constructed using the extracted kernels to evaluate the memory effects on prediction. The results show that the prediction errors are strongly affected by the choice of observables and can be further reduced by including the past history of the observables in the memory kernel.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Fluctuating hydrodynamics of chiral active fluids

Active materials are characterized by continuous injection of energy at the microscopic level and typically cannot be adequately described by equilibrium thermodynamics. Here we study a class of active fluids in which equilibrium-like properties emerge when fluctuating and activated degrees of freedom are statistically decoupled, such that their mutual information is negligible. We analyse three paradigmatic systems: chiral active fluids composed of spinning frictional particles that are free to translate, oscillating granular gases and active Brownian rollers. In all of these systems, a single effective temperature generated by activity parameterizes both the equation of state and the emergent Boltzmann statistics. The same effective temperature, renormalized by velocity correlations, relates viscosities to steady-state stress fluctuations via a Green–Kubo relation. To rationalize these observations, we develop a theory for the fluctuating hydrodynamics of these non-equilibrium fluids and validate it through large-scale molecular dynamics simulations. Our work sheds light on the microscopic origin of odd viscosities and stress fluctuations characteristic of parity-violating fluids, in which mirror symmetry and detailed balance are broken.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Disentangling the impact of quasiparticles and two-level systems on the statistics of superconducting-qubit lifetime

Temporal fluctuations in the superconducting qubit lifetime, T 1 , present additional challenges in the pursuit of fault-tolerant quantum computing. Although the exact mechanisms remain unclear, T 1 fluctuations are generally attributed to strong coupling between the qubit and a few near-resonant two-level systems (TLSs), which can exchange energy with an ensemble of thermally fluctuating two-level fluctuators (TLFs) at low frequencies. Here, we report T 1 measurements of qubits with varying geometrical footprints and surface dielectrics as a function of temperature. By analyzing the noise spectrum of the qubit depolarization rate, Γ 1 = 1 / T 1 , we disentangle the contributions of TLSs, nonequilibrium quasiparticles (QPs), and equilibrium (thermally excited) QPs to the variance in Γ 1 . We find that the Γ 1 variance in qubits with smaller footprints is more susceptible to QP and TLS fluctuations than that in larger-footprint qubits. Furthermore, the QP-induced variances in all qubits align with the theoretical framework of QP diffusion and fluctuation. These findings offer valuable insights for future qubit design and engineering optimization.

Zhu, Shaojiang [Fermilab] (ORCID:0000000293180092)↗
DOE OSTI2022

Linking Friction Scales from Nano to Macro via Avalanches

Abstract Steady-state fluctuations in the friction force of molybdenum disulfide (MoS 2 ), a prototypical lamellar solid, were analyzed experimentally for newton-scale forces and computationally via molecular dynamics simulations for nanonewton-scale forces. A mean field model links the statics and the dynamics of the friction behavior across these eight orders of magnitude in friction force and six orders of magnitude in friction force fluctuations (i.e., avalanches). Both the statistics and dynamics of the avalanches match model predictions, indicating that friction can be characterized as a series of avalanches with properties that are predictable over a wide range of scales.

Salners, Tyler (ORCID:0000000332647944)↗
DOE OSTI2026

Spectral Properties and Energy Injection in Mercury's Magnetotail Current Sheet

Mercury's magnetotail hosts a thin and highly dynamic current sheet (CS), where magnetic reconnection and strong fluctuations frequently occur. Here, we statistically analyze magnetic field power spectra across 370 magnetotail CSs observed by MESSENGER. About 20% of the events are quasi-laminar, showing single power-law spectra, whereas ∼80% are turbulent, exhibiting a spectral break separating inertial and kinetic ranges. A dawn–dusk asymmetry is identified: inertial-range slopes are systematically shallower on the dawnside, whereas kinetic-range slopes are steeper, indicating more developed turbulence there, consistent with the higher occurrence of reconnection-related processes on the dawnside. Component analysis shows that the transverse components, orthogonal to the tail-aligned principal field ( B X ), display shallow slopes near −1 in the inertial range, suggesting energy injection at ion scales rather than a classical inertial range. These results demonstrate that Mercury's unique plasma environment fundamentally reshapes the initiation of turbulence and the redistribution of energy in the magnetotail.

Mercury↗
DOE OSTI2023

Absolutely Stable Time Crystals at Finite Temperature

Here, we show that locally interacting, periodically driven (Floquet) Hamiltonian dynamics coupled to a Langevin bath support finite-temperature discrete time crystals (DTCs) with an infinite autocorrelation time. By contrast to both prethermal and many-body localized DTCs, the time crystalline order we uncover is stable to arbitrary perturbations, including those that break the time translation symmetry of the underlying drive. Our approach utilizes a general mapping from probabilistic cellular automata to open classical Floquet systems undergoing continuous-time Langevin dynamics. Applying this mapping to a variant of the Toom cellular automaton, which we dub the "π-Toom time crystal," leads to a 2D Floquet Hamiltonian with a finite-temperature DTC phase transition. We provide numerical evidence for the existence of this transition, and analyze the statistics of the finite temperature fluctuations. Finally, we discuss how general results from the field of probabilistic cellular automata imply the existence of discrete time crystals (with an infinite autocorrelation time) in all dimensions, d ≥ 1.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Landau Modes are Eigenmodes of Stellar Systems in the Limit of Zero Collisions

Abstract We consider the spectrum of eigenmodes in a stellar system dominated by gravitational forces in the limit of zero collisions. We show analytically and numerically using the Lenard–Bernstein collision operator that the Landau modes, which are not true eigenmodes in a strictly collisionless system (except for the Jeans unstable mode), become part of the true eigenmode spectrum in the limit of zero collisions. Under these conditions, the continuous spectrum of true eigenmodes in a collisionless system, also known as the Case–van Kampen modes, is eliminated. Furthermore, because the background distribution function in a weakly collisional system can exhibit significant deviations from a Maxwellian distribution function over long times, we show that the spectrum of Landau modes can change drastically even in the presence of slight deviations from a Maxwellian, primarily through the appearance of weakly damped modes that may be otherwise heavily damped for a Maxwellian distribution. Our results provide important insights for developing statistical theories to describe thermal fluctuations in a stellar system, which are currently a subject of great interest for N -body simulations as well as observations of gravitational systems.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

X-ray thermal diffuse scattering as a texture-robust temperature diagnostic for dynamically compressed solids

We present a model of x-ray thermal diffuse scattering (TDS) from a cubic polycrystal with an arbitrary crystallographic texture, based on the classic approach of Warren [B. E. Warren, Acta Crystallogr. 6, 803 (1953)]. We compare the predictions of our model with femtosecond x-ray diffraction patterns gathered from ambient and dynamically compressed rolled copper foils obtained at the High Energy Density instrument of the European X-Ray Free-Electron Laser facility and find that the texture-aware TDS model yields more accurate results than does the conventional powder model owed to Warren. Nevertheless, we further show: with sufficient angular detector coverage, the TDS signal is largely unchanged by sample orientation and in all cases strongly resembles the signal from a perfectly random powder; shot-to-shot fluctuations in the TDS signal resulting from grain-sampling statistics are at the percent level, in stark contrast to the fluctuations in the Bragg-peak intensities (which are over an order of magnitude greater); and TDS is largely unchanged even following texture evolution caused by compression-induced plastic deformation. We conclude that TDS is robust against texture variation, making it a flexible temperature diagnostic applicable just as well to off-the-shelf commercial foils as to ideal powders.

Crystal lattices↗
DOE OSTI2023

Beyond the Condon limit: Condensed phase optical spectra from atomistic simulations

While dark transitions made bright by molecular motions determine the optoelectronic properties of many materials, simulating such non-Condon effects in condensed phase spectroscopy remains a fundamental challenge. We derive a Gaussian theory to predict and analyze condensed phase optical spectra beyond the Condon limit. Our theory introduces novel quantities that encode how nuclear motions modulate the energy gap and transition dipole of electronic transitions in the form of spectral densities. By formulating the theory through a statistical framework of thermal averages and fluctuations, we circumvent the limitations of widely used microscopically harmonic theories, allowing us to tackle systems with generally anharmonic atomistic interactions and non-Condon fluctuations of arbitrary strength. We show how to calculate these spectral densities using first-principles simulations, capturing realistic molecular interactions and incorporating finite-temperature, disorder, and dynamical effects. Our theory accurately predicts the spectra of systems known to exhibit strong non-Condon effects (phenolate in various solvents) and reveals distinct mechanisms for electronic peak splitting: timescale separation of modes that tune non-Condon effects and spectral interference from correlated energy gap and transition dipole fluctuations. We further introduce analysis tools to identify how intramolecular vibrations, solute–solvent interactions, and environmental polarization effects impact dark transitions. Moreover, we prove an upper bound on the strength of cross correlated energy gap and transition dipole fluctuations, thereby elucidating a simple condition that a system must follow for our theory to accurately predict its spectrum.

Chemistry↗
DOE OSTI2024

Adaptive nonequilibrium design of actin-based metamaterials: Fundamental and practical limits of control

The adaptive and surprising emergent properties of biological materials self-assembled in far-from-equilibrium environments serve as an inspiration for efforts to design nanomaterials. In particular, controlling the conditions of self-assembly can modulate material properties, but there is no systematic understanding of either how to parameterize external control or how controllable a given material can be. Here, we demonstrate that branched actin networks can be encoded with metamaterial properties by dynamically controlling the applied force under which they grow and that the protocols can be selected using multi-task reinforcement learning. These actin networks have tunable responses over a large dynamic range depending on the chosen external protocol, providing a pathway to encoding “memory” within these structures. Interestingly, we obtain a bound that relates the dissipation rate and the rate of “encoding” that gives insight into the constraints on control—both physical and information theoretical. Taken together, these results emphasize the utility and necessity of nonequilibrium control for designing self-assembled nanostructures.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

The PDF perspective on the tracer-matter connection: Lagrangian bias and non-Poissonian shot noise

ABSTRACT We study the connection of matter density and its tracers from the probability density function (PDF) perspective. One aspect of this connection is the conditional expectation value 〈δtracer|δm〉 when averaging both tracer and matter density over some scale. We present a new way to incorporate a Lagrangian bias expansion of this expectation value into standard frameworks for modelling the PDF of density fluctuations and counts-in-cells statistics. Using N-body simulations and mock galaxy catalogues we confirm the accuracy of this expansion and compare it to the more commonly used Eulerian parametrization. For haloes hosting typical luminous red galaxies, the Lagrangian model provides a significantly better description of 〈δtracer|δm〉 at second order in perturbations. A second aspect of the matter-tracer connection is shot-noise, i.e. the scatter of tracer density around 〈δtracer|δm〉. It is well known that this noise can be significantly non-Poissonian and we validate the performance of a more general, two-parameter shot-noise model for different tracers and simulations. Both parts of our analysis are meant to pave the way for forthcoming applications to survey data.

Friedrich, Oliver↗