Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “correlation functions”

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 OSTI2023

Role of three-body dynamics in nucleon-deuteron correlation functions

Correlation functions of hadrons can be accessed in high-energy collisions of atomic nuclei, revealing information about the underlying interaction. This work complements experimental efforts to study nucleon-deuteron Nd—with N=p (proton) or N=n (neutron)—correlations with theory evaluations using different techniques. The correlation functions C nd and C pd are calculated based on a scattering wave function, extending previous benchmarks for the Nd scattering matrix to this new observable. We use hyperspherical harmonics and Faddeev techniques with one of the widely used nucleon-nucleon (NN) interactions, the Argonne v18 potential. Moreover, in the low-energy region we perform additional calculations in the framework of pionless effective field theory. The pd correlation function is computed in the large-energy region to make contact with a recent measurement by the ALICE Collaboration. We show that the scattering wave function has the proper dynamical input to describe an initial rise and subsequent oscillations of C pd as a function of the energy. Effects on the observables using different NN and three-nucleon potentials are evaluated with the conclusion that variations of around 2% are observed. Although these effects are small, future measurements can go beyond this accuracy, allowing for new detailed studies of strong interaction in light nuclear systems. Furthermore, the present study supports the current efforts devoted to the measurement of correlation functions in systems dominated by the strong interactions, such as pd, ppp, Λd, and ppΛ.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1992

Measuring Omega and the real correlation function from the redshift correlation function

Peculiar velocities distort the correlation function of galaxies in redshift space. In the linear regime, the distortion has a characteristic quadrupole plus hexadecapole form. The amplitude of the distortion depends on the cosmological density parameter Omega. Practical formulas are derived here which can be applied to redshift galaxy catalogs to measure Omega in the linear regime. The formulas also yield the real underlying correlation function in the linear regime, corrected for peculiar velocities.

Hamilton, A. J. S.↗
DOE OSTI2025

SILVERRUSH. XIV. Ly α Luminosity Functions and Angular Correlation Functions from 20,000 Ly α Emitters at z ∼ 2.2–7.3 from up to 24 deg 2 HSC-SSP and CHORUS Surveys: Linking the Postreionization Epoch to the Heart of Reionization

Abstract We present luminosity functions (LFs) and angular correlation functions (ACFs) derived from 18,960 Ly α emitters (LAEs) at z = 2.2−7.3 over a wide survey area of ≲24 deg 2 that are identified in the narrowband data of the HSC-SSP and CHORUS surveys. Confirming the large sample with 241 spectroscopically identified LAEs, we determine Ly α LFs and ACFs in the brighter luminosity range down to 0.5 L ⋆ , and confirm that our measurements are consistent with previous studies but offer significantly reduced statistical uncertainties. The improved precision of our ACFs allows us to clearly detect one-halo terms at some redshifts, and provides large-scale bias measurements that indicate host halo masses of ∼10 11 M ⊙ over z ≃ 2−7. By comparing our Ly α LF (ACF) measurements with reionization models, we estimate the neutral hydrogen fractions in the intergalactic medium to be x H i < 0.05 (= 0.06 − 0.03 + 0.12 ) at z = 5.7 and x H i = 0.1 5 − 0.08 + 0.10 ( 0.21 − 0.14 + 0.19 ), 0.1 8 − 0.12 + 0.14 , and 0.7 5 − 0.13 + 0.09 at z = 6.6, 7.0, and 7.3, respectively. Our findings suggest that the neutral hydrogen fraction remains relatively low, x H i ≲ 0.2, at z = 5−7, but increases sharply at z > 7, reaching x H i ∼ 0.9 by z ≃ 8−9, as indicated by recent JWST studies. The combination of our results from LAE observations with recent JWST observations suggests that the major epoch of reionization occurred at z ∼ 7−8, likely driven by the emergence of massive sources emitting significant ionizing photons.

Umeda, Hiroya (ORCID:0009000801675129)↗
NASA NTRS1978

The correlation function for density perturbations in an expanding universe. III The three-point and predictions of the four-point and higher order correlation functions

Higher-order correlation functions for the large-scale distribution of galaxies in space are investigated. It is demonstrated that the three-point correlation function observed by Peebles and Groth (1975) is not consistent with a distribution of perturbations that at present are randomly distributed in space. The two-point correlation function is shown to be independent of how the perturbations are distributed spatially, and a model of clustered perturbations is developed which incorporates a nonuniform perturbation distribution and which explains the three-point correlation function. A model with hierarchical perturbations incorporating the same nonuniform distribution is also constructed; it is found that this model also explains the three-point correlation function, but predicts different results for the four-point and higher-order correlation functions than does the model with clustered perturbations. It is suggested that the model of hierarchical perturbations might be explained by the single assumption of having density fluctuations or discrete objects all of the same mass randomly placed at some initial epoch.

Mcclelland, J.↗
NASA NTRS1979

The correlation function for density perturbations in an expanding universe. IV - The evolution of the correlation function

The evolution of the two-point correlation function for the large-scale distribution of galaxies in an expanding universe is studied on the assumption that the perturbation densities lie in a Gaussian distribution centered on any given mass scale. The perturbations are evolved according to the Friedmann equation, and the correlation function for the resulting distribution of perturbations at the present epoch is calculated. It is found that: (1) the computed correlation function gives a satisfactory fit to the observed function in cosmological models with a density parameter (Omega) of approximately unity, provided that a certain free parameter is suitably adjusted; (2) the power-law slope in the nonlinear regime reflects the initial fluctuation spectrum, provided that the density profile of individual perturbations declines more rapidly than the -2.4 power of distance; and (3) both positive and negative contributions to the correlation function are predicted for cosmological models with Omega less than unity.

Mcclelland, J.↗
NASA NTRS1989

Studies in astronomical time series analysis. III - Fourier transforms, autocorrelation functions, and cross-correlation functions of unevenly spaced data

This paper develops techniques to evaluate the discrete Fourier transform (DFT), the autocorrelation function (ACF), and the cross-correlation function (CCF) of time series which are not evenly sampled. The series may consist of quantized point data (e.g., yes/no processes such as photon arrival). The DFT, which can be inverted to recover the original data and the sampling, is used to compute correlation functions by means of a procedure which is effectively, but not explicitly, an interpolation. The CCF can be computed for two time series not even sampled at the same set of times. Techniques for removing the distortion of the correlation functions caused by the sampling, determining the value of a constant component to the data, and treating unequally weighted data are also discussed. FORTRAN code for the Fourier transform algorithm and numerical examples of the techniques are given.

Scargle, Jeffrey D.↗
DOE OSTI2025

Source function from two-particle correlation function through entropy-regularized Richardson-Lucy deblurring

Source functions are obtained from p – p and d – α correlation functions by applying the Richardson-Lucy (RL) deblurring to the Koonin-Pratt (KP) equation. To prevent fitting of noise in the correlation function, total-variation (TV) regularization is employed that has been effective in ordinary image restoration. TV alone cannot ensure normalization of the source functions. To ensure the latter, we propose a maximum-entropy regularized RL algorithm (MEM-RL). We outline the MEM-RL formalism and optimization strategy for the KP equation, demonstrating its effectiveness on both simulated and experimental data, including the p – p and d – α correlation functions.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
NASA NTRSDate not supplied

A Generalized, Compactly-Supported Correlation Function for Data Assimilation Applications

Correlation functions play an essential role in modern data assimilation, where they are used to model covariances given a set of tunable parameters or applied as tapering functions to localize covariances in ensemble-based schemes. One of the most widely-used correlation functions in data assimilation is the Gaspari and Cohn (1999) piecewise-rational, compactly-supported parametric correlation function (hereafter referred to as GC99). The GC99 correlation function is useful due to its tunable cut-off parameter c and Gaussian-like shape achieved when the parameter a is set to one-half. These properties are attractive for tapering functions in data assimilation applications. However, the GC99 correlation function is homogeneous over Euclidean 3-space and isotropic when restricted to the sphere, properties that may be less than ideal for some geophysical applications. GC99 is also compactly-supported on a sphere of fixed radius, which requires tuning of the cut-off parameter c that can depend on the specific application. This work presents a generalization of the GC99 correlation function that allows the cut-off parameter c and shape parameter a to vary over space to gain more flexibility in shape while maintaining its compact support property. The function, which we call the Generalized Gaspari Cohn (GenGC) correlation function, introduces inhomogeneity in Euclidean 3-space and anisotropy when restricted to the sphere by allowing both parameters c and a to vary, as functions, over the spatial domain. The GC99 correlation function is a special case of GenGC where the functions c and a are held constant, as fixed parameters rather than functions. The GenGC correlation function also generalizes the follow-on to the work of Gaspari and Cohn (1999) presented in Gaspari et al. (2006), which allowed a to vary while keeping c fixed. We illustrate through simple one- and two-dimensional examples the variety of inhomogeneous and anisotropic correlation functions GenGC can produce by varying c and a over space, and suggest applications where they may be useful in data assimilation, such as covariance modeling or localization. In particular, we describe how the GenGC correlation function can be used to construct covariances using correlation length and variance fields derived from dynamics. For example, the correlation length field for advective dynamics is governed by a partial differential equation (PDE) in N spatial dimensions, where N is the number of space dimensions of the state. Correlation length fields can be determined from this PDE and used with GenGC to construct the corresponding correlations. We can then approximate the full covariance by rescaling by the variance, which also satisfies a PDE in N spatial dimensions for advective dynamics. Thus we can approximate the full covariance without solving the covariance PDE, which is in 2N spatial dimensions, by solving just two PDEs each in only N spatial dimensions. This approach to evolving the correlation length and variance fields, then reconstructing the correlations using GenGC, is suggested as an alternative to current methods of covariance modeling in data assimilation algorithms.

GC99↗
DOE OSTI2022

Adversarial autoencoder ensemble for fast and probabilistic reconstructions of few-shot photon correlation functions for solid-state quantum emitters

Second-order photon correlation measurements [g (2) (τ) functions] are widely used to classify single-photon emission purity in quantum emitters or to measure the multiexciton quantum yield of emitters that can simultaneously host multiple excitations – such as quantum dots – by evaluating the value of g (2) (τ = 0). Accumulating enough photons to accurately calculate this value is time consuming and could be accelerated by fitting of few-shot photon correlations. Here, we develop an uncertainty-aware, deep adversarial autoencoder ensemble (AAE) that reconstructs noise-free g (2) (τ) functions from noise-dominated, few-shot inputs. The model is trained with simulated g (2) (τ) functions that are facilely generated by Poisson sampling time bins. The AAE reconstructions are performed orders-of-magnitude faster, with reconstruction errors and estimates of g (2) (τ = 0) that are lower in variance and similar in accuracy compared to Maximum likelihood estimation and Levenberg-Marquardt least-squares fitting approaches, for simulated and experimentally measured few-shot g (2) (τ) functions (~100 two-photon events) of InP/ZnS/ZnSe and CdS/CdSe/CdS quantum dots. The deep-ensemble model comprises eight individual autoencoders, allowing for probabilistic reconstructions of noise-free g (2) (τ) functions, and we show that the predicted variance scales inversely with number of shots, with comparable uncertainties to computationally intensive Markov chain Monte Carlo sampling. Furthermore, this work demonstrates the advantage of machine learning models to perform uncertainty-aware, fast, and accurate reconstructions of simple Poisson-distributed photon correlation functions, allowing for on-the-fly reconstructions and accelerated materials characterization of solid-state quantum emitters.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

MICCO: An Enhanced Multi-GPU Scheduling Framework for Many-Body Correlation Functions

Calculation of many-body correlation functions is one of the critical kernels utilized in many scientific computing areas, especially in Lattice Quantum Chromodynamics (Lattice QCD). It is formalized as a sum of a large number of contraction terms each of which can be represented by a graph consisting of vertices describing quarks inside a hadron node and edges designating quark propagations at specific time intervals. Due to its computation- and memory-intensive nature, real-world physics systems (e.g., multi-meson or multi-baryon systems) explored by Lattice QCD prefer to leverage multi-GPUs. Different from general graph processing, many-body correlation function calculations show two specific features: a large number of computation-/data-intensive kernels and frequently repeated appearances of original and intermediate data. The former results in expensive memory operations such as tensor movements and evictions. The latter offers data reuse opportunities to mitigate the data-intensive nature of many-body correlation function calculations. However, existing graph-based multi-GPU schedulers cannot capture these data-centric features, thus resulting in a sub-optimal performance for many-body correlation function calculations. To address this issue, this paper presents a multi-GPU scheduling framework, MICCO, to accelerate contractions for correlation functions particularly by taking the data dimension (e.g., data reuse and data eviction) into account. This work first performs a comprehensive study on the interplay of data reuse and load balance, and designs two new concepts: local reuse pattern and reuse bound to study the opportunity of achieving the optimal trade-off between them. Based on this study, MICCO proposes a heuristic scheduling algorithm and a machine-learning-based regression model to generate the optimal setting of reuse bounds. Specifically, MICCO is integrated into a real-world Lattice QCD system, Redstar, for the first time running on multiple GPUs. The evaluation demonstrates MICCO outperforms other state-of-art works, achieving up to 2.25× speedup in synthesized datasets, and 1.49× speedup in real-world correlation functions.

Wang, Qihan↗
DOE OSTI2023

Graph Contractions for Calculating Correlation Functions in Lattice QCD

Computing correlation functions for many-particle systems in Lattice QCD is vital to extract nuclear physics observables like the energy spectrum of hadrons such as protons. However, this type of calculation has long been considered to be very challenging and computing-resource intensive because of the complex nature of a hadron composed of quarks with many degrees of freedom. In particular, a correlation function can be calculated through a sum of all possible pairs of quark contractions, each of which is a batched tensor contraction, dictated by Wick's theorem. Because the number of terms of this sum can be very large for any hadronic system of interest, fast evaluation of the sum faces several challenges: an extremely large number of contractions, a huge memory footprint at runtime, and the speed of tensor contractions. In this paper, we present a Lattice QCD analysis software suite, Redstar, which addresses these challenges by utilizing novel algorithmic and software engineering methods targeting modern computing platforms such as many-core CPUs and GPUs. In particular, Redstar represents every term in the sum of a correlation function by a graph, applies efficient graph algorithms to reduce the number of contractions to lower the cost of computations, and minimizes the total memory footprint. Moreover, Redstar carries out the contractions on either CPUs or GPUs utilizing an internal and highly efficient Hadron contraction library. Specifically, we illustrate some important algorithmic optimizations of Redstar, show various key design features of Hadron library, and present the speedup values due to the optimizations along with performance figures for calculating six correlations functions on four computing platforms.

Chen, Jie↗
DOE OSTI2023

Graph Contractions for Calculating Correlation Functions in Lattice QCD

Computing correlation functions for many-particle systems in Lattice QCD is vital to extract nuclear physics observables like the energy spectrum of hadrons such as protons. However, this type of calculation has long been considered to be very challenging and computing-resource intensive because of the complex nature of a hadron composed of quarks with many degrees of freedom. In particular, a correlation function can be calculated through a sum of all possible pairs of quark contractions, each of which is a batched tensor contraction, dictated by Wick's theorem. Because the number of terms of this sum can be very large for any hadronic system of interest, fast evaluation of the sum faces several challenges: an extremely large number of contractions, a huge memory footprint at runtime, and the speed of tensor contractions. In this paper, we present a Lattice QCD analysis software suite, Redstar, which addresses these challenges by utilizing novel algorithmic and software engineering methods targeting modern computing platforms such as many-core CPUs and GPUs. In particular, Redstar represents every term in the sum of a correlation function by a graph, applies efficient graph algorithms to reduce the number of contractions to lower the cost of computations, and minimizes the total memory footprint. Moreover, Redstar carries out the contractions on either CPUs or GPUs utilizing an internal and highly efficient Hadron contraction library. Specifically, we illustrate some important algorithmic optimizations of Redstar, show various key design features of Hadron library, and present the speedup values due to the optimizations along with performance figures for calculating six correlations functions on four computing platforms.

Chen, Jie↗
DOE OSTI2024

Generalized quantum master equations can improve the accuracy of semiclassical predictions of multitime correlation functions

Multitime quantum correlation functions are central objects in physical science, offering a direct link between the experimental observables and the dynamics of an underlying model. While experiments such as 2D spectroscopy and quantum control can now measure such quantities, the accurate simulation of such responses remains computationally expensive and sometimes impossible, depending on the system’s complexity. A natural tool to employ is the generalized quantum master equation (GQME), which can offer computational savings by extending reference dynamics at a comparatively trivial cost. However, dynamical methods that can tackle chemical systems with atomistic resolution, such as those in the semiclassical hierarchy, often suffer from poor accuracy, limiting the credence one might lend to their results. By combining work on the accuracy-boosting formulation of semiclassical memory kernels with recent work on the multitime GQME, here we show for the first time that one can exploit a multitime semiclassical GQME to dramatically improve both the accuracy of coarse mean-field Ehrenfest dynamics and obtain orders of magnitude efficiency gains.

Chemistry↗
NASA NTRS1984

An analysis of the AVE-SESAME I period using statistical structure and correlation functions

Structure and correlation functions are used to describe atmospheric variability during the 10-11 April day of AVE-SESAME 1979 that coincided with the Red River Valley tornado outbreak. The special mesoscale rawinsonde data are employed in calculations involving temperature, geopotential height, horizontal wind speed and mixing ratio. Functional analyses are performed in both the lower and upper troposphere for the composite 24 h experiment period and at individual 3 h observation times. Results show that mesoscale features are prominent during the composite period. Fields of mixing ratio and horizontal wind speed exhibit the greatest amounts of small-scale variance, whereas temperature and geopotential height contain the least. Results for the nine individual times show that small-scale variance is greatest during the convective outbreak. The functions also are used to estimate random errors in the rawinsonde data. Finally, sensitivity analyses are presented to quantify confidence limits of the structure functions.

Fuelberg, H. E.↗
NASA NTRS2006

Spatial Correlation Function of the Chandra Selected Active Galactic Nuclei

We present the spatial correlation function analysis of non-stellar X-ray point sources in the Chandra Large Area Synoptic X-ray Survey of Lockman Hole Northwest (CLASXS). Our 9 ACIS-I fields cover a contiguous solid angle of 0.4 deg(exp 2) and reach a depth of 3 x 10(exp -15) erg/square cm/s in the 2-8 keV band. We supplement our analysis with data from the Chandra Deep Field North (CDFN). The addition of this field allows better probe of the correlation function at small scales. A total of 233 and 252 sources with spectroscopic information are used in the study of the CLASXS and CDFN fields respectively. We calculate both redshift-space and projected correlation functions in co-moving coordinates, averaged over the redshift range of 0.1 < z < 3.0, for both CLASXS and CDFN fields for a standard cosmology with Omega(sub Lambda) = 0.73,Omega(sub M) = 0.27, and h = 0.71 (H(sub 0) = 100h km/s Mpc(exp -1). The correlation function for the CLASXS field over scales of 3 Mpc< s < 200 Mpc can be modeled as a power-law of the form xi(s) = (S/SO)(exp - gamma), with gamma = 1.6(sup +0.4 sub -0.3) and S(sub o) = 8.0(sup +.14 sub -1.5) Mpc. The redshift-space correlation function for CDFN on scales of 1 Mpc< s < 100 Mpc is found to have a similar correlation length so = 8.55(sup +0.74 sub -0.74) Mpc, but a shallower slope (gamma = 1.3 +/- 0.1). The real-space correlation functions derived from the projected correlation functions, are found to be tau(sub 0 = 8.1(sup +1.2 sub -2.2) Mpc, and gamma = 2.1 +/- 0.5 for the CLASXS field, and tau(sub 0) = 5.8(sup +.1.0 sub -1.5) Mpc, gamma = 1.38(sup +0.12 sub -0.14 for the CDFN field. By comparing the real- and redshift-space correlation functions in the combined CLASXS and CDFN samples, we are able to estimate the redshift distortion parameter Beta = 0.4 +/- 0.2 at an effective redshift z = 0.94. We compare the correlation functions for hard and soft spectra sources in the CLASXS field and find no significant difference between the two groups. We have also found that the correlation between X-ray luminosity and clustering amplitude is weak, which, however, is fully consistent with the expectation using the simplest relations between X-ray luminosity, black hole mass, and dark halo mass. We study the evolution of the AGN clustering by dividing the samples into 4 redshift bins over 0.1 Mpc< z <3.0 Mpc. We find a very mild evolution in the clustering amplitude, which show the same evolution trend found in optically selected quasars in the 2dF survey. We estimate the evolution of the bias, and find that the bias increases rapidly with redshift (b(z = 0.45) = 0.95 +/- 0.15 and b(z = 2.07) = 3.03 +/- 0.83): The typical mass of the dark matter halo derived from the bias estimates show little change with redshift. The average halo mass is found to be log (M(sub halo)/M(sun))approximates 12.1. Subject headings: cosmology: observations - large-scale structure of the universe - x-rays: diffuse background - galaxies: nuclei

Yang, Y.↗