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 1,459 records · Page 81

Gas Transfer Across Air‐Water Interfaces in Inland Waters: From Micro‐Eddies to Super‐Statistics

In inland water covering lakes, reservoirs, and ponds, the gas exchange of slightly soluble gases such as carbon dioxide, dimethyl sulfide, methane, or oxygen across a clean and nearly flat air‐water interface is routinely described using a water‐side mean gas transfer velocity $\overline{k_{L}}$, where overline indicates time or ensemble averaging. The micro‐eddy surface renewal model predicts $\overline{k_{L}}$ = α o Sc -1/2 ($v\bar{ϵ}$) 1/4 , where Sc is the molecular Schmidt number, $v$ is the water kinematic viscosity, and $\bar{ϵ}$ is the waterside mean turbulent kinetic energy dissipation rate at or near the interface. While α o = 0.39 - 0.46 has been reported across a number of data sets, others report large scatter or variability around this value range. It is shown here that this scatter can be partly explained by high temporal variability in instantaneous ϵ around $\bar{ϵ}$, a mechanism that was not previously considered. As the coefficient of variation (CV e ) in ϵ increases, α o must be adjusted by a multiplier (1 = CV e 2 ) -3/32 that was derived from a log‐normal model for the probability density function of ϵ. Reported variations in α o with a macro‐scale Reynolds number can also be partly attributed to intermittency effects in ϵ. Such intermittency is characterized by the long‐range (i.e., power‐law decay) spatial auto‐correlation function of ϵ. That α o varies with a macro‐scale Reynolds number does not necessarily violate the micro‐eddy model. Instead, it points to a coordination between the macro‐ and micro‐scales arising from the transfer of energy across scales in the energy cascade.

Batchelor scale↗

Implementation and benchmark of a long-range corrected functional in the density functional based tight-binding method

Bridging the gap between first principles methods and empirical schemes, the density functional based tight-binding method (DFTB) has become a versatile tool in predictive atomistic simulations over the past years. One of the major restrictions of this method is the limitation to local or gradient corrected exchange-correlation functionals. This excludes the important class of hybrid or long-range corrected functionals, which are advantageous in thermochemistry, as well as in the computation of vibrational, photoelectron, and optical spectra. The present work provides a detailed account of the implementation of DFTB for a long-range corrected functional in generalized Kohn-Sham theory. We apply the method to a set of organic molecules and compare ionization potentials and electron affinities with the original DFTB method and higher level theory. The new scheme cures the significant overpolarization in electric fields found for local DFTB, which parallels the functional dependence in first principles density functional theory (DFT). At the same time, the computational savings with respect to full DFT calculations are not compromised as evidenced by numerical benchmark data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mock data sets for the Eboss and DESI Lyman-α forest surveys

We present a publicly-available code to generate sets of mock Lyman-α (Lyα) forest data that have realistic large-scale correlations including those due to the Baryonic Acoustic Oscillations (BAO). The primary purpose of these mocks is to test the analysis procedures of the Extended Baryon Oscillation Survey (eBOSS) and the Dark Energy Spectroscopy Instrument (DESI) surveys. The transmitted flux fraction, F(λ), of background quasars due to Lyα absorption in the intergalactic medium (IGM) is simulated using the Fluctuating Gunn-Petterson Approximation (FGPA) applied to Gaussian random fields produced through the use of fast Fourier transforms (FFT). The output includes the IGM-Lyα transmitted flux fraction along quasar lines of sight and a catalog of high-column-density systems appropriately placed at high-density regions of the IGM. This output serves as input to additional code that superimposes the IGM tranmission on realistic quasar spectra, adds absorption by high-column-density systems and metals, and simulates instrumental transmission and noise. Redshift space distortions (RSD) of the flux correlations are implemented by including the large-scale velocity-gradient field in the FGPA resulting in a correlation function of F(λ) that can be accurately predicted. One hundred realizations have been produced over the 14,000 deg 2 DESI survey footprint with 100 quasars per deg 2 . The analysis of these realizations shows that the correlations of F(λ) follows the prediction within the accuracy of eBOSS survey. Here, the most time-consuming part of the mock production occurs before application of the FGPA, and the existing pre-FGPA forests can be used to easily produce new mock sets with modified redshift-dependent bias parameters or observational conditions.

79 ASTRONOMY AND ASTROPHYSICS↗

Study of Azimuthal Correlations in the Target Fragmentation Region in p, d, He, C+C, Ta and C+Ne, Cu Collisions at a Momentum of 4.2, 4.5 and 10 AGeV/c

Azimuthal correlations between the same type of particles (protons or pions) in the target fragmentation region was studied in d, He, C + C, Ta (4.2 AGeV/c), C + Ne, Cu (4.5AGeV/c) and p + C, Ta (10 GeV/c) interactions. The data were obtained from the SKM-200-GIBS streamer chamber and from Propane Bubble Chamber (PBL-500) systems utilized at JINR. Study of multiparticle azimuthal correlations offers unique information about space-time evolution of the interactions. Azimuthal correlations were investigated by using correlation function C($\Delta\phi$)=dN/d($\Delta\phi$), where $\Delta\phi$ represents the angle between the sums of transverse momenta vectors for particles emitted in the forward and backward hemispheres. For protons a "back-to back" ("negative") azimuthal correlations were observed in the above mentioned interactions. The absolute values of the correlation coefficient $|\xi|$ -- the slope parameter of C($\Delta\phi$), strongly depend on the mass number of the target ($A_T$) nuclei in the nucleon-nucleus and nucleus-nucleus collisions. Namely, $|\xi|$ -- decreases with increase of $A_T$ in p+C and p+Ta collisions, while $|\xi|$ decreases from d+C up to C+Ne and then almost does not change with increase of $A_P$, $A_T$ in (d+He)Ta, C+Cu and C+Ta collisions. For pions a "back-to-back" correlations were obtained for a light targets (C, Ne), and a "side-by-side" ("positive") correlations for a heavy targets (Cu, Ta). The $|\xi|$ insignificantly changes with increase of the momenta per nucleon and almost does not change with increase of $A_P$ and $A_T$. Models, used for description of the data -- the Ultra relativistic Quantum Molecular Dynamic (UrQMD) and Quark-Gluon String Model (QGSM), satisfactorily describe the obtained experimental results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Statistics of density maxima and the large-scale matter distribution

High peaks in Gaussian noise display enhanced clustering. The enhancement takes two forms: on large scales one obtains a linear amplification of the correlation function which is independent of scale. On smaller scales, but larger than the mass scale of the peaks themselves, a nonlinear (exponential) enhancement of the number density of high peaks in overdense regions arises. The large-scale correlations of Abell's rich clusters can be understood as a manifestation of this phenomenon. If the formation of bright galaxies favors the high overdensity peaks then the number of galaxies (per unit mass) in clusters and groups may be considerably enhanced. Consequences of these ideas for the density parameter and the large-scale matter distribution are discussed.

Kaiser, N.↗

Path Integral Monte Carlo Simulations of Warm Dense Matter and Plasmas

New path integral Monte Carlo simulation (PIMC) techniques will be developed and applied to derive the equation of state (EOS) for the regime of warm dense matter and dense plasmas where existing first-principles methods cannot be applied. While standard density functional theory has been used to accurately predict the structure of many solids and liquids up to temperatures on the order of 10,000 K, this method is not applicable at much higher temperature where electronic excitations become important because the number of partially occupied electronic orbitals reaches intractably large numbers and, more importantly, the use of zero-temperature exchange-correlation functionals introduces an uncontrolled approximation. Here we focus on PIMC methods that become more and more efficient with increasing temperatures and still include all electronic correlation effects. In this approach, electronic excitations increase the efficiency rather than reduce it. While it has commonly been assumed such methods can only be applied to elements without core electrons like hydrogen and helium, we recently showed how to extend PIMC to heavier elements by performing the first PIMC simulations of carbon and water plasmas [Driver, Militzer, Phys. Rev. Lett. 108 (2012) 115502]. Here we propose to continue this important development to extend the reach of PIMC simulations to yet heavier elements and also lower temperatures. The goal is to provide a robust first-principles simulation method that can accurately and efficiently study materials with excited electrons at solid-state densities in order to access parts of the phase diagram such the regime of warm dense matter and plasmas where so far only more approximate, semi-analytical methods could be applied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Notes on resonances and unitarity from celestial amplitudes

We study the celestial description of the O(N) sigma model in the large N limit as introduced by Coleman, Jackiw and Politzer. Focusing on three dimensions, we analyze the implications of a UV complete, all-loop order 4-point amplitude of pions in terms of correlation functions defined on the celestial circle. We find these retain many key features from the previously studied tree-level case, such as their relation to Generalized Free Field theories and crossing-symmetry, but also incorporate new properties such as IR/UV softness and S-matrix metastable states. In particular, to understand unitarity, we propose a form of the optical theorem that controls the imaginary part of the correlator based solely on the presence of these resonances. We also explicitly analyze the conformal block expansions and factorization of four-point functions into three-point functions. We find that summing over resonances is key for these factorization properties to hold. We end with some topics for future study.

1/N expansion↗

Imaging the wakes of jets with energy-energy-energy correlators

As the partons in a high energy jet propagate through the droplet of quark-gluon plasma (QGP) produced in a heavy-ion collision they lose energy to, kick, and are kicked by the medium. The resulting modifications to the parton shower encode information about the microscopic nature of QGP. A direct consequence, however, is that the momentum and energy lost by the parton shower are gained by the medium and, since QGP is a strongly coupled liquid, this means that the jet excites a wake in the droplet of QGP. After freezeout, this wake becomes soft hadrons with net momentum in the jet direction meaning that what an experimentalist later reconstructs as a jet includes hadrons originating from both the modified parton shower and its wake. This has made it challenging to find experimental observables that provide an unambiguous view of the dynamical response of a droplet of QGP to a jet shooting through it. Recent years have seen significant substantial advances in the theoretical and experimental understanding of the substructure of jets, in particular, using correlation functions, $\langle$ ($\vec{n}$ 1 )•••$\mathcal{E}$($\vec{n}$ k ) $\rangle$, of the energy flux operator in proton-proton collisions and, recently, in heavy-ion collisions. So far, such studies have focused primarily on the two-point correlator, which allows for the identification of the angular scale of the underlying dynamics. Higher-point correlators hold the promise of mapping out the dynamics themselves. In this paper we perform the first study of the shape-dependent three-point energy-energy-energy correlator in heavy-ion collisions. Using the Hybrid Model to simulate the interactions of high energy jets with the QGP medium, we show that the three-point correlator presents us with a striking new opportunity. We find that hadrons originating from wakes are the dominant contribution to the three-point correlator in the kinematic regime in which the three points are well-separated in angle, forming a roughly equilateral triangle. This equilateral region of the correlator is far from the region populated by collinear vacuum emissions, making it a canvas on which jet wakes are laid out, where experimentalists can map their shapes. Our work provides a key step towards the systematic use of energy correlators to image and unravel the dynamical response of a droplet of QGP that has been probed by a passing jet, and motivates numerous experimental and theoretical studies.

Jets and Jet Substructure↗

Scattering from Colloid-Polymer Conjugates with Excluded Volume Effect

This work presents scattering functions of conjugates consisting of a colloid particle and a self-avoiding polymer chain as a model for protein-polymer conjugates and nanoparticle-polymer conjugates in solution. The model is directly derived from the two-point correlation function with the inclusion of excluded volume effects. The dependence of the calculated scattering function on the geometric shape of the colloid and polymer stiffness is investigated. The model is able to describe the experimental scattering signature of the solutions of suspending hard particle-polymer conjugates and provide additional conformational information. This model explicitly elucidates the link between the global conformation of a conjugate and the microstructure of its constituent components.

36 MATERIALS SCIENCE↗

Can Kohn–Sham density functional theory predict accurate charge distributions for both single-reference and multi-reference molecules?

Dipole moments are the first moment of electron density and are fundamental quantities that are often available from experiments. An exchange–correlation functional that leads to an accurate representation of the charge distribution of a molecule should accurately predict the dipole moments of the molecule. It is well known that Kohn–Sham density functional theory (DFT) is more accurate for the energetics of single-reference systems than for the energetics of multi-reference ones, but there has been less study of charge distributions. Here, in this work, we benchmark 48 density functionals chosen with various combinations of ingredients, against accurate experimental data for dipole moments of 78 molecules, in particular 55 single-reference molecules and 23 multi-reference ones. We chose both organic and inorganic molecules, and within the category of inorganic molecules there are both main-group and transition-metal-containing molecules, with some of them being multi-reference. As one would expect, the multi-reference molecules are not as well described by single-reference DFT, and the functionals tested in this work do show larger mean unsigned errors (MUEs) for the 23 multi-reference molecules than the single-reference ones. Five of the 78 molecules have relatively large experimental error bars and were therefore not included in calculating the overall MUEs. For the 73 molecules not excluded, we find that three of the hybrid functionals, B97-1, PBE0, and TPSSh (each with less than or equal to 25% Hartree–Fock (HF) exchange), the range-separated hybrid functional, HSE06 (with HF exchange decreasing from 25% to 0 as interelectronic distance increases), and the hybrid functional, PW6B95 (with 28% HF exchange) are the best performing functionals with each yielding an MUE of 0.18 D. Perhaps the most significant finding of this study is that there exists great similarity among the success rate of various functionals in predicting dipole moments. In particular, of 39 functionals designed as general-purpose functionals and that do not have a global value of 100% HF exchange, the average MUE is 0.23 D, with a standard deviation of only 0.04 D. Among gradient approximations, which are especially interesting because of their speed and portability, the best overall performance is by PBE, HCTH/407, OLYP, OreLYP, and GAM, each with MUE of 0.22 D.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Temperature dependence of structural, dynamical, and electronic properties of amorphous Bi 2 Te 3 : an ab initio study

Bismuth telluride (Bi 2 Te 3 ) has garnered significant interest in thermoelectric applications and three-dimensional topological insulators due to its unique electronic, transport, and thermal properties. Bi 2 Te 3 and Sb 2 Te 3 chalcogenide compounds have the same crystal structure. While Sb 2 Te 3 has been shown to be a prototypical phase change memory (PCM) compound along the pseudobinary tie-line of Ge-Sb-Te alloys, whether Bi 2 Te 3 can also exhibit PCM functionality is still not well established. In this work, a systematic study on the structural, dynamical, and electronic properties of amorphous Bi 2 Te 3 during the quenching process has been performed by using ab initio molecular dynamics simulations. Pair correlation function, coordination number, bond-angle distribution functions, and a novel atomistic cluster alignment method are used to explore the structural characteristics of Bi 2 Te 3 as a function of temperature. Our study shows that there are many distorted octahedral clusters in amorphous Bi 2 Te 3 . In comparison with the local structures in Sb 2 Te 3 , we found that the degree of distortion of the octahedrons in the Bi 2 Te 3 system is smaller than that in Sb 2 Te 3 system. Moreover, the changes in the dynamical properties of Bi 2 Te 3 from liquid to glassy state are also explored. The approximate range of liquid-to-glass transition temperature is determined to be between 673 and 723 K. The electronic properties of Bi2Te 3 and Sb 2 Te 3 are also analysed by density-of-states and Bader charge calculations, both of them in glass state are semiconductors. Our studies provide useful insights into the local structure and dynamical properties of Bi 2 Te 3 at the atomistic level during the fast cooling process, and suggest that the compound can be a candidate for PCM materials.

36 MATERIALS SCIENCE↗

Cosmic shear power spectra in practice

Cosmic shear is one of the most powerful probes of Dark Energy, targeted by several current and future galaxy surveys. Lensing shear, however, is only sampled at the positions of galaxies with measured shapes in the catalog, making its associated sky window function one of the most complicated amongst all projected cosmological probes of inhomogeneities, as well as giving rise to inhomogeneous noise. Partly for this reason, cosmic shear analyses have been mostly carried out in real-space, making use of correlation functions, as opposed to Fourier-space power spectra. Since the use of power spectra can yield complementary information and has numerical advantages over real-space pipelines, it is important to develop a complete formalism describing the standard unbiased power spectrum estimators as well as their associated uncertainties. Building on previous work, this paper contains a study of the main complications associated with estimating and interpreting shear power spectra, and presents fast and accurate methods to estimate two key quantities needed for their practical usage: the noise bias and the Gaussian covariance matrix, fully accounting for survey geometry, with some of these results also applicable to other cosmological probes. Here, we demonstrate the performance of these methods by applying them to the latest public data releases of the Hyper Suprime-Cam and the Dark Energy Survey collaborations, quantifying the presence of systematics in our measurements and the validity of the covariance matrix estimate. We make the resulting power spectra, covariance matrices, null tests and all associated data necessary for a full cosmological analysis publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗

A redshift survey of IRAS galaxies. III - Reconstruction of the velocity and density fields in N-body model universes

N-body simulations of a 'cold dark matter' universe are presently used to calibrate the accuracy, and assess the limitations, of the procedure previously employed to predict the velocity field within 8000 km/sec of the Local Group through the application of linear gravitational theory to a full-sky, flux-limited sample of IRAS galaxies. The rms difference between the one-dimensional acceleration and velocity of field particles is an increasing function of local density; linear theory can in this way account for all but one-sixth of kinetic energy. A series of artificial IRAS catalogs closely matching the real sample in space density and clustering amplitude is constructed. Velocity correlation functions are used to demonstrate that the predicted velocity fields are in good agreement with the true velocity fields on large scales.

Davis, Marc↗

$DK$ and $D^* K$ scattering near threshold

We study the three $D_s$ quantum channels $J^P = 0^+$, $1^+$ and $2^+$ where experiments have identified the charm-strange states $D^*_{s0} (2317)$, $D_{s1}(2460)$, $D_{s1}(2536)$ near the $DK$ and $D^*K$ thresholds, and $D^*_{s2}(2573)$. We consider correlation functions for sets of $\overline q q$ operators and, for $J^P = 0^+$, $1^+$, also the $DK$ and $D^*K$ meson-meson interpolators and determine for these cases values of the elastic scattering amplitude. Constructing the full set of correlators requires propagators which connect any pair of lattice sites. For one ensemble of gauge configurations ($32^3\times 64$, $m_\pi\approx 156$ MeV) a stochastic distillation variant is employed and for another ensemble ($16^3\times 32$, $m_\pi\approx 266$ MeV) we use the full distillation method. Both, $D^*_{s0} (2317)$ and $D_{s1}(2460)$, are found as bound states below threshold, whereas $D_{s1}(2536)$, and $D^*_{s2}(2573)$ are identified as narrow resonances close to the experimental masses.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Transverse single-spin asymmetries: Challenges and recent progress

In this study, transverse single-spin asymmetries are among the most intriguing observables in hadronic physics. Though such asymmetries were already measured for the first time about four decades ago, their origin is still under debate. Here we consider transverse single-spin asymmetries in semi-inclusive lepton–nucleon scattering, in nucleon–nucleon scattering, and in inclusive lepton–nucleon scattering. It is argued that, according to recent work, the single-spin asymmetries for those three processes may be simultaneously described in perturbative QCD, where the re-scattering of the active partons plays a crucial role. A comparison of single-spin asymmetries in different reactions can also shed light on the universality of transverse momentum dependent parton correlation functions. In particular, we discuss what existing data may tell us about the predicted process dependence of the Sivers function.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Isotope effects in molecular structures and electronic properties of liquid water via deep potential molecular dynamics based on the SCAN functional

Here, Feynman path-integral deep potential molecular dynamics (PI-DPMD) calculations have been employed to study both light (H 2 O) and heavy water (D 2 O) within the isothermalisobaric ensemble. In particular, the deep neural network is trained based on ab initio data obtained from the strongly constrained and appropriately normed (SCAN) exchange-correlation functional. Because of the lighter mass of hydrogen than deuteron, the properties of light water are more influenced by nuclear quantum effect than those of heavy water. Clear isotope effects are observed and analyzed in terms of hydrogen-bond structure and electronic properties of water that are closely associated with experimental observables. The molecular structures of both liquid H 2 O and D 2 O agree well with the data extracted from scattering experiments. The delicate isotope effects on radial distribution functions and angular distribution functions are well reproduced as well. Our approach demonstrates that deep neural network combined with SCAN functional based ab initio molecular dynamics provides an accurate theoretical tool for modeling water and its isotope effects.

36 MATERIALS SCIENCE↗

Sequestration of Radionuclides in Metal–Organic Frameworks from Density Functional Theory Calculations

The high porosity, modularity, and synthetic diversity of metal-organic frameworks (MOFs) make them attractive candidate materials for selective sensing, separation, and sequestration of the radionuclides present in the nuclear wastes. Here, we use Density Functional Theory (DFT) calculations to determine the favorability of ion-exchange at the metal node of Tc, U, Th, Am, and Cm in Zr-, U-, and Th-based MOF clusters. A range of DFT methods, including several exchange-correlation functionals, DFT+U technique, relativistic effects, and magnetic effects are employed to establish robustness in the results. We also explore various reference states of ions including the vacuum, continuum water medium and explicit water molecules surrounding ions. Substitution of Tc, Am, Cm, and Th is believed to be energetically favorable in all the MOFs using ions surrounded by explicit water molecules as a reference.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Delocalization error poisons the density-functional many-body expansion

The many-body expansion is a fragment-based approach to large-scale quantum chemistry that partitions a single monolithic calculation into manageable subsystems. This technique is increasingly being used as a basis for fitting classical force fields to electronic structure data, especially for water and aqueous ions, and for machine learning. Here, we show that the many-body expansion based on semilocal density functional theory affords wild oscillations and runaway error accumulation for ion–water interactions, typified by F − (H 2 O) N with N ≳ 15. We attribute these oscillations to self-interaction error in the density-functional approximation. The effect is minor or negligible in small water clusters, explaining why it has not been noticed previously, but grows to catastrophic proportion in clusters that are only moderately larger. This behavior can be counteracted with hybrid functionals but only if the fraction of exact exchange is ≳50%, whereas modern meta-generalized gradient approximations including ωB97X-V, SCAN, and SCAN0 are insufficient to eliminate divergent behavior. Other mitigation strategies including counterpoise correction, density correction (i.e., exchange–correlation functionals evaluated atop Hartree–Fock densities), and dielectric continuum boundary conditions do little to curtail the problematic oscillations. In contrast, energy-based screening to cull unimportant subsystems can successfully forestall divergent behavior. These results suggest that extreme caution is warranted when the many-body expansion is combined with density functional theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗