Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CORRELATION FUNCTION”

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 523 records · Page 29

A new class of ensemble conserving algorithms for approximate quantum dynamics: Theoretical formulation and model problems

We develop two classes of quasi-classical dynamics that are shown to conserve the initial quantum ensemble when used in combination with the Feynman-Kleinert approximation of the density operator. These dynamics are used to improve the Feynman-Kleinert implementation of the classical Wigner approximation for the evaluation of quantum time correlation functions known as Feynman-Kleinert linearized path-integral. As shown, both classes of dynamics are able to recover the exact classical and high temperature limits of the quantum time correlation function, while a subset is able to recover the exact harmonic limit. A comparison of the approximate quantum time correlation functions obtained from both classes of dynamics is made with the exact results for the challenging model problems of the quartic and double-well potentials. It is found that these dynamics provide a great improvement over the classical Wigner approximation, in which purely classical dynamics are used. In a special case, our first method becomes identical to centroid molecular dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Rotational memory function of SPC/E water

Memory effects are essential for the dynamics of condensed materials and are responsible for non-exponential relaxation of correlation functions of dynamic variables through the memory function. Memory functions of dipole rotations for water have never been calculated directly from molecular dynamics simulations. We present here calculations of memory functions for single-dipole rotations and for the overall dipole moment of the sample for SPC/E water. The normalized memory functions for single-particle and collective dipole dynamics turn out to be nearly identical. This result validates theories of dielectric spectroscopy in terms of single-particle time correlation functions and the connection between the collective and single-particle relaxation times through the Kirkwood factor. The dielectric function in this formalism contains no new dynamic information that does not exist in the single-dipole correlation function. A short memory time, ≲1 fs, justifies the use of the mathematics of rotational diffusion to describe the dynamics of a single molecular dipole moment in bulk water. Here, an analytical equation for the rotational memory time is derived.

Asthagiri, Dilipkumar N. [Oak Ridge National Labor↗

UV background fluctuations and three-point correlations in the large-scale clustering of the Lyman α forest

Using the Lyman α (Lyα) Mass Association Scheme, we make theoretical predictions for the three-dimensional three-point correlation function (3PCF) of the Lyα forest at redshift z = 2.3. We bootstrap results from the (100 h -1 Mpc) 3 Horizon hydrodynamic simulation to a (1 h -1 Gpc) 3 N -body simulation, considering both a uniform ultraviolet background (UVB) and a fluctuating UVB sourced by quasars with a comoving n q ≈ 10 -5 h 3 Mpc -3 placed either in massive haloes or randomly. On scales of 10–30 h -1 Mpc, the flux 3PCF displays hierarchical scaling with the square of the two-point correlation function (2PCF), but with an unusual value of Q ≡ ζ 123 /(ξ 12 ξ 13 + ξ 12 ξ 23 + ξ 13 ξ 23 ) ≈ -4.5 that reflects the low bias of the Lyα forest and the anticorrelation between mass density and transmitted flux. For halo-based quasars and an ionizing photon mean free path of λ = 300 h -1 Mpc comoving, UVB fluctuations moderately depress the 2PCF and 3PCF, with cancelling effects on Q . For λ = 100 or 50 h -1 Mpc, UVB fluctuations substantially boost the 2PCF and 3PCF on large scales, shifting the hierarchical ratio to Q ≈ -3. We scale our simulation results to derive rough estimate of the detectability of the 3PCF in current and future observational data sets for the redshift range z = 2.1–2.6. At r = 10 and 20 h -1 Mpc, we predict a signal-to-noise ratio (SNR) of ~9 and ~7, respectively, for both Baryon Oscillation Spectroscopic Survey (BOSS) and extended BOSS (eBOSS), and ~37 and ~25 for Dark Energy Spectroscopic Instrument (DESI). At r = 40 h -1 Mpc the predicted SNR is lower by a factor of ~3–5. Measuring the flux 3PCF would provide a novel test of the conventional paradigm of the Lyα forest and help separate the contributions of UVB fluctuations and density fluctuations to Lyα forest clustering, thereby solidifying its foundation as a tool of precision cosmology.

Astronomy & Astrophysics↗

Accuracy of cosmological parameters using the baryon acoustic scale

Percent-level measurements of the comoving baryon acoustic scale standard ruler can be used to break degeneracies in parameter constraints from the CMB alone. The sound horizon at the epoch of baryon drag is often used as a proxy for the scale of the peak in the matter density correlation function, and can conveniently be calculated quickly for different cosmological models. However, the measurements are not directly constraining this scale, but rather a measurement of the full correlation function, which depends on the detailed evolution through decoupling. We assess the level of reliability of parameter constraints based on a simple approximation of the acoustic scale compared to a more direct determination from the full numerical two-point correlation function. Using a five-parameter fitting technique similar to recent BAO data analyses, we find that for standard ΛCDM models and extensions with massive neutrinos and additional relativistic degrees of freedom, the approximation is at better than 0.15% for most parameter combinations varying over reasonable ranges.

79 ASTRONOMY AND ASTROPHYSICS↗

Accurate potential energy curve of the LiH{sup +} molecule calculated with explicitly correlated Gaussian functions

Very accurate calculations of the ground-state potential energy curve (PEC) of the LiH{sup +} ion performed with all-electron explicitly correlated Gaussian functions with shifted centers are presented. The variational method is employed. The calculations involve optimization of nonlinear exponential parameters of the Gaussians performed with the aid of the analytical first derivatives of the energy determined with respect to the parameters. The diagonal adiabatic correction is also calculated for each PEC point. The PEC is then used to calculate the vibrational energies of the system. In that calculation, the non-adiabatic effects are accounted for by using an effective vibrational mass obtained by the minimization of the difference between the vibrational energies obtained from the calculations where the Born-Oppenheimer approximation was not assumed and the results of the present calculations.

74 ATOMIC AND MOLECULAR PHYSICS↗

Redshift space clustering of galaxies and cold dark matter model

The distorting effect of peculiar velocities on the power speturm and correlation function of IRAS and optical galaxies is studied. The observed redshift space power spectra and correlation functions of IRAS and optical the galaxies over the entire range of scales are directly compared with the corresponding redshift space distributions using large-scale computer simulations of cold dark matter (CDM) models in order to study the distortion effect of peculiar velocities on the power spectrum and correlation function of the galaxies. It is found that the observed power spectrum of IRAS and optical galaxies is consistent with the spectrum of an Omega = 1 CDM model. The problems that such a model currently faces may be related more to the high value of Omega in the model than to the shape of the spectrum. A low-density CDM model is also investigated and found to be consistent with the data.

Bahcall, Neta A.↗

Range-separated thermal hybrid exchange-correlation density functional for accurate band-gap calculations of warm dense matter

Here, in this study, we present a range-separated thermal hybrid exchange-correlation (XC) functional, termed RS-KDT0, which provides highly accurate band-gap predictions at elevated temperature while reducing the computational cost relative to thermal global hybrid functionals. At zero temperature, RS-KDT0 recovers the widely used ground-state range-separated hybrid functional HSE06. To assess its performance across temperature regimes, we conducted static calculations of electronic pressure, band gaps, and band structures. RS-KDT0 exhibits the same level of accuracy as HSE06 in the low-temperature regime and demonstrates improved agreement with finite-temperature GW calculations at elevated temperatures when compared to other thermal XC functionals such as the generalized gradient approximation (GGA)-level KDT16 and the global hybrid KDT0.

density functional theory↗

Game-Theoretic strategies for systems of components using product-form utilities

Many critical infrastructures are composed of multiple systems of components which are correlated so that disruptions to one may propagate to others. We consider such infrastructures with correlations characterized in two ways: (i) an aggregate failure correlation function specifies the conditional failure probability of the infrastructure given the failure of an individual system, and (ii) a pairwise correlation function between two systems specifies the failure probability of one system given the failure of the other. We formulate a game for ensuring the resilience of the infrastructure, wherein the utility functions of the provider and attacker are products of an infrastructure survival probability term and a cost term, both expressed in terms of the numbers of system components attacked and reinforced. The survival probabilities of individual systems satisfy first-order differential conditions that lead to simple Nash Equilibrium conditions. We then derive sensitivity functions that highlight the dependence of infrastructure resilience on the cost terms, correlation functions, and individual system survival probabilities. We apply these results to simplified models of distributed cloud computing and energy grid infrastructures.

Rao, Nageswara S↗

Defense strategies for infrastructures with multiple systems of components

In several critical infrastructures correlations between the constituent systems represent certain vulnerabilities: disruptions to one may propagate to others and possibly to the entire infrastructure. The correlations between the systems are characterized in two ways in this paper: (i) the aggregate failure correlation function specifies the conditional failure probability of the infrastructure given the failure of an individual system, and (ii) the pairwise correlation function between two systems specifies the failure probability of one system given the failure of the other. The survival probabilities of individual systems satisfy firstorder differential conditions that generalize the contest success functions and statistical independence conditions. We formulate a problem of ensuring the resilience of an infrastructure as a game between the provider and attacker; their utility functions are sums of infrastructure survival probability terms and cost terms, both expressed in terms of the numbers of system components attacked and reinforced. We derive Nash Equilibrium conditions and sensitivity functions that highlight the dependence of infrastructure resilience on the cost terms, correlation functions, and individual system survival probabilities. We apply these results to models of distributed cloud computing and energy grid infrastructures.

Rao, Nageswara S↗

Diffuse X-ray scattering from polished silicon: application of the distorted wave Born approximation

Measured diffuse X-ray scattering data for a `smooth' as well as for a `rough' silicon sample were fit to theoretical expressions within the distorted wave Born approximation (DWBA). Data for the power spectral density (PSD) for both samples were also obtained by means of atomic force microscopy and optical interferometry. The Fourier transforms of trial correlation functions were fit to the PSD data and then applied to the DWBA formalism. The net correlation functions needed to fit the PSD data for each sample comprised the sum of two terms with different cutoff lengths and different self-affine fractal exponents. At zero distance these correlation functions added up to yield net values of σ 2 = (2) 2 and (71) 2 Å 2 for the smooth and rough samples, respectively. X-ray scattering data were obtained at beamline 1-BM of the Advanced Photon Source. Data and fits at values of q z = 0.05 and 0.10 Å −1 for the smooth sample are reported. Good fits for the smooth sample were obtained at both q z values simultaneously, that is, identical fitting parameters were applied at both values of q z . The smooth sample also exhibited weak Yoneda wings and a clear distinction between the strong specular scattering and the weak diffuse scattering. Data for the rough sample were qualitatively different and exhibited very weak scattering at the specular condition in contrast to extremely large Yoneda wings. Fits for the rough sample are reported for q z = 0.04, 0.05, and 0.06 Å −1 . Although the large Yoneda wings could be fit quite well in both position and amplitude, scattering near the specular condition could not be equally well fit by applying the same fitting parameters at all values of q z . Albeit imperfect, best-fitting results at the specular condition were obtained by invoking only diffuse scattering, that is, without including a separate theoretical expression for specular scattering.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Machine-learning-informed scattering correlation analysis of sheared colloids

We have carried out theoretical analysis, Monte Carlo simulations and machine-learning analysis to quantify microscopic rearrangements of dilute dispersions of spherical colloidal particles from coherent scattering intensity. Both monodisperse and polydisperse dispersions of colloids were created and underwent a rearrangement consisting of an affine simple shear and non-affine rearrangement using the Monte Carlo method. We calculated the coherent scattering intensity of the dispersions and the correlation function of intensity before and after the rearrangement and generated a large data set of angular correlation functions for varying system parameters, including number density, polydispersity, shear strain and non-affine rearrangement. Singular value decomposition of the data set shows the feasibility of machine-learning inversion from the correlation function for the polydispersity, shear strain and non-affine rearrangement using only three parameters. A Gaussian process regressor is then trained on the data set and can retrieve the affine shear strain, non-affine rearrangement and polydispersity with relative errors of 3%, 1% and 6%, respectively. Altogether, our model provides a framework for quantitative studies of both steady and non-steady microscopic dynamics of colloidal dispersions using coherent scattering methods.

Gaussian process regression↗

On the replica structure of Sachdev-Ye-Kitaev model

We investigate existence of replica off-diagonal solutions in the field-theoretical description of Sachdev-Ye-Kitaev model. To this end we evaluate a set of local and non-local dynamic correlation functions in the long time limit. We argue that the structure of the soft-mode Schwarzian action is qualitatively different in replica-diagonal vs. replica-off-diagonal scenarios, leading to distinct long-time predictions for the correlation functions. We then evaluate the corresponding correlation functions numerically and compare the simulations with analytical predictions of replica-diagonal and replica-off-diagonal calculations. We conclude that all our numerical results are in a quantitative agreement with the theory based on the replica-diagonal saddle point plus Schwarzian and massive Gaussian fluctuations (the latter do contain replica off-diagonal components). This seems to exclude any contributions from replica-off-diagonal saddle points, at least on the time scales shorter than the inverse many-body level spacing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Grid-generated isotropic homogeneous turbulence at high Reynolds numbers

Consideration is given to an empirical formula for the longitudinal correlation function for grid-generated incompressible fluid turbulence at Reynolds numbers above 12,800. The formula, which relates the longitudinal correlation function to the inverse cube of a dimensionless geometrical ratio, is shown to minimize the global correlation integrals into which the two-point velocity correlation tensor has been substituted subject to a global constraint on the Sobolev concomitent of the longitudinal correlation function. Furthermore, the energy spectrum function associated with the empirical formula is shown to satisfy a tertiary Helmholtz-type linear condition throughout the initial period of decay.

Rosen, G.↗

Reconstructing the primordial spectrum of fluctuations of the universe from the observed nonlinear clustering of galaxies

It was discovered that the nonlinear evolution of the two point correlation function in N-body experiments of galaxy clustering with Omega = 1 appears to be described to good approximation by a simple general formula. The underlying form of the formula is physically motivated, but its detailed representation is obtained empirically by fitting to N-body experiments. In this paper, the formula is presented along with an inverse formula which converts a final, nonlinear correlation function into the initial linear correlation function. The inverse formula is applied to observational data from the CfA, IRAs, and APM galaxy surveys, and the initial spectrum of fluctuations of the universe, if Omega = 1.

Hamilton, A. J. S.↗

Beyond Generalized Eigenvalues

Two analysis techniques, the generalized eigenvalue method (GEM) or Prony's method (PM), are commonly used to analyze statistical estimates of correlation functions produced in lattice quantum field theory calculations. GEM takes full advantage of the matrix structure of correlation functions but only considers individual pairs of time separations when much more data exists. PM can be applied to many time separations and many individual matrix elements simultaneously but does not fully exploit the matrix structure of the correlation function. We combine both these methods into a single framework based on matrix polynomials which we call block Prony method (BPM).

Fleming, George T.↗

Beyond Generalized Eigenvalues in Lattice Quantum Field Theory

Two analysis techniques, the generalized eigenvalue method (GEM) or Prony's (or related) method (PM), are commonly used to analyze statistical estimates of correlation functions produced in lattice quantum field theory calculations. GEM takes full advantage of the matrix structure of correlation functions but only considers individual pairs of time separations when much more data exists. PM can be applied to many time separations and many individual matrix elements simultaneously but does not fully exploit the matrix structure of the correlation function. We combine both these methods into a single framework based on matrix polynomials. As these algebraic methods are well known for producing extensive spectral information about statistically-noisy data, the method should be paired with some information criteria, like the recently proposed Bayesean model averaging.

Fleming, George T.↗

Monopole Fluctuations in Galaxy Surveys

Galaxy clustering provides a powerful way to probe cosmology. This requires understanding of the background mean density of galaxy samples, which is estimated from the survey itself by averaging the observed galaxy number density over the angular position. The angle average includes not only the background mean density but also the monopole fluctuation at each redshift. Here for the first time we compute the monopole fluctuations in galaxy surveys and investigate their impact on galaxy clustering. The monopole fluctuations vary as a function of redshift, and it is correlated with other fluctuations, affecting the two-point correlation function measurements. In an idealized all-sky survey, the rms fluctuation at z = 0.5 can be as large as 7% of the two-point correlation function in amplitude at the baryonic acoustic oscillation scale, and it becomes smaller than 1% at z > 2. The monopole fluctuations are unavoidable, but they can be modeled. We discuss its relation to the integral constraint and the implications for the galaxy clustering analysis.

Baryon acoustic oscillations↗

Fluorescence Correlation Spectroscopy with Photobleaching Correction in Slowly Diffusing Systems

Fluorescence correlation spectroscopy (FCS) is a powerful tool to quantitatively study the diffusion of fluorescently labeled molecules. It allows in principle important questions of macromolecular transport and supramolecular aggregation in living cells to be addressed. However, the crowded environment inside the cells slows diffusion and limits the reservoir of labeled molecules, causing artifacts that arise especially from photobleaching and limit the utility of FCS in these applications. We present a method to compute the time correlation function from weighted photon arrival times, which compensates computationally during the data analysis for the effect of photobleaching. We demonstrate the performance of this method using numerical simulations and experimental data from model solutions. Using this technique, we obtain correlation functions in which the effect of photobleaching has been removed and in turn recover quantitatively accurate mean-square displacements of the fluorophores, especially when deviations from an ideal Gaussian excitation volume are accounted for by using a reference calibration correlation function. This allows quantitative FCS studies of transport processes in challenging environments with substantial photobleaching like in living cells in the future.

74 ATOMIC AND MOLECULAR PHYSICS↗