What do the two times in two-time correlation functions mean for interpreting tr-ARPES?
Not provided.
Spectroscopy↗
SEARCH · Search NASA
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.
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↗
Correlation Functions as Nests of Self-Avoiding Paths
We discuss a connection between the XXZ Heisenberg spin chain in the limiting case of zero anisotropy and some aspects of enumerative combinatorics. The representation of the Bethe wave functions in terms of Schur functions allows us to apply the theory of symmetric functions to calculating correlation functions. We provide a combinatorial derivation of the dynamical correlation functions of the projection operator in terms of nests of self-avoiding lattice paths.
97 MATHEMATICS AND COMPUTING↗
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↗
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↗
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↗
Information content of higher order galaxy correlation functions
The shapes of galaxy N-point correlation functions can be used as standard rulers to constrain the distance–redshift relationship. The cosmological density fields traced by late-time galaxy formation are initially nearly Gaussian, and hence, all the cosmological information can be extracted from their two-point correlation function. Subsequent non-linear evolution under gravity, as well as halo and then galaxy formation, generates higher order correlation functions. Since the mapping of the initial to the final density field is, on large scales, invertible, it is often claimed that the information content of the initial field’s power spectrum is equal to that of all the higher order functions of the final, non-linear field. This claim implies that reconstruction of the initial density field from the non-linear field renders analysis of higher order correlation functions of the latter superfluous. We show that this claim is false when the N-point functions are used as standard rulers. Constraints available from joint analysis of the two and three-point correlation functions can, in some cases, exceed those offered by the initial power spectrum. We provide a mathematical justification for this claim and demonstrate it using a large suite of N-body simulations. In particular, we show that for the z = 0 real-space matter field in the limit of vanishing shot-noise, taking modes up to k max = 0.2 h Mpc −1 , using the bispectrum alone offers a factor of 2 reduction in the variance on the cosmic distance scale relative to that available from the linear power spectrum.
79 ASTRONOMY AND ASTROPHYSICS↗
Two-point correlation functions in inhomogeneous and anisotropic cosmologies
Two-point correlation functions are ubiquitous tools of modern cosmology, appearing in disparate topics ranging from cosmological inflation to late-time astrophysics. When the background spacetime is maximally symmetric, invariance arguments can be used to fix the functional dependence of this function as the invariant distance between any two points. In this paper we introduce a novel formalism which fixes this functional dependence directly from the isometries of the background metric, thus allowing one to quickly assess the overall features of Gaussian correlators without resorting to the full machinery of perturbation theory. As an application we construct the CMB temperature correlation function in one inhomogeneous (namely, an off-center LTB model) and two spatially flat and anisotropic (Bianchi) universes, and derive their covariance matrices in the limit of almost Friedmannian symmetry. We show how the method can be extended to arbitrary N -point correlation functions and illustrate its use by constructing three-point correlation functions in some simple geometries.
79 ASTRONOMY AND ASTROPHYSICS↗
Empirical investigation of nuclear correlation function distributions in lattice QCD
Two-point correlation functions of systems with baryon number 𝐵 ∈ {1,2,3,4} are investigated using lattice quantum chromodynamics (QCD). In particular, the empirical distributions of importance-sampling Monte-Carlo samples of these correlation functions are examined as a function of the spacetime separation between the two points and the baryon number. While the exact forms of these distributions are not known for QCD, recent work has determined asymptotic expressions for analogous correlation function distributions in simpler theories such as scalar field theory and the disordered phase of the 𝑂(𝑁) model. The theoretical 𝑂(𝑁) model distributions are found to provide an accurate description of the empirical QCD distributions at zero momentum over a wide range of temporal separations for each baryon number when assessed with a range of different statistical tests. In particular, the behavior of the baryon number 𝐵 QCD correlation function at large temporal separation is well reproduced by the 𝑂(𝑁 ∼ 2/𝐵) model distribution.
Lattice field theory↗
Correlation function of the luminosity distances
We present the correlation function of the luminosity distances in a flat ΛCDM universe. Decomposing the luminosity distance fluctuation into the velocity, the gravitational potential, and the lensing contributions in linear perturbation theory, we study their individual contributions to the correlation function. The lensing contribution is important at large redshift ( z ∼> 0.5) but only for small angular separation (θ ∼< 3°), while the velocity contribution dominates over the other contributions at low redshift or at larger separation. However, the gravitational potential contribution is always subdominant at all scale, if the correct gauge-invariant expression is used. The correlation function of the luminosity distances depends significantly on the matter content, especially for the lensing contribution, thus providing a novel tool of estimating cosmological parameters.
79 ASTRONOMY AND ASTROPHYSICS↗
Euclidean correlation functions in quantum gravity
We calculate Euclidean correlation functions through next-to-leading order in the low-energy effective theory of gravity. We focus on correlation functions of curvature and volume operators, calculating these functions through one-loop order. We show that quantum fluctuations of the background spacetime must be taken into account in order to obtain gauge invariant expressions, and we point out a subtlety associated with the analytic continuation of the conformal mode. Our final expressions for the correlation functions involve only Newton’s constant and the source-sink separation, and they are a universal prediction of the low-energy effective theory. Thus, they serve as a useful point of comparison for nonperturbative lattice formulations of gravity.
Effective field theory↗
Correlation functions from tensor network influence functionals: The case of the spin-boson model
We investigate the application of matrix product state (MPS) representations of the influence functionals (IFs) for the calculation of real-time equilibrium correlation functions in open quantum systems. Focusing specifically on the unbiased spin-boson model, we explore the use of IF-MPSs for complex time propagation, as well as IF-MPSs for constructing correlation functions in the steady state. We examine three different IF approaches: one based on the Kadanoff–Baym contour targeting correlation functions at all times, one based on a complex contour targeting the correlation function at a single time, and a steady state formulation, which avoids imaginary or complex times, while providing access to correlation functions at all times. We show that within the IF language, the steady state formulation provides a powerful approach to evaluate equilibrium correlation functions.
Chemistry↗
Cumulants and correlation functions versus the QCD phase diagram
Here, we discuss the relation of particle number cumulants and correlation functions. It is argued that measuring couplings of the genuine multiparticle correlation functions could provide cleaner information on possible nontrivial dynamics in heavy-ion collisions. We also extract integrated multiproton correlation functions from the presently available experimental data on proton cumulants. We find that the STAR data contain significant four-proton correlations, at least at the lower energies, with indication of changing dynamics in central collisions. We also find that these correlations are rather long ranged in rapidity. Finally, using the Ising model, we demonstrate how the signs of the multiproton correlation functions may be used to exclude certain regions of the phase diagram close to the critical point.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Statistics of baryon correlation functions in lattice QCD
A systematic analysis of the structure of single-baryon correlation functions calculated with lattice QCD is performed, with a particular focus on characterizing the structure of the noise associated with quantum fluctuations. The signal-to-noise problem in these correlation functions is shown, as long suspected, to result from a sign problem. The log-magnitude and complex phase are found to be approximately described by normal and wrapped normal distributions respectively. Properties of circular statistics are used to understand the emergence of a large time noise region where standard energy measurements are unreliable. Power-law tails in the distribution of baryon correlation functions, associated with stable distributions and “Lévy flights,” are found to play a central role in their time evolution. A new method of analyzing correlation functions is considered for which the signal-to-noise ratio of energy measurements is constant, rather than exponentially degrading, with increasing source-sink separation time. This new method includes an additional systematic uncertainty that can be removed by performing an extrapolation, and the signal-to-noise problem reemerges in the statistics of this extrapolation. It is demonstrated that this new method allows accurate results for the nucleon mass to be extracted from the large-time noise region inaccessible to standard methods. The observations presented here are expected to apply to quantum Monte Carlo calculations more generally. Similar methods to those introduced here may lead to practical improvements in analysis of noisier systems.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
The double-soft limit in cosmological correlation functions and graviton exchange effects
The graviton exchange effect on cosmological correlation functions is examined by employing the double-soft limit technique. A new relation among correlation functions that contain the effects due to graviton exchange diagrams in addition to those due to scalar-exchange and scalar-contact-interaction, is derived by using the background field method and independently by the method of Ward identities associated with dilatation symmetry. We compare these three terms, putting small values for the slow-roll parameters and (1− n {sub s} ) ≈ 0.042, where n {sub s} is the scalar spectral index. It is argued that the graviton exchange effects are more dominant than the other two and could be observed in the trispectrum in the double-soft limit. Our observation strengthens the previous work by Seery, Sloth and Vernizzi, in which it has been argued that the graviton exchange dominates in the counter-collinear limit for single field slow-roll inflation.
79 ASTRONOMY AND ASTROPHYSICS↗
A perturbative approach to the redshift space correlation function: beyond the Standard Model
We extend our previous redshift space power spectrum code to the redshift space correlation function. Here we focus on the Gaussian Streaming Model (GSM). Again, the code accommodates a wide range of modified gravity and dark energy models. For the non-linear real space correlation function used in the GSM we use the Fourier transform of the RegPT 1-loop matter power spectrum. We compare predictions of the GSM for a Vainshtein screened and Chameleon screened model as well as GR. These predictions are compared to the Fourier transform of the Taruya, Nishimichi and Saito (TNS) redshift space power spectrum model which is fit to N-body data. We find very good agreement between the Fourier transform of the TNS model and the GSM predictions, with ≤ 6% deviations in the first two correlation function multipoles for all models for redshift space separations in 50Mpc h ≤ s ≤ 180Mpc/ h . Excellent agreement is found in the differences between the modified gravity and GR multipole predictions for both approaches to the redshift space correlation function, highlighting their matched ability in picking up deviations from GR. We elucidate the timeliness of such non-standard templates at the dawn of stage-IV surveys and discuss necessary preparations and extensions needed for upcoming high quality data.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Localizing Holes as Polarons and Predicting Band Gaps, Defect Levels, and Delithiation Energies of Solid-State Materials with a Local Exchange-Correlation Functional
Here, we assess the performance of seven exchange–correlation functionals (some with and some without a Hubbard U correction) for their ability (i) to predict band gaps of silicon, diamond, and Li-ion battery cathode materials, (ii) to localize hole polarons and predict delithiation energies in Li-ion battery cathode materials, and (iii) to predict transition levels of charge carriers of doped silicon and diamond. Both local and hybrid exchange–correlation functionals were tested. The local functionals tend to underestimate band gaps and delocalize polarons. The hybrid functionals very often give a good description of both properties, but they may not be practical for calculations involving large unit cells, large ensembles, or dynamics, and therefore a local functional with a Hubbard U correction is often used (giving the method called DFT+U), where the value of a parameter U is adjusted according to the system and the property being investigated. Keeping in mind the importance of computational cost and the undesirability of having to adjust an empirical parameter for each system or property of interest, we recently developed a local functional, namely HLE17, to try to accurately predict band gaps and excitation energies, and we validated it using mostly main-group solids and molecules. Here we test the performance of HLE17 for its ability to predict band gaps and localize polarons in other solid-state materials, and we compare its performance to that of popular local functionals (PBE, PBEsol, and TPSS), a range-separated hybrid functional with screened exchange (HSE06), and DFT+U. We find that HLE17 predicts more accurate band gaps than other local functionals and can localize holes as polarons, which other local functionals usually fail to do, and for a number of cases it is comparable in performance quality to Hubbard-corrected functionals without the need for system-specific parametrization and to hybrid functionals without the high cost. Because HLE17 does not predict accurate lattice constants, we use the single-point method of quantum chemistry, where the geometry is optimized with one functional and the band gap is calculated with HLE17, or we perform calculations with the lattice constants obtained by TPSS and both the fractional intracell coordinates and the electronic structure obtained by HLE17 (a new method denoted HLE17\\TPSS). In particular, we performed calculations by HLE17//TPSS, HLE17//HSE06, HLE17//DFT+U, and HLE17\\TPSS, and these methods usually agree well with each other and give values similar to experiment.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Correlation functions on the Half-BPS Wilson loop: perturbation and hexagonalization
We compute correlation functions of protected primaries on the 1/2-BPS Wilson loop in N=4 super Yang-Mills theory at weak coupling. We first perform direct perturbative computation at one loop in the planar limit and present explicit formulae for general two-, three- and four-point functions. The results for two- and three-point functions as well as four-point functions in special kinematics are in perfect agreement with the localization computation performed in arXiv:1802.05201 . We then analyze the results in view of the integrability-based approach called “hexagonalization”, which was introduced previously to study the correlation functions in the absence of the Wilson loop. In this approach, one decomposes the correlator into fundamental building blocks called “hexagons”, and glues them back together by summing over the intermediate states. Through the comparison, we conjecture that the correlation functions on the Wilson loop can be computed by contracting hexagons with boundary states, where each boundary state represents a segment of the Wilson loop. As a byproduct, we make predictions for the large-charge asymptotics of the structure constants on the Wilson loop. Along the way, we refine the conjecture for the integrability-based approach to the general non-BPS structure constants on the Wilson loop, proposed originally in arXiv:1706.02989 .
1/N Expansion↗