Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “AUTOCORRELATION”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Autocorrelation and Multifractal Detrended Fluctuation Analyses Reveal Superdiffusive Mass Transport in Solvent-Filled Nanoporous Media

Fluorescence fluctuation spectroscopy experiments were conducted to better understand the complex mass transport dynamics of organic molecules in liquid-filled nanoporous media. Anodic aluminum oxide (AAO) membranes incorporating 10 and 20 nm diameter cylindrical pores were employed as model materials. Nile red (NR) dye was used as a fluorescent tracer. The dye was dissolved separately in ethanol and toluene at a concentration of 20 nM and used to fill the membrane nanopores. Confocal fluorescence microscopy was employed to capture photon intensity time series data reflecting apparent diffusion of the dye within the pores. Autocorrelation of these data revealed that NR diffusion within the membranes occurred over a broad range of time scales. The autocorrelation decays were fit to a model for one-dimensional diffusion incorporating both fast and slow components having apparent diffusion coefficients, D f and D s , differing by a factor of ∼100. The fast mechanism was attributed to hindered bulk-like diffusion in the central pore cavity, while slow diffusion likely involved absorption of the dye to the pore surfaces. Unfortunately, important evidence of diffusion anomalies is lost in the broad autocorrelation decays obtained. The method of multifractal detrended fluctuation analysis (MF-DFA) was applied to the same data as a means to overcome this limitation. MF-DFA revealed that time series acquired from within the nanopores were multifractal and exhibited evidence of anomalous superdiffusion, likely resulting from the participation of a desorption-mediated diffusion mechanism. Monte Carlo simulations of time series modeling desorption-mediated diffusion in cylindrical nanopores provided support for this assignment. Here, the new knowledge gained affords an improved understanding of hydrocarbon dynamics within nanoporous oil and gas shales.

Diffusion↗
DOE OSTI2025

Forecasting constraints on the high-z IGM thermal state from the Lyman-α forest flux autocorrelation function

ABSTRACT The autocorrelation function of the Lyman-$\alpha$ (Ly $\alpha$) forest flux from high-z quasars probes the small-scale structure of the intergalactic medium (IGM). The thermal state of the IGM, determined by the physics of reionization, sets the small-scale power observed in the Ly $\alpha$ forest. To explore the sensitivity of the autocorrelation function to the IGM’s thermal state, we compute the autocorrelation function from a cosmological hydrodynamical simulation with an instantaneous reionization model and 135 post-processed thermal states. Using mock data sets of 20 quasars, we forecast constraints on $T_0$ and $\gamma$, which characterize the post-processed IGM thermal state, at $5.4 \le z \le 6$. While this model simplifies the IGM’s thermal state, it serves as a key first step in assessing future observational prospects. We also perform an inference test on mocks and re-weight out posterior distributions to guarantee that they exhibit statistically correct behaviour. At $z = 5.4$, we find that an idealized data set constrains $T_0$ to 59 per cent and $\gamma$ to 16 per cent at the 1$\sigma$ equivalent confidence level. To explore more realistic, non-instantaneous reionization scenarios, we analyse four models combining temperature and ultraviolet background (UVB) fluctuations at $z = 5.8$. We find that mock data generated from a model with both temperature and UVB fluctuations can rule out a model with only temperature fluctuations at the $> 1\sigma$ level 73.9 per cent of the time.

Wolfson, Molly↗
DOE OSTI2025

Operator dynamics in Floquet many-body systems

We study operator dynamics in many-body quantum systems, focusing on generic features of systems that are ergodic, spatially extended, and lack conserved densities. Quantum circuits of various types provide simple models for such systems. We focus on Floquet quantum circuits, comparing their behavior with what has been found previously for circuits that are random in time. Floquet circuits, which have discrete time-translation symmetry, represent an intermediate case between circuits that are random in time and lack any symmetry, and systems with a time-independent Hamiltonian and continuous time-translation invariance. By making this comparison, one of our aims is to identify signatures of time-translation symmetry in Floquet operator dynamics. To characterize behavior we examine a variety of quantities in solvable models and numerically: operator autocorrelation functions; the partial spectral form factor; the out-of-time-order correlator (OTOC); and the paths in operator space that make the dominant contributions to the ensemble-averaged autocorrelation functions. Our most striking result is that ensemble-averaged autocorrelation functions show behavior that is distinctively different in Floquet systems compared to systems in which successive time-steps are independent. Specifically, while average autocorrelation functions decay on a microscopic timescale for circuits that are random in time, in Floquet systems they have a late-time tail with a duration that grows parametrically with the size of the operator support. In the simplest models this tail is separated from the initial decay by a minimum, so that the average autocorrelation function has an intermediate-time peak. The existence of these tails provides a way to understand deviations of the spectral form factor from random matrix behavior at times shorter than the Thouless time. In contrast to this feature in autocorrelation functions, we find no new aspects to the behavior of OTOCs for Floquet models compared to random-in-time circuits. We show that this difference between averaged autocorrelation functions and OTOCs can be understood in terms of the paths in operator space that contribute to the two quantities: paths for the former retain a limited support at late times, while paths for the latter are dominated by operator spreading. Published by the American Physical Society 2025

Yoshimura, Takato (ORCID:0000000309159846)↗
DOE OSTI2026

The biogeography of soil and airborne fungi in the Southwestern USA in relation to climate and vegetation

To assess how fungal dispersal might respond to climate change, we examined how climate and geography influence the regional distribution of fungi in soil and air. Specifically, we hypothesized that neighboring fungal communities should be more similar than distant communities (i.e. spatially autocorrelated) and that fungal dispersal should be more limited in soil than in air. We collected soil and air samples from 60 sites across five states in the Southwestern USA. Then, we sequenced the ITS2 region to identify fungal taxa in each sample. Next, we used distance-based redundancy analysis to partition variation in fungal community composition between climate variables and spatial structure. Fungi were indeed spatially autocorrelated. Moreover, precipitation, maximum vapor pressure deficit, and soil moisture were significantly related to fungal community composition in soils. In comparison, only precipitation was significantly related to community composition in the air. After accounting for climate, the strength of spatial autocorrelation did not differ significantly in soilborne versus airborne fungi. Dispersal limitation was evident in soilborne fungi at short distances (<100 km) and was not observed at any distance in airborne fungi. Altogether, climate may influence which fungal taxa are present in soil and air, and fungi could feasibly wind disperse over regional scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Constraints on primordial non-Gaussianity from the cross-correlation of DESI luminous red galaxies and Planck CMB lensing

We use the angular cross-correlation between a luminous red galaxy (LRG) sample from the Dark Energy Spectroscopic Instrument (DESI) Legacy Survey data release DR9 and the Planck cosmic microwave background (CMB) lensing maps to constrain the local primordial non-Gaussianity parameter, f NL , using the scale-dependent galaxy bias effect. The galaxy sample covers approximately 40% of the sky, contains galaxies up to redshift z ∼ 1.4, and is calibrated with the LRG spectra that have been observed for DESI Year 1 (Y1). We apply a nonlinear imaging systematics treatment based on neural networks to remove observational effects that could potentially bias the f NL measurement. Our measurement is performed without blinding, but the full analysis pipeline is tested with simulations including systematics. Using the two-point angular cross-correlation between LRG and CMB lensing only, we find f NL = 39 −38 +40 at the 68% confidence level, and our result is robust in terms of systematics and cosmological assumptions. If we combine this information with the autocorrelation of LRG, applying a scale cut to limit the impact of systematics, we find f NL = 24 −21 +20 at the 68% confidence level. Our results motivate the use of CMB lensing cross-correlations to measure f NL with future datasets, given its stability in terms of observational systematics compared to the angular autocorrelation. Furthermore, performing accurate systematics mitigation is crucially important in order to achieve competitive constraints on f NL from CMB lensing cross-correlation in combination with the tracers’ autocorrelation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Moment method and continued fraction expansion in Floquet operator Krylov space

Recursion methods such as Krylov techniques map complex dynamics to an effective noninteracting problem in one dimension. For example, the operator Krylov space for Floquet dynamics can be mapped to the dynamics of an edge operator of the one-dimensional Floquet inhomogeneous transverse field Ising model (ITFIM), where the latter, after a Jordan-Wigner transformation, is a Floquet model of noninteracting Majorana fermions and the couplings correspond to Krylov angles. We present an application of this showing that a moment method exists where given an autocorrelation function, one can construct the corresponding Krylov angles and from that the corresponding Floquet ITFIM. Consequently, when no solutions for the Krylov angles are obtained, it indicates that the autocorrelation is not generated by unitary dynamics. We highlight this by studying certain special cases: stable m-period dynamics derived using the method of continued fractions, exponentially decaying, and power-law decaying stroboscopic dynamics. Remarkably, our examples of stable m-period dynamics correspond to m-period edge modes for the Floquet ITFIM where, deep in the chain, the couplings correspond to a critical phase. Furthermore our results pave the way to engineer Floquet systems with desired properties of edge modes and also provide examples of persistent edge modes in gapless Floquet systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Shear-Induced Conformations of Salt-Free Polyelectrolytes in Semidilute Solutions

Here, we present coarse-grained molecular dynamics simulations of salt-free polyelectrolyte chains in semidilute solutions under simple shear flow, with full hydrodynamic interactions and explicit dipolar solvent. At equilibrium, chain orientation statistics follow a pseudo-Voigt distribution, and the structural correlation length and chain end-to-end vector autocorrelation function exhibit scaling behavior consistent with theoretical predictions for polyelectrolytes. Under shear, chains transition from coiled to stretched states and the end-to-end vector autocorrelation function reveals oscillatory dynamics at high Weissenberg numbers. Analysis of the gyration tensor and shear strain distributions identified three distinct chain populations with directional alignment along and against the shear gradient. Compared with their neutral polymer counterparts, polyelectrolytes exhibit stronger shear thinning and enhanced chain alignment under the same Weissenberg number, which is attributed to electrostatic interactions and shear-induced counterion release. These findings provide molecular insight into the distinct flow response of charged polymers and have implications for tailoring the rheological properties of polyelectrolyte-based materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Measuring and unbiasing the BAO shift in the Ly α forest with AbacusSummit

ABSTRACT The Dark Energy Spectroscopic Instrument (DESI) places sub- per cent constraints on measurements of the Baryon Acoustic Oscillation (BAO) scaling parameters from the Ly $\alpha$ forest. However, no systematic error budget stemming from non-linearities in the three-dimensional clustering of the Ly $\alpha$ forest is included in the DESI-Ly $\alpha$ analysis. In this work, we measure the size of the shift of the BAO peak using large Ly $\alpha$ forest mocks produced on the N-body simulation suite AbacusSummit, which adopt the Fluctuating–Gunn–Peterson Approximation (FGPA). Specifically, we measure the Ly $\alpha$ autocorrelation and the Ly $\alpha$-quasar cross-correlation functions. To mitigate the noise, we adopt a linear control variates technique, reducing the error bars by a factor of up to $\sim \sqrt{50}$ on large scales. From the autocorrelation, we detect a small positive shift in radial direction of $\Delta \alpha _{\parallel }= 0.35~{{\ \rm per\ cent}}$ at the 3$\sigma$ level and virtually no shift in the transverse direction, $\alpha _\perp$. From the cross-correlation, we see a similar shift to $\Delta \alpha _\parallel$, albeit with larger error bars, and a small negative shift, $\Delta \alpha _{\perp }=\sim$0.25 per cent, at the 2$\sigma$ level. We also make a connection with the Ly $\alpha$ forest effective field theory (EFT) framework and find that the one-loop EFT power spectrum yields unbiased measurements of the BAO shift parameters in radial and transverse direction for Ly $\alpha$ auto- and the Ly $\alpha$-quasar cross-correlation measurements. When using the one-loop EFT framework, we find that we can recover the BAO parameters without a shift, which has important implications for future Ly $\alpha$ forest analyses based on EFT. This work paves the way for novel full-shape analyses of the currently observing DESI and future surveys such as the PFS, WEAVE-QSO, and 4MOST.

Hadzhiyska, Boryana↗
DOE OSTI2026

Floquet operator dynamics and orthogonal polynomials on the unit circle

Operator spreading under stroboscopic time evolution under a unitary is studied. An operator Krylov space is constructed and mapped to orthogonal polynomials on a unit circle (OPUC), as well as to the Krylov space of the edge operator of the Floquet transverse field Ising model with inhomogeneous couplings (ITFIM). The Verblunsky coefficients in the OPUC representation are related to the Krylov angles parametrizing the ITFIM. The relations between the OPUC and spectral functions are summarized and several applications are presented. These include derivation of analytic expressions for the OPUC under persistent m-periodic dynamics, and the numerical construction of the OPUC for autocorrelations of the homogeneous Floquet-Ising model as well as the Z 3 clock model. The numerically obtained Krylov angles of the Z 3 clock model with long-lived period tripled autocorrelations show a spatial periodicity of six, and this observation is used to develop an analytically solvable model for the ITFIM that mimics this behavior.

Floquet systems↗
DOE OSTI2024

PYSIMFRAC: A Python library for synthetic fracture generation and analysis

In this paper, we introduce PYSIMFRAC, an open-source python library for generating 3-D synthetic fracture realizations, integrating with fluid simulators, and performing analysis. PYSIMFRAC allows the user to specify one of three fracture generation techniques (Box, Gaussian, or Spectral) and perform statistical analysis including the autocorrelation, moments, and probability density functions of the fracture surfaces and aperture. This analysis and accessibility of a python library allows the user to create realistic fracture realizations and vary properties of interest. In addition, PYSIMFRAC includes integration examples to two different pore-scale simulators and the discrete fracture network simulator, dfnWorks. The capabilities developed in this work provides opportunity for quick and smooth adoption and implementation by the wider scientific community for accurate characterization of fluid transport in geologic media. We present PYSIMFRAC along with integration examples and discuss the ability to extend PYSIMFRAC from a single complex fracture to complex fracture networks.

58 GEOSCIENCES↗
DOE OSTI2026

Extending SLUSCHI for Automated Diffusion Calculations

We present an extension of the SLUSCHI package (Solid and Liquid in Ultra Small Coexistence with Hovering Interfaces) to enable automated diffusion calculations from first-principles molecular dynamics. While the original SLUSCHI workflow was designed for melting temperature estimation via solid-liquid coexistence, we adapt its input and output handling to isolate the volume search stage and generate one production trajectory suitable for diffusion analysis. Post-processing tools parse VASP outputs, compute mean-square displacements (MSD), and extract tracer diffusivities using the Einstein relation with robust error estimates through block averaging. Diagnostic plots, including MSD curves, running slopes, and velocity autocorrelations, are produced automatically to help identify diffusive regimes. The method has been validated through representative case studies: self-diffusion in Al-Cu liquid alloys, sublattice melting in Li7La3Zr2O12 and Er2O3, interstitial oxygen transport in bcc and fcc Fe, and oxygen diffusivity in Fe-O liquids with variable Si and Al contents. Viscosity and diffusivity are linked through the Stokes-Einstein relation, with composition dependence assessed via simple linear mixing. This capability broadens SLUSCHI from melting-point predictions to transport property evaluation, enabling high-throughput, fully first-principles datasets of diffusion coefficients and viscosities across metals and oxides.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Ensemble Monte Carlo calculations with five novel moves

We introduce five novel types of Monte Carlo (MC) moves that brings the number of moves of ensemble MC calculations from three to eight. So far such calculations have relied on affine invariant stretch moves that were originally introduced by Christen (2007), walk moves by Goodman and Weare (2010) and quadratic moves by Militzer (2023). Ensemble MC methods have been very popular because they harness information about the fitness landscape from a population of walkers rather than relying on expert knowledge. Here we modified the affine method and employed a simplex of points to set the stretch direction. We adopt the simplex concept to quadratic moves. We also generalize quadratic moves to arbitrary order. Finally, we introduce directed moves that employ the values of the probability density while all other types of moves rely solely on the location of the walkers. We apply all algorithms to the Rosenbrock density in 2 and 20 dimensions and to the ring potential in 12 and 24 dimensions. We evaluate their efficiency by comparing error bars, autocorrelation time, travel time, and the level of cohesion that measures whether any walkers were left behind. Our code is open source.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Urban morphology from a landscape perspective: How building morphology distribution land models (BMDLM) emulate pattern and process

Urban form (e.g., building morphology such as height or footprint) can be used to predict environmental footprints, such as energy/water consumption and carbon emissions. Although progress has been made in predicting building characteristics to fill gaps in observation or derive 3-D representations, the relationships between morphology and other variables such as land use and population are poorly understood. Understanding these relationships may enable projections for how cities will evolve with landscapes in the future. A suite of random forest models, the Building Morphology Distribution Land Models (BMDLM), was developed to determine how well building morphology for two distinct statistical measures (central tendency and frequency) can be predicted using land use (e.g., zoning) and population at different resolutions. Clark County, Nevada and Los Angeles County, California are explored as case studies. Generally, 1-km models outperformed 30-m models. Frequency distribution models had the best performance, especially in LA County. Frequency models significantly outperformed spatial autocorrelative models using inverse distance weighting (IDW). BMDLM offers a new take on modeling urban form in which generalized landscape patterns are characterized to understand the influence of population and zoning on urban development, as described by urban scaling theory.

Sturtevant, Jillian [Baylor Univ., Waco, TX (Unite↗
DOE OSTI2024

Stochastic modal velocity field in rough-wall turbulence

Stochastically generated instantaneous velocity profiles are used to reproduce the outer region of rough-wall turbulent boundary layers in a range of Reynolds numbers extending from the wind tunnel to field conditions. Each profile consists in a sequence of steps, defined by the modal velocities and representing uniform momentum zones (UMZs), separated by velocity jumps representing the internal shear layers. Height-dependent UMZ is described by a minimal set of attributes: thickness, mid-height elevation, and streamwise (modal) and vertical velocities. These are informed by experimental observations and reproducing the statistical behaviour of rough-wall turbulence and attached eddy scaling, consistent with the corresponding experimental datasets. Sets of independently generated profiles are reorganized in the streamwise direction to form a spatially consistent modal velocity field, starting from any randomly selected profile. The operation allows one to stretch or compress the velocity field in space, increases the size of the domain and adjusts the size of the largest emerging structures to the Reynolds number of the simulated flow. By imposing the autocorrelation function of the modal velocity field to be anchored on the experimental measurements, we obtain a physically based spatial resolution, which is employed in the computation of the velocity spectrum, and second-order structure functions. The results reproduce the Kolmogorov inertial range extending from the UMZ and their attached-eddy vertical organization to the very-large-scale motions (VLSMs) introduced with the reordering process. The dynamic role of VLSM is confirmed in the –u'w' co-spectra and in their vertical derivative, representing a scale-dependent pressure gradient contribution.

42 ENGINEERING↗
DOE OSTI2024

Adsorption and Solvation Modulate Rhodamine B Diffusion in Ethanol/Water-Filled Anodic Alumina Nanopores

Confinement of solvents and solutes within nanoporous materials frequently leads to the emergence of unique mass transport behaviors that, once fully understood, may lead to improved chemical separations. Here, the diffusion of Rhodamine B (RhB) dye within 10 and 20 nm diameter anodic aluminum oxide (AAO) nanopores filled with binary ethanol/water mixtures is investigated. Mixture compositions spanning from pure ethanol to pure water are employed. The results of confocal fluorescence correlation spectroscopy studies reveal that RhB diffusion occurs by a two-component mechanism comprising composition-dependent fast and slow motions, characterized by diffusion coefficients D f and D s . The results are consistent with those of previous studies performed under more limited conditions [J. Phys. Chem. C, 2023, 127, 411-420]. The fast component scales with mixture viscosity and is assigned to hindered bulk-like diffusion in central pore regions. Slow diffusion likely involves adsorption of RhB to the pore surface and may be described by a desorption mediated mechanism. The occurrence of RhB adsorption to the AAO surface is verified at the single-molecule level by wide-field fluorescence imaging of membrane cross-sectional surfaces. Unique composition dependent trends in the autocorrelation amplitude and in D s that mimic bulk RhB solubility are revealed. D s is found to be smallest in pure ethanol and pure water and largest in intermediate mixtures. These results suggest that RhB surface adsorption is strongest in the pure liquids and weakest in mixtures of intermediate composition, where the dye is least soluble, and most soluble, respectively. As a result, molecular dynamics simulations reveal that a water layer appears on the pore surface under most conditions, while RhB is solvated primarily by ethanol. The composition dependence of RhB diffusion is concluded to reflect its solvation dependent interactions with the pore walls.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Molecular-Level Insights into the NMR Relaxivity of Gadobutrol Using Quantum and Classical Molecular Simulations

MRI is an indispensable diagnostic tool in modern medicine; however, understanding the molecular-level processes governing NMR relaxation of water in the presence of MRI contrast agents remains a challenge, hindering the molecular-guided development of more effective contrast agents. By using quantum-based polarizable force fields, the first-of-its-kind molecular dynamics (MD) simulations of Gadobutrol are reported where the 1 H NMR longitudinal relaxivity r 1 of the aqueous phase is determined without any adjustable parameters. The MD simulations of r 1 dispersion (i.e., frequency dependence) show good agreement with measurements at frequencies of interest in clinical MRI. Importantly, the simulations reveal key insights into the molecular level processes leading to r 1 dispersion by decomposing the NMR dipole–dipole autocorrelation function G(t) into a discrete set of molecular modes, analogous to the eigenmodes of a quantum harmonic oscillator. The molecular modes reveal important aspects of the underlying mechanisms governing r 1 , such as its multiexponential nature and the importance of the second eigenmodal decay. By simply analyzing the MD trajectories on a parameter-free approach, the Gadobutrol simulations show that the outer-shell water contributes ∼50% of the total relaxivity r 1 compared to the inner-shell water, in contrast to simulations of (nonchelated) gadolinium-aqua where the outer shell contributes only ∼15% of r 1 . The deviation between simulations and measurements of r 1 below clinical MRI frequencies is used to determine the low-frequency electron-spin relaxation time for Gadobutrol, in good agreement with independent studies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Light scattering by V 4 O 7 film across the metal–insulator transition

The experimental study of the angle-resolved hemispherical light scattering by V 4 O 7 film within a broad temperature range across metal–insulator transition reveals complex structural reorganization of the film deposited on the c-cut sapphire crystal. The bidirectional scattering distribution function and the surface autocorrelation function were obtained from scattering data to visualize statistics of the spatially resolved contributions of optical inhomogeneities in normal and lateral directions to the surface. The measurements reveal an anisotropic surface roughness distribution due to the twinned domain structure, with significant anisotropy changes across the phase transition. The V 4 O 7 film deposited on sapphire leads to a polydomain structure, minimizing elastic strain energy with distinct multiscale distributions of surface domains. Near T c , the material shows the lowest roughness but the highest lateral disorder of the surface.

36 MATERIALS SCIENCE↗